Hello,
we cannot give you a self-contained C++ function for the element matrices, you need the ngsolve fem library at runtime.
We do (optional) code generation for the integration point, but shape functions are in the runtime.
The element function your are looking for is in
github.com/NGSolve/ngsolve/blob/master/f...mbolicintegrator.hpp
,
line 668:
NGS_DLL_HEADER virtual void
CalcElementMatrix (const FiniteElement & fel,
const ElementTransformation & trafo,
FlatMatrix<double> elmat,
LocalHeap & lh) const override;
It could be called from an extern fem code.
Joachim