\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 -4.458146329554135732988892774226101337842 \cdot 10^{62}:\\
\;\;\;\;\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\frac{\ell}{\sqrt[3]{Om}} \cdot \frac{\ell}{\sqrt[3]{Om}}}{\sqrt[3]{Om}}\right) - \left(U - U*\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right)\right)}\\
\mathbf{elif}\;U \le 1.43075486013443336421549580762136466865 \cdot 10^{-67}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(\left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(U - U*\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot n\right)\right) \cdot U\right)}\\
\mathbf{elif}\;U \le 3.638498443778453206704785881478597925308 \cdot 10^{-26}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot U} \cdot \sqrt{\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(U - U*\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot n\right)}\\
\mathbf{elif}\;U \le 1.219991701569961735518499105488225569099 \cdot 10^{291}:\\
\;\;\;\;\sqrt{\sqrt{\left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(U - U*\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot n\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}} \cdot \sqrt{\sqrt{\left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(U - U*\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot n\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot U} \cdot \sqrt{\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(U - U*\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot n\right)}\\
\end{array}double f(double n, double U, double t, double l, double Om, double U_) {
double r2930154 = 2.0;
double r2930155 = n;
double r2930156 = r2930154 * r2930155;
double r2930157 = U;
double r2930158 = r2930156 * r2930157;
double r2930159 = t;
double r2930160 = l;
double r2930161 = r2930160 * r2930160;
double r2930162 = Om;
double r2930163 = r2930161 / r2930162;
double r2930164 = r2930154 * r2930163;
double r2930165 = r2930159 - r2930164;
double r2930166 = r2930160 / r2930162;
double r2930167 = pow(r2930166, r2930154);
double r2930168 = r2930155 * r2930167;
double r2930169 = U_;
double r2930170 = r2930157 - r2930169;
double r2930171 = r2930168 * r2930170;
double r2930172 = r2930165 - r2930171;
double r2930173 = r2930158 * r2930172;
double r2930174 = sqrt(r2930173);
return r2930174;
}
double f(double n, double U, double t, double l, double Om, double U_) {
double r2930175 = U;
double r2930176 = -4.458146329554136e+62;
bool r2930177 = r2930175 <= r2930176;
double r2930178 = 2.0;
double r2930179 = n;
double r2930180 = r2930178 * r2930179;
double r2930181 = r2930180 * r2930175;
double r2930182 = t;
double r2930183 = l;
double r2930184 = Om;
double r2930185 = cbrt(r2930184);
double r2930186 = r2930183 / r2930185;
double r2930187 = r2930186 * r2930186;
double r2930188 = r2930187 / r2930185;
double r2930189 = r2930178 * r2930188;
double r2930190 = r2930182 - r2930189;
double r2930191 = U_;
double r2930192 = r2930175 - r2930191;
double r2930193 = r2930183 / r2930184;
double r2930194 = pow(r2930193, r2930178);
double r2930195 = r2930179 * r2930194;
double r2930196 = r2930192 * r2930195;
double r2930197 = r2930190 - r2930196;
double r2930198 = r2930181 * r2930197;
double r2930199 = sqrt(r2930198);
double r2930200 = 1.4307548601344334e-67;
bool r2930201 = r2930175 <= r2930200;
double r2930202 = r2930183 * r2930193;
double r2930203 = r2930178 * r2930202;
double r2930204 = r2930182 - r2930203;
double r2930205 = 2.0;
double r2930206 = r2930178 / r2930205;
double r2930207 = pow(r2930193, r2930206);
double r2930208 = r2930192 * r2930207;
double r2930209 = r2930207 * r2930179;
double r2930210 = r2930208 * r2930209;
double r2930211 = r2930204 - r2930210;
double r2930212 = r2930211 * r2930175;
double r2930213 = r2930180 * r2930212;
double r2930214 = sqrt(r2930213);
double r2930215 = 3.638498443778453e-26;
bool r2930216 = r2930175 <= r2930215;
double r2930217 = sqrt(r2930181);
double r2930218 = sqrt(r2930211);
double r2930219 = r2930217 * r2930218;
double r2930220 = 1.2199917015699617e+291;
bool r2930221 = r2930175 <= r2930220;
double r2930222 = r2930211 * r2930181;
double r2930223 = sqrt(r2930222);
double r2930224 = sqrt(r2930223);
double r2930225 = r2930224 * r2930224;
double r2930226 = r2930221 ? r2930225 : r2930219;
double r2930227 = r2930216 ? r2930219 : r2930226;
double r2930228 = r2930201 ? r2930214 : r2930227;
double r2930229 = r2930177 ? r2930199 : r2930228;
return r2930229;
}



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 < -4.458146329554136e+62Initial program 30.0
rmApplied add-cube-cbrt30.1
Applied associate-/r*30.1
Simplified27.8
if -4.458146329554136e+62 < U < 1.4307548601344334e-67Initial program 37.2
rmApplied *-un-lft-identity37.2
Applied times-frac34.4
Simplified34.4
rmApplied sqr-pow34.4
Applied associate-*r*33.3
rmApplied associate-*l*32.9
rmApplied associate-*l*29.3
if 1.4307548601344334e-67 < U < 3.638498443778453e-26 or 1.2199917015699617e+291 < U Initial program 30.9
rmApplied *-un-lft-identity30.9
Applied times-frac28.6
Simplified28.6
rmApplied sqr-pow28.6
Applied associate-*r*27.6
rmApplied associate-*l*28.2
rmApplied sqrt-prod41.8
if 3.638498443778453e-26 < U < 1.2199917015699617e+291Initial program 29.6
rmApplied *-un-lft-identity29.6
Applied times-frac26.1
Simplified26.1
rmApplied sqr-pow26.1
Applied associate-*r*25.4
rmApplied associate-*l*26.0
rmApplied add-sqr-sqrt26.2
Final simplification29.0
herbie shell --seed 2019168
(FPCore (n U t l Om U*)
:name "Toniolo and Linder, Equation (13)"
(sqrt (* (* (* 2.0 n) U) (- (- t (* 2.0 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2.0)) (- U U*))))))