solver_interface.c

Go to the documentation of this file.

#include "solver_interface.h"
#include <math.h>

#ifdef GENERATE_DOCS
namespace generic {
#endif

void accelerInt_initialize(int num_threads) {
    initialize_solver(num_threads);
}


void accelerInt_integrate(const int NUM, const double t_start, const double t_end, const double stepsize,
                          double * __restrict__ y_host, const double * __restrict__ var_host)
{
    double t = t_start;
    double step = stepsize < 0 ? t_end - t : stepsize;
    double t_next = fmin(end_time, t + step);
    int numSteps = 0;

    // time integration loop
    while (t + EPS < t_end)
    {
        numSteps++;
        intDriver(NUM, t, t_end, var_host, y_host);
        t = t_next;
        t_next = fmin(t_end, (numSteps + 1) * step);
    }
}


void accelerInt_cleanup(int num_threads) {
    cleanup_solver(num_threads);
}




#ifdef GENERATE_DOCS
}
#endif