12 #include <cuComplex.h> 27 #define STRIDE (M_MAX + P) 29 #define MAX_STEPS (10000) 31 #define MAX_CONSECUTIVE_ERRORS (5) 86 #define EC_success (0) 88 #define EC_consecutive_steps (1) 90 #define EC_max_steps_exceeded (2) 92 #define EC_h_plus_t_equals_h (3) double * Vm
the Arnoldi basis array
double * k6
the stage 6 results
double * work2
a work array
double * k3
the stage 3 results
double * phiHm
the exponential Krylov subspace array for EXP4
double * k4
the stage 4 results
double * k5
the stage 5 results
int * ipiv
the pivot indicies
Structure containing memory needed for EXP4 algorithm.
double * Hm
The Hessenberg Kyrlov subspace array for EXP4, to take the exponential action on. ...
double * k1
the stage 1 results
double * work3
a work array
cuDoubleComplex * work4
a (complex) work array
double * work1
a work array
double * k2
the stage 2 results
double * k7
the stage 7 results
double * sc
the scaled error coefficients
int * result
an array of integration results for the various threads
cuDoubleComplex * invA
the inverse of the Hessenberg Krylov subspace