def pressure_projection(self, state):
    'Enforce the incompressibility constraint.'
    raise NotImplementedError
Enforce the incompressibility constraint.
def pressure_projection(self, state):
    
    raise NotImplementedError