<div dir="ltr">Dear PETSc --<div><br></div><div>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.</div><div><br></div><div>It suffices to show the problem with</div><div><br></div><div>src/ts/examples/tutorials/ex2.c</div><div><br></div><div>(Note ex2.c does not use DMDA so ..._color is not available for ex2.c.)</div><div><br></div><div>Error as follows:<br></div><div><br></div><div><div>~/petsc/src/ts/examples/tutorials[master*]$ ./ex2 -ts_dt 0.001 -ts_final_time 0.001 -snes_monitor </div><div> 0 SNES Function norm 1.059316422854e+01 </div><div> 1 SNES Function norm 1.035505461114e-05 </div><div> 2 SNES Function norm 5.498223366328e-12 </div><div>~/petsc/src/ts/examples/tutorials[master*]$ ./ex2 -ts_dt 0.001 -ts_final_time 0.001 -snes_monitor -snes_fd</div><div> 0 SNES Function norm 1.059316422854e+01 </div><div> 1 SNES Function norm 1.208245988550e-05 </div><div> 2 SNES Function norm 6.022374930788e-12 </div><div>~/petsc/src/ts/examples/tutorials[master*]$ ./ex2 -ts_dt 0.001 -ts_final_time 0.001 -snes_monitor -snes_mf</div><div> 0 SNES Function norm 1.059316422854e+01 </div><div> 1 SNES Function norm 6.136984336801e-05 </div><div> 2 SNES Function norm 5.355730806625e-10 </div><div>~/petsc/src/ts/examples/tutorials[master*]$ ./ex2 -ts_dt 0.001 -ts_final_time 0.001 -snes_monitor -snes_mf_operator</div><div> 0 SNES Function norm 1.059316422854e+01 </div><div>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div><div>[0]PETSC ERROR: No support for this operation for this object type</div><div>[0]PETSC ERROR: Mat type mffd</div><div>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</div><div>[0]PETSC ERROR: Petsc Development GIT revision: v3.7.5-3426-g0c7851c GIT Date: 2017-04-01 18:40:06 -0600</div><div>[0]PETSC ERROR: ./ex2 on a linux-c-dbg named ed-lemur by ed Sun Apr 2 17:02:47 2017</div><div>[0]PETSC ERROR: Configure options --download-mpich --with-debugging=1</div><div>[0]PETSC ERROR: #1 MatZeroEntries() line 5471 in /home/ed/petsc/src/mat/interface/matrix.c</div><div>[0]PETSC ERROR: #2 TSComputeIJacobian() line 965 in /home/ed/petsc/src/ts/interface/ts.c</div><div>[0]PETSC ERROR: #3 SNESTSFormJacobian_Theta() line 515 in /home/ed/petsc/src/ts/impls/implicit/theta/theta.c</div><div>[0]PETSC ERROR: #4 SNESTSFormJacobian() line 5078 in /home/ed/petsc/src/ts/interface/ts.c</div><div>[0]PETSC ERROR: #5 SNESComputeJacobian() line 2276 in /home/ed/petsc/src/snes/interface/snes.c</div><div>[0]PETSC ERROR: #6 SNESSolve_NEWTONLS() line 222 in /home/ed/petsc/src/snes/impls/ls/ls.c</div><div>[0]PETSC ERROR: #7 SNESSolve() line 3967 in /home/ed/petsc/src/snes/interface/snes.c</div><div>[0]PETSC ERROR: #8 TS_SNESSolve() line 171 in /home/ed/petsc/src/ts/impls/implicit/theta/theta.c</div><div>[0]PETSC ERROR: #9 TSStep_Theta() line 211 in /home/ed/petsc/src/ts/impls/implicit/theta/theta.c</div><div>[0]PETSC ERROR: #10 TSStep() line 3843 in /home/ed/petsc/src/ts/interface/ts.c</div><div>[0]PETSC ERROR: #11 TSSolve() line 4088 in /home/ed/petsc/src/ts/interface/ts.c</div><div>[0]PETSC ERROR: #12 main() line 194 in /home/ed/petsc/src/ts/examples/tutorials/ex2.c</div><div>[0]PETSC ERROR: PETSc Option Table entries:</div><div>[0]PETSC ERROR: -snes_mf_operator</div><div>[0]PETSC ERROR: -snes_monitor</div><div>[0]PETSC ERROR: -ts_dt 0.001</div><div>[0]PETSC ERROR: -ts_final_time 0.001</div><div>[0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov----------</div><div>application called MPI_Abort(MPI_COMM_WORLD, 56) - process 0</div><div>[unset]: aborting job:</div><div>application called MPI_Abort(MPI_COMM_WORLD, 56) - process 0</div></div><div><br></div><div>Is this intended or a bug? If it is intended, what is the issue?</div><div><br></div><div>I am using current master branch (commit 0c7851c55cba8e40da5083f79ba1ff846acd45b2).</div><div><br></div><div>Thanks for your help and awesome library!</div><div><br></div><div>Ed</div><div><br></div><div><br></div><div><br></div><div><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div>Ed Bueler<br>Dept of Math and Stat and Geophysical Institute<br>University of Alaska Fairbanks<br>Fairbanks, AK 99775-6660<br>301C Chapman<br></div></div></div>
</div></div>