accelerInt/sparse__multiplier_8cu_source.html

 #include "sparse_multiplier.cuh"

 #ifdef GENERATE_DOCS
 //put this in the van der Pol namespace for documentation
 namespace van_der_pol_cu {
 #endif


 __device__
 void sparse_multiplier(const double * A, const double * Vm, double* w) {
   w[INDEX(0)] =  A[INDEX(0)] * Vm[INDEX(0)] +  A[INDEX(NSP)] * Vm[INDEX(1)];
   w[INDEX(1)] =  A[INDEX(1)] * Vm[INDEX(0)] +  A[INDEX(NSP + 1)] * Vm[INDEX(1)];
 }


 #ifdef GENERATE_DOCS
 }
 #endif
NSP
#define NSP
The IVP system size.
Definition: header.cuh:20

sparse_multiplier.cuh
Header definition for CUDA Jacobian vector multiplier, used in exponential integrators.

van_der_pol_cu
Definition: dydt.cu:14

van_der_pol_cu::sparse_multiplier
__device__ void sparse_multiplier(const double *A, const double *Vm, double *w)
Implements Jacobian \ vector multiplication in sparse (or unrolled) form.
Definition: sparse_multiplier.cu:22

INDEX
#define INDEX(i)
Convenience macro to get the value of a vector at index i, calculated as i * GRID_DIM + T_ID...
Definition: gpu_macros.cuh:24