Hello,
I have a question regarding the transformation between reference and physical elements in NGSolve. Say I have the following spaces:
V = HDiv(mesh, order=order)
Q = L2(mesh, order=order)
M = FacetFESpace(mesh, order=order)
fes = FESpace([V, Q, M, M])
u, p, uxhat, uyhat = fes.TrialFunction()
then I assume that a Piola map is used for u when mapping between reference and physical element.
But what if we have
V = L2(mesh, order=order)
Q = L2(mesh, order=order)
M = FacetFESpace(mesh, order=order)
fes = FESpace([V, V, Q, M, M])
ux, uy, p, uxhat, uyhat = fes.TrialFunction()
u = CoefficientFunction((ux, uy))
What mapping is used now for u? If F is an affine mapping from reference to physical element, is the mapping just \hat{ux} = ux \circ F and \hat{uy} = uy \circ F ?
Thanks,
Sander