void* nrn_ms_bksub(NrnThread* nth) { CTBEGIN nrn_multisplit_bksub(nth); second_order_cur(nth); update(nth); CTADD /* see above comment in nrn_fixed_step_thread */ if (!nrnthread_v_transfer_) { nrn_fixed_step_lastpart(nth); } return (void*)0; }
int Cvode::solvex_thread_part3(double* b, NrnThread* nt){ nrn_multisplit_bksub(nt); //for (i=0; i < v_node_count; ++i) { // printf("%d rhs %d %g t=%g\n", nrnmpi_myid, i, VEC_RHS(i), t); //} if (ncv_->stiff() == 2) { solvemem(nt); }else{ // bug here should multiply by gam } gather_ydot(b, nt->id); //printf("\texit b\n"); //for (i=0; i < neq_; ++i) { printf("\t\t%d %g\n", i, b[i]);} return 0; }