hp-adaptive loop using NGSpy


I’ve been using the Netgen/NGSolve python interface now for a while to solve 3D PDE’s using an adaptive mesh refinement loop.

Now I was wondering if there is some built-in functionality that would allow me to extend my programs to support a hp-adaptive loop. And if there is, are there any code examples that demonstrate its usage?

Greetings Patrik