[petsc-users] block matrix without MatCreateNest

[Klaij, Christiaan]

Jed,

Just a reminder, you haven't responded to this thread, notably
Matt's question below whether you can fix the L2G mapping.

Chris

> > From: Matthew Knepley <knepley at gmail.com>
> > Sent: Tuesday, August 02, 2016 12:28 AM
> > To: Klaij, Christiaan
> > Cc: petsc-users at mcs.anl.gov; Jed Brown
> > Subject: Re: [petsc-users] block matrix without MatCreateNest
> >
> > On Mon, Aug 1, 2016 at 9:36 AM, Klaij, Christiaan <C.Klaij at marin.nl>
> wrote:
> >
> >     Matt,
> >
> >
> >     1) great!
> >
> >
> >     2) ??? that's precisely why I paste the output of "cat
mattry.F90"
> in the emails, so you have a small example that produces the errors I
> mention. Now I'm also attaching it to this email.
> >
> > Okay, I have gone through it. You are correct that it is completely
> broken.
> >
> > The way that MatNest currently works is that it trys to use L2G
mappings
> from individual blocks
> > and then builds a composite L2G map for the whole matrix. This is
> obviously incompatible with
> > the primary use case, and should be changed to break up the full L2G
> into one for each block.
> >
> > Jed, can you fix this? I am not sure I know enough about how Nest
works.
> >
> >    Matt
> >
> >     Thanks,
> >
> >     Chris


dr. ir. Christiaan Klaij  | CFD Researcher | Research & Development
MARIN | T +31 317 49 33 44 | mailto:C.Klaij at marin.nl | http://www.marin.nl

MARIN news: http://www.marin.nl/web/News/News-items/MARIN-at-Monaco-Yacht-Show-September-28October-1.htm