[petsc-users] -snes_mf_operator yields "No support for this operation for this object type" in TS codes?

[Ed Bueler]

Dear PETSc --

I have a TS-using and DMDA-using code in which I want to set a RHSJacobian
which is only approximate.  (The Jacobian uses first-order upwinding MOL
while the RHSFunction uses a flux-limited MOL.)  While it works with the
analytical Jacobian, and -snes_fd, and -snes_fd_color, and -snes_mf, I get
a "No support ..." message for -snes_mf_operator.

It suffices to show the problem with

src/ts/examples/tutorials/ex2.c

(Note ex2.c does not use DMDA so ..._color is not available for ex2.c.)

Error as follows:

~/petsc/src/ts/examples/tutorials[master*]$ ./ex2 -ts_dt 0.001
-ts_final_time 0.001 -snes_monitor
    0 SNES Function norm 1.059316422854e+01
    1 SNES Function norm 1.035505461114e-05
    2 SNES Function norm 5.498223366328e-12
~/petsc/src/ts/examples/tutorials[master*]$ ./ex2 -ts_dt 0.001
-ts_final_time 0.001 -snes_monitor -snes_fd
    0 SNES Function norm 1.059316422854e+01
    1 SNES Function norm 1.208245988550e-05
    2 SNES Function norm 6.022374930788e-12
~/petsc/src/ts/examples/tutorials[master*]$ ./ex2 -ts_dt 0.001
-ts_final_time 0.001 -snes_monitor -snes_mf
    0 SNES Function norm 1.059316422854e+01
    1 SNES Function norm 6.136984336801e-05
    2 SNES Function norm 5.355730806625e-10
~/petsc/src/ts/examples/tutorials[master*]$ ./ex2 -ts_dt 0.001
-ts_final_time 0.001 -snes_monitor -snes_mf_operator
    0 SNES Function norm 1.059316422854e+01
[0]PETSC ERROR: --------------------- Error Message
--------------------------------------------------------------
[0]PETSC ERROR: No support for this operation for this object type
[0]PETSC ERROR: Mat type mffd
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for
trouble shooting.
[0]PETSC ERROR: Petsc Development GIT revision: v3.7.5-3426-g0c7851c  GIT
Date: 2017-04-01 18:40:06 -0600
[0]PETSC ERROR: ./ex2 on a linux-c-dbg named ed-lemur by ed Sun Apr  2
17:02:47 2017
[0]PETSC ERROR: Configure options --download-mpich --with-debugging=1
[0]PETSC ERROR: #1 MatZeroEntries() line 5471 in
/home/ed/petsc/src/mat/interface/matrix.c
[0]PETSC ERROR: #2 TSComputeIJacobian() line 965 in
/home/ed/petsc/src/ts/interface/ts.c
[0]PETSC ERROR: #3 SNESTSFormJacobian_Theta() line 515 in
/home/ed/petsc/src/ts/impls/implicit/theta/theta.c
[0]PETSC ERROR: #4 SNESTSFormJacobian() line 5078 in
/home/ed/petsc/src/ts/interface/ts.c
[0]PETSC ERROR: #5 SNESComputeJacobian() line 2276 in
/home/ed/petsc/src/snes/interface/snes.c
[0]PETSC ERROR: #6 SNESSolve_NEWTONLS() line 222 in
/home/ed/petsc/src/snes/impls/ls/ls.c
[0]PETSC ERROR: #7 SNESSolve() line 3967 in
/home/ed/petsc/src/snes/interface/snes.c
[0]PETSC ERROR: #8 TS_SNESSolve() line 171 in
/home/ed/petsc/src/ts/impls/implicit/theta/theta.c
[0]PETSC ERROR: #9 TSStep_Theta() line 211 in
/home/ed/petsc/src/ts/impls/implicit/theta/theta.c
[0]PETSC ERROR: #10 TSStep() line 3843 in
/home/ed/petsc/src/ts/interface/ts.c
[0]PETSC ERROR: #11 TSSolve() line 4088 in
/home/ed/petsc/src/ts/interface/ts.c
[0]PETSC ERROR: #12 main() line 194 in
/home/ed/petsc/src/ts/examples/tutorials/ex2.c
[0]PETSC ERROR: PETSc Option Table entries:
[0]PETSC ERROR: -snes_mf_operator
[0]PETSC ERROR: -snes_monitor
[0]PETSC ERROR: -ts_dt 0.001
[0]PETSC ERROR: -ts_final_time 0.001
[0]PETSC ERROR: ----------------End of Error Message -------send entire
error message to petsc-maint at mcs.anl.gov----------
application called MPI_Abort(MPI_COMM_WORLD, 56) - process 0
[unset]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 56) - process 0

Is this intended or a bug?  If it is intended, what is the issue?

I am using current master branch (commit
0c7851c55cba8e40da5083f79ba1ff846acd45b2).

Thanks for your help and awesome library!

Ed





-- 
Ed Bueler
Dept of Math and Stat and Geophysical Institute
University of Alaska Fairbanks
Fairbanks, AK 99775-6660
301C Chapman
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170402/a635077e/attachment.html>