Barry’s comment resolved my issue. -gideon
> On Aug 11, 2015, at 5:40 PM, Barry Smith <[email protected]> wrote: > > You also have to KEEP the -snes_check_jacobian option > > Barry > >> On Aug 11, 2015, at 1:09 PM, Gideon Simpson <[email protected]> wrote: >> >> I don’t see it listed in -help, but I do get >> >> ./blowup -xmax 50 -nx 1000 -snes_check_jacobian >> Testing hand-coded Jacobian, if the ratio is O(1.e-8), the hand-coded >> Jacobian is probably correct. >> Run with -snes_check_jacobian_view [viewer][:filename][:format] to show >> difference of hand-coded and finite difference Jacobian. >> 63386.1 = ||J - Jfd||//J|| 63386.1 = ||J - Jfd|| >> >> -gideon >> >>> On Aug 11, 2015, at 2:07 PM, Matthew Knepley <[email protected]> wrote: >>> >>> On Tue, Aug 11, 2015 at 1:04 PM, Gideon Simpson <[email protected]> >>> wrote: >>> Macports installation of 3.5.3. >>> >>> Use -help to find the option name. Maybe its -snes_test. >>> >>> Thanks, >>> >>> Matt >>> >>> -gideon >>> >>>> On Aug 11, 2015, at 2:03 PM, Matthew Knepley <[email protected]> wrote: >>>> >>>> On Tue, Aug 11, 2015 at 12:51 PM, Gideon Simpson >>>> <[email protected]> wrote: >>>> #End of PETSc Option Table entries >>>> There is one unused database option. It is: >>>> Option left: name:-snes_check_jacobian_view (no value) >>>> >>>> This is the option for the newest release. What are you using? >>>> >>>> Matt >>>> >>>> -gideon >>>> >>>>> On Aug 11, 2015, at 1:50 PM, Matthew Knepley <[email protected]> wrote: >>>>> >>>>> On Tue, Aug 11, 2015 at 12:49 PM, Gideon Simpson >>>>> <[email protected]> wrote: >>>>> Maybe it’s a quirk of the macports installation of petsc, but nothing >>>>> seems to be getting generated. >>>>> >>>>> Run with -options_left. Is it reading the option? >>>>> >>>>> Matt >>>>> >>>>> -gideon >>>>> >>>>>> On Aug 11, 2015, at 1:40 PM, Jed Brown <[email protected]> wrote: >>>>>> >>>>>> Gideon Simpson <[email protected]> writes: >>>>>> >>>>>>> I’m a bit confused by the following options: >>>>>>> >>>>>>> Run with -snes_check_jacobian_view [viewer][:filename][:format] to show >>>>>>> difference of hand-coded and finite difference Jacobian. >>>>>>> >>>>>>> What flags do I pass it to get some output to diagnose my Jacobian >>>>>>> error? >>>>>> >>>>>> Nothing to display ASCII to the screen. You might use >>>>>> "binary:thematrix" if you want to read it in with Python or MATLAB, for >>>>>> example. >>>>> >>>>> >>>>> >>>>> >>>>> -- >>>>> 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 >>>> >>>> >>>> >>>> >>>> -- >>>> 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 >>> >>> >>> >>> >>> -- >>> 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 >> >
