PETSc's solver for block matrices PCFieldsplit http://www.mcs.anl.gov/petsc/petsc-3.3/docs/manualpages/PC/PCFIELDSPLIT.html#PCFIELDSPLIT is what you need. More information on fieldsplit is given in section 4.5 of the manual.
----- Original Message ----- > Dear all, > > I have a linear system in blocks that I would like to solve. > > The system is written as > > [ A C^T ] ? b_1 > ? ? ? ? ? ? ? ?= > [ C ? ? ? B] ? b_2 > > where C matrix results from a coupling condition in fluid-structure > interaction and rather sparse depending on the size of the coupling > surfaces. A and B are also sparse symmetric matrices. > > I was wondering if I could find a solution to this problem without > factorizing the complete operator matrix but only using, somehow, A > and > B and using some algebraic tricks and the sparsity of C. > > Any pointers/ideas are appreciated. > > Best regards, > > Umut >
