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.
Thanks,
Dave
On Wed, 8 Feb 2017 at 22:51, Manav Bhatia wrote:
> aha.. that might be it.
>
> Does that need to be called for the global matrix after each assembly of
aha.. that might be it.
Does that need to be called for the global matrix after each assembly of the
Jacobian blocks, or just once for the whole matrix?
-Manav
> On Feb 8, 2017, at 3:47 PM, Barry Smith wrote:
>
>
>> On Feb 8, 2017, at 3:40 PM, Manav Bhatia wrote:
>>
>> Hi,
>>
>> I ha
> On Feb 8, 2017, at 3:40 PM, Manav Bhatia wrote:
>
> Hi,
>
>I have a nested matrix with 2x2 blocks. The blocks (1,1) and (2,2) are AIJ
> matrices and blocks (1,2) and (2,1) and shell matrices. I am calling the code
> with the following arguments: -pc_type fieldsplit , and get the error
Hi,
I have a nested matrix with 2x2 blocks. The blocks (1,1) and (2,2) are AIJ
matrices and blocks (1,2) and (2,1) and shell matrices. I am calling the code
with the following arguments: -pc_type fieldsplit , and get the error shown
below.
I see that the error is complaining about the