A numerical integration scheme for special finite elements for the Helmholtz equation

The theory for integrating the element matrices for rectangular, triangular and quadrilateral finite elements for the solution of the Helmholtz equation for very short waves is presented. A numerical integration scheme is developed. Samples of Maple and Fortran code for the evaluation of integration abscissae and weights are made available. The results are compared with those obtained using large numbers of Gauss-Legendre integration points for a range of testing wave problems. The results demonstrate that the method gives correct results, which gives confidence in the procedures, and show that large savings in computation time can be achieved. Copyright (C) 2002 John Wiley Sons, Ltd.