12 #include <cuComplex.h> 28 #define STRIDE (M_MAX + P) 30 #define MAX_STEPS (100000) 32 #define MAX_CONSECUTIVE_ERRORS (5) 74 #define EC_success (0) 76 #define EC_consecutive_steps (1) 78 #define EC_max_steps_exceeded (2) 80 #define EC_h_plus_t_equals_h (3) double * Vm
the Arnoldi basis array
double * work3
a work array
double * gy
The difference between RHS function and the Jacobian state vector product.
double * work1
a work array
double * work2
a work array
cuDoubleComplex * invA
the inverse of the Hessenberg Krylov subspace
double * sc
the scaled error coefficients
int * result
an array of integration results for the various threads
double * savedActions
Saved stage results.
cuDoubleComplex * work4
a (complex) work array
double * phiHm
the exponential Krylov subspace array for EXP4
int * ipiv
the pivot indicies
double * Hm
The Hessenberg Kyrlov subspace array for EXP4, to take the exponential action on. ...