[petsc-users] Use boomeramg to solve system of PDEs

[Jed Brown]

On Thu, Feb 2, 2012 at 15:43, Mark F. Adams <mark.adams at columbia.edu> wrote:

> Use MatSetBlockSize(mat,ndof); and that info will get passed down to HYPRE.


This assumes that you have ordered your equations so that all components
are next to each other.

This is usually practical for systems in H^1 (e.g. compressible
elasticity), but for systems with constraints (mixed methods for Stokes) or
systems in H(div) or H(curl), you need more specialized methods. Hypre has
some support for H(curl), but it's for a specific discretization and you
would have to use their custom interface. Unfortunately, those systems
always tend to be a lot of work, so the alternative is to use fieldsplit
(also "auxilliary space") preconditioning to transform those systems into
ones that can be solved with standard methods.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120202/45622080/attachment.htm>