jacob.c

Go to the documentation of this file.

#include "header.h"

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

void eval_jacob (const double t, const double mu, const double * __restrict__ y, double * __restrict__ jac)
{
    
    jac[0 * NSP + 0] = 0;
    jac[1 * NSP + 0] = -2 * mu * y[0] * y[1] - 1;
    jac[1 * NSP + 0] = 1;
    jac[1 * NSP + 1] = mu * (1 - y[0] * y[0]);
}

#ifdef GENERATE_DOCS
}
#endif