Hi all,
I have two overlapping meshes, mesh_coarse
(with quadratic basis functions) and mesh_fine
(with linear basis functions):
mesh_fine = refine(mesh_coarse)
V1 = FunctionSpace(mesh_fine, "CG", 1)
V2 = FunctionSpace(mesh_coarse, "CG", 2)
I now want to map vectors between both meshes (in parallel). Right now, I’m using LagrangeInterpolator.interpolate
, which works. This is, however, very slow and inefficient since both meshes have the same number of DOFs and the mapping actually stays the same.
Does anybody how to map between two vectors in V1
and V2
in parallel?
Thank you very much!