In this paper, we describe a framework to compute expected convergence rates for residuals based on the Calder\'on identities for general second order differential operators for which fundamental solutions are known. The idea is that these rates could be used to validate implementations of boundary integral operators and allow to test operators separately by choosing solutions where parts of the Calder\'on identities vanish. Our estimates rely on simple vector norms, and thus avoid the use of hard-to-compute norms and the residual computation can be easily implemented in existing boundary element codes. We test the proposed Calder\'on residuals as debugging tool by introducing artificial errors into the Galerkin matrices of some of the boundary integral operators for the Laplacian and time-harmonic Maxwell's equations. From this, we learn that our estimates are not sharp enough to always detect errors, but still provide a simple and useful debugging tool in many situations.
翻译:暂无翻译