[petsc-users] Error with KSPSetUp and MatNest

[Matthew Knepley]

On Thu, Apr 11, 2019 at 9:58 AM Klaij, Christiaan <C.Klaij at marin.nl> wrote:

> 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.
>
> Okay, I will make an example once I get back from my trip and fix anything
that does not work.

  Thanks,

     Matt

> Chris
>
> dr. ir. Christiaan Klaij | Senior Researcher | Research & Development
> MARIN | T +31 317 49 33 44 | C.Klaij at marin.nl | www.marin.nl
>
> [image: LinkedIn] <https://www.linkedin.com/company/marin> [image:
> YouTube] <http://www.youtube.com/marinmultimedia> [image: Twitter]
> <https://twitter.com/MARIN_nieuws> [image: 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 <knepley at gmail.com>
> *Sent:* Thursday, April 11, 2019 2:16 PM
> *To:* Klaij, Christiaan
> *Cc:* petsc-users at 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 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/>
>
>
> *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*
>
>

-- 
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/b473f4ae/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: imagea896fb.PNG
Type: image/png
Size: 293 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190411/b473f4ae/attachment-0004.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image0539c8.PNG
Type: image/png
Size: 331 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190411/b473f4ae/attachment-0005.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: imagefee82c.PNG
Type: image/png
Size: 333 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190411/b473f4ae/attachment-0006.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image188b81.PNG
Type: image/png
Size: 253 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190411/b473f4ae/attachment-0007.png>