void KCa::dynamics(const FDouble & v, const double & t, tuple<const FDouble&> & cIntTuple, tuple<const double&> & cExtTuple) const { (void)t; (void)cExtTuple; const FDouble & caInt = get<0>(cIntTuple); mInf = scaleLogistic(caInt / (cHalfActive + caInt), (v + mOffsetF) / mScaleF); mTau = mTauAF - scaleLogistic(mTauBF, (v + mTauOffsetF) / mTauVScaleF); }
void ACS::dynamics(const FDouble & v, const double & t) const { (void)t; mInf = pow(scaleLogistic(mCoef * exp((v + mOffsetF) / mScaleF), -(v + mOffset2F)/ mScale2F), 1.0/3.0); if(!(mInf < 1.0)) mInf = 1.0; mTau = mTauAF + scaleLogistic(mTauBF, -(v + mTauOffsetF) / mTauVScaleF); hInf = power<4>( logistic(-(v + hOffsetF) / hScaleF) ); hTau = hTauAF + scaleLogistic(hTauBF, -(v + hTauOffsetF) / hTauVScaleF); }