Accurate field computation with the boundary element method

Interface problems in magnetostatics are formulated as boundary integral equations of the second kind involving the appropriate scalar (no current sources) or vector potentials. The boundary element method (BEM), which employs parametric representation of surfaces and sources, is used to solve some two-dimensional examples by way of illustration. A novel approach, automated to address Green´s functions singularities over arbitrarily-shaped geometries, is introduced.