Toniolo and Linder, Equation (13)

Average Error: 34.9 → 29.5

Time: 38.4s

Precision: 64

Math

\[\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}\;t \le -7.152776844878304065152644239470045087319 \cdot 10^{-89}:\\ \;\;\;\;\sqrt{\left(\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(2 \cdot n\right)\right) \cdot U}\\ \mathbf{elif}\;t \le 6.298121475305343130041003071629157124028 \cdot 10^{116}:\\ \;\;\;\;\sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \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 \left(U - U*\right)\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\\ \mathbf{else}:\\ \;\;\;\;\sqrt{t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)} \cdot \sqrt{\left(2 \cdot n\right) \cdot U}\\ \end{array}\]

C

double f(double n, double U, double t, double l, double Om, double U_) {
        double r99286 = 2.0;
        double r99287 = n;
        double r99288 = r99286 * r99287;
        double r99289 = U;
        double r99290 = r99288 * r99289;
        double r99291 = t;
        double r99292 = l;
        double r99293 = r99292 * r99292;
        double r99294 = Om;
        double r99295 = r99293 / r99294;
        double r99296 = r99286 * r99295;
        double r99297 = r99291 - r99296;
        double r99298 = r99292 / r99294;
        double r99299 = pow(r99298, r99286);
        double r99300 = r99287 * r99299;
        double r99301 = U_;
        double r99302 = r99289 - r99301;
        double r99303 = r99300 * r99302;
        double r99304 = r99297 - r99303;
        double r99305 = r99290 * r99304;
        double r99306 = sqrt(r99305);
        return r99306;
}

↓

double f(double n, double U, double t, double l, double Om, double U_) {
        double r99307 = t;
        double r99308 = -7.152776844878304e-89;
        bool r99309 = r99307 <= r99308;
        double r99310 = 2.0;
        double r99311 = l;
        double r99312 = Om;
        double r99313 = r99311 / r99312;
        double r99314 = r99311 * r99313;
        double r99315 = n;
        double r99316 = pow(r99313, r99310);
        double r99317 = r99315 * r99316;
        double r99318 = U;
        double r99319 = U_;
        double r99320 = r99318 - r99319;
        double r99321 = r99317 * r99320;
        double r99322 = fma(r99310, r99314, r99321);
        double r99323 = r99307 - r99322;
        double r99324 = r99310 * r99315;
        double r99325 = r99323 * r99324;
        double r99326 = r99325 * r99318;
        double r99327 = sqrt(r99326);
        double r99328 = 6.298121475305343e+116;
        bool r99329 = r99307 <= r99328;
        double r99330 = 2.0;
        double r99331 = r99310 / r99330;
        double r99332 = pow(r99313, r99331);
        double r99333 = r99315 * r99332;
        double r99334 = r99332 * r99320;
        double r99335 = r99333 * r99334;
        double r99336 = fma(r99310, r99314, r99335);
        double r99337 = r99307 - r99336;
        double r99338 = r99324 * r99318;
        double r99339 = r99337 * r99338;
        double r99340 = sqrt(r99339);
        double r99341 = sqrt(r99323);
        double r99342 = sqrt(r99338);
        double r99343 = r99341 * r99342;
        double r99344 = r99329 ? r99340 : r99343;
        double r99345 = r99309 ? r99327 : r99344;
        return r99345;
}

Error

Bits error versus n

Bits error versus U

Bits error versus t

Bits error versus l

Bits error versus Om

Bits error versus U*

Derivation

1. Split input into 3 regimes

2. if t < -7.152776844878304e-89

   1. Initial program 33.3

      \[\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)}\]

   2. Simplified 33.3

      \[\leadsto \color{blue}{\sqrt{\left(t - \mathsf{fma}\left(2, \frac{\ell \cdot \ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}}\]
   3. Using strategy rm

   4. Applied *-un-lft-identity 33.3

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \frac{\ell \cdot \ell}{\color{blue}{1 \cdot Om}}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   5. Applied times-frac 31.3

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \color{blue}{\frac{\ell}{1} \cdot \frac{\ell}{Om}}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   6. Simplified 31.3

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \color{blue}{\ell} \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]
   7. Using strategy rm

   8. Applied associate-*r* 30.3

      \[\leadsto \sqrt{\color{blue}{\left(\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(2 \cdot n\right)\right) \cdot U}}\]

   if -7.152776844878304e-89 < t < 6.298121475305343e+116

   1. Initial program 35.0

      \[\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)}\]

   2. Simplified 35.0

      \[\leadsto \color{blue}{\sqrt{\left(t - \mathsf{fma}\left(2, \frac{\ell \cdot \ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}}\]
   3. Using strategy rm

   4. Applied *-un-lft-identity 35.0

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \frac{\ell \cdot \ell}{\color{blue}{1 \cdot Om}}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   5. Applied times-frac 32.0

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \color{blue}{\frac{\ell}{1} \cdot \frac{\ell}{Om}}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   6. Simplified 32.0

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \color{blue}{\ell} \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]
   7. Using strategy rm

   8. Applied sqr-pow 32.0

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot \color{blue}{\left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right)}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   9. Applied associate-*r* 30.8

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \color{blue}{\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)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]
   10. Using strategy rm

   11. Applied associate-*l* 30.3

       \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \color{blue}{\left(n \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 \left(U - U*\right)\right)}\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   if 6.298121475305343e+116 < t

   1. Initial program 37.8

      \[\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)}\]

   2. Simplified 37.8

      \[\leadsto \color{blue}{\sqrt{\left(t - \mathsf{fma}\left(2, \frac{\ell \cdot \ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}}\]
   3. Using strategy rm

   4. Applied *-un-lft-identity 37.8

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \frac{\ell \cdot \ell}{\color{blue}{1 \cdot Om}}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   5. Applied times-frac 35.3

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \color{blue}{\frac{\ell}{1} \cdot \frac{\ell}{Om}}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]

   6. Simplified 35.3

      \[\leadsto \sqrt{\left(t - \mathsf{fma}\left(2, \color{blue}{\ell} \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\]
   7. Using strategy rm

   8. Applied sqrt-prod 25.5

      \[\leadsto \color{blue}{\sqrt{t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)} \cdot \sqrt{\left(2 \cdot n\right) \cdot U}}\]
3. Recombined 3 regimes into one program.

4. Final simplification 29.5

   \[\leadsto \begin{array}{l} \mathbf{if}\;t \le -7.152776844878304065152644239470045087319 \cdot 10^{-89}:\\ \;\;\;\;\sqrt{\left(\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)\right) \cdot \left(2 \cdot n\right)\right) \cdot U}\\ \mathbf{elif}\;t \le 6.298121475305343130041003071629157124028 \cdot 10^{116}:\\ \;\;\;\;\sqrt{\left(t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \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 \left(U - U*\right)\right)\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}\\ \mathbf{else}:\\ \;\;\;\;\sqrt{t - \mathsf{fma}\left(2, \ell \cdot \frac{\ell}{Om}, \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)} \cdot \sqrt{\left(2 \cdot n\right) \cdot U}\\ \end{array}\]

Reproduce

herbie shell --seed 2019347 +o rules:numerics
(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*))))))