Method of Lines Transpose: An Implicit Solution to the Wave Equation

We present a new method for solving the wave equation implic-itly. Our approach is to discretize the wave equation in time, following the method of lines transpose, sometimes referred to as the transverse method of lines, or Rothe’s method. We differ from conventional meth-ods that follow this approach, in that we solve the resulting system of partial differential equations using boundary integral methods. Our algorithm extends to higher spatial dimensions using an alter-nating direction implicit (ADI) framework. Thus we develop a bound-ary integral solution, that is competitive with explicit finite differ-ence methods, both in terms of accuracy and speed. However, it pro-vides more flexibility in the treatment of source functions, and complex boundaries. We provide the analytical details of our one-dimensional method herein, along with a proof of the convergence of our schemes, in free space and on a bounded domain. We find that the method is uncon-ditionally stable, and achieves second order accuracy. A caveat of the analysis is the derivation of a unique and novel optimal quadrature method, which can be viewed as a Lax-type correction