Numerical Integration Two-Dimensional Complex-Valued Functions GSL C wrapper
Institute of Mathematics and Informatics Bulgarian Academy of Sciences
Serdica Journal of Computing, Vol. 14, No 1-2, (2020), 1p-12p
Numerical integration is a key problem with numerous applications, including the boundary integral equation method (BIEM). GSL  provides powerful implementations of a broad variety of well-known numerical methods and algorithms. However, its routines have certain limitations: the numerical integration based on QUADPACK library  is implemented only in the case of one-dimensional functions. In this paper we present an extension of the GSL numerical integration routines for the special case of two-dimensional complex-valued functions. The presented approach is part of our effort to build a BIEM software system for solving a 3D elastodynamics problem for wave propagation in a continuously inhomogeneous half-space.