calc_meshdata

Calculates the affine transformations needed for integration in 2D/3D via the unit reference triangle/tetrahedron. In 2D the function takes as its input the matrices p,e,t given by initmesh/poimesh (PDEtool toolbox), and calculates the affine transformations for each triangle simultaneously in a vectorized manner. In 3D we use Tetgen meshes.
With a 2009 MacBook Pro laptop, calculating the affine transformations for a mesh of 2 million triangles takes about 1.2 seconds.

We also demonstrate how to use this data to perform fast vectorized calculation of integral values over 2D and 3D meshes.