[petsc-users] What do with singular blocks in block matrix preconditioning?
Thomas Witkowski
thomas.witkowski at tu-dresden.de
Thu Feb 16 05:49:28 CST 2012
I consider a 2x2 block matrix (saddle point) with the left upper block
being singular due to Neumann boundary conditions. The whole block
matrix is still non-singular. I worked on some ideas for block
preconditioning, but there is always some problem with the singular
block. All publications I know assume the block to be definite. There is
also some work on highly singular blocks, but this is here not the
case. Does some of you know papers about block preconditioners for some
class of 2x2 saddle point problems, where the left upper block is
assumed to be positive semi-definite?
From a more practical point of view, I have the problem that,
independently of a special kind of block preconditioner, one has always
to solve (or to approximate the solution) a system with the singular
block with an arbitrary right hand side. But in general the right hand
side does not fulfill the compatibility condition of having zero mean.
Is there a way out of this problem?
Thomas
More information about the petsc-users
mailing list