Are you providing a shell matrix to the solver? Barry
> On Jun 20, 2016, at 11:32 AM, Xiangdong <[email protected]> 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 >
