\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell \cdot \ell}{Om}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}\begin{array}{l}
\mathbf{if}\;U* \le -20628.158740343424:\\
\;\;\;\;\sqrt{{\left(\left(\left(2 \cdot n\right) \cdot \left(t - \left(2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right) - \left(-n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(2 \cdot \frac{2}{2}\right)}\right) \cdot \left(U - U*\right)\right)\right)\right) \cdot U\right)}^{1}}\\
\mathbf{elif}\;U* \le -3.38886052338754529 \cdot 10^{-282}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(t - \left(2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right) - \left(-n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(2 \cdot \frac{2}{2}\right)}\right) \cdot \left(U - U*\right)\right)\right)\right)}\\
\mathbf{elif}\;U* \le 2.66582174134472983 \cdot 10^{-120}:\\
\;\;\;\;\sqrt{{\left(\left(\left(2 \cdot n\right) \cdot \left(t - \left(2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right) - \left(-n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(2 \cdot \frac{2}{2}\right)}\right) \cdot \left(U - U*\right)\right)\right)\right) \cdot U\right)}^{1}}\\
\mathbf{elif}\;U* \le 6.37844133114075821 \cdot 10^{-28}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot U} \cdot \sqrt{\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left(U - U*\right)}\\
\mathbf{elif}\;U* \le 1.3559215670195436 \cdot 10^{141}:\\
\;\;\;\;\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(\left(n \cdot \ell\right) \cdot {\left(\frac{1}{{Om}^{1}}\right)}^{1}\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left(U - U*\right)\right)}\\
\mathbf{elif}\;U* \le 7.30429266906966537 \cdot 10^{231}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(t - \left(2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right) - \left(-n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(2 \cdot \frac{2}{2}\right)}\right) \cdot \left(U - U*\right)\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(\left(n \cdot \ell\right) \cdot {\left(\frac{1}{{Om}^{1}}\right)}^{1}\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left(U - U*\right)\right)}\\
\end{array}double f(double n, double U, double t, double l, double Om, double U_) {
double r176072 = 2.0;
double r176073 = n;
double r176074 = r176072 * r176073;
double r176075 = U;
double r176076 = r176074 * r176075;
double r176077 = t;
double r176078 = l;
double r176079 = r176078 * r176078;
double r176080 = Om;
double r176081 = r176079 / r176080;
double r176082 = r176072 * r176081;
double r176083 = r176077 - r176082;
double r176084 = r176078 / r176080;
double r176085 = pow(r176084, r176072);
double r176086 = r176073 * r176085;
double r176087 = U_;
double r176088 = r176075 - r176087;
double r176089 = r176086 * r176088;
double r176090 = r176083 - r176089;
double r176091 = r176076 * r176090;
double r176092 = sqrt(r176091);
return r176092;
}
double f(double n, double U, double t, double l, double Om, double U_) {
double r176093 = U_;
double r176094 = -20628.158740343424;
bool r176095 = r176093 <= r176094;
double r176096 = 2.0;
double r176097 = n;
double r176098 = r176096 * r176097;
double r176099 = t;
double r176100 = l;
double r176101 = Om;
double r176102 = r176100 / r176101;
double r176103 = r176100 * r176102;
double r176104 = r176096 * r176103;
double r176105 = 2.0;
double r176106 = r176096 / r176105;
double r176107 = r176105 * r176106;
double r176108 = pow(r176102, r176107);
double r176109 = r176097 * r176108;
double r176110 = -r176109;
double r176111 = U;
double r176112 = r176111 - r176093;
double r176113 = r176110 * r176112;
double r176114 = r176104 - r176113;
double r176115 = r176099 - r176114;
double r176116 = r176098 * r176115;
double r176117 = r176116 * r176111;
double r176118 = 1.0;
double r176119 = pow(r176117, r176118);
double r176120 = sqrt(r176119);
double r176121 = -3.3888605233875453e-282;
bool r176122 = r176093 <= r176121;
double r176123 = r176111 * r176115;
double r176124 = r176098 * r176123;
double r176125 = sqrt(r176124);
double r176126 = 2.66582174134473e-120;
bool r176127 = r176093 <= r176126;
double r176128 = 6.378441331140758e-28;
bool r176129 = r176093 <= r176128;
double r176130 = r176098 * r176111;
double r176131 = sqrt(r176130);
double r176132 = r176099 - r176104;
double r176133 = pow(r176102, r176106);
double r176134 = r176097 * r176133;
double r176135 = r176134 * r176133;
double r176136 = r176135 * r176112;
double r176137 = r176132 - r176136;
double r176138 = sqrt(r176137);
double r176139 = r176131 * r176138;
double r176140 = 1.3559215670195436e+141;
bool r176141 = r176093 <= r176140;
double r176142 = r176097 * r176100;
double r176143 = 1.0;
double r176144 = pow(r176101, r176143);
double r176145 = r176118 / r176144;
double r176146 = pow(r176145, r176143);
double r176147 = r176142 * r176146;
double r176148 = r176147 * r176133;
double r176149 = r176148 * r176112;
double r176150 = r176132 - r176149;
double r176151 = r176130 * r176150;
double r176152 = sqrt(r176151);
double r176153 = 7.304292669069665e+231;
bool r176154 = r176093 <= r176153;
double r176155 = r176154 ? r176125 : r176152;
double r176156 = r176141 ? r176152 : r176155;
double r176157 = r176129 ? r176139 : r176156;
double r176158 = r176127 ? r176120 : r176157;
double r176159 = r176122 ? r176125 : r176158;
double r176160 = r176095 ? r176120 : r176159;
return r176160;
}



Bits error versus n



Bits error versus U



Bits error versus t



Bits error versus l



Bits error versus Om



Bits error versus U*
Results
if U* < -20628.158740343424 or -3.3888605233875453e-282 < U* < 2.66582174134473e-120Initial program 35.7
rmApplied *-un-lft-identity35.7
Applied times-frac33.4
Simplified33.4
rmApplied sqr-pow33.4
Applied associate-*r*32.4
rmApplied pow132.4
Applied pow132.4
Applied pow132.4
Applied pow132.4
Applied pow-prod-down32.4
Applied pow-prod-down32.4
Applied pow-prod-down32.4
Simplified33.2
if -20628.158740343424 < U* < -3.3888605233875453e-282 or 1.3559215670195436e+141 < U* < 7.304292669069665e+231Initial program 33.3
rmApplied *-un-lft-identity33.3
Applied times-frac30.7
Simplified30.7
rmApplied sqr-pow30.7
Applied associate-*r*29.8
rmApplied associate-*l*30.6
Simplified31.8
if 2.66582174134473e-120 < U* < 6.378441331140758e-28Initial program 34.9
rmApplied *-un-lft-identity34.9
Applied times-frac31.7
Simplified31.7
rmApplied sqr-pow31.7
Applied associate-*r*30.8
rmApplied sqrt-prod42.6
if 6.378441331140758e-28 < U* < 1.3559215670195436e+141 or 7.304292669069665e+231 < U* Initial program 35.4
rmApplied *-un-lft-identity35.4
Applied times-frac32.8
Simplified32.8
rmApplied sqr-pow32.8
Applied associate-*r*31.6
Taylor expanded around 0 32.2
Final simplification33.2
herbie shell --seed 2020049
(FPCore (n U t l Om U*)
:name "Toniolo and Linder, Equation (13)"
:precision binary64
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))))