Ejemplo n.º 1
0
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;
}
Ejemplo n.º 2
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;
}