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