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 <[email protected]> > > Sent: Tuesday, August 02, 2016 12:28 AM > > To: Klaij, Christiaan > > Cc: [email protected]; Jed Brown > > Subject: Re: [petsc-users] block matrix without MatCreateNest > > > > On Mon, Aug 1, 2016 at 9:36 AM, Klaij, Christiaan <[email protected]> > 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:[email protected] | http://www.marin.nl MARIN news: http://www.marin.nl/web/News/News-items/MARIN-at-Monaco-Yacht-Show-September-28October-1.htm
