TSComputeIJacobian

Evaluates the Jacobian of the DAE

Synopsis

#include "petscts.h"  
PetscErrorCode TSComputeIJacobian(TS ts, PetscReal t, Vec U, Vec Udot, PetscReal shift, Mat A, Mat B, PetscBool imex)

Collective

Input Parameters

• ts - the TS context

• t - current timestep

• U - state vector

• Udot - time derivative of state vector

• shift - shift to apply, see note below

• imex - flag indicates if the method is TSARKIMEX so that the RHSJacobian should be kept separate

Output Parameters

• A - Jacobian matrix

• B - matrix from which the preconditioner is constructed; often the same as A

Notes

If F(t,U,Udot)=0 is the DAE, the required Jacobian is

   dF/dU + shift*dF/dUdot

Most users should not need to explicitly call this routine, as it is used internally within the nonlinear solvers.

See Also

TS: Scalable ODE and DAE Solvers, TS, TSSetIJacobian()

Level

developer

Location

src/ts/interface/ts.c

---

Index of all TS routines
Table of Contents for all manual pages
Index of all manual pages