HBE: HELMHOLTZ BOUNDARY ELEMENTS

by boundary-element-method.com


Click on routine to view source. Key [File] / [Save As] to download.

SubroutineTest Program
2-dimensional H2LC.FOR H2LC_T.FOR
3-dimensional H3LC.FOR H3LC_T.FOR
axisymmetric H3ALC.FOR H3ALC_T.FOR

The Subroutines
The Fortran subroutines H2LC.FOR, H3LC.FOR, H3ALC.FOR evaluate the discrete form of the Helmholtz integral operators L_k, M_k, M_k^t and N_k for two- dimensional, three-dimensional and three-dimensional axisymmetric problems.The subroutines have been designed to be easy to use, reliable and efficient. The subroutines are also flexible in that the quadrature rule is defined as a parameter and the library functions (such as the Hankel, exponential and square root functions) are called from external routines.

Test Problems
The subroutines are demonstrated on test problems arising from the solution of the Neumann problem exterior to a closed boundary via the Burton and Miller equation. The test problems are enacted through files H2LC_T.FOR, H3LC_T.FOR, H3ALC_T.FOR.

Some of the test programs call NAG routines.
For an alternative GLRULE routine, see GLRULE2.FOR
For an alternative FNHANK routine, see FNHANK2.FOR
A code for complex k is available, written by Nick Yas'ko: FNHANKNY.FOR

An on-line Manual is available.

The manual is an on-line version of the research paper describing the codes:
Fortran codes for computing the discrete Helmholtz integral operators
by Stephen Kirkup
Advances in Computational Mathematics 9 (1998) 391-409.