SUBROUTINE MBEM{2 *or* 3 *or* A(}

real wavenumber,

description of boundary *and* set of interior solution points,

homogeneous boundary condition,

validation parameters,

characteristic wavenumbers and their corresponding
mode shape at boundary points and at interior solution points
*(solution)*,

working space )

The subroutines require input a geometrical description of the boundary of the domain (as covered in Chapter 2) and a list of the points in the interior (from which the values of the f at those points allow us to construct the mode shapes), the homogeneous boundary condition, the subroutine returns the characteristic wavenumbers and their corresponding mode shape at boundary points and at interior solution points. The use of the subroutines are demonstrated by the programs HMBEM2_T, HMBEM3_T and HMBEMA_T in the next Section.

In the subroutines HMBEM2, HMBEM3 and HMBEMA,
the technique
employed for deriving the polynomial approximation
(6.12)
involves computing A_{k} at the m+1 Chebyshev (¥ norm)
interpolation points for any selected wavenumber range [k_{A}, k_{B}]. The
coefficient matrices
A_{[0]}, A_{[1]}, ..., A_{[m]} in (6.12) are obtained through
Newton's divided differences using the value of A_{k} at the selected values of
k in [k_{A}, k_{B}].

The subroutine parameters that specify the interior
two-dimensional Helmholtz
problem must be set up in the main program. Let this be called
*MAIN*.FOR.
The following files must be linked together to construct the
complete program:

*MAIN*.FOR (and files containing any user-defined sub-programs),

HIBEM2.FOR,

H2LC.FOR, the file for computing the discrete operators - see Chapter 3,

FNHANK.FOR, the Hankel function - see Appendix 5,

INTEIG.FOR, the file for computing the solution to a linear
system - see Appendix 4,

GEOM2D.FOR, the file for 2D geometry - see Appendix 6.

The subroutine parameters that specify the exterior
three-dimensional Helmholtz
problem must be set up in the main program. Let this be called
*MAIN*.FOR.
The following files must be linked together to construct the
complete program:

*MAIN*.FOR (and files containing any user-defined sub-programs),

HEBEM3.FOR,

H3LC.FOR, the file for computing the discrete operators - see Chapter 3,

INTEIG.FOR, the file for computing the solution to a linear
system - see Appendix 4,

GEOM3D.FOR, the file for 3D geometry - see Appendix 6.

The subroutine parameters that specify the exterior
three-dimensional axisymmetric Helmholtz
problem must be set up in the main program. Let this be called
*MAIN*.FOR.
The following files must be linked together to construct the
complete program:

*MAIN*.FOR (and files containing any user-defined sub-programs),

HIBEMA.FOR,

H3ALC.FOR, the file for computing the discrete operators - see Chapter 3,

INTEIG.FOR, the file for computing the solution to a linear
system - see Appendix 4,

GEOM2D.FOR, the file for 2D geometry - see Appendix 6,

GEOM3D.FOR, the file for 3D geometry - see Appendix 6.

Return to boundary-element-method.com/helmholtz

Return to boundary-element-method.com