13 int eval_jacob_cvodes(
long int N,
double t, N_Vector y, N_Vector ydot, DlsMat jac,
void* f, N_Vector tmp1, N_Vector tmp2, N_Vector tmp3)
15 double* local_y = NV_DATA_S(y);
16 eval_jacob((
double)t, *(
double*)f, local_y, (
double*)jac->data);
void eval_jacob(const double t, const double pres, const double *cy, double *jac)
Computes a finite difference Jacobian of order FD_ORD of the RHS function dydt at the given pressure ...
int eval_jacob_cvodes(long int N, double t, N_Vector y, N_Vector ydot, DlsMat jac, void *f, N_Vector tmp1, N_Vector tmp2, N_Vector tmp3)
The CVODEs Jacobian interface for a direct dense Jacobian.
A simple wrapper, allowing for use of the analytical jacobian w/ CVODES.