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.
@article{arxiv.2410.12026,
title = {A Generic MATLAB Toolbox to Approximate PDEs Using Computational Geometry},
author = {Kiefer Green and Harbir Antil},
journal= {arXiv preprint arXiv:2410.12026},
year = {2024}
}