[petsc-users] advice on a linear system solution for a matrix in blocks

[Shri]

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
>