RESTRICTION/PROLONGATION operators

Question

Hello,
Is there any way how to obtain restriction/prelongation operators in fenics? I am building solver, where I need to partition my domain into non-overlaping subdomains, solve problem on every subdomain and from solution corrections build preconditioner to obtian global solution.

Therefore, I need easy way how to switch between subdomains and global solution.

Thanks for your help!

Answer 1

Hello,

I'll chime in here to say that I am also interested in this topic; specifically if there is a natural way to do this in fenics, and if the same is possible on a distributed mesh

I'll add that I have had some luck using a vertex_to_dof_map, restricting on the nodes using the classical boolean restriction operator, and then using a dof_to_vertex_map to get back into the finite element space.

Ben