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 _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : dolfin@lists.launchpad.net Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
