[petsc-users] A question on Schur complement

saturday luis luis.saturday at gmail.com
Fri Jan 26 13:18:16 CST 2018

Hi PETSc team:

I want to consult for a good way of doing an approximated Schur complement
matrix. So I just want to handle a matrix in the format C + B^t B, or
sometimes B^t B.

An easy way is to run with the MatCreateSchurComplement to create C + B^t
inv A B, wherein A is an identity matrix. The downside is that the matrix
free style does not allow for algebraic preconditioners.

Another way is to generate the matrix B^t B. For my case, I checked the
sparsity of the resulting matrix from MatMatMult and it is pretty sparse
still. I found a similar approach in
where the B^t B matrix is created explicitly (line 367-388). However, there
is a line of comment saying  " in real life this matrix would be build
directly without MatMatMult". So I just want to ask the best way of
assemblying and solving the B^t B matrix.


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180126/2027a343/attachment.html>

More information about the petsc-users mailing list