[petsc-users] A question on Schur complement

saturday luis luis.saturday at gmail.com
Fri Jan 26 16:28:14 CST 2018

Thanks, Jed.

It seems that MatCreateLRC requires two dense matrices for U and V
according to the page

Not sure if it is applicable to sparse matrices.

2018-01-26 12:09 GMT-08:00 Jed Brown <jed at jedbrown.org>:

> saturday luis <luis.saturday at gmail.com> writes:
> > 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.
> MatCreateLRC is more natural in this context.
> Even if you solve with LSQR, you'll likely need some explicit
> representation for preconditioning.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180126/0013a331/attachment.html>

More information about the petsc-users mailing list