A Generic MATLAB Toolbox to Approximate PDEs Using Computational Geometry

Abstract

This article introduces a general purpose framework and software to approximate partial differential equations (PDEs). The sparsity patterns of finite element discretized operators is identified automatically using the tools from computational geometry. They may enable experimentation with novel mesh generation techniques and could simplify the implementation of methods such as multigrid. We also implement quadrature methods following the work of Grundmann and Moller. These methods have been overlooked in the past but are more efficient than traditional tensor product methods. The proposed framework is applied to several standard examples.