Methods for singularity treatment in Bempp?

Dear Bempp users,

I am a graduate student specializing in boundary integral equations. While reading the Bempp book, I found that the explanation regarding the treatment of singularities in boundary integral operators is only partial: “The Green’s function will have a singularity when x=y, so internally Bempp will use appropriate singular quadrature rules to handle this.” I would like to understand the specific treatment methods in Bempp for these cases or the literature on these methods, including the handling of singularities in the Green’s functions for Laplace, Helmholtz, and especially Maxwell. I would be very grateful if anyone could provide information.

Thank you in advance for your time, and I am looking forward to your advice.