Numerical conformal mapping based on the generalised conjugation operator

An iterative procedure for numerical conformal mapping is presented which imposes no restriction on the boundary complexity. The formulation involves two analytically equivalent boundary integral equations established by applying the conjugation operator to the real and the imaginary parts of an analytical function. The conventional approach is to use only one and ignore the other equation. However, the discrete version of the operator using the boundary element method (BEM) leads to two non-equivalent sets of linear equations forming an over-determined system. The generalised conjugation operator is introduced so that both sets of equations can be utilised and their least-square solution determined without any additional computational cost, a strategy largely responsible for the stability and efficiency of the proposed method. Numerical tests on various samples including problems with cracked domains suggest global convergence, although this cannot be proved theoretically. The computational efficiency appears significantly higher than that reported earlier by other investigators.