[petsc-users] snes_type aspin error

[Barry Smith]

   Are you providing a shell matrix to the solver? 

 Barry

> On Jun 20, 2016, at 11:32 AM, Xiangdong <epscodes at gmail.com> wrote:
> 
> Hello everyone,
> 
> 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.
> 
> Any hint for me to fix this? Thanks.
> 
> Best,
> Xiangdong
> 
> [0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
> [0]PETSC ERROR: No support for this operation for this object type
> [0]PETSC ERROR: Mat type shell
> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
> [0]PETSC ERROR: Petsc Release Version 3.6.3, unknown
> [0]PETSC ERROR: #1 MatZeroEntries() line 5404 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/mat/interface/matrix.c
> [0]PETSC ERROR: #2 FormJacobianLocal() line 37 in /home/MyCodes/epscprplus/Jacobian.c
> [0]PETSC ERROR: #3 SNESComputeJacobian_DMDA() line 175 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/utils/dmdasnes.c
> [0]PETSC ERROR: #4 SNESComputeJacobian() line 2232 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/interface/snes.c
> [0]PETSC ERROR: #5 SNESNASMComputeFinalJacobian_Private() line 770 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/impls/nasm/nasm.c
> [0]PETSC ERROR: #6 SNESSolve_NASM() line 888 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/impls/nasm/nasm.c
> [0]PETSC ERROR: #7 SNESSolve() line 3906 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/interface/snes.c
> [0]PETSC ERROR: #8 SNESApplyNPC() line 48 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/interface/snespc.c
> [0]PETSC ERROR: #9 SNESSolve_NEWTONLS() line 173 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/impls/ls/ls.c
> [0]PETSC ERROR: #10 SNESSolve() line 3906 in /home/shared/MyLocal_stor/petsc/petsc-dev/src/snes/interface/snes.c
> [0]PETSC ERROR: #11 main() line 794 in /home/MyCodes/epscprplus/twophaseow.c
>