<div dir="ltr">Hi all,<div><br></div><div>So I am running into some issues with the SNESVINEWTONRSLS solver. I originally had this done in firedrake, but have ported it to petsc4py. Attached is the code as well as the two required input files.</div><div><br></div><div>First issue, when I run the program like this:</div><div><br></div><div>python testVI.py psiA_1 psiB_1 2 1</div><div><br></div><div>I get this error:</div><div><br>Traceback (most recent call last):<br>  File "testVI.py", line 103, in <module></div><div><p style="margin:0px;line-height:normal">    snes.solve(None,X)<br>  File "PETSc/SNES.pyx", line 520, in petsc4py.PETSc.SNES.solve (src/petsc4py.PETSc.c:169677)<br>petsc4py.PETSc.Error: error code 60<br>[0] SNESSolve() line 3992 in /Users/justin/Software/petsc/src/snes/interface/snes.c<br>[0] SNESSolve_VINEWTONRSLS() line 507 in /Users/justin/Software/petsc/src/snes/impls/vi/rs/virs.c<br>[0] KSPSetOperators() line 544 in /Users/justin/Software/petsc/src/ksp/ksp/interface/itcreate.c<br>[0] PCSetOperators() line 1170 in /Users/justin/Software/petsc/src/ksp/pc/interface/precon.c<br>[0] Nonconforming object sizes<br>[0] Cannot change local size of Amat after use old sizes 2 2 new sizes 3 3</p><p style="margin:0px;line-height:normal"><br></p><p style="margin:0px;line-height:normal">No such error occurred when this was ran in firedrake, though I did notice that -snes_view did indicate that some of the iterations had 2x2 instead of 3x3 matrices. Why is this happening? I thought I should be solving a 3x3 system each time so why would a 2x2 ever arise?</p><p style="margin:0px;line-height:normal"><br></p><p style="margin:0px;line-height:normal">More issues to come :)</p><p style="margin:0px;line-height:normal"><br></p><p style="margin:0px;line-height:normal">Thanks,</p><p style="margin:0px;line-height:normal">Justin</p></div></div>