Reuse LU decomposition for pressure correction.

Only with linear divergence
Add optional domain argument
Update examples