Just like Manuel Colera Rico, I would like to begin with existing (sub)matrices and put them directly into a matnest.
You seem to have understood that problem from the archive 2.5 years ago... If my memory is correct, it was an attempt to create a mat and switch between -mat_type nest and aij on the command line. Chris dr. ir. Christiaan Klaij | Senior Researcher | Research & Development MARIN | T +31 317 49 33 44 | c.kl...@marin.nl<mailto:c.kl...@marin.nl> | www.marin.nl<http://www.marin.nl> [LinkedIn]<https://www.linkedin.com/company/marin> [YouTube] <http://www.youtube.com/marinmultimedia> [Twitter] <https://twitter.com/MARIN_nieuws> [Facebook] <https://www.facebook.com/marin.wageningen> MARIN news: First autonomous manoeuvring vessel trials held on North Sea<http://www.marin.nl/web/News/News-items/First-autonomous-manoeuvring-vessel-trials-held-on-North-Sea.htm> ________________________________ From: Matthew Knepley <knep...@gmail.com> Sent: Thursday, April 11, 2019 2:16 PM To: Klaij, Christiaan Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] Error with KSPSetUp and MatNest On Thu, Apr 11, 2019 at 7:51 AM Klaij, Christiaan via petsc-users <petsc-users@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov>> wrote: Matt, As a happy MATNEST user, I got a bit worried when you wrote "we should remove MatCreateNest() completely". This would not remove any of the Nest functionality, just the direct interface to it, which is the problem. This happened last time I tried to use AIJ instead: https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2016-August/029973.html Has this problem been fixed in the meantime? This looks like a problem with the Nest implementation, not AIJ? Maybe I am misunderstanding it. Thanks, Matt Chris > I think we should remove MatCreateNest() completely. We should have 1 > mechanism for getting submatrices for creation, not 2, so we retain > MatCreateLocalRef(). Then if -mat_type nest was specified, it gives you > back an actually submatrix, otherwise it gives a view. This would do > everything that we currently do without this horrible MatNest interface > bubbling to the top. > > Matt > dr. ir. Christiaan Klaij | Senior Researcher | Research & Development MARIN | T +31 317 49 33 44 | mailto:c.kl...@marin.nl<mailto:c.kl...@marin.nl> | http://www.marin.nl MARIN news: http://www.marin.nl/web/News/News-items/42nd-FPSO-JIP-Week-April-812-Singapore.htm -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener https://www.cse.buffalo.edu/~knepley/<http://www.cse.buffalo.edu/~knepley/> Help us improve the spam filter. If this message contains SPAM, click here<https://www.mailcontrol.com/sr/VaxGuPX0yKPGX2PQPOmvUqEZpHzWxmiI7TcWSqBM5JDxEDemOSATW7gNYv05LkIR8iuercSaQ0xkxH4nYzTQYQ==> to report. Thank you, MARIN Support Group