Hello,

I would like to use FreeFEM to solve a PDE with an infinte number of solutions, and recover a unique solution by imposing a constraint.

I found some other post on the forum linking to this page: http://jolivet.perso.enseeiht.fr/FreeFem-tutorial/section_8/example5.edp where an example is given for the parallel version. I tried, but not managed to make the example work in plain FreeFem. I keep getting some weird “segmentation faults” and I don’t see where is the problem.

In particular, the command changeNumbering does not seem to work (I don’t know if this is important). Anyway, if it is not too much trouble, would it be possible to have a working code for the simple case:

-\Delta u = f, u \in H^1(\omega), \int_\omega u=0, \partial_n u = 0.

Thank you!