dimplic.cpp

/*
# =============================================================================
# Copyright (c) 2016 - 2021 Blue Brain Project/EPFL
#
# See top-level LICENSE file for details.
# =============================================================================
*/

/** When we use solve methods like euler, newton or kinetic schemes,
 *  the state/current updates function need to call solver methods
 *  defined in coreneuron. This is typically done via function pointers.
 *  But for GPU implementation using OpenACC, we can't pass function
 *  pointers.The temporary "workaround" for this was to generate
 *  switch case to select the proper callback function. This is implemented
 *  using python script that look into translated file and generate
 *  _kinderiv.h which has cases for steer functions defined below.
 *  This allows OpenACC to select gpu implementations at compile
 *  time.
 *  \todo: eulerfun/difun are legacy macros and can be replaced with
 *         actual steer function for euler/derivimplicit methods.
 */

#include "coreneuron/mechanism/mech/cfile/scoplib.h"
#include "coreneuron/mechanism/mech/mod2c_core_thread.hpp"
#include "_kinderiv.h"
namespace coreneuron {
nrn_pragma_omp(declare target)
int derivimplicit_thread(int n, int* slist, int* dlist, DIFUN fun, _threadargsproto_) {
    difun(fun);
    return 0;
}

int nrn_derivimplicit_steer(int fun, _threadargsproto_) {
    switch (fun) { _NRN_DERIVIMPLICIT_CASES }
    return 0;
}

int nrn_euler_steer(int fun, _threadargsproto_) {
    switch (fun) { _NRN_EULER_CASES }
    return 0;
}

int nrn_newton_steer(int fun, _threadargsproto_) {
    switch (fun) { _NRN_DERIVIMPLICIT_NEWTON_CASES }
    return 0;
}

int nrn_kinetic_steer(int fun, SparseObj* so, double* rhs, _threadargsproto_) {
    switch (fun) { _NRN_KINETIC_CASES }
    return 0;
}
nrn_pragma_omp(end declare target)

}  // namespace coreneuron