[petsc-users] Issues with the Variational Inequality solver

[Justin Chang]

Hi all,

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.

First issue, when I run the program like this:

python testVI.py psiA_1 psiB_1 2 1

I get this error:

Traceback (most recent call last):
  File "testVI.py", line 103, in <module>

    snes.solve(None,X)
  File "PETSc/SNES.pyx", line 520, in petsc4py.PETSc.SNES.solve
(src/petsc4py.PETSc.c:169677)
petsc4py.PETSc.Error: error code 60
[0] SNESSolve() line 3992 in
/Users/justin/Software/petsc/src/snes/interface/snes.c
[0] SNESSolve_VINEWTONRSLS() line 507 in
/Users/justin/Software/petsc/src/snes/impls/vi/rs/virs.c
[0] KSPSetOperators() line 544 in
/Users/justin/Software/petsc/src/ksp/ksp/interface/itcreate.c
[0] PCSetOperators() line 1170 in
/Users/justin/Software/petsc/src/ksp/pc/interface/precon.c
[0] Nonconforming object sizes
[0] Cannot change local size of Amat after use old sizes 2 2 new sizes 3 3


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?


More issues to come :)


Thanks,

Justin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151201/3cfad916/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: testVI.py
Type: text/x-python-script
Size: 3088 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151201/3cfad916/attachment-0001.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: psiA_1
Type: application/octet-stream
Size: 69448 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151201/3cfad916/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: psiB_1
Type: application/octet-stream
Size: 68176 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151201/3cfad916/attachment-0003.obj>