<div dir="ltr"><div dir="ltr">On Thu, Apr 11, 2019 at 7:51 AM Klaij, Christiaan via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Matt,<br>
<br>
As a happy MATNEST user, I got a bit worried when you wrote "we<br>
should remove MatCreateNest() completely". </blockquote><div><br></div><div>This would not remove any of the Nest functionality, just the direct interface to it, which is the problem.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">This happened last<br>
time I tried to use AIJ instead:<br>
<br>
<a href="https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2016-August/029973.html" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2016-August/029973.html</a><br>
<br>
Has this problem been fixed in the meantime?<br></blockquote><div><br></div><div>This looks like a problem with the Nest implementation, not AIJ? Maybe I am misunderstanding it.</div><div><br></div><div>   Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Chris<br>
<br>
> I think we should remove MatCreateNest() completely. We should have 1<br>
> mechanism for getting submatrices for creation, not 2, so we retain<br>
> MatCreateLocalRef(). Then if -mat_type nest was specified, it gives you<br>
> back an actually submatrix, otherwise it gives a view. This would do<br>
> everything that we currently do without this horrible MatNest interface<br>
> bubbling to the top.<br>
><br>
>    Matt<br>
><br>
<br>
<br>
dr. ir. Christiaan Klaij  | Senior Researcher | Research & Development<br>
MARIN | T +31 317 49 33 44 | mailto:<a href="mailto:C.Klaij@marin.nl" target="_blank">C.Klaij@marin.nl</a> | <a href="http://www.marin.nl" rel="noreferrer" target="_blank">http://www.marin.nl</a><br>
<br>
MARIN news: <a href="http://www.marin.nl/web/News/News-items/42nd-FPSO-JIP-Week-April-812-Singapore.htm" rel="noreferrer" target="_blank">http://www.marin.nl/web/News/News-items/42nd-FPSO-JIP-Week-April-812-Singapore.htm</a><br>
<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>