petsc-3.7.1 2016-05-15
Report Typos and Errors

TSComputeIJacobianConstant

Reuses a time-independent for a semi-implicit DAE or ODE

Synopsis

#include "petscts.h"  
PetscErrorCode TSComputeIJacobianConstant(TS ts,PetscReal t,Vec U,Vec Udot,PetscReal shift,Mat A,Mat B,void *ctx)
Collective on TS Many br

Input Arguments

ts - time stepping context Many br
t - time at which to evaluate Many br
U - state at which to evaluate Many br
Udot - time derivative of state vector Many br
shift - shift to apply Many br
ctx - context Many br

Output Arguments

A - pointer to operator Many br
B - pointer to preconditioning matrix Many br
flg - matrix structure flag Many br

Many br

Notes

This function is intended to be passed to TSSetIJacobian() to evaluate the Jacobian for linear time-independent problems. Many br

It is only appropriate for problems of the form Many br

    M Udot = F(U,t)

where M is constant and F is non-stiff. The user must pass M to TSSetIJacobian(). The current implementation only Many brworks with IMEX time integration methods such as TSROSW and TSARKIMEX, since there is no support for de-constructing Many bran implicit operator of the form Many br

   shift*M + J

where J is the Jacobian of -F(U). Support may be added in a future version of PETSc, but for now, the user must store Many bra copy of M or reassemble it when requested. Many br

See Also

TSSetIFunction(), TSSetIJacobian(), TSComputeIFunctionLinear()

Level:advanced
Location:
src/ts/interface/ts.c
Index of all TS routines
Table of Contents for all manual pages
Index of all manual pages