[petsc-users] Solving/creating SPD systems

Jed Brown jed at jedbrown.org
Fri Dec 11 16:02:49 CST 2015

Justin Chang <jychang48 at gmail.com> writes:

> Pardon me for my apparent lack of understanding over what may be simple
> concepts, but why is div[u]*div[v] singular in the context of LSFEM?

The corresponding strong form is grad(div(u)).  If u has n*d entries (in
d dimensions -- it is a vector), div(u) has only n entries.  That means
grad(div(u)) must have a null space of dimension at least n*(d-1).

If you a matrix A and a singular matrix B, then

  A + c*B

cannot be well-conditioned for sufficiently large c.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151211/90e4f721/attachment.pgp>

More information about the petsc-users mailing list