[petsc-users] A number of questions about DMDA with SNES and Quasi-Newton methods

[zakaryah .]

OK - I ran with -snes_monitor -snes_converged_reason
-snes_linesearch_monitor -ksp_monitor -ksp_monitor_true_residual
-ksp_converged_reason -pc_type svd -pc_svd_monitor -snes_type newtonls
-snes_compare_explicit

and here is the full error message, output immediately after

Finite difference Jacobian
Mat Object: 24 MPI processes
  type: mpiaij

[0]PETSC ERROR: --------------------- Error Message
--------------------------------------------------------------

[0]PETSC ERROR: Invalid argument

[0]PETSC ERROR: Matrix not generated from a DMDA

[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for
trouble shooting.

[0]PETSC ERROR: Petsc Release Version 3.7.6, Apr, 24, 2017

[0]PETSC ERROR: ./CalculateOpticalFlow on a arch-linux2-c-opt named
node046.hpc.rockefeller.internal by zfrentz Sat Oct  7 13:44:44 2017

[0]PETSC ERROR: Configure options
--prefix=/ru-auth/local/home/zfrentz/PETSc-3.7.6 --download-fblaslapack
-with-debugging=0

[0]PETSC ERROR: #1 MatView_MPI_DA() line 551 in
/rugpfs/fs0/home/zfrentz/PETSc/build/petsc-3.7.6/src/dm/impls/da/fdda.c

[0]PETSC ERROR: #2 MatView() line 901 in
/rugpfs/fs0/home/zfrentz/PETSc/build/petsc-3.7.6/src/mat/interface/matrix.c

[0]PETSC ERROR: #3 SNESComputeJacobian() line 2371 in
/rugpfs/fs0/home/zfrentz/PETSc/build/petsc-3.7.6/src/snes/interface/snes.c

[0]PETSC ERROR: #4 SNESSolve_NEWTONLS() line 228 in
/rugpfs/fs0/home/zfrentz/PETSc/build/petsc-3.7.6/src/snes/impls/ls/ls.c

[0]PETSC ERROR: #5 SNESSolve() line 4005 in
/rugpfs/fs0/home/zfrentz/PETSc/build/petsc-3.7.6/src/snes/interface/snes.c

[0]PETSC ERROR: #6 solveWarp3D() line 659 in
/ru-auth/local/home/zfrentz/Code/OpticalFlow/working/October6_2017/mshs.c

On Tue, Oct 3, 2017 at 5:37 PM, Jed Brown <jed at jedbrown.org> wrote:

> Always always always send the whole error message.
>
> "zakaryah ." <zakaryah at gmail.com> writes:
>
> > I tried -snes_compare_explicit, and got the following error:
> >
> > [0]PETSC ERROR: Invalid argument
> >
> > [0]PETSC ERROR: Matrix not generated from a DMDA
> >
> > What am I doing wrong?
> >
> > On Tue, Oct 3, 2017 at 10:08 AM, Jed Brown <jed at jedbrown.org> wrote:
> >
> >> Barry Smith <bsmith at mcs.anl.gov> writes:
> >>
> >> >> On Oct 3, 2017, at 5:54 AM, zakaryah . <zakaryah at gmail.com> wrote:
> >> >>
> >> >> I'm still working on this.  I've made some progress, and it looks
> like
> >> the issue is with the KSP, at least for now.  The Jacobian may be
> >> ill-conditioned.  Is it possible to use -snes_test_display during an
> >> intermediate step of the analysis?  I would like to inspect the Jacobian
> >> after several solves have already completed,
> >> >
> >> >    No, our currently code for testing Jacobians is poor quality and
> >> poorly organized. Needs a major refactoring to do things properly. Sorry
> >>
> >> You can use -snes_compare_explicit or -snes_compare_coloring to output
> >> differences on each Newton step.
> >>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20171007/dc9afe67/attachment.html>