Field interpolation in p-refinement
How can I transfer a p
-refined solution to Gmsh for visualization?
For a give p
, I understand and successfully used
gmsh.view.setInterpolationMatrices(..., coef, exp, ...)
My problem is that in a p
-refined mesh, a certain element type ("Quadrangle", "Tetrahedron", etc.) will have different polynomial bases, based on the value of p
. However, the setInterpolationMatrices
method assumes the same change-of-basis matrices coef
and exp
for every element of type "Quadrangle", "Tetrahedron", etc.
Is there a way to tell Gmsh to use (coef_1, exp_1)
for element_set_1
, ..., (coef_n, exp_n)
for element_set_n
within the same element type? That would allow me to use the very powerful visualization mesh adaptation strategy of Gmsh in a p
-refined context (relevant in CFD, for instance).