Re: [petsc-users] ksp solve error with nested matrix

2017-02-11 Thread Matthew Knepley
On Sat, Feb 11, 2017 at 11:27 AM, Manav Bhatia wrote: > I am now getting this error, which is complaining for an operation for the > shell matrix. I have only specified the following operation for the matrix: > “MATOP_MULT”. > > I am calling the solver with -pc_type

Re: [petsc-users] ksp solve error with nested matrix

2017-02-11 Thread Manav Bhatia
I am now getting this error, which is complaining for an operation for the shell matrix. I have only specified the following operation for the matrix: “MATOP_MULT”. I am calling the solver with -pc_type fieldsplit. What can I do to move forward with the solution? Thanks, Manav

Re: [petsc-users] ksp solve error with nested matrix

2017-02-09 Thread Lawrence Mitchell
On 09/02/17 17:29, Barry Smith wrote: > >> On Feb 8, 2017, at 11:39 PM, Dave May 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

Re: [petsc-users] ksp solve error with nested matrix

2017-02-09 Thread Barry Smith
> On Feb 8, 2017, at 11:39 PM, Dave May 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 > >

Re: [petsc-users] ksp solve error with nested matrix

2017-02-08 Thread Dave May
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

Re: [petsc-users] ksp solve error with nested matrix

2017-02-08 Thread Manav Bhatia
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

Re: [petsc-users] ksp solve error with nested matrix

2017-02-08 Thread Barry Smith
> 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

[petsc-users] ksp solve error with nested matrix

2017-02-08 Thread Manav Bhatia
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