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

[Umut Tabak]

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