I’ve been playing around with FF++ to solve some problems quite successfully. We are using a hybrid discontinuous Galerkin formulation. Currently, we are solving a time-dependent Cahn-Hilliard-like problem.
However, we were not able to implement the static condensation (Schur complement) approach to reduce the size of the global linear system.
I’ve looked around and in the documentation, but I haven’t found many tips about this.
Also I’ve seen papers citing FreeFEM++ and which seem to use this approach.
I wonder if there is an example somewhere so that I could learn how to do it on FF++ efficiently.
Thanks for your attention.