3-column integer matrix with 0-based vertex indices for each triangle
fem_order_max
integer; the highest operator order to compute
aniso
If non-NULL, a list(gamma, v). Calculates anisotropic structure
matrices (in addition to the regular) for \(\gamma\) and \(v\) for
an anisotropic operator \(\nabla\cdot H \nabla\), where
\(H=\gamma I + v v^\top\).
Currently (2023-08-05) the fields need to be given per vertex.