<div dir="ltr">Hi Matt,<div><br></div><div>Here is the output of snesview before snessolve.</div><div><br></div><div><div>SNES Object: 2 MPI processes</div><div> type: newtonls</div><div> SNES has not been set up so information may be incomplete</div><div> maximum iterations=10, maximum function evaluations=10000</div><div> tolerances: relative=0.001, absolute=1e-50, solution=1e-08</div><div> total number of linear solver iterations=0</div><div> total number of function evaluations=0</div><div> norm schedule ALWAYS</div><div> SNESLineSearch Object: 2 MPI processes</div><div> type: basic</div><div> maxstep=1.000000e+08, minlambda=1.000000e-12</div><div> tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08</div><div> maximum iterations=40</div><div> SNES Object: (npc_) 2 MPI processes</div><div> type: nasm</div><div> SNES has not been set up so information may be incomplete</div><div> Nonlinear Additive Schwarz: total subdomain blocks = -2</div><div> maximum iterations=10000, maximum function evaluations=10000</div><div> tolerances: relative=1e-08, absolute=1e-50, solution=1e-08</div><div> total number of linear solver iterations=0</div><div> total number of function evaluations=0</div><div> norm schedule ALWAYS</div><div>use the results from previous step as initial guess.</div><div>the option for switch linesearch method if Newton fails is OFF;</div><div>SNES Object: 2 MPI processes</div><div> type: newtonls</div><div> SNES has not been set up so information may be incomplete</div><div> maximum iterations=10, maximum function evaluations=10000</div><div> tolerances: relative=0.001, absolute=1e-50, solution=1e-08</div><div> total number of linear solver iterations=0</div><div> total number of function evaluations=0</div><div> norm schedule ALWAYS</div><div> SNESLineSearch Object: 2 MPI processes</div><div> type: basic</div><div> maxstep=1.000000e+08, minlambda=1.000000e-12</div><div> tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08</div><div> maximum iterations=40</div><div> SNES Object: (npc_) 2 MPI processes</div><div> type: nasm</div><div> SNES has not been set up so information may be incomplete</div><div> Nonlinear Additive Schwarz: total subdomain blocks = -2</div><div> maximum iterations=10000, maximum function evaluations=10000</div><div> tolerances: relative=1e-08, absolute=1e-50, solution=1e-08</div><div> total number of linear solver iterations=0</div><div> total number of function evaluations=0</div><div> norm schedule ALWAYS</div></div><div><br></div><div>Thanks.</div><div><br></div><div>Xiangdong</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 20, 2016 at 12:54 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Can you give us the output of -snes_view? It would help to see how you are setting everything up, or maybe<div>run SNESView() before SNESSolve()?</div><div><br></div><div> Matt</div></div><div class="gmail_extra"><div><div class="h5"><br><div class="gmail_quote">On Mon, Jun 20, 2016 at 9:32 AM, Xiangdong <span dir="ltr"><<a href="mailto:epscodes@gmail.com" target="_blank">epscodes@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello everyone,<div><br></div><div>I am trying different snes_type in petsc. The newtonls and nasm all works well. However, when I switched to aspin, I got the error like "No support for this operation for this object type" "Mat type shell". I am using the dmda 3d in petsc 3.6.3. The full error messages are below.</div><div><br></div><div>Any hint for me to fix this? Thanks.</div><div><br></div><div>Best,</div><div>Xiangdong</div><div><br></div><div><div>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div><div>[0]PETSC ERROR: No support for this operation for this object type</div><div>[0]PETSC ERROR: Mat type shell</div><div>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</div><div>[0]PETSC ERROR: Petsc Release Version 3.6.3, unknown</div><div>[0]PETSC ERROR: #1 MatZeroEntries() line 5404 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/mat/interface/matrix.c</div><div>[0]PETSC ERROR: #2 FormJacobianLocal() line 37 in /home/MyCodes/epscprplus/Jacobian.c</div><div>[0]PETSC ERROR: #3 SNESComputeJacobian_DMDA() line 175 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/utils/dmdasnes.c</div><div>[0]PETSC ERROR: #4 SNESComputeJacobian() line 2232 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/interface/snes.c</div><div>[0]PETSC ERROR: #5 SNESNASMComputeFinalJacobian_Private() line 770 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/impls/nasm/nasm.c</div><div>[0]PETSC ERROR: #6 SNESSolve_NASM() line 888 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/impls/nasm/nasm.c</div><div>[0]PETSC ERROR: #7 SNESSolve() line 3906 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/interface/snes.c</div><div>[0]PETSC ERROR: #8 SNESApplyNPC() line 48 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/interface/snespc.c</div><div>[0]PETSC ERROR: #9 SNESSolve_NEWTONLS() line 173 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/impls/ls/ls.c</div><div>[0]PETSC ERROR: #10 SNESSolve() line 3906 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/interface/snes.c</div><div>[0]PETSC ERROR: #11 main() line 794 in /home/MyCodes/epscprplus/twophaseow.c</div></div><div><br></div></div>
</blockquote></div><br><br clear="all"><div><br></div></div></div><span class="HOEnZb"><font color="#888888">-- <br><div data-smartmail="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</font></span></div>
</blockquote></div><br></div>