Comment 3 for bug 831050

Revision history for this message
Garth Wells (garth-wells) wrote :

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.