P1ISOP2-P1

Works like a charm !
Thank you very much,

Dear Christopher,
I'm working on the C++ side, and I was wondering is there a way to get the coordinates of the vertex of the element inside the function CalcShape and CalcDShape.
Kind regards,
Umberto

Hi Umberto,

if you have a MappedIntegrationPoint, you can get the ElementTransformation from it, and from the ElementTransformation you can get the Mesh and the ElementId.

A short example is attached.

It uses the brand-new feature of including C++ code in Python scripts.

Joachim

Hi Umberto,

we can now do it by refining the mesh, and keeping the pressure in the range of the prolongation operator.

I added the method to the Stokes-tutorial, available in the latest nightly docu:

docu.ngsolve.org/nightly/i-tutorials/uni....html#(P1-iso-P2)-P1

Joachim

Dear Joachim,

thank you very much for the example on the prolongation operator and for the example on the class MappedIntegrationPoint !!

Kind regards,
Umberto

Dear Joachim,
is there a way to implement the prolongation operator also for the custom elements that I've implemented using mylittlengsolve ?
Kind regards,
Umberto