which doesn't not work in parallel. I will add an error to dolfin::refine.
Second problem looks to be that the demo uses ILU preconditioning. This can be problematic with the zero pressure block in the matrix. There are some options to use a non-zero shift in the preconditioner.
First problem is that the demo make a call to
mesh = refine(mesh)
which doesn't not work in parallel. I will add an error to dolfin::refine.
Second problem looks to be that the demo uses ILU preconditioning. This can be problematic with the zero pressure block in the matrix. There are some options to use a non-zero shift in the preconditioner.