Accurate C++ BEM Library for accurate computation of integrals involving 3D Helmholtz and Laplace type kernels arising in BEM with constant elements. The aim of the library is to provide the proven set of methods that deliver guaranteed machine accuracy of BEM integrals.


The main idea is to convert surface integrals arising from BEM to contour integrals. This way, adaptive algorithms for numerical computation of integrals can be applied more effectively because it is easier to control the accuracy if we compute contour integrals only. In the case of BEM integrals involving Laplace kernels, completely analytical solution over constant element is found and used for computation of such integrals. This analytical solution for integrals involving Laplace kernels is not the sum of terms obtained by converting integrands to series.