[petsc-users] Test Jacobian

[Barry Smith]

On Apr 15, 2014, at 11:09 AM, Peter Brune <prbrune at gmail.com> wrote:

> Perhaps there might be a clue somewhere else in the message:
> 
> [0]PETSC ERROR: SNESTest aborts after Jacobian test!
> 
> All SNESTest does is compute the difference between the finite-difference Jacobian and the user Jacobian (congratulations! yours isn't very different!) and then aborts with the error that you have experienced here.
> 
> It may finally be time to excise SNESTest entirely as we have way better alternatives now.  Anyone mind if I do that?  I'll make sure to update the docs appropriately.

   Peter,

   This has not been done because there are two sets of alternatives and both have problems. We need to merge/fix up the other testing infrastructures before removing the -snes_type test

   Barry

> 
> - Peter
> 
> 
> On Tue, Apr 15, 2014 at 11:01 AM, Que Cat <quecat001 at gmail.com> wrote:
> Hello Petsc-users,
> 
> I tested the hand-coded Jacobian, and received the follow error:
> 
> 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 7.23832e-10 difference 1.39168e-07 (user-defined state)
> Norm of matrix ratio 1.15746e-08 difference 2.2254e-06 (constant state -1.0)
> Norm of matrix ratio 7.23832e-10 difference 1.39168e-07 (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!
> [0]PETSC ERROR: ------------------------------------------------------------------------
> [0]PETSC ERROR: Petsc Release Version 3.4.3, Oct, 15, 2013 
> [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.
> 
> What does "Object is in wrong state!" mean? What are the possible sources of this error? Thank you for your time.
> 
> Que
>