A package program for corner issue in BEM

The boundary element method (BEM) is the most important alternative to the finite element method (FEM) due in part to its high precision. However, the high precision will be compromised if corners are not properly treated. The discontinuous element method (DEM), most commonly used in corner treatment, is satisfactory in numerical properties but hard to couple itself with FEM based codes because the FE-domain and the BE-domain do not share the interface nodes. A package program, designated by PP, is proposed to deal with corners on the boundary or the material interface, where, in each element connecting to a corner, another source point instead of the corner is set. The source point has no its own independent variables but is interpolated with the element nodes. With half of the number of degrees of freedom used by DEM, PP has the same precision as DEM, and exceeds any other existing methods in numerical properties.