Statistics.Math.RootFinding:ridders from math-functions-0.1.5.2

Average Error: 25.2 → 6.3

Time: 52.6s

Precision: 64

Math

\[\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\]

↓

\[\begin{array}{l} \mathbf{if}\;z \le -4.611098198946549 \cdot 10^{108}:\\ \;\;\;\;-1 \cdot \left(x \cdot y\right)\\ \mathbf{elif}\;z \le 3.1007095082743632 \cdot 10^{95}:\\ \;\;\;\;\left(\frac{\sqrt[3]{y}}{\frac{\frac{\frac{\left|\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}\right|}{\sqrt[3]{z}}}{\sqrt[3]{z}}}{\sqrt[3]{y}}} \cdot \frac{x}{\sqrt{\sqrt{z \cdot z - t \cdot a}}}\right) \cdot \frac{\sqrt[3]{y}}{\frac{\sqrt{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}}{\sqrt[3]{z}}}\\ \mathbf{else}:\\ \;\;\;\;x \cdot y\\ \end{array}\]

C

double f(double x, double y, double z, double t, double a) {
        double r492310 = x;
        double r492311 = y;
        double r492312 = r492310 * r492311;
        double r492313 = z;
        double r492314 = r492312 * r492313;
        double r492315 = r492313 * r492313;
        double r492316 = t;
        double r492317 = a;
        double r492318 = r492316 * r492317;
        double r492319 = r492315 - r492318;
        double r492320 = sqrt(r492319);
        double r492321 = r492314 / r492320;
        return r492321;
}

↓

double f(double x, double y, double z, double t, double a) {
        double r492322 = z;
        double r492323 = -4.611098198946549e+108;
        bool r492324 = r492322 <= r492323;
        double r492325 = -1.0;
        double r492326 = x;
        double r492327 = y;
        double r492328 = r492326 * r492327;
        double r492329 = r492325 * r492328;
        double r492330 = 3.100709508274363e+95;
        bool r492331 = r492322 <= r492330;
        double r492332 = cbrt(r492327);
        double r492333 = r492322 * r492322;
        double r492334 = t;
        double r492335 = a;
        double r492336 = r492334 * r492335;
        double r492337 = r492333 - r492336;
        double r492338 = sqrt(r492337);
        double r492339 = cbrt(r492338);
        double r492340 = fabs(r492339);
        double r492341 = cbrt(r492322);
        double r492342 = r492340 / r492341;
        double r492343 = r492342 / r492341;
        double r492344 = r492343 / r492332;
        double r492345 = r492332 / r492344;
        double r492346 = sqrt(r492338);
        double r492347 = r492326 / r492346;
        double r492348 = r492345 * r492347;
        double r492349 = sqrt(r492339);
        double r492350 = r492349 / r492341;
        double r492351 = r492332 / r492350;
        double r492352 = r492348 * r492351;
        double r492353 = r492331 ? r492352 : r492328;
        double r492354 = r492324 ? r492329 : r492353;
        return r492354;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Bits error versus t

Bits error versus a

Target

Original	25.2
Target	7.9
Herbie	6.3

\[\begin{array}{l} \mathbf{if}\;z \lt -3.1921305903852764 \cdot 10^{46}:\\ \;\;\;\;-y \cdot x\\ \mathbf{elif}\;z \lt 5.9762681209208942 \cdot 10^{90}:\\ \;\;\;\;\frac{x \cdot z}{\frac{\sqrt{z \cdot z - a \cdot t}}{y}}\\ \mathbf{else}:\\ \;\;\;\;y \cdot x\\ \end{array}\]

Derivation

1. Split input into 3 regimes

2. if z < -4.611098198946549e+108

   1. Initial program 44.9

      \[\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\]

   2. Taylor expanded around -inf 2.4

      \[\leadsto \color{blue}{-1 \cdot \left(x \cdot y\right)}\]

   if -4.611098198946549e+108 < z < 3.100709508274363e+95

   1. Initial program 11.9

      \[\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\]
   2. Using strategy rm

   3. Applied associate-/l* 10.0

      \[\leadsto \color{blue}{\frac{x \cdot y}{\frac{\sqrt{z \cdot z - t \cdot a}}{z}}}\]
   4. Using strategy rm

   5. Applied *-un-lft-identity 10.0

      \[\leadsto \frac{x \cdot y}{\frac{\sqrt{z \cdot z - t \cdot a}}{\color{blue}{1 \cdot z}}}\]

   6. Applied add-sqr-sqrt 10.0

      \[\leadsto \frac{x \cdot y}{\frac{\sqrt{\color{blue}{\sqrt{z \cdot z - t \cdot a} \cdot \sqrt{z \cdot z - t \cdot a}}}}{1 \cdot z}}\]

   7. Applied sqrt-prod 10.3

      \[\leadsto \frac{x \cdot y}{\frac{\color{blue}{\sqrt{\sqrt{z \cdot z - t \cdot a}} \cdot \sqrt{\sqrt{z \cdot z - t \cdot a}}}}{1 \cdot z}}\]

   8. Applied times-frac 10.2

      \[\leadsto \frac{x \cdot y}{\color{blue}{\frac{\sqrt{\sqrt{z \cdot z - t \cdot a}}}{1} \cdot \frac{\sqrt{\sqrt{z \cdot z - t \cdot a}}}{z}}}\]

   9. Applied times-frac 11.2

      \[\leadsto \color{blue}{\frac{x}{\frac{\sqrt{\sqrt{z \cdot z - t \cdot a}}}{1}} \cdot \frac{y}{\frac{\sqrt{\sqrt{z \cdot z - t \cdot a}}}{z}}}\]

   10. Simplified 11.2

       \[\leadsto \color{blue}{\frac{x}{\sqrt{\sqrt{z \cdot z - t \cdot a}}}} \cdot \frac{y}{\frac{\sqrt{\sqrt{z \cdot z - t \cdot a}}}{z}}\]
   11. Using strategy rm

   12. Applied add-cube-cbrt 11.7

       \[\leadsto \frac{x}{\sqrt{\sqrt{z \cdot z - t \cdot a}}} \cdot \frac{y}{\frac{\sqrt{\sqrt{z \cdot z - t \cdot a}}}{\color{blue}{\left(\sqrt[3]{z} \cdot \sqrt[3]{z}\right) \cdot \sqrt[3]{z}}}}\]

   13. Applied add-cube-cbrt 11.5

       \[\leadsto \frac{x}{\sqrt{\sqrt{z \cdot z - t \cdot a}}} \cdot \frac{y}{\frac{\sqrt{\color{blue}{\left(\sqrt[3]{\sqrt{z \cdot z - t \cdot a}} \cdot \sqrt[3]{\sqrt{z \cdot z - t \cdot a}}\right) \cdot \sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}}}{\left(\sqrt[3]{z} \cdot \sqrt[3]{z}\right) \cdot \sqrt[3]{z}}}\]

   14. Applied sqrt-prod 11.6

       \[\leadsto \frac{x}{\sqrt{\sqrt{z \cdot z - t \cdot a}}} \cdot \frac{y}{\frac{\color{blue}{\sqrt{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}} \cdot \sqrt[3]{\sqrt{z \cdot z - t \cdot a}}} \cdot \sqrt{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}}}{\left(\sqrt[3]{z} \cdot \sqrt[3]{z}\right) \cdot \sqrt[3]{z}}}\]

   15. Applied times-frac 11.5

       \[\leadsto \frac{x}{\sqrt{\sqrt{z \cdot z - t \cdot a}}} \cdot \frac{y}{\color{blue}{\frac{\sqrt{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}} \cdot \sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}}{\sqrt[3]{z} \cdot \sqrt[3]{z}} \cdot \frac{\sqrt{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}}{\sqrt[3]{z}}}}\]

   16. Applied add-cube-cbrt 11.8

       \[\leadsto \frac{x}{\sqrt{\sqrt{z \cdot z - t \cdot a}}} \cdot \frac{\color{blue}{\left(\sqrt[3]{y} \cdot \sqrt[3]{y}\right) \cdot \sqrt[3]{y}}}{\frac{\sqrt{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}} \cdot \sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}}{\sqrt[3]{z} \cdot \sqrt[3]{z}} \cdot \frac{\sqrt{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}}{\sqrt[3]{z}}}\]

   17. Applied times-frac 11.3

       \[\leadsto \frac{x}{\sqrt{\sqrt{z \cdot z - t \cdot a}}} \cdot \color{blue}{\left(\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\frac{\sqrt{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}} \cdot \sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}}{\sqrt[3]{z} \cdot \sqrt[3]{z}}} \cdot \frac{\sqrt[3]{y}}{\frac{\sqrt{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}}{\sqrt[3]{z}}}\right)}\]

   18. Applied associate-*r* 9.0

       \[\leadsto \color{blue}{\left(\frac{x}{\sqrt{\sqrt{z \cdot z - t \cdot a}}} \cdot \frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\frac{\sqrt{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}} \cdot \sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}}{\sqrt[3]{z} \cdot \sqrt[3]{z}}}\right) \cdot \frac{\sqrt[3]{y}}{\frac{\sqrt{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}}{\sqrt[3]{z}}}}\]

   19. Simplified 9.0

       \[\leadsto \color{blue}{\left(\frac{\sqrt[3]{y}}{\frac{\frac{\frac{\left|\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}\right|}{\sqrt[3]{z}}}{\sqrt[3]{z}}}{\sqrt[3]{y}}} \cdot \frac{x}{\sqrt{\sqrt{z \cdot z - t \cdot a}}}\right)} \cdot \frac{\sqrt[3]{y}}{\frac{\sqrt{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}}{\sqrt[3]{z}}}\]

   if 3.100709508274363e+95 < z

   1. Initial program 42.8

      \[\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\]

   2. Taylor expanded around inf 2.5

      \[\leadsto \color{blue}{x \cdot y}\]
3. Recombined 3 regimes into one program.

4. Final simplification 6.3

   \[\leadsto \begin{array}{l} \mathbf{if}\;z \le -4.611098198946549 \cdot 10^{108}:\\ \;\;\;\;-1 \cdot \left(x \cdot y\right)\\ \mathbf{elif}\;z \le 3.1007095082743632 \cdot 10^{95}:\\ \;\;\;\;\left(\frac{\sqrt[3]{y}}{\frac{\frac{\frac{\left|\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}\right|}{\sqrt[3]{z}}}{\sqrt[3]{z}}}{\sqrt[3]{y}}} \cdot \frac{x}{\sqrt{\sqrt{z \cdot z - t \cdot a}}}\right) \cdot \frac{\sqrt[3]{y}}{\frac{\sqrt{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}}{\sqrt[3]{z}}}\\ \mathbf{else}:\\ \;\;\;\;x \cdot y\\ \end{array}\]

Reproduce

herbie shell --seed 2020060 
(FPCore (x y z t a)
  :name "Statistics.Math.RootFinding:ridders from math-functions-0.1.5.2"
  :precision binary64

  :herbie-target
  (if (< z -3.1921305903852764e+46) (- (* y x)) (if (< z 5.976268120920894e+90) (/ (* x z) (/ (sqrt (- (* z z) (* a t))) y)) (* y x)))

  (/ (* (* x y) z) (sqrt (- (* z z) (* t a)))))