jacob.cu

Go to the documentation of this file.

#include "header.cuh"

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

__device__
void eval_jacob (const double t, const double mu, const double * __restrict__ y, double * __restrict__ jac, const mechanism_memory * __restrict__ d_mem)
{
    //Note, to reach index [i, j] of the Jacobian, we multiply `i` by NSP, the size of the first dimension of Jacobian and add j, i.e.:
    //jac[i, j] -> jac[i * NSP + j]
    jac[INDEX(0 * NSP + 0)] = 0;
    jac[INDEX(0 * NSP + 1)] = -2 * mu * y[INDEX(0)] * y[INDEX(1)] - 1;
    jac[INDEX(1 * NSP + 0)] = 1;
    jac[INDEX(1 * NSP + 1)] = mu * (1 - y[INDEX(0)] * y[INDEX(0)]);
}

#ifdef GENERATE_DOCS
}
#endif