[petsc-users] snes test fails

[Jed Brown]

Dharmendar Reddy <dharmareddy84 at gmail.com> writes:

> Hello,
>          I am solving a one dimensional (linear) Poisson equation. I have
> setup a snes problem using DM object. I am using Dirichlet boundary
> conditons at both ends of the one dimensional domain. For the test case i
> use the same value for potential at both ends.
>
> I am using DMplex object for the mesh and dof lay out.
>
> If i run the solver, i get a constant potential profile as expected.
> However, If i run the solver with snes_type test
> It is giving an error. what could i have done wrong ?
> Testing hand-coded Jacobian, if the ratio is
> O(1.e-8), the hand-coded Jacobian is probably correct.
> Run with -snes_test_display to show difference
> of hand-coded and finite difference Jacobian.
> Norm of matrix ratio 1.0913e-10 difference 5.85042e-10 (user-defined state)
> Norm of matrix ratio 0.5 difference 5.36098 (constant state -1.0)
> Norm of matrix ratio 0.666667 difference 10.722 (constant state 1.0)
> [0]PETSC ERROR: --------------------- Error Message
> -----------------------------------
> -
> [0]PETSC ERROR: Object is in wrong state!
> [0]PETSC ERROR: SNESTest aborts after Jacobian test!

This is what '-snes_type test' is supposed to do.  So you're fine and
your Jacobian is fine at the initial state, but not at the constant
value 1.0 or -1.0.  (That's okay if those are non-physical states,
otherwise your Jacobian evaluation is incorrect.)

> [0]PETSC ERROR:
> -----------------------------------------------------------------------
> -
> [0]PETSC ERROR: Petsc Development GIT revision:
> e0030536e6573667cee5340eb367e8213e67d68
> 9  GIT Date: 2013-04-16 21:48:15 -0500
> [0]PETSC ERROR: See docs/changes/index.html for recent updates.
> [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> [0]PETSC ERROR: See docs/index.html for manual pages.
> [0]PETSC ERROR:
> -----------------------------------------------------------------------
> -
> [0]PETSC ERROR: ./PoisTest on a mpi_rScalar_Debug named
> login2.stampede.tacc.utexas.edu
>  by Reddy135 Wed Apr 17 17:22:02 2013
> [0]PETSC ERROR: Libraries linked from
> /home1/00924/Reddy135/LocalApps/petsc/mpi_rScalar
> _Debug/lib
> [0]PETSC ERROR: Configure run at Tue Apr 16 22:30:58 2013
> [0]PETSC ERROR: Configure options   --download-blacs=1 --download-ctetgen=1
> --download-
> metis=1 --download-mumps=1 --download-parmetis=1 --download-scalapack=1
> --download-supe
> rlu_dist=1 --download-triangle=1 --download-umfpack=1
> --with-blas-lapack-dir=/opt/apps/
> intel/13/composer_xe_2013.2.146/mkl/lib/intel64/ --with-debugging=1
> --with-mpi-dir=/opt
> /apps/intel13/mvapich2/1.9/ --with-petsc-arch=mpi_rScalar_Debug
> --with-petsc-dir=/home1
> /00924/Reddy135/LocalApps/petsc PETSC_ARCH=mpi_rScalar_Debug
> [0]PETSC ERROR:
> -----------------------------------------------------------------------
> -
> [0]PETSC ERROR: SNESSolve_Test() line 127 in
> /home1/00924/Reddy135/LocalApps/petsc/src/
> snes/impls/test/snestest.c
> [0]PETSC ERROR: SNESSolve() line 3755 in
> /home1/00924/Reddy135/LocalApps/petsc/src/snes
> /interface/snes.c
>  SNESDivergedReason            0
>  Exiting solve
> -- 
> -----------------------------------------------------
> Dharmendar Reddy Palle
> Graduate Student
> Microelectronics Research center,
> University of Texas at Austin,
> 10100 Burnet Road, Bldg. 160
> MER 2.608F, TX 78758-4445
> e-mail: dharmareddy84 at gmail.com
> Phone: +1-512-350-9082
> United States of America.
> Homepage: https://webspace.utexas.edu/~dpr342