Factors affecting the accuracy of the boundary element method in the forward problem. I. Calculating surface potentials

A comprehensive review of factors affecting the accuracy of the boundary element method (BEM) for calculating surface potentials is presented. A relative-error statistic is developed which is only sensitive to calculation errors that could affect the inverse solution for source position, and insensitive to errors that only affect the solution for source strength. The factors considered in this paper are numerical approximations intrinsic to the BEM, such as constant-potential versus linear-potential basis functions and sharp-edged versus smooth-surfaced volumes; aspects of the volume conductor including the volume shape, density of surface elements, and element shape; source position and orientation; and effects of "refinements" in the numerical methods. The effects of these factors are considered in both smooth-shaped (spheres and spheroids) and sharp-edged (cubes) volume conductors. This represents the first attempt to assess the effects of many of these factors pertaining to the numerical methods commonly used in fields such as electrocardiography (ECG) and electroencephalography (EEG). Strategies for obtaining the most accurate solutions are presented.