Hi All, In 'SNESSetJacobian' I could pass 'SNESDefaultComputeJacobian' to let PETSc do the finite differencing for me to calculate Jacobian, i.e.,
SNESSetJacobian(snes, J, J, SNESDefaultComputeJacobian, PETSC_NULL); In TSSetRHSJacobian, this options seems to be invalid, i.e., it requires a function to evaluate the Jacobian, i.e., TSSetRHSJacobian(ts, J, J, FormRHSJacobian, PETSC_NULL); Is it possible that TS could also use PETSc finite differencing Jacobian as SNES does? Best, Ling
