← Back to team overview

dolfin team mailing list archive

Solving singular systems

 

I've added a demo for solving singular systems using a Krylov solver via PETSc:

   demo/undocumented/singular

This can provide a good alternative to, for example, pinning the
pressure in a Stokes problem or using a global Lagrange multiplier. In
contrast to pinning or Lagrange multipliers, it doesn't mess up
preconditioners or introduce a global row to the matrix.

Does anyone know if Trilinos/AztecOO has similar functionality to
PETSc for solving singular problems?

Garth


Follow ups