Finite Difference Jacobian implementation based on CVODEs. More...
#include "fd_jacob.cuh"
Include dependency graph for fd_jacob.cu:
Go to the source code of this file.
Macros | |
| #define | FD_ORD 1 |
| The finite difference order [Default: 1]. More... | |
Functions | |
| __device__ void | eval_jacob (const double t, const double pres, const double *__restrict__ cy, double *__restrict__ jac, const mechanism_memory *__restrict__ d_mem, double *__restrict__ y_temp, double *__restrict__ ewt) |
| Computes a finite difference Jacobian of order FD_ORD of the RHS function dydt at the given pressure and state. More... | |
Detailed Description
Finite Difference Jacobian implementation based on CVODEs.
Definition in file fd_jacob.cu.
Macro Definition Documentation
◆ FD_ORD
| #define FD_ORD 1 |
The finite difference order [Default: 1].
Definition at line 9 of file fd_jacob.cu.
Function Documentation
◆ eval_jacob()
| __device__ void eval_jacob | ( | const double | t, |
| const double | pres, | ||
| const double *__restrict__ | cy, | ||
| double *__restrict__ | jac, | ||
| const mechanism_memory *__restrict__ | d_mem, | ||
| double *__restrict__ | y_temp, | ||
| double *__restrict__ | ewt | ||
| ) |
Computes a finite difference Jacobian of order FD_ORD of the RHS function dydt at the given pressure and state.
- Parameters
-
[in] t the current system time [in] pres the current system pressure [in] cy the system state vector [out] jac the resulting Jacobian [in] d_mem the mechanism_memory object used in computing dydt [in] y_temp a work array for the state vector [in] ewt a storage for the error weights in computing the Jacobian perturbation factor
Definition at line 35 of file fd_jacob.cu.
