TSSetRHSJacobianP#

Sets the function that computes the Jacobian of G w.r.t. the parameters P where U_t = G(U,P,t), as well as the location to store the matrix.

Synopsis#

#include "petscts.h"  
PetscErrorCode TSSetRHSJacobianP(TS ts, Mat Amat, PetscErrorCode (*func)(TS, PetscReal, Vec, Mat, void *), void *ctx)

Logically Collective

Input Parameters#

  • ts - TS context obtained from TSCreate()

  • Amat - JacobianP matrix

  • func - function

  • ctx - [optional] user-defined function context

Calling sequence of func#

PetscErrorCode func(TS ts, PetscReal t, Vec y, Mat A, void *ctx)

  • t - current timestep

  • U - input vector (current ODE solution)

  • A - output matrix

  • ctx - [optional] user-defined function context

Note#

Amat has the same number of rows and the same row parallel layout as u, Amat has the same number of columns and parallel layout as p

See Also#

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

Level#

intermediate

Location#

src/ts/interface/sensitivity/tssen.c

Examples#

src/ts/tutorials/ex16fwd.c
src/ts/tutorials/ex20adj.c
src/ts/tutorials/ex20fwd.c
src/ts/tutorials/ex20opt_p.c
src/ts/tutorials/ex20td.c

Edit on GitLab

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