The Performance of the Collocation and Galerkin Methods with Hermite Bi-Cubics

This paper presents a study of the performance of the collocation and Galerkin methods using Hermite bi-cubic basis functions. It is a sequel to the studies of Houstis et al. [6] and Weiser et al. [15]. The two methods have the linear systems solved by direct methods, band Gauss elimination or Cholesky factorization. The problem domain consists of linear, self-adjoint elliptic equations on two-dimensional rectangular domains. The measures of performance are computer time and memory needed to achieve moderate accuracy. The earlier study comparing finite element and finite difference methods observes that collocation uses less computer time than Galerkin. The second study gave detailed operation counts which support this observation, but also gave substantial experimental evidence to the contrary. We use a new implementation of the collocation method by E. N. Houstis which is tailored for rectangular domains (the one used in Houstis et al. [6] was designed for general domains). We use the same Galerkin impl...