[petsc-users] Error with KSPSetUp and MatNest

[Matthew Knepley]

On Thu, Apr 11, 2019 at 7:51 AM Klaij, Christiaan via petsc-users <
petsc-users at 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.Klaij at 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/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190411/8d989374/attachment-0001.html>