<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 20, 2016 at 1:02 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
   Are you providing a shell matrix to the solver?<br></blockquote><div><br></div><div>No. The matrix (global Jacobian) is setup by MatSetValuesStencil. It works for the standard newtonls as well as nasm preconditioner. I think the aspin will assemble the new Jacobian for me based on the information of old Jacobian.<br><br></div><div>Do I need to provide something special to use the aspin?<br><br></div><div>Thanks.<br><br></div><div>Best,<br></div><div>Xiangdong <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<span class=""><font color="#888888"><br>
 Barry<br>
</font></span><div class=""><div class="h5"><br>
> On Jun 20, 2016, at 11:32 AM, Xiangdong <<a href="mailto:epscodes@gmail.com">epscodes@gmail.com</a>> wrote:<br>
><br>
> Hello everyone,<br>
><br>
> 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.<br>
><br>
> Any hint for me to fix this? Thanks.<br>
><br>
> Best,<br>
> Xiangdong<br>
><br>
> [0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>
> [0]PETSC ERROR: No support for this operation for this object type<br>
> [0]PETSC ERROR: Mat type shell<br>
> [0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" rel="noreferrer" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.<br>
> [0]PETSC ERROR: Petsc Release Version 3.6.3, unknown<br>
> [0]PETSC ERROR: #1 MatZeroEntries() line 5404 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/mat/interface/matrix.c<br>
> [0]PETSC ERROR: #2 FormJacobianLocal() line 37 in /home/MyCodes/epscprplus/Jacobian.c<br>
> [0]PETSC ERROR: #3 SNESComputeJacobian_DMDA() line 175 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/utils/dmdasnes.c<br>
> [0]PETSC ERROR: #4 SNESComputeJacobian() line 2232 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/interface/snes.c<br>
> [0]PETSC ERROR: #5 SNESNASMComputeFinalJacobian_Private() line 770 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/impls/nasm/nasm.c<br>
> [0]PETSC ERROR: #6 SNESSolve_NASM() line 888 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/impls/nasm/nasm.c<br>
> [0]PETSC ERROR: #7 SNESSolve() line 3906 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/interface/snes.c<br>
> [0]PETSC ERROR: #8 SNESApplyNPC() line 48 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/interface/snespc.c<br>
> [0]PETSC ERROR: #9 SNESSolve_NEWTONLS() line 173 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/impls/ls/ls.c<br>
> [0]PETSC ERROR: #10 SNESSolve() line 3906 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/interface/snes.c<br>
> [0]PETSC ERROR: #11 main() line 794 in /home/MyCodes/epscprplus/twophaseow.c<br>
><br>
<br>
</div></div></blockquote></div><br></div></div>