Parallel coordinate free implementation of local meshless method

This paper presents an implementation of a Meshless Local Strong Form Method that allows users to write elegant code expressed in terms of abstract mathematical objects, such as operators and fields, consequently avoiding working directly with matrix and array indices, which is tedious and error prone. This is achieved by using object oriented programming techniques for definition of abstract concepts and leveraging C++'s powerful templating mechanism. It is demonstrated that code written this way has little-to-no performance overhead compared to classical numerical code while being more expressive and readable, which gravely shortens model development and testing phases. The overall functionality of presented implementation is illustrated on numerical examples from classical thermodynamics, linear elasticity and fluid dynamics in one, two and three dimensions.