void NetParEvent::deliver(double tt, NetCvode* nc, NrnThread* nt){
int seq;
if (nrn_use_selfqueue_) { //first handle pending flag=1 self events
nrn_pending_selfqueue(tt, nt);
}
// has to be the last event at this time in order to avoid a race
// condition with HocEvent that may call things such as pc.barrier
// actually allthread HocEvent (cvode.event(tev) and cvode.event(tev,stmt)
// will be executed last after a thread join when nrn_allthread_handle
// is called.
net_cvode_instance->deliver_events(tt, nt);
nt->_stop_stepping = 1;
nt->_t = tt;
#if NRNMPI
if (nrnmpi_numprocs > 0) {
MUTLOCK
seq = ++seqcnt_;
MUTUNLOCK
if (seq == nrn_nthread) {
last_nt_ = nt;
#if BGPDMA
if (use_bgpdma_) {
bgp_dma_receive(nt);
}else{
nrn_spike_exchange(nt);
}
#else
nrn_spike_exchange(nt);
#endif
wx_ += wt_;
ws_ += wt1_;
seqcnt_ = 0;
}
}
void nrn_finitialize(int setv, double v) {
int i;
NrnThread* _nt;
t = 0.;
dt2thread(-1.);
nrn_thread_table_check();
clear_event_queue();
nrn_spike_exchange_init();
#if VECTORIZE
nrn_play_init(); /* Vector.play */
///Play events should be executed before initializing events
for (i=0; i < nrn_nthread; ++i) {
nrn_deliver_events(nrn_threads + i); /* The play events at t=0 */
}
if (setv) {
for (_nt = nrn_threads; _nt < nrn_threads + nrn_nthread; ++_nt) {
for (i=0; i < _nt->end; ++i) {
VEC_V(i) = v;
}
}
}
for (i=0; i < nrn_nthread; ++i) {
nrn_ba(nrn_threads + i, BEFORE_INITIAL);
}
/* the INITIAL blocks are ordered so that mechanisms that write
concentrations are after ions and before mechanisms that read
concentrations.
*/
/* the memblist list in NrnThread is already so ordered */
for (i=0; i < nrn_nthread; ++i) {
NrnThread* nt = nrn_threads + i;
NrnThreadMembList* tml;
for (tml = nt->tml; tml; tml = tml->next) {
mod_f_t s = memb_func[tml->index].initialize;
if (s) {
(*s)(nt, tml->ml, tml->index);
}
}
}
#endif
init_net_events();
for (i = 0; i < nrn_nthread; ++i) {
nrn_ba(nrn_threads + i, AFTER_INITIAL);
}
for (i=0; i < nrn_nthread; ++i) {
nrn_deliver_events(nrn_threads + i); /* The INITIAL sent events at t=0 */
}
for (i=0; i < nrn_nthread; ++i) {
setup_tree_matrix_minimal(nrn_threads + i);
}
for (i=0; i < nrn_nthread; ++i) {
nrn_deliver_events(nrn_threads + i); /* The record events at t=0 */
}
#if NRNMPI
nrn_spike_exchange(nrn_threads);
#endif
}
