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.
Examples of vectorized calculation of integrals.
Added a function giving 2d integration quadratures.
The same in 3D too.
Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.