This is for ngsxfem.

I would like to define my level set fully numerically, i.e. there will be no analytic function for it. For instance, I will start with an analytically defined level set, then perturb it with a FE function that is in the same finite element space as my interpolated level set function.

If I use P1 for the level set function, I can do this no problem. But for higher order, I’m not sure. I looked at the example file “unf_interf_prob.py”, but it isn’t clear.