solver_generic.c

Go to the documentation of this file.

#include "header.h"
#include "solver.h"

#ifdef GENERATE_DOCS
 namespace generic {
#endif

void intDriver (const int NUM, const double t, const double t_end,
                const double *pr_global, double *y_global)
{
    int tid;
    #pragma omp parallel for shared(y_global, pr_global) private(tid)
    for (tid = 0; tid < NUM; ++tid) {

        // local array with initial values
        double y_local[NSP];
        double pr_local = pr_global[tid];

        // load local array with initial values from global array

        for (int i = 0; i < NSP; i++)
        {
            y_local[i] = y_global[tid + i * NUM];
        }

        // call integrator for one time step
        check_error(tid, integrate (t, t_end, pr_local, y_local));

        // update global array with integrated values

        for (int i = 0; i < NSP; i++)
        {
            y_global[tid + i * NUM] = y_local[i];
        }

    } //end tid loop

} // end intDriver

#ifdef GENERATE_DOCS
 }
#endif