Owe Axelsson ¹, Evgeny Glushkov ², Natalia Glushkova ²

Previous theoretical and computational investigations have shown high efficiency of the local Green's function method for the numerical solution of singularly perturbed problems with sharp boundary layers. However, in several space variables those functions, used as projectors in the Petrov-Galerkin scheme, cannot be derived in a closed analytical form. This is an obstacle for the application of the method when applied to multi-dimensional problems. The present work proposes a semi-analytical approach to calculate the local Green's function, which opens a way to effective practical application of the method. Besides very accurate approximation, the matrix stencils obtained with these functions allow the use of fast and stable iterative solution of the large sparse algebraic systems that arise from the grid-discretization. The advantages of the method are illustrated by numerical examples.