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 <[email protected]> wrote: > Always always always send the whole error message. > > "zakaryah ." <[email protected]> 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 <[email protected]> wrote: > > > >> Barry Smith <[email protected]> writes: > >> > >> >> On Oct 3, 2017, at 5:54 AM, zakaryah . <[email protected]> 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. > >> >
