[petsc-users] ksp solve error with nested matrix

[Barry Smith]

  Hmm, we might need to overload the MatAssemblyEnd() on the inner matrices to have them update the state of the outer matrix.

  Barry

> On Feb 9, 2017, at 11:33 AM, Lawrence Mitchell <lawrence.mitchell at imperial.ac.uk> wrote:
> 
> 
> 
> On 09/02/17 17:29, Barry Smith wrote:
>> 
>>> On Feb 8, 2017, at 11:39 PM, Dave May <dave.mayhem23 at gmail.com> wrote:
>>> 
>>> Any time you modify one of the submats, you need to call assembly begin/end on that sub matrix AND on the outer matnest object.
>> 
>>   Weird, and prone to errors it seems to me. Perhaps this needs to be rethought
> 
> If you just call assembly begin/end on the outer object, that takes
> care of calling it on the nested matrices.  But that's still not
> ideal, I think.
> 
> Lawrence
>