[petsc-users] snes_type aspin error

[Xiangdong]

Hi Matt,

Here is the output of snesview before snessolve.

SNES Object: 2 MPI processes
  type: newtonls
  SNES has not been set up so information may be incomplete
  maximum iterations=10, maximum function evaluations=10000
  tolerances: relative=0.001, absolute=1e-50, solution=1e-08
  total number of linear solver iterations=0
  total number of function evaluations=0
  norm schedule ALWAYS
  SNESLineSearch Object:   2 MPI processes
    type: basic
    maxstep=1.000000e+08, minlambda=1.000000e-12
    tolerances: relative=1.000000e-08, absolute=1.000000e-15,
lambda=1.000000e-08
    maximum iterations=40
  SNES Object:  (npc_)   2 MPI processes
    type: nasm
    SNES has not been set up so information may be incomplete
      Nonlinear Additive Schwarz: total subdomain blocks = -2
    maximum iterations=10000, maximum function evaluations=10000
    tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
    total number of linear solver iterations=0
    total number of function evaluations=0
    norm schedule ALWAYS
use the results from previous step as initial guess.
the option for switch linesearch method if Newton fails is OFF;
SNES Object: 2 MPI processes
  type: newtonls
  SNES has not been set up so information may be incomplete
  maximum iterations=10, maximum function evaluations=10000
  tolerances: relative=0.001, absolute=1e-50, solution=1e-08
  total number of linear solver iterations=0
  total number of function evaluations=0
  norm schedule ALWAYS
  SNESLineSearch Object:   2 MPI processes
    type: basic
    maxstep=1.000000e+08, minlambda=1.000000e-12
    tolerances: relative=1.000000e-08, absolute=1.000000e-15,
lambda=1.000000e-08
    maximum iterations=40
  SNES Object:  (npc_)   2 MPI processes
    type: nasm
    SNES has not been set up so information may be incomplete
      Nonlinear Additive Schwarz: total subdomain blocks = -2
    maximum iterations=10000, maximum function evaluations=10000
    tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
    total number of linear solver iterations=0
    total number of function evaluations=0
    norm schedule ALWAYS

Thanks.

Xiangdong

On Mon, Jun 20, 2016 at 12:54 PM, Matthew Knepley <knepley at gmail.com> wrote:

> Can you give us the output of -snes_view? It would help to see how you are
> setting everything up, or maybe
> run SNESView() before SNESSolve()?
>
>    Matt
>
> On Mon, Jun 20, 2016 at 9: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
>>
>>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160620/bd556652/attachment.html>