Re: [petsc-users] pcfieldsplitting not working for VI

2016-11-29 Thread Lawrence Mitchell
I suspect this is not a petsc issue per se. Shall we migrate this to the firedrake list? Lawrence > On 29 Nov 2016, at 07:40, Justin Chang wrote: > > FWIW I ended up doing this: > > rs_solver.setDM(W._dm()) > > and was left with this error: > > Traceback (most recent

Re: [petsc-users] pcfieldsplitting not working for VI

2016-11-28 Thread Justin Chang
FWIW I ended up doing this: rs_solver.setDM(W._dm()) and was left with this error: Traceback (most recent call last): File "sphere.py", line 180, in virs() File "sphere.py", line 163, in virs rs_solver.solve(None,sol_vec) File "PETSc/SNES.pyx", line 537, in

Re: [petsc-users] pcfieldsplitting not working for VI

2016-11-28 Thread Justin Chang
Lawrence, I added the following line: rs_solver.setDM(a._dm) And it gives this error: Traceback (most recent call last): File "sphere.py", line 128, in rs_solver.setDM(a._dm) AttributeError: 'Form' object has no attribute '_dm' Thanks, Justin On Tue, Nov 29, 2016 at 1:12 AM, Lawrence

Re: [petsc-users] pcfieldsplitting not working for VI

2016-11-28 Thread Lawrence Mitchell
> On 29 Nov 2016, at 02:02, Justin Chang wrote: > > Why is this happening? It seems to me VINEWTONRSLS is not understanding that > I have a two-field formulation. For the built in solvers, we set up a dmshell on the snes that knows how to do field splitting. Because