Forum Message

 

 

We have moved the forum to https://forum.ngsolve.org . This is an archived version of the topics until 05/05/23. All the topics were moved to the new forum and conversations can be continued there. This forum is just kept as legacy to not invalidate old links. If you want to continue a conversation just look for the topic in the new forum.

Notice

The forum is in read only mode.

Adaptive mesh refinement: Interpolation operator

More
4 years 5 months ago #2146 by joachim
Hi Anand,

yes, HDiv is missing. I have two recommendations:

1. If you are in 2D, you can use HCurl instead of HDiv, rotate your field by 90 deg, and replace div by curl. For lowest order HCurl we have the prolongation. This is a pure Python solution.

2. Copy the C++ code from HCurlSpace to HDivFESpace. Also prolongations for lowest order RT tetrahedral elements can be done similarly. Have a look into

ngsolve/multigrid/prolongation.hpp, line 183, class EdgeProlongation
and
comp/hcurlhdivfes.cpp, line 140 for building the edge hierarchy

We only assume from the mesh interface that it provides the hierarchy of vertices and elements, no edges or faces, this is more troublesome for p-version prolongations.
I have heard from people who have added at least second order H1 prolongations.

Joachim
More
4 years 5 months ago #2147 by anandes
Hi Joachim,

Thank you for your clarifications and your pointers on code modifications.

I'm working on problems in 3D domains and so will have to look into option 2 from your recommendations.

I'll write back to you when I have more questions in this regard.

Thank you,
Anand
Time to create page: 0.114 seconds