nrnste.cpp

Go to the documentation of this file.

#include <../../nrnconf.h>
#include <stdio.h>
#include <InterViews/resource.h>
#include <string.h>
#include <nrnoc2iv.h>
#include <nrniv_mf.h>
#include <classreg.h>
#include <objcmd.h>
#include <nrnste.h>
#include <netcon.h>
 
extern int hoc_return_type_code;
 
static double ste_transition(void* v) {
    StateTransitionEvent* ste = (StateTransitionEvent*) v;
    int src = (int) chkarg(1, 0, ste->nstate() - 1);
    int dest = (int) chkarg(2, 0, ste->nstate() - 1);
    double* var1 = hoc_pgetarg(3);
    double* var2 = hoc_pgetarg(4);
    HocCommand* hc = NULL;
    if (ifarg(5)) {
        Object* obj = NULL;
        if (hoc_is_str_arg(5)) {
            char* stmt = NULL;
            stmt = gargstr(5);
            if (ifarg(6)) {
                obj = *hoc_objgetarg(6);
            }
            hc = new HocCommand(stmt, obj);
        } else {
            obj = *hoc_objgetarg(5);
            hc = new HocCommand(obj);
        }
    }
    ste->transition(src, dest, var1, var2, hc);
    return 1.;
}
 
static double ste_state(void* v) {
    StateTransitionEvent* ste = (StateTransitionEvent*) v;
    hoc_return_type_code = 1;  // integer
    int state = ste->state();
    if (ifarg(1)) {
        ste->state((int) chkarg(1, 0, ste->nstate() - 1));
    }
    return (double) state;
}
 
static Member_func members[] = {"transition", ste_transition, "state", ste_state, 0, 0};
 
static void* ste_cons(Object*) {
    int nstate = (int) chkarg(1, 1, 1e6);
    Point_process* pnt = NULL;
    if (ifarg(2)) {
        Object* obj = *hoc_objgetarg(2);
        pnt = ob2pntproc(obj);
    }
    StateTransitionEvent* ste = new StateTransitionEvent(nstate, pnt);
    return ste;
}
 
static void ste_destruct(void* v) {
    StateTransitionEvent* ste = (StateTransitionEvent*) v;
    delete ste;
}
 
void StateTransitionEvent_reg() {
    class2oc("StateTransitionEvent", ste_cons, ste_destruct, members, NULL, NULL, NULL);
}
 
StateTransitionEvent::StateTransitionEvent(int nstate, Point_process* pnt) {
    nstate_ = nstate;
    states_ = new STEState[nstate_];
    istate_ = 0;
    pnt_ = pnt;
    activated_ = -1;
}
 
StateTransitionEvent::~StateTransitionEvent() {
    deactivate();
    delete[] states_;
}
 
void StateTransitionEvent::deactivate() {
    if (activated_ < 0) {
        return;
    }
    STEState& s = states_[activated_];
    for (int i = 0; i < s.ntrans_; ++i) {
        s.transitions_[i].deactivate();
    }
    activated_ = -1;
}
 
void StateTransitionEvent::activate() {
    if (activated_ >= 0) {
        deactivate();
    }
    STEState& s = states_[istate_];
    for (int i = 0; i < s.ntrans_; ++i) {
        s.transitions_[i].activate();
    }
    activated_ = istate_;
}
 
void StateTransitionEvent::state(int ist) {
    assert(ist >= 0 && ist < nstate_);
    deactivate();
    istate_ = ist;
    activate();
}
 
STEState::STEState() {
    ntrans_ = 0;
    transitions_ = NULL;
}
 
STEState::~STEState() {
    if (transitions_) {
        delete[] transitions_;
    }
}
 
STETransition* STEState::add_transition() {
    ++ntrans_;
    STETransition* st = transitions_;
    transitions_ = new STETransition[ntrans_];
    if (st) {
        int i;
        for (i = 0; i < ntrans_ - 1; ++i) {
            transitions_[i].hc_ = st[i].hc_;
            st[i].hc_ = NULL;
            transitions_[i].ste_ = st[i].ste_;
            st[i].ste_ = NULL;
            transitions_[i].stec_ = st[i].stec_;
            st[i].stec_ = NULL;
            transitions_[i].stec_->stet_ = transitions_ + i;
            transitions_[i].var1_ = st[i].var1_;
            transitions_[i].var2_ = st[i].var2_;
            transitions_[i].dest_ = st[i].dest_;
            transitions_[i].var1_is_time_ = st[i].var1_is_time_;
        }
        delete[] st;
    }
    return transitions_ + ntrans_ - 1;
}
 
STETransition::STETransition() {
    hc_ = NULL;
    ste_ = NULL;
    stec_ = NULL;
    var1_is_time_ = false;
}
 
STETransition::~STETransition() {
    if (hc_) {
        delete hc_;
    }
    if (stec_) {
        delete stec_;
    }
}
 
void STETransition::event() {
    ste_->deactivate();
    ste_->istate_ = dest_;
    if (hc_) {
        nrn_hoc_lock();
        hc_->execute();
        nrn_hoc_unlock();
    }
    ste_->activate();
}