
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept) :precision binary64 (+ (/ NdChar (+ 1.0 (exp (/ (- (- (- (- Ec Vef) EDonor) mu)) KbT)))) (/ NaChar (+ 1.0 (exp (/ (+ (+ (+ Ev Vef) EAccept) (- mu)) KbT))))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
return (NdChar / (1.0 + exp((-(((Ec - Vef) - EDonor) - mu) / KbT)))) + (NaChar / (1.0 + exp(((((Ev + Vef) + EAccept) + -mu) / KbT))));
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
code = (ndchar / (1.0d0 + exp((-(((ec - vef) - edonor) - mu) / kbt)))) + (nachar / (1.0d0 + exp(((((ev + vef) + eaccept) + -mu) / kbt))))
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
return (NdChar / (1.0 + Math.exp((-(((Ec - Vef) - EDonor) - mu) / KbT)))) + (NaChar / (1.0 + Math.exp(((((Ev + Vef) + EAccept) + -mu) / KbT))));
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): return (NdChar / (1.0 + math.exp((-(((Ec - Vef) - EDonor) - mu) / KbT)))) + (NaChar / (1.0 + math.exp(((((Ev + Vef) + EAccept) + -mu) / KbT))))
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) return Float64(Float64(NdChar / Float64(1.0 + exp(Float64(Float64(-Float64(Float64(Float64(Ec - Vef) - EDonor) - mu)) / KbT)))) + Float64(NaChar / Float64(1.0 + exp(Float64(Float64(Float64(Float64(Ev + Vef) + EAccept) + Float64(-mu)) / KbT))))) end
function tmp = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) tmp = (NdChar / (1.0 + exp((-(((Ec - Vef) - EDonor) - mu) / KbT)))) + (NaChar / (1.0 + exp(((((Ev + Vef) + EAccept) + -mu) / KbT)))); end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := N[(N[(NdChar / N[(1.0 + N[Exp[N[((-N[(N[(N[(Ec - Vef), $MachinePrecision] - EDonor), $MachinePrecision] - mu), $MachinePrecision]) / KbT), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(NaChar / N[(1.0 + N[Exp[N[(N[(N[(N[(Ev + Vef), $MachinePrecision] + EAccept), $MachinePrecision] + (-mu)), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{NdChar}{1 + e^{\frac{-\left(\left(\left(Ec - Vef\right) - EDonor\right) - mu\right)}{KbT}}} + \frac{NaChar}{1 + e^{\frac{\left(\left(Ev + Vef\right) + EAccept\right) + \left(-mu\right)}{KbT}}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 18 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept) :precision binary64 (+ (/ NdChar (+ 1.0 (exp (/ (- (- (- (- Ec Vef) EDonor) mu)) KbT)))) (/ NaChar (+ 1.0 (exp (/ (+ (+ (+ Ev Vef) EAccept) (- mu)) KbT))))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
return (NdChar / (1.0 + exp((-(((Ec - Vef) - EDonor) - mu) / KbT)))) + (NaChar / (1.0 + exp(((((Ev + Vef) + EAccept) + -mu) / KbT))));
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
code = (ndchar / (1.0d0 + exp((-(((ec - vef) - edonor) - mu) / kbt)))) + (nachar / (1.0d0 + exp(((((ev + vef) + eaccept) + -mu) / kbt))))
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
return (NdChar / (1.0 + Math.exp((-(((Ec - Vef) - EDonor) - mu) / KbT)))) + (NaChar / (1.0 + Math.exp(((((Ev + Vef) + EAccept) + -mu) / KbT))));
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): return (NdChar / (1.0 + math.exp((-(((Ec - Vef) - EDonor) - mu) / KbT)))) + (NaChar / (1.0 + math.exp(((((Ev + Vef) + EAccept) + -mu) / KbT))))
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) return Float64(Float64(NdChar / Float64(1.0 + exp(Float64(Float64(-Float64(Float64(Float64(Ec - Vef) - EDonor) - mu)) / KbT)))) + Float64(NaChar / Float64(1.0 + exp(Float64(Float64(Float64(Float64(Ev + Vef) + EAccept) + Float64(-mu)) / KbT))))) end
function tmp = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) tmp = (NdChar / (1.0 + exp((-(((Ec - Vef) - EDonor) - mu) / KbT)))) + (NaChar / (1.0 + exp(((((Ev + Vef) + EAccept) + -mu) / KbT)))); end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := N[(N[(NdChar / N[(1.0 + N[Exp[N[((-N[(N[(N[(Ec - Vef), $MachinePrecision] - EDonor), $MachinePrecision] - mu), $MachinePrecision]) / KbT), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(NaChar / N[(1.0 + N[Exp[N[(N[(N[(N[(Ev + Vef), $MachinePrecision] + EAccept), $MachinePrecision] + (-mu)), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{NdChar}{1 + e^{\frac{-\left(\left(\left(Ec - Vef\right) - EDonor\right) - mu\right)}{KbT}}} + \frac{NaChar}{1 + e^{\frac{\left(\left(Ev + Vef\right) + EAccept\right) + \left(-mu\right)}{KbT}}}
\end{array}
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept) :precision binary64 (+ (/ NdChar (+ (exp (/ (+ mu (+ EDonor (- Vef Ec))) KbT)) 1.0)) (/ NaChar (+ (exp (/ (- (+ (+ Vef Ev) EAccept) mu) KbT)) 1.0))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
return (NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
code = (ndchar / (exp(((mu + (edonor + (vef - ec))) / kbt)) + 1.0d0)) + (nachar / (exp(((((vef + ev) + eaccept) - mu) / kbt)) + 1.0d0))
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
return (NdChar / (Math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (Math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): return (NdChar / (math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0))
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) return Float64(Float64(NdChar / Float64(exp(Float64(Float64(mu + Float64(EDonor + Float64(Vef - Ec))) / KbT)) + 1.0)) + Float64(NaChar / Float64(exp(Float64(Float64(Float64(Float64(Vef + Ev) + EAccept) - mu) / KbT)) + 1.0))) end
function tmp = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) tmp = (NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0)); end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := N[(N[(NdChar / N[(N[Exp[N[(N[(mu + N[(EDonor + N[(Vef - Ec), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + N[(NaChar / N[(N[Exp[N[(N[(N[(N[(Vef + Ev), $MachinePrecision] + EAccept), $MachinePrecision] - mu), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{NdChar}{e^{\frac{mu + \left(EDonor + \left(Vef - Ec\right)\right)}{KbT}} + 1} + \frac{NaChar}{e^{\frac{\left(\left(Vef + Ev\right) + EAccept\right) - mu}{KbT}} + 1}
\end{array}
Initial program 100.0%
Final simplification100.0%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(let* ((t_0 (/ NaChar (+ (exp (/ (- (+ (+ Vef Ev) EAccept) mu) KbT)) 1.0)))
(t_1 (/ NdChar (+ (exp (/ (+ mu (+ EDonor (- Vef Ec))) KbT)) 1.0)))
(t_2 (+ t_1 (* NaChar 0.5)))
(t_3 (+ t_1 t_0)))
(if (<= t_3 -4e+164)
(+
(/ NaChar (+ (exp (/ EAccept KbT)) 1.0))
(/ NdChar (+ (exp (/ Ec (- KbT))) 1.0)))
(if (<= t_3 -1e-65)
t_2
(if (<= t_3 4e-258)
(/ NdChar (+ (exp (/ (+ (+ Vef EDonor) (- mu Ec)) KbT)) 1.0))
(if (<= t_3 1e+50)
(/ NaChar (+ (exp (/ (+ EAccept (+ Ev (- Vef mu))) KbT)) 1.0))
(if (<= t_3 2e+231) t_2 (+ t_0 (* NdChar 0.5)))))))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0);
double t_1 = NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0);
double t_2 = t_1 + (NaChar * 0.5);
double t_3 = t_1 + t_0;
double tmp;
if (t_3 <= -4e+164) {
tmp = (NaChar / (exp((EAccept / KbT)) + 1.0)) + (NdChar / (exp((Ec / -KbT)) + 1.0));
} else if (t_3 <= -1e-65) {
tmp = t_2;
} else if (t_3 <= 4e-258) {
tmp = NdChar / (exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0);
} else if (t_3 <= 1e+50) {
tmp = NaChar / (exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0);
} else if (t_3 <= 2e+231) {
tmp = t_2;
} else {
tmp = t_0 + (NdChar * 0.5);
}
return tmp;
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_0 = nachar / (exp(((((vef + ev) + eaccept) - mu) / kbt)) + 1.0d0)
t_1 = ndchar / (exp(((mu + (edonor + (vef - ec))) / kbt)) + 1.0d0)
t_2 = t_1 + (nachar * 0.5d0)
t_3 = t_1 + t_0
if (t_3 <= (-4d+164)) then
tmp = (nachar / (exp((eaccept / kbt)) + 1.0d0)) + (ndchar / (exp((ec / -kbt)) + 1.0d0))
else if (t_3 <= (-1d-65)) then
tmp = t_2
else if (t_3 <= 4d-258) then
tmp = ndchar / (exp((((vef + edonor) + (mu - ec)) / kbt)) + 1.0d0)
else if (t_3 <= 1d+50) then
tmp = nachar / (exp(((eaccept + (ev + (vef - mu))) / kbt)) + 1.0d0)
else if (t_3 <= 2d+231) then
tmp = t_2
else
tmp = t_0 + (ndchar * 0.5d0)
end if
code = tmp
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = NaChar / (Math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0);
double t_1 = NdChar / (Math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0);
double t_2 = t_1 + (NaChar * 0.5);
double t_3 = t_1 + t_0;
double tmp;
if (t_3 <= -4e+164) {
tmp = (NaChar / (Math.exp((EAccept / KbT)) + 1.0)) + (NdChar / (Math.exp((Ec / -KbT)) + 1.0));
} else if (t_3 <= -1e-65) {
tmp = t_2;
} else if (t_3 <= 4e-258) {
tmp = NdChar / (Math.exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0);
} else if (t_3 <= 1e+50) {
tmp = NaChar / (Math.exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0);
} else if (t_3 <= 2e+231) {
tmp = t_2;
} else {
tmp = t_0 + (NdChar * 0.5);
}
return tmp;
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): t_0 = NaChar / (math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0) t_1 = NdChar / (math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0) t_2 = t_1 + (NaChar * 0.5) t_3 = t_1 + t_0 tmp = 0 if t_3 <= -4e+164: tmp = (NaChar / (math.exp((EAccept / KbT)) + 1.0)) + (NdChar / (math.exp((Ec / -KbT)) + 1.0)) elif t_3 <= -1e-65: tmp = t_2 elif t_3 <= 4e-258: tmp = NdChar / (math.exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0) elif t_3 <= 1e+50: tmp = NaChar / (math.exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0) elif t_3 <= 2e+231: tmp = t_2 else: tmp = t_0 + (NdChar * 0.5) return tmp
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = Float64(NaChar / Float64(exp(Float64(Float64(Float64(Float64(Vef + Ev) + EAccept) - mu) / KbT)) + 1.0)) t_1 = Float64(NdChar / Float64(exp(Float64(Float64(mu + Float64(EDonor + Float64(Vef - Ec))) / KbT)) + 1.0)) t_2 = Float64(t_1 + Float64(NaChar * 0.5)) t_3 = Float64(t_1 + t_0) tmp = 0.0 if (t_3 <= -4e+164) tmp = Float64(Float64(NaChar / Float64(exp(Float64(EAccept / KbT)) + 1.0)) + Float64(NdChar / Float64(exp(Float64(Ec / Float64(-KbT))) + 1.0))); elseif (t_3 <= -1e-65) tmp = t_2; elseif (t_3 <= 4e-258) tmp = Float64(NdChar / Float64(exp(Float64(Float64(Float64(Vef + EDonor) + Float64(mu - Ec)) / KbT)) + 1.0)); elseif (t_3 <= 1e+50) tmp = Float64(NaChar / Float64(exp(Float64(Float64(EAccept + Float64(Ev + Float64(Vef - mu))) / KbT)) + 1.0)); elseif (t_3 <= 2e+231) tmp = t_2; else tmp = Float64(t_0 + Float64(NdChar * 0.5)); end return tmp end
function tmp_2 = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0); t_1 = NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0); t_2 = t_1 + (NaChar * 0.5); t_3 = t_1 + t_0; tmp = 0.0; if (t_3 <= -4e+164) tmp = (NaChar / (exp((EAccept / KbT)) + 1.0)) + (NdChar / (exp((Ec / -KbT)) + 1.0)); elseif (t_3 <= -1e-65) tmp = t_2; elseif (t_3 <= 4e-258) tmp = NdChar / (exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0); elseif (t_3 <= 1e+50) tmp = NaChar / (exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0); elseif (t_3 <= 2e+231) tmp = t_2; else tmp = t_0 + (NdChar * 0.5); end tmp_2 = tmp; end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := Block[{t$95$0 = N[(NaChar / N[(N[Exp[N[(N[(N[(N[(Vef + Ev), $MachinePrecision] + EAccept), $MachinePrecision] - mu), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(NdChar / N[(N[Exp[N[(N[(mu + N[(EDonor + N[(Vef - Ec), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 + N[(NaChar * 0.5), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(t$95$1 + t$95$0), $MachinePrecision]}, If[LessEqual[t$95$3, -4e+164], N[(N[(NaChar / N[(N[Exp[N[(EAccept / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + N[(NdChar / N[(N[Exp[N[(Ec / (-KbT)), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, -1e-65], t$95$2, If[LessEqual[t$95$3, 4e-258], N[(NdChar / N[(N[Exp[N[(N[(N[(Vef + EDonor), $MachinePrecision] + N[(mu - Ec), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 1e+50], N[(NaChar / N[(N[Exp[N[(N[(EAccept + N[(Ev + N[(Vef - mu), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 2e+231], t$95$2, N[(t$95$0 + N[(NdChar * 0.5), $MachinePrecision]), $MachinePrecision]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{NaChar}{e^{\frac{\left(\left(Vef + Ev\right) + EAccept\right) - mu}{KbT}} + 1}\\
t_1 := \frac{NdChar}{e^{\frac{mu + \left(EDonor + \left(Vef - Ec\right)\right)}{KbT}} + 1}\\
t_2 := t\_1 + NaChar \cdot 0.5\\
t_3 := t\_1 + t\_0\\
\mathbf{if}\;t\_3 \leq -4 \cdot 10^{+164}:\\
\;\;\;\;\frac{NaChar}{e^{\frac{EAccept}{KbT}} + 1} + \frac{NdChar}{e^{\frac{Ec}{-KbT}} + 1}\\
\mathbf{elif}\;t\_3 \leq -1 \cdot 10^{-65}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_3 \leq 4 \cdot 10^{-258}:\\
\;\;\;\;\frac{NdChar}{e^{\frac{\left(Vef + EDonor\right) + \left(mu - Ec\right)}{KbT}} + 1}\\
\mathbf{elif}\;t\_3 \leq 10^{+50}:\\
\;\;\;\;\frac{NaChar}{e^{\frac{EAccept + \left(Ev + \left(Vef - mu\right)\right)}{KbT}} + 1}\\
\mathbf{elif}\;t\_3 \leq 2 \cdot 10^{+231}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_0 + NdChar \cdot 0.5\\
\end{array}
\end{array}
if (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -4e164Initial program 100.0%
Taylor expanded in EAccept around inf
lower-/.f6485.7
Applied rewrites85.7%
Taylor expanded in Ec around inf
mul-1-negN/A
lower-neg.f6482.4
Applied rewrites82.4%
if -4e164 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -9.99999999999999923e-66 or 1.0000000000000001e50 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 2.0000000000000001e231Initial program 99.9%
Taylor expanded in KbT around inf
*-commutativeN/A
lower-*.f6469.5
Applied rewrites69.5%
if -9.99999999999999923e-66 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 3.99999999999999982e-258Initial program 100.0%
Taylor expanded in NdChar around inf
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6482.5
Applied rewrites82.5%
if 3.99999999999999982e-258 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 1.0000000000000001e50Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6466.7
Applied rewrites66.7%
if 2.0000000000000001e231 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) Initial program 99.9%
Taylor expanded in KbT around inf
*-commutativeN/A
lower-*.f6477.9
Applied rewrites77.9%
Final simplification75.4%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(let* ((t_0 (/ NdChar (+ (exp (/ (+ mu (+ EDonor (- Vef Ec))) KbT)) 1.0)))
(t_1 (+ t_0 (/ NaChar (+ (exp (/ EAccept KbT)) 1.0))))
(t_2
(+
t_0
(/ NaChar (+ (exp (/ (- (+ (+ Vef Ev) EAccept) mu) KbT)) 1.0)))))
(if (<= t_2 -5e-81)
t_1
(if (<= t_2 -1e-129)
(/ NaChar (+ (exp (/ mu (- KbT))) 1.0))
(if (<= t_2 4e-258)
(/ NdChar (+ (exp (/ (+ (+ Vef EDonor) (- mu Ec)) KbT)) 1.0))
(if (<= t_2 5e-163)
(/ NaChar (+ (exp (/ (+ EAccept (+ Ev (- Vef mu))) KbT)) 1.0))
t_1))))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0);
double t_1 = t_0 + (NaChar / (exp((EAccept / KbT)) + 1.0));
double t_2 = t_0 + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
double tmp;
if (t_2 <= -5e-81) {
tmp = t_1;
} else if (t_2 <= -1e-129) {
tmp = NaChar / (exp((mu / -KbT)) + 1.0);
} else if (t_2 <= 4e-258) {
tmp = NdChar / (exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0);
} else if (t_2 <= 5e-163) {
tmp = NaChar / (exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = ndchar / (exp(((mu + (edonor + (vef - ec))) / kbt)) + 1.0d0)
t_1 = t_0 + (nachar / (exp((eaccept / kbt)) + 1.0d0))
t_2 = t_0 + (nachar / (exp(((((vef + ev) + eaccept) - mu) / kbt)) + 1.0d0))
if (t_2 <= (-5d-81)) then
tmp = t_1
else if (t_2 <= (-1d-129)) then
tmp = nachar / (exp((mu / -kbt)) + 1.0d0)
else if (t_2 <= 4d-258) then
tmp = ndchar / (exp((((vef + edonor) + (mu - ec)) / kbt)) + 1.0d0)
else if (t_2 <= 5d-163) then
tmp = nachar / (exp(((eaccept + (ev + (vef - mu))) / kbt)) + 1.0d0)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = NdChar / (Math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0);
double t_1 = t_0 + (NaChar / (Math.exp((EAccept / KbT)) + 1.0));
double t_2 = t_0 + (NaChar / (Math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
double tmp;
if (t_2 <= -5e-81) {
tmp = t_1;
} else if (t_2 <= -1e-129) {
tmp = NaChar / (Math.exp((mu / -KbT)) + 1.0);
} else if (t_2 <= 4e-258) {
tmp = NdChar / (Math.exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0);
} else if (t_2 <= 5e-163) {
tmp = NaChar / (Math.exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0);
} else {
tmp = t_1;
}
return tmp;
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): t_0 = NdChar / (math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0) t_1 = t_0 + (NaChar / (math.exp((EAccept / KbT)) + 1.0)) t_2 = t_0 + (NaChar / (math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0)) tmp = 0 if t_2 <= -5e-81: tmp = t_1 elif t_2 <= -1e-129: tmp = NaChar / (math.exp((mu / -KbT)) + 1.0) elif t_2 <= 4e-258: tmp = NdChar / (math.exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0) elif t_2 <= 5e-163: tmp = NaChar / (math.exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0) else: tmp = t_1 return tmp
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = Float64(NdChar / Float64(exp(Float64(Float64(mu + Float64(EDonor + Float64(Vef - Ec))) / KbT)) + 1.0)) t_1 = Float64(t_0 + Float64(NaChar / Float64(exp(Float64(EAccept / KbT)) + 1.0))) t_2 = Float64(t_0 + Float64(NaChar / Float64(exp(Float64(Float64(Float64(Float64(Vef + Ev) + EAccept) - mu) / KbT)) + 1.0))) tmp = 0.0 if (t_2 <= -5e-81) tmp = t_1; elseif (t_2 <= -1e-129) tmp = Float64(NaChar / Float64(exp(Float64(mu / Float64(-KbT))) + 1.0)); elseif (t_2 <= 4e-258) tmp = Float64(NdChar / Float64(exp(Float64(Float64(Float64(Vef + EDonor) + Float64(mu - Ec)) / KbT)) + 1.0)); elseif (t_2 <= 5e-163) tmp = Float64(NaChar / Float64(exp(Float64(Float64(EAccept + Float64(Ev + Float64(Vef - mu))) / KbT)) + 1.0)); else tmp = t_1; end return tmp end
function tmp_2 = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0); t_1 = t_0 + (NaChar / (exp((EAccept / KbT)) + 1.0)); t_2 = t_0 + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0)); tmp = 0.0; if (t_2 <= -5e-81) tmp = t_1; elseif (t_2 <= -1e-129) tmp = NaChar / (exp((mu / -KbT)) + 1.0); elseif (t_2 <= 4e-258) tmp = NdChar / (exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0); elseif (t_2 <= 5e-163) tmp = NaChar / (exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0); else tmp = t_1; end tmp_2 = tmp; end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := Block[{t$95$0 = N[(NdChar / N[(N[Exp[N[(N[(mu + N[(EDonor + N[(Vef - Ec), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 + N[(NaChar / N[(N[Exp[N[(EAccept / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$0 + N[(NaChar / N[(N[Exp[N[(N[(N[(N[(Vef + Ev), $MachinePrecision] + EAccept), $MachinePrecision] - mu), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -5e-81], t$95$1, If[LessEqual[t$95$2, -1e-129], N[(NaChar / N[(N[Exp[N[(mu / (-KbT)), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 4e-258], N[(NdChar / N[(N[Exp[N[(N[(N[(Vef + EDonor), $MachinePrecision] + N[(mu - Ec), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 5e-163], N[(NaChar / N[(N[Exp[N[(N[(EAccept + N[(Ev + N[(Vef - mu), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{NdChar}{e^{\frac{mu + \left(EDonor + \left(Vef - Ec\right)\right)}{KbT}} + 1}\\
t_1 := t\_0 + \frac{NaChar}{e^{\frac{EAccept}{KbT}} + 1}\\
t_2 := t\_0 + \frac{NaChar}{e^{\frac{\left(\left(Vef + Ev\right) + EAccept\right) - mu}{KbT}} + 1}\\
\mathbf{if}\;t\_2 \leq -5 \cdot 10^{-81}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq -1 \cdot 10^{-129}:\\
\;\;\;\;\frac{NaChar}{e^{\frac{mu}{-KbT}} + 1}\\
\mathbf{elif}\;t\_2 \leq 4 \cdot 10^{-258}:\\
\;\;\;\;\frac{NdChar}{e^{\frac{\left(Vef + EDonor\right) + \left(mu - Ec\right)}{KbT}} + 1}\\
\mathbf{elif}\;t\_2 \leq 5 \cdot 10^{-163}:\\
\;\;\;\;\frac{NaChar}{e^{\frac{EAccept + \left(Ev + \left(Vef - mu\right)\right)}{KbT}} + 1}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -4.99999999999999981e-81 or 4.99999999999999977e-163 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) Initial program 99.9%
Taylor expanded in EAccept around inf
lower-/.f6476.7
Applied rewrites76.7%
if -4.99999999999999981e-81 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -9.9999999999999993e-130Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f64100.0
Applied rewrites100.0%
Taylor expanded in mu around inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
if -9.9999999999999993e-130 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 3.99999999999999982e-258Initial program 100.0%
Taylor expanded in NdChar around inf
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6486.8
Applied rewrites86.8%
if 3.99999999999999982e-258 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 4.99999999999999977e-163Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6468.8
Applied rewrites68.8%
Final simplification79.8%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(let* ((t_0 (/ NaChar (+ (exp (/ (- (+ (+ Vef Ev) EAccept) mu) KbT)) 1.0)))
(t_1 (+ t_0 (* NdChar 0.5)))
(t_2 (/ NdChar (+ (exp (/ (+ mu (+ EDonor (- Vef Ec))) KbT)) 1.0)))
(t_3 (+ t_2 t_0)))
(if (<= t_3 -5e+40)
t_1
(if (<= t_3 4e-258)
(/ NdChar (+ (exp (/ (+ (+ Vef EDonor) (- mu Ec)) KbT)) 1.0))
(if (<= t_3 1e+50)
(/ NaChar (+ (exp (/ (+ EAccept (+ Ev (- Vef mu))) KbT)) 1.0))
(if (<= t_3 2e+231) (+ t_2 (* NaChar 0.5)) t_1))))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0);
double t_1 = t_0 + (NdChar * 0.5);
double t_2 = NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0);
double t_3 = t_2 + t_0;
double tmp;
if (t_3 <= -5e+40) {
tmp = t_1;
} else if (t_3 <= 4e-258) {
tmp = NdChar / (exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0);
} else if (t_3 <= 1e+50) {
tmp = NaChar / (exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0);
} else if (t_3 <= 2e+231) {
tmp = t_2 + (NaChar * 0.5);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_0 = nachar / (exp(((((vef + ev) + eaccept) - mu) / kbt)) + 1.0d0)
t_1 = t_0 + (ndchar * 0.5d0)
t_2 = ndchar / (exp(((mu + (edonor + (vef - ec))) / kbt)) + 1.0d0)
t_3 = t_2 + t_0
if (t_3 <= (-5d+40)) then
tmp = t_1
else if (t_3 <= 4d-258) then
tmp = ndchar / (exp((((vef + edonor) + (mu - ec)) / kbt)) + 1.0d0)
else if (t_3 <= 1d+50) then
tmp = nachar / (exp(((eaccept + (ev + (vef - mu))) / kbt)) + 1.0d0)
else if (t_3 <= 2d+231) then
tmp = t_2 + (nachar * 0.5d0)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = NaChar / (Math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0);
double t_1 = t_0 + (NdChar * 0.5);
double t_2 = NdChar / (Math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0);
double t_3 = t_2 + t_0;
double tmp;
if (t_3 <= -5e+40) {
tmp = t_1;
} else if (t_3 <= 4e-258) {
tmp = NdChar / (Math.exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0);
} else if (t_3 <= 1e+50) {
tmp = NaChar / (Math.exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0);
} else if (t_3 <= 2e+231) {
tmp = t_2 + (NaChar * 0.5);
} else {
tmp = t_1;
}
return tmp;
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): t_0 = NaChar / (math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0) t_1 = t_0 + (NdChar * 0.5) t_2 = NdChar / (math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0) t_3 = t_2 + t_0 tmp = 0 if t_3 <= -5e+40: tmp = t_1 elif t_3 <= 4e-258: tmp = NdChar / (math.exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0) elif t_3 <= 1e+50: tmp = NaChar / (math.exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0) elif t_3 <= 2e+231: tmp = t_2 + (NaChar * 0.5) else: tmp = t_1 return tmp
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = Float64(NaChar / Float64(exp(Float64(Float64(Float64(Float64(Vef + Ev) + EAccept) - mu) / KbT)) + 1.0)) t_1 = Float64(t_0 + Float64(NdChar * 0.5)) t_2 = Float64(NdChar / Float64(exp(Float64(Float64(mu + Float64(EDonor + Float64(Vef - Ec))) / KbT)) + 1.0)) t_3 = Float64(t_2 + t_0) tmp = 0.0 if (t_3 <= -5e+40) tmp = t_1; elseif (t_3 <= 4e-258) tmp = Float64(NdChar / Float64(exp(Float64(Float64(Float64(Vef + EDonor) + Float64(mu - Ec)) / KbT)) + 1.0)); elseif (t_3 <= 1e+50) tmp = Float64(NaChar / Float64(exp(Float64(Float64(EAccept + Float64(Ev + Float64(Vef - mu))) / KbT)) + 1.0)); elseif (t_3 <= 2e+231) tmp = Float64(t_2 + Float64(NaChar * 0.5)); else tmp = t_1; end return tmp end
function tmp_2 = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0); t_1 = t_0 + (NdChar * 0.5); t_2 = NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0); t_3 = t_2 + t_0; tmp = 0.0; if (t_3 <= -5e+40) tmp = t_1; elseif (t_3 <= 4e-258) tmp = NdChar / (exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0); elseif (t_3 <= 1e+50) tmp = NaChar / (exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0); elseif (t_3 <= 2e+231) tmp = t_2 + (NaChar * 0.5); else tmp = t_1; end tmp_2 = tmp; end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := Block[{t$95$0 = N[(NaChar / N[(N[Exp[N[(N[(N[(N[(Vef + Ev), $MachinePrecision] + EAccept), $MachinePrecision] - mu), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 + N[(NdChar * 0.5), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(NdChar / N[(N[Exp[N[(N[(mu + N[(EDonor + N[(Vef - Ec), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(t$95$2 + t$95$0), $MachinePrecision]}, If[LessEqual[t$95$3, -5e+40], t$95$1, If[LessEqual[t$95$3, 4e-258], N[(NdChar / N[(N[Exp[N[(N[(N[(Vef + EDonor), $MachinePrecision] + N[(mu - Ec), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 1e+50], N[(NaChar / N[(N[Exp[N[(N[(EAccept + N[(Ev + N[(Vef - mu), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 2e+231], N[(t$95$2 + N[(NaChar * 0.5), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{NaChar}{e^{\frac{\left(\left(Vef + Ev\right) + EAccept\right) - mu}{KbT}} + 1}\\
t_1 := t\_0 + NdChar \cdot 0.5\\
t_2 := \frac{NdChar}{e^{\frac{mu + \left(EDonor + \left(Vef - Ec\right)\right)}{KbT}} + 1}\\
t_3 := t\_2 + t\_0\\
\mathbf{if}\;t\_3 \leq -5 \cdot 10^{+40}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_3 \leq 4 \cdot 10^{-258}:\\
\;\;\;\;\frac{NdChar}{e^{\frac{\left(Vef + EDonor\right) + \left(mu - Ec\right)}{KbT}} + 1}\\
\mathbf{elif}\;t\_3 \leq 10^{+50}:\\
\;\;\;\;\frac{NaChar}{e^{\frac{EAccept + \left(Ev + \left(Vef - mu\right)\right)}{KbT}} + 1}\\
\mathbf{elif}\;t\_3 \leq 2 \cdot 10^{+231}:\\
\;\;\;\;t\_2 + NaChar \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -5.00000000000000003e40 or 2.0000000000000001e231 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) Initial program 99.9%
Taylor expanded in KbT around inf
*-commutativeN/A
lower-*.f6471.8
Applied rewrites71.8%
if -5.00000000000000003e40 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 3.99999999999999982e-258Initial program 100.0%
Taylor expanded in NdChar around inf
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6477.5
Applied rewrites77.5%
if 3.99999999999999982e-258 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 1.0000000000000001e50Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6466.7
Applied rewrites66.7%
if 1.0000000000000001e50 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 2.0000000000000001e231Initial program 100.0%
Taylor expanded in KbT around inf
*-commutativeN/A
lower-*.f6473.2
Applied rewrites73.2%
Final simplification73.6%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(let* ((t_0 (+ (* NaChar 0.5) (/ NdChar (+ (exp (/ EDonor KbT)) 1.0))))
(t_1
(+
(/ NdChar (+ (exp (/ (+ mu (+ EDonor (- Vef Ec))) KbT)) 1.0))
(/ NaChar (+ (exp (/ (- (+ (+ Vef Ev) EAccept) mu) KbT)) 1.0))))
(t_2 (/ NaChar (+ (exp (/ mu (- KbT))) 1.0)))
(t_3 (+ (- Vef mu) (+ Ev EAccept))))
(if (<= t_1 -4e-62)
t_0
(if (<= t_1 -1e-282)
t_2
(if (<= t_1 0.0)
(/
NaChar
(-
2.0
(/
(fma -0.5 (/ (* t_3 t_3) KbT) (- (- mu Vef) (+ Ev EAccept)))
KbT)))
(if (<= t_1 2e-114) t_2 t_0))))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = (NaChar * 0.5) + (NdChar / (exp((EDonor / KbT)) + 1.0));
double t_1 = (NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
double t_2 = NaChar / (exp((mu / -KbT)) + 1.0);
double t_3 = (Vef - mu) + (Ev + EAccept);
double tmp;
if (t_1 <= -4e-62) {
tmp = t_0;
} else if (t_1 <= -1e-282) {
tmp = t_2;
} else if (t_1 <= 0.0) {
tmp = NaChar / (2.0 - (fma(-0.5, ((t_3 * t_3) / KbT), ((mu - Vef) - (Ev + EAccept))) / KbT));
} else if (t_1 <= 2e-114) {
tmp = t_2;
} else {
tmp = t_0;
}
return tmp;
}
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = Float64(Float64(NaChar * 0.5) + Float64(NdChar / Float64(exp(Float64(EDonor / KbT)) + 1.0))) t_1 = Float64(Float64(NdChar / Float64(exp(Float64(Float64(mu + Float64(EDonor + Float64(Vef - Ec))) / KbT)) + 1.0)) + Float64(NaChar / Float64(exp(Float64(Float64(Float64(Float64(Vef + Ev) + EAccept) - mu) / KbT)) + 1.0))) t_2 = Float64(NaChar / Float64(exp(Float64(mu / Float64(-KbT))) + 1.0)) t_3 = Float64(Float64(Vef - mu) + Float64(Ev + EAccept)) tmp = 0.0 if (t_1 <= -4e-62) tmp = t_0; elseif (t_1 <= -1e-282) tmp = t_2; elseif (t_1 <= 0.0) tmp = Float64(NaChar / Float64(2.0 - Float64(fma(-0.5, Float64(Float64(t_3 * t_3) / KbT), Float64(Float64(mu - Vef) - Float64(Ev + EAccept))) / KbT))); elseif (t_1 <= 2e-114) tmp = t_2; else tmp = t_0; end return tmp end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := Block[{t$95$0 = N[(N[(NaChar * 0.5), $MachinePrecision] + N[(NdChar / N[(N[Exp[N[(EDonor / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(NdChar / N[(N[Exp[N[(N[(mu + N[(EDonor + N[(Vef - Ec), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + N[(NaChar / N[(N[Exp[N[(N[(N[(N[(Vef + Ev), $MachinePrecision] + EAccept), $MachinePrecision] - mu), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(NaChar / N[(N[Exp[N[(mu / (-KbT)), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(Vef - mu), $MachinePrecision] + N[(Ev + EAccept), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -4e-62], t$95$0, If[LessEqual[t$95$1, -1e-282], t$95$2, If[LessEqual[t$95$1, 0.0], N[(NaChar / N[(2.0 - N[(N[(-0.5 * N[(N[(t$95$3 * t$95$3), $MachinePrecision] / KbT), $MachinePrecision] + N[(N[(mu - Vef), $MachinePrecision] - N[(Ev + EAccept), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 2e-114], t$95$2, t$95$0]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := NaChar \cdot 0.5 + \frac{NdChar}{e^{\frac{EDonor}{KbT}} + 1}\\
t_1 := \frac{NdChar}{e^{\frac{mu + \left(EDonor + \left(Vef - Ec\right)\right)}{KbT}} + 1} + \frac{NaChar}{e^{\frac{\left(\left(Vef + Ev\right) + EAccept\right) - mu}{KbT}} + 1}\\
t_2 := \frac{NaChar}{e^{\frac{mu}{-KbT}} + 1}\\
t_3 := \left(Vef - mu\right) + \left(Ev + EAccept\right)\\
\mathbf{if}\;t\_1 \leq -4 \cdot 10^{-62}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{-282}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;\frac{NaChar}{2 - \frac{\mathsf{fma}\left(-0.5, \frac{t\_3 \cdot t\_3}{KbT}, \left(mu - Vef\right) - \left(Ev + EAccept\right)\right)}{KbT}}\\
\mathbf{elif}\;t\_1 \leq 2 \cdot 10^{-114}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -4.0000000000000002e-62 or 2.0000000000000001e-114 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) Initial program 99.9%
Taylor expanded in KbT around inf
*-commutativeN/A
lower-*.f6462.6
Applied rewrites62.6%
Taylor expanded in EDonor around inf
lower-/.f6448.6
Applied rewrites48.6%
if -4.0000000000000002e-62 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -1e-282 or 0.0 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 2.0000000000000001e-114Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6453.9
Applied rewrites53.9%
Taylor expanded in mu around inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6441.9
Applied rewrites41.9%
if -1e-282 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 0.0Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6497.9
Applied rewrites97.9%
Taylor expanded in KbT around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
Applied rewrites90.6%
Final simplification54.0%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(let* ((t_0 (* 0.5 (+ NdChar NaChar)))
(t_1 (+ (- Vef mu) (+ Ev EAccept)))
(t_2
(+
(/ NdChar (+ (exp (/ (+ mu (+ EDonor (- Vef Ec))) KbT)) 1.0))
(/ NaChar (+ (exp (/ (- (+ (+ Vef Ev) EAccept) mu) KbT)) 1.0)))))
(if (<= t_2 -1e-65)
t_0
(if (<= t_2 -2e-285)
(* NdChar 0.5)
(if (<= t_2 5e-240)
(/
NaChar
(-
2.0
(/
(fma -0.5 (/ (* t_1 t_1) KbT) (- (- mu Vef) (+ Ev EAccept)))
KbT)))
t_0)))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = 0.5 * (NdChar + NaChar);
double t_1 = (Vef - mu) + (Ev + EAccept);
double t_2 = (NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
double tmp;
if (t_2 <= -1e-65) {
tmp = t_0;
} else if (t_2 <= -2e-285) {
tmp = NdChar * 0.5;
} else if (t_2 <= 5e-240) {
tmp = NaChar / (2.0 - (fma(-0.5, ((t_1 * t_1) / KbT), ((mu - Vef) - (Ev + EAccept))) / KbT));
} else {
tmp = t_0;
}
return tmp;
}
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = Float64(0.5 * Float64(NdChar + NaChar)) t_1 = Float64(Float64(Vef - mu) + Float64(Ev + EAccept)) t_2 = Float64(Float64(NdChar / Float64(exp(Float64(Float64(mu + Float64(EDonor + Float64(Vef - Ec))) / KbT)) + 1.0)) + Float64(NaChar / Float64(exp(Float64(Float64(Float64(Float64(Vef + Ev) + EAccept) - mu) / KbT)) + 1.0))) tmp = 0.0 if (t_2 <= -1e-65) tmp = t_0; elseif (t_2 <= -2e-285) tmp = Float64(NdChar * 0.5); elseif (t_2 <= 5e-240) tmp = Float64(NaChar / Float64(2.0 - Float64(fma(-0.5, Float64(Float64(t_1 * t_1) / KbT), Float64(Float64(mu - Vef) - Float64(Ev + EAccept))) / KbT))); else tmp = t_0; end return tmp end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := Block[{t$95$0 = N[(0.5 * N[(NdChar + NaChar), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(Vef - mu), $MachinePrecision] + N[(Ev + EAccept), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(NdChar / N[(N[Exp[N[(N[(mu + N[(EDonor + N[(Vef - Ec), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + N[(NaChar / N[(N[Exp[N[(N[(N[(N[(Vef + Ev), $MachinePrecision] + EAccept), $MachinePrecision] - mu), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -1e-65], t$95$0, If[LessEqual[t$95$2, -2e-285], N[(NdChar * 0.5), $MachinePrecision], If[LessEqual[t$95$2, 5e-240], N[(NaChar / N[(2.0 - N[(N[(-0.5 * N[(N[(t$95$1 * t$95$1), $MachinePrecision] / KbT), $MachinePrecision] + N[(N[(mu - Vef), $MachinePrecision] - N[(Ev + EAccept), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.5 \cdot \left(NdChar + NaChar\right)\\
t_1 := \left(Vef - mu\right) + \left(Ev + EAccept\right)\\
t_2 := \frac{NdChar}{e^{\frac{mu + \left(EDonor + \left(Vef - Ec\right)\right)}{KbT}} + 1} + \frac{NaChar}{e^{\frac{\left(\left(Vef + Ev\right) + EAccept\right) - mu}{KbT}} + 1}\\
\mathbf{if}\;t\_2 \leq -1 \cdot 10^{-65}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_2 \leq -2 \cdot 10^{-285}:\\
\;\;\;\;NdChar \cdot 0.5\\
\mathbf{elif}\;t\_2 \leq 5 \cdot 10^{-240}:\\
\;\;\;\;\frac{NaChar}{2 - \frac{\mathsf{fma}\left(-0.5, \frac{t\_1 \cdot t\_1}{KbT}, \left(mu - Vef\right) - \left(Ev + EAccept\right)\right)}{KbT}}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -9.99999999999999923e-66 or 5.0000000000000004e-240 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) Initial program 99.9%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6431.0
Applied rewrites31.0%
if -9.99999999999999923e-66 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -2.00000000000000015e-285Initial program 100.0%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6410.2
Applied rewrites10.2%
Taylor expanded in NaChar around 0
*-commutativeN/A
lower-*.f6422.6
Applied rewrites22.6%
if -2.00000000000000015e-285 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 5.0000000000000004e-240Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6496.1
Applied rewrites96.1%
Taylor expanded in KbT around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
Applied rewrites85.6%
Final simplification39.9%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(let* ((t_0 (* 0.5 (+ NdChar NaChar)))
(t_1
(+
(/ NdChar (+ (exp (/ (+ mu (+ EDonor (- Vef Ec))) KbT)) 1.0))
(/ NaChar (+ (exp (/ (- (+ (+ Vef Ev) EAccept) mu) KbT)) 1.0)))))
(if (<= t_1 -1e-65)
t_0
(if (<= t_1 -2e-277)
(* NdChar 0.5)
(if (<= t_1 4e-141)
(/ NaChar (- 2.0 (/ (/ (* -0.5 (* Vef Vef)) KbT) KbT)))
t_0)))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = 0.5 * (NdChar + NaChar);
double t_1 = (NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
double tmp;
if (t_1 <= -1e-65) {
tmp = t_0;
} else if (t_1 <= -2e-277) {
tmp = NdChar * 0.5;
} else if (t_1 <= 4e-141) {
tmp = NaChar / (2.0 - (((-0.5 * (Vef * Vef)) / KbT) / KbT));
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = 0.5d0 * (ndchar + nachar)
t_1 = (ndchar / (exp(((mu + (edonor + (vef - ec))) / kbt)) + 1.0d0)) + (nachar / (exp(((((vef + ev) + eaccept) - mu) / kbt)) + 1.0d0))
if (t_1 <= (-1d-65)) then
tmp = t_0
else if (t_1 <= (-2d-277)) then
tmp = ndchar * 0.5d0
else if (t_1 <= 4d-141) then
tmp = nachar / (2.0d0 - ((((-0.5d0) * (vef * vef)) / kbt) / kbt))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = 0.5 * (NdChar + NaChar);
double t_1 = (NdChar / (Math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (Math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
double tmp;
if (t_1 <= -1e-65) {
tmp = t_0;
} else if (t_1 <= -2e-277) {
tmp = NdChar * 0.5;
} else if (t_1 <= 4e-141) {
tmp = NaChar / (2.0 - (((-0.5 * (Vef * Vef)) / KbT) / KbT));
} else {
tmp = t_0;
}
return tmp;
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): t_0 = 0.5 * (NdChar + NaChar) t_1 = (NdChar / (math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0)) tmp = 0 if t_1 <= -1e-65: tmp = t_0 elif t_1 <= -2e-277: tmp = NdChar * 0.5 elif t_1 <= 4e-141: tmp = NaChar / (2.0 - (((-0.5 * (Vef * Vef)) / KbT) / KbT)) else: tmp = t_0 return tmp
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = Float64(0.5 * Float64(NdChar + NaChar)) t_1 = Float64(Float64(NdChar / Float64(exp(Float64(Float64(mu + Float64(EDonor + Float64(Vef - Ec))) / KbT)) + 1.0)) + Float64(NaChar / Float64(exp(Float64(Float64(Float64(Float64(Vef + Ev) + EAccept) - mu) / KbT)) + 1.0))) tmp = 0.0 if (t_1 <= -1e-65) tmp = t_0; elseif (t_1 <= -2e-277) tmp = Float64(NdChar * 0.5); elseif (t_1 <= 4e-141) tmp = Float64(NaChar / Float64(2.0 - Float64(Float64(Float64(-0.5 * Float64(Vef * Vef)) / KbT) / KbT))); else tmp = t_0; end return tmp end
function tmp_2 = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = 0.5 * (NdChar + NaChar); t_1 = (NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0)); tmp = 0.0; if (t_1 <= -1e-65) tmp = t_0; elseif (t_1 <= -2e-277) tmp = NdChar * 0.5; elseif (t_1 <= 4e-141) tmp = NaChar / (2.0 - (((-0.5 * (Vef * Vef)) / KbT) / KbT)); else tmp = t_0; end tmp_2 = tmp; end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := Block[{t$95$0 = N[(0.5 * N[(NdChar + NaChar), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(NdChar / N[(N[Exp[N[(N[(mu + N[(EDonor + N[(Vef - Ec), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + N[(NaChar / N[(N[Exp[N[(N[(N[(N[(Vef + Ev), $MachinePrecision] + EAccept), $MachinePrecision] - mu), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -1e-65], t$95$0, If[LessEqual[t$95$1, -2e-277], N[(NdChar * 0.5), $MachinePrecision], If[LessEqual[t$95$1, 4e-141], N[(NaChar / N[(2.0 - N[(N[(N[(-0.5 * N[(Vef * Vef), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision] / KbT), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.5 \cdot \left(NdChar + NaChar\right)\\
t_1 := \frac{NdChar}{e^{\frac{mu + \left(EDonor + \left(Vef - Ec\right)\right)}{KbT}} + 1} + \frac{NaChar}{e^{\frac{\left(\left(Vef + Ev\right) + EAccept\right) - mu}{KbT}} + 1}\\
\mathbf{if}\;t\_1 \leq -1 \cdot 10^{-65}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-277}:\\
\;\;\;\;NdChar \cdot 0.5\\
\mathbf{elif}\;t\_1 \leq 4 \cdot 10^{-141}:\\
\;\;\;\;\frac{NaChar}{2 - \frac{\frac{-0.5 \cdot \left(Vef \cdot Vef\right)}{KbT}}{KbT}}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -9.99999999999999923e-66 or 4.0000000000000002e-141 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) Initial program 99.9%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6432.2
Applied rewrites32.2%
if -9.99999999999999923e-66 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -1.99999999999999994e-277Initial program 100.0%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6410.5
Applied rewrites10.5%
Taylor expanded in NaChar around 0
*-commutativeN/A
lower-*.f6423.2
Applied rewrites23.2%
if -1.99999999999999994e-277 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 4.0000000000000002e-141Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6489.1
Applied rewrites89.1%
Taylor expanded in KbT around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
Applied rewrites67.2%
Taylor expanded in Vef around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6445.5
Applied rewrites45.5%
Final simplification34.3%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(let* ((t_0 (* 0.5 (+ NdChar NaChar)))
(t_1
(+
(/ NdChar (+ (exp (/ (+ mu (+ EDonor (- Vef Ec))) KbT)) 1.0))
(/ NaChar (+ (exp (/ (- (+ (+ Vef Ev) EAccept) mu) KbT)) 1.0)))))
(if (<= t_1 -1e-65)
t_0
(if (<= t_1 -2e-277)
(* NdChar 0.5)
(if (<= t_1 2e-140)
(/ NaChar (- 2.0 (/ (/ (* -0.5 (* EAccept EAccept)) KbT) KbT)))
t_0)))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = 0.5 * (NdChar + NaChar);
double t_1 = (NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
double tmp;
if (t_1 <= -1e-65) {
tmp = t_0;
} else if (t_1 <= -2e-277) {
tmp = NdChar * 0.5;
} else if (t_1 <= 2e-140) {
tmp = NaChar / (2.0 - (((-0.5 * (EAccept * EAccept)) / KbT) / KbT));
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = 0.5d0 * (ndchar + nachar)
t_1 = (ndchar / (exp(((mu + (edonor + (vef - ec))) / kbt)) + 1.0d0)) + (nachar / (exp(((((vef + ev) + eaccept) - mu) / kbt)) + 1.0d0))
if (t_1 <= (-1d-65)) then
tmp = t_0
else if (t_1 <= (-2d-277)) then
tmp = ndchar * 0.5d0
else if (t_1 <= 2d-140) then
tmp = nachar / (2.0d0 - ((((-0.5d0) * (eaccept * eaccept)) / kbt) / kbt))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = 0.5 * (NdChar + NaChar);
double t_1 = (NdChar / (Math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (Math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
double tmp;
if (t_1 <= -1e-65) {
tmp = t_0;
} else if (t_1 <= -2e-277) {
tmp = NdChar * 0.5;
} else if (t_1 <= 2e-140) {
tmp = NaChar / (2.0 - (((-0.5 * (EAccept * EAccept)) / KbT) / KbT));
} else {
tmp = t_0;
}
return tmp;
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): t_0 = 0.5 * (NdChar + NaChar) t_1 = (NdChar / (math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0)) tmp = 0 if t_1 <= -1e-65: tmp = t_0 elif t_1 <= -2e-277: tmp = NdChar * 0.5 elif t_1 <= 2e-140: tmp = NaChar / (2.0 - (((-0.5 * (EAccept * EAccept)) / KbT) / KbT)) else: tmp = t_0 return tmp
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = Float64(0.5 * Float64(NdChar + NaChar)) t_1 = Float64(Float64(NdChar / Float64(exp(Float64(Float64(mu + Float64(EDonor + Float64(Vef - Ec))) / KbT)) + 1.0)) + Float64(NaChar / Float64(exp(Float64(Float64(Float64(Float64(Vef + Ev) + EAccept) - mu) / KbT)) + 1.0))) tmp = 0.0 if (t_1 <= -1e-65) tmp = t_0; elseif (t_1 <= -2e-277) tmp = Float64(NdChar * 0.5); elseif (t_1 <= 2e-140) tmp = Float64(NaChar / Float64(2.0 - Float64(Float64(Float64(-0.5 * Float64(EAccept * EAccept)) / KbT) / KbT))); else tmp = t_0; end return tmp end
function tmp_2 = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = 0.5 * (NdChar + NaChar); t_1 = (NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0)); tmp = 0.0; if (t_1 <= -1e-65) tmp = t_0; elseif (t_1 <= -2e-277) tmp = NdChar * 0.5; elseif (t_1 <= 2e-140) tmp = NaChar / (2.0 - (((-0.5 * (EAccept * EAccept)) / KbT) / KbT)); else tmp = t_0; end tmp_2 = tmp; end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := Block[{t$95$0 = N[(0.5 * N[(NdChar + NaChar), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(NdChar / N[(N[Exp[N[(N[(mu + N[(EDonor + N[(Vef - Ec), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + N[(NaChar / N[(N[Exp[N[(N[(N[(N[(Vef + Ev), $MachinePrecision] + EAccept), $MachinePrecision] - mu), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -1e-65], t$95$0, If[LessEqual[t$95$1, -2e-277], N[(NdChar * 0.5), $MachinePrecision], If[LessEqual[t$95$1, 2e-140], N[(NaChar / N[(2.0 - N[(N[(N[(-0.5 * N[(EAccept * EAccept), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision] / KbT), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.5 \cdot \left(NdChar + NaChar\right)\\
t_1 := \frac{NdChar}{e^{\frac{mu + \left(EDonor + \left(Vef - Ec\right)\right)}{KbT}} + 1} + \frac{NaChar}{e^{\frac{\left(\left(Vef + Ev\right) + EAccept\right) - mu}{KbT}} + 1}\\
\mathbf{if}\;t\_1 \leq -1 \cdot 10^{-65}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-277}:\\
\;\;\;\;NdChar \cdot 0.5\\
\mathbf{elif}\;t\_1 \leq 2 \cdot 10^{-140}:\\
\;\;\;\;\frac{NaChar}{2 - \frac{\frac{-0.5 \cdot \left(EAccept \cdot EAccept\right)}{KbT}}{KbT}}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -9.99999999999999923e-66 or 2e-140 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) Initial program 99.9%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6432.4
Applied rewrites32.4%
if -9.99999999999999923e-66 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -1.99999999999999994e-277Initial program 100.0%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6410.5
Applied rewrites10.5%
Taylor expanded in NaChar around 0
*-commutativeN/A
lower-*.f6423.2
Applied rewrites23.2%
if -1.99999999999999994e-277 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 2e-140Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6487.9
Applied rewrites87.9%
Taylor expanded in KbT around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
Applied rewrites66.2%
Taylor expanded in EAccept around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6433.1
Applied rewrites33.1%
Final simplification31.2%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(let* ((t_0 (* 0.5 (+ NdChar NaChar)))
(t_1
(+
(/ NdChar (+ (exp (/ (+ mu (+ EDonor (- Vef Ec))) KbT)) 1.0))
(/ NaChar (+ (exp (/ (- (+ (+ Vef Ev) EAccept) mu) KbT)) 1.0)))))
(if (<= t_1 -1e-65)
t_0
(if (<= t_1 -2e-285)
(* NdChar 0.5)
(if (<= t_1 5e-295)
(/ (* 2.0 (* NaChar (* KbT KbT))) (* EAccept EAccept))
t_0)))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = 0.5 * (NdChar + NaChar);
double t_1 = (NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
double tmp;
if (t_1 <= -1e-65) {
tmp = t_0;
} else if (t_1 <= -2e-285) {
tmp = NdChar * 0.5;
} else if (t_1 <= 5e-295) {
tmp = (2.0 * (NaChar * (KbT * KbT))) / (EAccept * EAccept);
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = 0.5d0 * (ndchar + nachar)
t_1 = (ndchar / (exp(((mu + (edonor + (vef - ec))) / kbt)) + 1.0d0)) + (nachar / (exp(((((vef + ev) + eaccept) - mu) / kbt)) + 1.0d0))
if (t_1 <= (-1d-65)) then
tmp = t_0
else if (t_1 <= (-2d-285)) then
tmp = ndchar * 0.5d0
else if (t_1 <= 5d-295) then
tmp = (2.0d0 * (nachar * (kbt * kbt))) / (eaccept * eaccept)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = 0.5 * (NdChar + NaChar);
double t_1 = (NdChar / (Math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (Math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
double tmp;
if (t_1 <= -1e-65) {
tmp = t_0;
} else if (t_1 <= -2e-285) {
tmp = NdChar * 0.5;
} else if (t_1 <= 5e-295) {
tmp = (2.0 * (NaChar * (KbT * KbT))) / (EAccept * EAccept);
} else {
tmp = t_0;
}
return tmp;
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): t_0 = 0.5 * (NdChar + NaChar) t_1 = (NdChar / (math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0)) tmp = 0 if t_1 <= -1e-65: tmp = t_0 elif t_1 <= -2e-285: tmp = NdChar * 0.5 elif t_1 <= 5e-295: tmp = (2.0 * (NaChar * (KbT * KbT))) / (EAccept * EAccept) else: tmp = t_0 return tmp
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = Float64(0.5 * Float64(NdChar + NaChar)) t_1 = Float64(Float64(NdChar / Float64(exp(Float64(Float64(mu + Float64(EDonor + Float64(Vef - Ec))) / KbT)) + 1.0)) + Float64(NaChar / Float64(exp(Float64(Float64(Float64(Float64(Vef + Ev) + EAccept) - mu) / KbT)) + 1.0))) tmp = 0.0 if (t_1 <= -1e-65) tmp = t_0; elseif (t_1 <= -2e-285) tmp = Float64(NdChar * 0.5); elseif (t_1 <= 5e-295) tmp = Float64(Float64(2.0 * Float64(NaChar * Float64(KbT * KbT))) / Float64(EAccept * EAccept)); else tmp = t_0; end return tmp end
function tmp_2 = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = 0.5 * (NdChar + NaChar); t_1 = (NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0)); tmp = 0.0; if (t_1 <= -1e-65) tmp = t_0; elseif (t_1 <= -2e-285) tmp = NdChar * 0.5; elseif (t_1 <= 5e-295) tmp = (2.0 * (NaChar * (KbT * KbT))) / (EAccept * EAccept); else tmp = t_0; end tmp_2 = tmp; end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := Block[{t$95$0 = N[(0.5 * N[(NdChar + NaChar), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(NdChar / N[(N[Exp[N[(N[(mu + N[(EDonor + N[(Vef - Ec), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + N[(NaChar / N[(N[Exp[N[(N[(N[(N[(Vef + Ev), $MachinePrecision] + EAccept), $MachinePrecision] - mu), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -1e-65], t$95$0, If[LessEqual[t$95$1, -2e-285], N[(NdChar * 0.5), $MachinePrecision], If[LessEqual[t$95$1, 5e-295], N[(N[(2.0 * N[(NaChar * N[(KbT * KbT), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(EAccept * EAccept), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.5 \cdot \left(NdChar + NaChar\right)\\
t_1 := \frac{NdChar}{e^{\frac{mu + \left(EDonor + \left(Vef - Ec\right)\right)}{KbT}} + 1} + \frac{NaChar}{e^{\frac{\left(\left(Vef + Ev\right) + EAccept\right) - mu}{KbT}} + 1}\\
\mathbf{if}\;t\_1 \leq -1 \cdot 10^{-65}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-285}:\\
\;\;\;\;NdChar \cdot 0.5\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{-295}:\\
\;\;\;\;\frac{2 \cdot \left(NaChar \cdot \left(KbT \cdot KbT\right)\right)}{EAccept \cdot EAccept}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -9.99999999999999923e-66 or 5.00000000000000008e-295 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) Initial program 99.9%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6430.4
Applied rewrites30.4%
if -9.99999999999999923e-66 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -2.00000000000000015e-285Initial program 100.0%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6410.2
Applied rewrites10.2%
Taylor expanded in NaChar around 0
*-commutativeN/A
lower-*.f6422.6
Applied rewrites22.6%
if -2.00000000000000015e-285 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 5.00000000000000008e-295Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f64100.0
Applied rewrites100.0%
Taylor expanded in KbT around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
Applied rewrites90.6%
Taylor expanded in EAccept around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6435.4
Applied rewrites35.4%
Final simplification30.0%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(let* ((t_0 (* 0.5 (+ NdChar NaChar)))
(t_1
(+
(/ NdChar (+ (exp (/ (+ mu (+ EDonor (- Vef Ec))) KbT)) 1.0))
(/ NaChar (+ (exp (/ (- (+ (+ Vef Ev) EAccept) mu) KbT)) 1.0)))))
(if (<= t_1 -1e-65)
t_0
(if (<= t_1 -2e-285)
(* NdChar 0.5)
(if (<= t_1 4e-258) (/ (* -0.25 (* NaChar Ev)) KbT) t_0)))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = 0.5 * (NdChar + NaChar);
double t_1 = (NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
double tmp;
if (t_1 <= -1e-65) {
tmp = t_0;
} else if (t_1 <= -2e-285) {
tmp = NdChar * 0.5;
} else if (t_1 <= 4e-258) {
tmp = (-0.25 * (NaChar * Ev)) / KbT;
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = 0.5d0 * (ndchar + nachar)
t_1 = (ndchar / (exp(((mu + (edonor + (vef - ec))) / kbt)) + 1.0d0)) + (nachar / (exp(((((vef + ev) + eaccept) - mu) / kbt)) + 1.0d0))
if (t_1 <= (-1d-65)) then
tmp = t_0
else if (t_1 <= (-2d-285)) then
tmp = ndchar * 0.5d0
else if (t_1 <= 4d-258) then
tmp = ((-0.25d0) * (nachar * ev)) / kbt
else
tmp = t_0
end if
code = tmp
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = 0.5 * (NdChar + NaChar);
double t_1 = (NdChar / (Math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (Math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
double tmp;
if (t_1 <= -1e-65) {
tmp = t_0;
} else if (t_1 <= -2e-285) {
tmp = NdChar * 0.5;
} else if (t_1 <= 4e-258) {
tmp = (-0.25 * (NaChar * Ev)) / KbT;
} else {
tmp = t_0;
}
return tmp;
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): t_0 = 0.5 * (NdChar + NaChar) t_1 = (NdChar / (math.exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (math.exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0)) tmp = 0 if t_1 <= -1e-65: tmp = t_0 elif t_1 <= -2e-285: tmp = NdChar * 0.5 elif t_1 <= 4e-258: tmp = (-0.25 * (NaChar * Ev)) / KbT else: tmp = t_0 return tmp
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = Float64(0.5 * Float64(NdChar + NaChar)) t_1 = Float64(Float64(NdChar / Float64(exp(Float64(Float64(mu + Float64(EDonor + Float64(Vef - Ec))) / KbT)) + 1.0)) + Float64(NaChar / Float64(exp(Float64(Float64(Float64(Float64(Vef + Ev) + EAccept) - mu) / KbT)) + 1.0))) tmp = 0.0 if (t_1 <= -1e-65) tmp = t_0; elseif (t_1 <= -2e-285) tmp = Float64(NdChar * 0.5); elseif (t_1 <= 4e-258) tmp = Float64(Float64(-0.25 * Float64(NaChar * Ev)) / KbT); else tmp = t_0; end return tmp end
function tmp_2 = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = 0.5 * (NdChar + NaChar); t_1 = (NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0)); tmp = 0.0; if (t_1 <= -1e-65) tmp = t_0; elseif (t_1 <= -2e-285) tmp = NdChar * 0.5; elseif (t_1 <= 4e-258) tmp = (-0.25 * (NaChar * Ev)) / KbT; else tmp = t_0; end tmp_2 = tmp; end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := Block[{t$95$0 = N[(0.5 * N[(NdChar + NaChar), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(NdChar / N[(N[Exp[N[(N[(mu + N[(EDonor + N[(Vef - Ec), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + N[(NaChar / N[(N[Exp[N[(N[(N[(N[(Vef + Ev), $MachinePrecision] + EAccept), $MachinePrecision] - mu), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -1e-65], t$95$0, If[LessEqual[t$95$1, -2e-285], N[(NdChar * 0.5), $MachinePrecision], If[LessEqual[t$95$1, 4e-258], N[(N[(-0.25 * N[(NaChar * Ev), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.5 \cdot \left(NdChar + NaChar\right)\\
t_1 := \frac{NdChar}{e^{\frac{mu + \left(EDonor + \left(Vef - Ec\right)\right)}{KbT}} + 1} + \frac{NaChar}{e^{\frac{\left(\left(Vef + Ev\right) + EAccept\right) - mu}{KbT}} + 1}\\
\mathbf{if}\;t\_1 \leq -1 \cdot 10^{-65}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-285}:\\
\;\;\;\;NdChar \cdot 0.5\\
\mathbf{elif}\;t\_1 \leq 4 \cdot 10^{-258}:\\
\;\;\;\;\frac{-0.25 \cdot \left(NaChar \cdot Ev\right)}{KbT}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -9.99999999999999923e-66 or 3.99999999999999982e-258 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) Initial program 99.9%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6430.6
Applied rewrites30.6%
if -9.99999999999999923e-66 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -2.00000000000000015e-285Initial program 100.0%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6410.2
Applied rewrites10.2%
Taylor expanded in NaChar around 0
*-commutativeN/A
lower-*.f6422.6
Applied rewrites22.6%
if -2.00000000000000015e-285 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 3.99999999999999982e-258Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6497.9
Applied rewrites97.9%
Taylor expanded in KbT around inf
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
associate-*r/N/A
lower-/.f64N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
associate--l+N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f642.3
Applied rewrites2.3%
Taylor expanded in Ev around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6424.5
Applied rewrites24.5%
Final simplification28.2%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(let* ((t_0 (+ (- Vef mu) (+ Ev EAccept)))
(t_1 (/ NaChar (+ (exp (/ EAccept KbT)) 1.0)))
(t_2
(+
(/ NdChar (+ (exp (/ (+ mu (+ EDonor (- Vef Ec))) KbT)) 1.0))
(/ NaChar (+ (exp (/ (- (+ (+ Vef Ev) EAccept) mu) KbT)) 1.0)))))
(if (<= t_2 -1e-282)
t_1
(if (<= t_2 0.0)
(/
NaChar
(-
2.0
(/ (fma -0.5 (/ (* t_0 t_0) KbT) (- (- mu Vef) (+ Ev EAccept))) KbT)))
t_1))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = (Vef - mu) + (Ev + EAccept);
double t_1 = NaChar / (exp((EAccept / KbT)) + 1.0);
double t_2 = (NdChar / (exp(((mu + (EDonor + (Vef - Ec))) / KbT)) + 1.0)) + (NaChar / (exp(((((Vef + Ev) + EAccept) - mu) / KbT)) + 1.0));
double tmp;
if (t_2 <= -1e-282) {
tmp = t_1;
} else if (t_2 <= 0.0) {
tmp = NaChar / (2.0 - (fma(-0.5, ((t_0 * t_0) / KbT), ((mu - Vef) - (Ev + EAccept))) / KbT));
} else {
tmp = t_1;
}
return tmp;
}
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = Float64(Float64(Vef - mu) + Float64(Ev + EAccept)) t_1 = Float64(NaChar / Float64(exp(Float64(EAccept / KbT)) + 1.0)) t_2 = Float64(Float64(NdChar / Float64(exp(Float64(Float64(mu + Float64(EDonor + Float64(Vef - Ec))) / KbT)) + 1.0)) + Float64(NaChar / Float64(exp(Float64(Float64(Float64(Float64(Vef + Ev) + EAccept) - mu) / KbT)) + 1.0))) tmp = 0.0 if (t_2 <= -1e-282) tmp = t_1; elseif (t_2 <= 0.0) tmp = Float64(NaChar / Float64(2.0 - Float64(fma(-0.5, Float64(Float64(t_0 * t_0) / KbT), Float64(Float64(mu - Vef) - Float64(Ev + EAccept))) / KbT))); else tmp = t_1; end return tmp end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := Block[{t$95$0 = N[(N[(Vef - mu), $MachinePrecision] + N[(Ev + EAccept), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(NaChar / N[(N[Exp[N[(EAccept / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(NdChar / N[(N[Exp[N[(N[(mu + N[(EDonor + N[(Vef - Ec), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + N[(NaChar / N[(N[Exp[N[(N[(N[(N[(Vef + Ev), $MachinePrecision] + EAccept), $MachinePrecision] - mu), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -1e-282], t$95$1, If[LessEqual[t$95$2, 0.0], N[(NaChar / N[(2.0 - N[(N[(-0.5 * N[(N[(t$95$0 * t$95$0), $MachinePrecision] / KbT), $MachinePrecision] + N[(N[(mu - Vef), $MachinePrecision] - N[(Ev + EAccept), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(Vef - mu\right) + \left(Ev + EAccept\right)\\
t_1 := \frac{NaChar}{e^{\frac{EAccept}{KbT}} + 1}\\
t_2 := \frac{NdChar}{e^{\frac{mu + \left(EDonor + \left(Vef - Ec\right)\right)}{KbT}} + 1} + \frac{NaChar}{e^{\frac{\left(\left(Vef + Ev\right) + EAccept\right) - mu}{KbT}} + 1}\\
\mathbf{if}\;t\_2 \leq -1 \cdot 10^{-282}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq 0:\\
\;\;\;\;\frac{NaChar}{2 - \frac{\mathsf{fma}\left(-0.5, \frac{t\_0 \cdot t\_0}{KbT}, \left(mu - Vef\right) - \left(Ev + EAccept\right)\right)}{KbT}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < -1e-282 or 0.0 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) Initial program 99.9%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6453.8
Applied rewrites53.8%
Taylor expanded in EAccept around inf
lower-/.f6433.9
Applied rewrites33.9%
if -1e-282 < (+.f64 (/.f64 NdChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (neg.f64 (-.f64 (-.f64 (-.f64 Ec Vef) EDonor) mu)) KbT)))) (/.f64 NaChar (+.f64 #s(literal 1 binary64) (exp.f64 (/.f64 (+.f64 (+.f64 (+.f64 Ev Vef) EAccept) (neg.f64 mu)) KbT))))) < 0.0Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6497.9
Applied rewrites97.9%
Taylor expanded in KbT around -inf
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-/.f64N/A
Applied rewrites90.6%
Final simplification43.6%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(let* ((t_0 (/ NaChar (+ (exp (/ (+ EAccept (+ Ev (- Vef mu))) KbT)) 1.0))))
(if (<= NaChar -4e-104)
t_0
(if (<= NaChar 6.2e-72)
(/ NdChar (+ (exp (/ (+ (+ Vef EDonor) (- mu Ec)) KbT)) 1.0))
t_0))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = NaChar / (exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0);
double tmp;
if (NaChar <= -4e-104) {
tmp = t_0;
} else if (NaChar <= 6.2e-72) {
tmp = NdChar / (exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0);
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
real(8) :: t_0
real(8) :: tmp
t_0 = nachar / (exp(((eaccept + (ev + (vef - mu))) / kbt)) + 1.0d0)
if (nachar <= (-4d-104)) then
tmp = t_0
else if (nachar <= 6.2d-72) then
tmp = ndchar / (exp((((vef + edonor) + (mu - ec)) / kbt)) + 1.0d0)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = NaChar / (Math.exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0);
double tmp;
if (NaChar <= -4e-104) {
tmp = t_0;
} else if (NaChar <= 6.2e-72) {
tmp = NdChar / (Math.exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0);
} else {
tmp = t_0;
}
return tmp;
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): t_0 = NaChar / (math.exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0) tmp = 0 if NaChar <= -4e-104: tmp = t_0 elif NaChar <= 6.2e-72: tmp = NdChar / (math.exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0) else: tmp = t_0 return tmp
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = Float64(NaChar / Float64(exp(Float64(Float64(EAccept + Float64(Ev + Float64(Vef - mu))) / KbT)) + 1.0)) tmp = 0.0 if (NaChar <= -4e-104) tmp = t_0; elseif (NaChar <= 6.2e-72) tmp = Float64(NdChar / Float64(exp(Float64(Float64(Float64(Vef + EDonor) + Float64(mu - Ec)) / KbT)) + 1.0)); else tmp = t_0; end return tmp end
function tmp_2 = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = NaChar / (exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0); tmp = 0.0; if (NaChar <= -4e-104) tmp = t_0; elseif (NaChar <= 6.2e-72) tmp = NdChar / (exp((((Vef + EDonor) + (mu - Ec)) / KbT)) + 1.0); else tmp = t_0; end tmp_2 = tmp; end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := Block[{t$95$0 = N[(NaChar / N[(N[Exp[N[(N[(EAccept + N[(Ev + N[(Vef - mu), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[NaChar, -4e-104], t$95$0, If[LessEqual[NaChar, 6.2e-72], N[(NdChar / N[(N[Exp[N[(N[(N[(Vef + EDonor), $MachinePrecision] + N[(mu - Ec), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{NaChar}{e^{\frac{EAccept + \left(Ev + \left(Vef - mu\right)\right)}{KbT}} + 1}\\
\mathbf{if}\;NaChar \leq -4 \cdot 10^{-104}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;NaChar \leq 6.2 \cdot 10^{-72}:\\
\;\;\;\;\frac{NdChar}{e^{\frac{\left(Vef + EDonor\right) + \left(mu - Ec\right)}{KbT}} + 1}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if NaChar < -3.99999999999999971e-104 or 6.1999999999999996e-72 < NaChar Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6470.0
Applied rewrites70.0%
if -3.99999999999999971e-104 < NaChar < 6.1999999999999996e-72Initial program 100.0%
Taylor expanded in NdChar around inf
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6470.4
Applied rewrites70.4%
Final simplification70.2%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(if (<= KbT -1.85e+229)
(+ (* NaChar 0.5) (/ NdChar (+ (exp (/ EDonor KbT)) 1.0)))
(if (<= KbT 2.45e+241)
(/ NaChar (+ (exp (/ (+ EAccept (+ Ev (- Vef mu))) KbT)) 1.0))
(fma
-0.25
(* Ev (/ NaChar KbT))
(fma
0.25
(* NdChar (/ (- (- Ec EDonor) (+ Vef mu)) KbT))
(* 0.5 (+ NdChar NaChar)))))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double tmp;
if (KbT <= -1.85e+229) {
tmp = (NaChar * 0.5) + (NdChar / (exp((EDonor / KbT)) + 1.0));
} else if (KbT <= 2.45e+241) {
tmp = NaChar / (exp(((EAccept + (Ev + (Vef - mu))) / KbT)) + 1.0);
} else {
tmp = fma(-0.25, (Ev * (NaChar / KbT)), fma(0.25, (NdChar * (((Ec - EDonor) - (Vef + mu)) / KbT)), (0.5 * (NdChar + NaChar))));
}
return tmp;
}
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) tmp = 0.0 if (KbT <= -1.85e+229) tmp = Float64(Float64(NaChar * 0.5) + Float64(NdChar / Float64(exp(Float64(EDonor / KbT)) + 1.0))); elseif (KbT <= 2.45e+241) tmp = Float64(NaChar / Float64(exp(Float64(Float64(EAccept + Float64(Ev + Float64(Vef - mu))) / KbT)) + 1.0)); else tmp = fma(-0.25, Float64(Ev * Float64(NaChar / KbT)), fma(0.25, Float64(NdChar * Float64(Float64(Float64(Ec - EDonor) - Float64(Vef + mu)) / KbT)), Float64(0.5 * Float64(NdChar + NaChar)))); end return tmp end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := If[LessEqual[KbT, -1.85e+229], N[(N[(NaChar * 0.5), $MachinePrecision] + N[(NdChar / N[(N[Exp[N[(EDonor / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[KbT, 2.45e+241], N[(NaChar / N[(N[Exp[N[(N[(EAccept + N[(Ev + N[(Vef - mu), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], N[(-0.25 * N[(Ev * N[(NaChar / KbT), $MachinePrecision]), $MachinePrecision] + N[(0.25 * N[(NdChar * N[(N[(N[(Ec - EDonor), $MachinePrecision] - N[(Vef + mu), $MachinePrecision]), $MachinePrecision] / KbT), $MachinePrecision]), $MachinePrecision] + N[(0.5 * N[(NdChar + NaChar), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;KbT \leq -1.85 \cdot 10^{+229}:\\
\;\;\;\;NaChar \cdot 0.5 + \frac{NdChar}{e^{\frac{EDonor}{KbT}} + 1}\\
\mathbf{elif}\;KbT \leq 2.45 \cdot 10^{+241}:\\
\;\;\;\;\frac{NaChar}{e^{\frac{EAccept + \left(Ev + \left(Vef - mu\right)\right)}{KbT}} + 1}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-0.25, Ev \cdot \frac{NaChar}{KbT}, \mathsf{fma}\left(0.25, NdChar \cdot \frac{\left(Ec - EDonor\right) - \left(Vef + mu\right)}{KbT}, 0.5 \cdot \left(NdChar + NaChar\right)\right)\right)\\
\end{array}
\end{array}
if KbT < -1.85000000000000001e229Initial program 100.0%
Taylor expanded in KbT around inf
*-commutativeN/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in EDonor around inf
lower-/.f64100.0
Applied rewrites100.0%
if -1.85000000000000001e229 < KbT < 2.44999999999999986e241Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6463.4
Applied rewrites63.4%
if 2.44999999999999986e241 < KbT Initial program 99.6%
Applied rewrites99.6%
Taylor expanded in KbT around -inf
lower-fma.f64N/A
lower-/.f64N/A
lower-*.f64N/A
associate--l+N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f64N/A
lower-fma.f64N/A
Applied rewrites61.7%
Taylor expanded in Ev around inf
associate-/l*N/A
lower-*.f64N/A
lower-/.f6487.7
Applied rewrites87.7%
Final simplification66.0%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(let* ((t_0 (/ NaChar (+ (exp (/ Vef KbT)) 1.0))))
(if (<= Vef -5.5e+35)
t_0
(if (<= Vef -4.2e-263)
(/ NaChar (+ (exp (/ mu (- KbT))) 1.0))
(if (<= Vef 5e+66) (/ NaChar (+ (exp (/ EAccept KbT)) 1.0)) t_0)))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = NaChar / (exp((Vef / KbT)) + 1.0);
double tmp;
if (Vef <= -5.5e+35) {
tmp = t_0;
} else if (Vef <= -4.2e-263) {
tmp = NaChar / (exp((mu / -KbT)) + 1.0);
} else if (Vef <= 5e+66) {
tmp = NaChar / (exp((EAccept / KbT)) + 1.0);
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
real(8) :: t_0
real(8) :: tmp
t_0 = nachar / (exp((vef / kbt)) + 1.0d0)
if (vef <= (-5.5d+35)) then
tmp = t_0
else if (vef <= (-4.2d-263)) then
tmp = nachar / (exp((mu / -kbt)) + 1.0d0)
else if (vef <= 5d+66) then
tmp = nachar / (exp((eaccept / kbt)) + 1.0d0)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = NaChar / (Math.exp((Vef / KbT)) + 1.0);
double tmp;
if (Vef <= -5.5e+35) {
tmp = t_0;
} else if (Vef <= -4.2e-263) {
tmp = NaChar / (Math.exp((mu / -KbT)) + 1.0);
} else if (Vef <= 5e+66) {
tmp = NaChar / (Math.exp((EAccept / KbT)) + 1.0);
} else {
tmp = t_0;
}
return tmp;
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): t_0 = NaChar / (math.exp((Vef / KbT)) + 1.0) tmp = 0 if Vef <= -5.5e+35: tmp = t_0 elif Vef <= -4.2e-263: tmp = NaChar / (math.exp((mu / -KbT)) + 1.0) elif Vef <= 5e+66: tmp = NaChar / (math.exp((EAccept / KbT)) + 1.0) else: tmp = t_0 return tmp
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = Float64(NaChar / Float64(exp(Float64(Vef / KbT)) + 1.0)) tmp = 0.0 if (Vef <= -5.5e+35) tmp = t_0; elseif (Vef <= -4.2e-263) tmp = Float64(NaChar / Float64(exp(Float64(mu / Float64(-KbT))) + 1.0)); elseif (Vef <= 5e+66) tmp = Float64(NaChar / Float64(exp(Float64(EAccept / KbT)) + 1.0)); else tmp = t_0; end return tmp end
function tmp_2 = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = NaChar / (exp((Vef / KbT)) + 1.0); tmp = 0.0; if (Vef <= -5.5e+35) tmp = t_0; elseif (Vef <= -4.2e-263) tmp = NaChar / (exp((mu / -KbT)) + 1.0); elseif (Vef <= 5e+66) tmp = NaChar / (exp((EAccept / KbT)) + 1.0); else tmp = t_0; end tmp_2 = tmp; end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := Block[{t$95$0 = N[(NaChar / N[(N[Exp[N[(Vef / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[Vef, -5.5e+35], t$95$0, If[LessEqual[Vef, -4.2e-263], N[(NaChar / N[(N[Exp[N[(mu / (-KbT)), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[Vef, 5e+66], N[(NaChar / N[(N[Exp[N[(EAccept / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{NaChar}{e^{\frac{Vef}{KbT}} + 1}\\
\mathbf{if}\;Vef \leq -5.5 \cdot 10^{+35}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;Vef \leq -4.2 \cdot 10^{-263}:\\
\;\;\;\;\frac{NaChar}{e^{\frac{mu}{-KbT}} + 1}\\
\mathbf{elif}\;Vef \leq 5 \cdot 10^{+66}:\\
\;\;\;\;\frac{NaChar}{e^{\frac{EAccept}{KbT}} + 1}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if Vef < -5.50000000000000001e35 or 4.99999999999999991e66 < Vef Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6465.0
Applied rewrites65.0%
Taylor expanded in Vef around inf
lower-/.f6455.6
Applied rewrites55.6%
if -5.50000000000000001e35 < Vef < -4.20000000000000005e-263Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6458.0
Applied rewrites58.0%
Taylor expanded in mu around inf
associate-*r/N/A
mul-1-negN/A
lower-/.f64N/A
lower-neg.f6445.1
Applied rewrites45.1%
if -4.20000000000000005e-263 < Vef < 4.99999999999999991e66Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6460.0
Applied rewrites60.0%
Taylor expanded in EAccept around inf
lower-/.f6441.9
Applied rewrites41.9%
Final simplification48.1%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:precision binary64
(let* ((t_0 (/ NaChar (+ (exp (/ Vef KbT)) 1.0))))
(if (<= Vef -3.4e+19)
t_0
(if (<= Vef 5e+66) (/ NaChar (+ (exp (/ EAccept KbT)) 1.0)) t_0))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = NaChar / (exp((Vef / KbT)) + 1.0);
double tmp;
if (Vef <= -3.4e+19) {
tmp = t_0;
} else if (Vef <= 5e+66) {
tmp = NaChar / (exp((EAccept / KbT)) + 1.0);
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
real(8) :: t_0
real(8) :: tmp
t_0 = nachar / (exp((vef / kbt)) + 1.0d0)
if (vef <= (-3.4d+19)) then
tmp = t_0
else if (vef <= 5d+66) then
tmp = nachar / (exp((eaccept / kbt)) + 1.0d0)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double t_0 = NaChar / (Math.exp((Vef / KbT)) + 1.0);
double tmp;
if (Vef <= -3.4e+19) {
tmp = t_0;
} else if (Vef <= 5e+66) {
tmp = NaChar / (Math.exp((EAccept / KbT)) + 1.0);
} else {
tmp = t_0;
}
return tmp;
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): t_0 = NaChar / (math.exp((Vef / KbT)) + 1.0) tmp = 0 if Vef <= -3.4e+19: tmp = t_0 elif Vef <= 5e+66: tmp = NaChar / (math.exp((EAccept / KbT)) + 1.0) else: tmp = t_0 return tmp
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = Float64(NaChar / Float64(exp(Float64(Vef / KbT)) + 1.0)) tmp = 0.0 if (Vef <= -3.4e+19) tmp = t_0; elseif (Vef <= 5e+66) tmp = Float64(NaChar / Float64(exp(Float64(EAccept / KbT)) + 1.0)); else tmp = t_0; end return tmp end
function tmp_2 = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) t_0 = NaChar / (exp((Vef / KbT)) + 1.0); tmp = 0.0; if (Vef <= -3.4e+19) tmp = t_0; elseif (Vef <= 5e+66) tmp = NaChar / (exp((EAccept / KbT)) + 1.0); else tmp = t_0; end tmp_2 = tmp; end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := Block[{t$95$0 = N[(NaChar / N[(N[Exp[N[(Vef / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[Vef, -3.4e+19], t$95$0, If[LessEqual[Vef, 5e+66], N[(NaChar / N[(N[Exp[N[(EAccept / KbT), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{NaChar}{e^{\frac{Vef}{KbT}} + 1}\\
\mathbf{if}\;Vef \leq -3.4 \cdot 10^{+19}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;Vef \leq 5 \cdot 10^{+66}:\\
\;\;\;\;\frac{NaChar}{e^{\frac{EAccept}{KbT}} + 1}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if Vef < -3.4e19 or 4.99999999999999991e66 < Vef Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6463.8
Applied rewrites63.8%
Taylor expanded in Vef around inf
lower-/.f6454.6
Applied rewrites54.6%
if -3.4e19 < Vef < 4.99999999999999991e66Initial program 100.0%
Taylor expanded in NdChar around 0
lower-/.f64N/A
lower-+.f64N/A
lower-exp.f64N/A
lower-/.f64N/A
associate--l+N/A
lower-+.f64N/A
sub-negN/A
associate-+r+N/A
mul-1-negN/A
lower-+.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6459.8
Applied rewrites59.8%
Taylor expanded in EAccept around inf
lower-/.f6440.6
Applied rewrites40.6%
Final simplification46.1%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept) :precision binary64 (if (<= NaChar -2.9) (* NaChar 0.5) (if (<= NaChar 9e-38) (* NdChar 0.5) (* NaChar 0.5))))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double tmp;
if (NaChar <= -2.9) {
tmp = NaChar * 0.5;
} else if (NaChar <= 9e-38) {
tmp = NdChar * 0.5;
} else {
tmp = NaChar * 0.5;
}
return tmp;
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
real(8) :: tmp
if (nachar <= (-2.9d0)) then
tmp = nachar * 0.5d0
else if (nachar <= 9d-38) then
tmp = ndchar * 0.5d0
else
tmp = nachar * 0.5d0
end if
code = tmp
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
double tmp;
if (NaChar <= -2.9) {
tmp = NaChar * 0.5;
} else if (NaChar <= 9e-38) {
tmp = NdChar * 0.5;
} else {
tmp = NaChar * 0.5;
}
return tmp;
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): tmp = 0 if NaChar <= -2.9: tmp = NaChar * 0.5 elif NaChar <= 9e-38: tmp = NdChar * 0.5 else: tmp = NaChar * 0.5 return tmp
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) tmp = 0.0 if (NaChar <= -2.9) tmp = Float64(NaChar * 0.5); elseif (NaChar <= 9e-38) tmp = Float64(NdChar * 0.5); else tmp = Float64(NaChar * 0.5); end return tmp end
function tmp_2 = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) tmp = 0.0; if (NaChar <= -2.9) tmp = NaChar * 0.5; elseif (NaChar <= 9e-38) tmp = NdChar * 0.5; else tmp = NaChar * 0.5; end tmp_2 = tmp; end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := If[LessEqual[NaChar, -2.9], N[(NaChar * 0.5), $MachinePrecision], If[LessEqual[NaChar, 9e-38], N[(NdChar * 0.5), $MachinePrecision], N[(NaChar * 0.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;NaChar \leq -2.9:\\
\;\;\;\;NaChar \cdot 0.5\\
\mathbf{elif}\;NaChar \leq 9 \cdot 10^{-38}:\\
\;\;\;\;NdChar \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;NaChar \cdot 0.5\\
\end{array}
\end{array}
if NaChar < -2.89999999999999991 or 9.00000000000000018e-38 < NaChar Initial program 100.0%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6421.5
Applied rewrites21.5%
Taylor expanded in NaChar around inf
*-commutativeN/A
lower-*.f6421.5
Applied rewrites21.5%
if -2.89999999999999991 < NaChar < 9.00000000000000018e-38Initial program 99.9%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6423.5
Applied rewrites23.5%
Taylor expanded in NaChar around 0
*-commutativeN/A
lower-*.f6422.2
Applied rewrites22.2%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept) :precision binary64 (* 0.5 (+ NdChar NaChar)))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
return 0.5 * (NdChar + NaChar);
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
code = 0.5d0 * (ndchar + nachar)
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
return 0.5 * (NdChar + NaChar);
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): return 0.5 * (NdChar + NaChar)
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) return Float64(0.5 * Float64(NdChar + NaChar)) end
function tmp = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) tmp = 0.5 * (NdChar + NaChar); end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := N[(0.5 * N[(NdChar + NaChar), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.5 \cdot \left(NdChar + NaChar\right)
\end{array}
Initial program 100.0%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6422.4
Applied rewrites22.4%
Final simplification22.4%
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept) :precision binary64 (* NaChar 0.5))
double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
return NaChar * 0.5;
}
real(8) function code(ndchar, ec, vef, edonor, mu, kbt, nachar, ev, eaccept)
real(8), intent (in) :: ndchar
real(8), intent (in) :: ec
real(8), intent (in) :: vef
real(8), intent (in) :: edonor
real(8), intent (in) :: mu
real(8), intent (in) :: kbt
real(8), intent (in) :: nachar
real(8), intent (in) :: ev
real(8), intent (in) :: eaccept
code = nachar * 0.5d0
end function
public static double code(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
return NaChar * 0.5;
}
def code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept): return NaChar * 0.5
function code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) return Float64(NaChar * 0.5) end
function tmp = code(NdChar, Ec, Vef, EDonor, mu, KbT, NaChar, Ev, EAccept) tmp = NaChar * 0.5; end
code[NdChar_, Ec_, Vef_, EDonor_, mu_, KbT_, NaChar_, Ev_, EAccept_] := N[(NaChar * 0.5), $MachinePrecision]
\begin{array}{l}
\\
NaChar \cdot 0.5
\end{array}
Initial program 100.0%
Taylor expanded in KbT around inf
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f6422.4
Applied rewrites22.4%
Taylor expanded in NaChar around inf
*-commutativeN/A
lower-*.f6416.8
Applied rewrites16.8%
herbie shell --seed 2024219
(FPCore (NdChar Ec Vef EDonor mu KbT NaChar Ev EAccept)
:name "Bulmash initializePoisson"
:precision binary64
(+ (/ NdChar (+ 1.0 (exp (/ (- (- (- (- Ec Vef) EDonor) mu)) KbT)))) (/ NaChar (+ 1.0 (exp (/ (+ (+ (+ Ev Vef) EAccept) (- mu)) KbT))))))