Report error when 1 processor but 2 or more is OK

Dear all,
Here is my parallel code
try_movemesh.edp (4.0 KB)
movemesh.mesh (82.9 KB)

I find when I run it with 2 or more processors(ff-mpirun -n 2 try_movemesh.edp -wg -v 0
), it can work well.
However, when I choose only one processor(ff-mpirun -n 1 try_movemesh.edp -wg -v 0
), it report an error as


Does anyone know why?

I think you should set penal=-1.e-8 in order to prevent having a singular matrix (pressure not unique because of addition of a constant). When there are several processes the linear system is solved by approximations that eventually do not “see” that the matrix is singular.

Yes, it works. Many thanks, professor.