Define the trace operator

[hgauss]

Hi !

Is there a way or a predefined function to define the trace operator ? For a given fem function u, I would like to define its restriction to the boundary.

Thank you.

[prj]

If your function is u, on a fespace Vh, then you can do:

Vh tr;
varf vTrace(unused, v) = on(boundary, unused = u);
tr[] = vTrace(0, Vh, tgv = -1);