A method for solving a linear system of equations is required by both the interior and exterior boundary value problems. In this work the subroutine for this purpose is GLS and it is called by subroutines LIBEM2, LIBEM3, LIBEMA, LEBEM3 and LEBEMA. The subroutine is located in file GLS.FOR .
The use of the direct boundary element methods and the general nature of the boundary condition generally results in a linear system of equations of the following form


This is not a standard matrixvector problem but it can be turned into one by rearranging the rows. This solution procedure is carried out in subroutine GLS and the subroutine LINSL is called for solving the resulting standard problem by an LU factorisation method that is supplied without calling an external routine. The GLS code is currently being revised. Further information on the following link .