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

Average Error: 25.2 → 6.3

Time: 52.8s

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 r436311 = x;
        double r436312 = y;
        double r436313 = r436311 * r436312;
        double r436314 = z;
        double r436315 = r436313 * r436314;
        double r436316 = r436314 * r436314;
        double r436317 = t;
        double r436318 = a;
        double r436319 = r436317 * r436318;
        double r436320 = r436316 - r436319;
        double r436321 = sqrt(r436320);
        double r436322 = r436315 / r436321;
        return r436322;
}

↓

double f(double x, double y, double z, double t, double a) {
        double r436323 = z;
        double r436324 = -4.611098198946549e+108;
        bool r436325 = r436323 <= r436324;
        double r436326 = -1.0;
        double r436327 = x;
        double r436328 = y;
        double r436329 = r436327 * r436328;
        double r436330 = r436326 * r436329;
        double r436331 = 3.100709508274363e+95;
        bool r436332 = r436323 <= r436331;
        double r436333 = cbrt(r436328);
        double r436334 = r436323 * r436323;
        double r436335 = t;
        double r436336 = a;
        double r436337 = r436335 * r436336;
        double r436338 = r436334 - r436337;
        double r436339 = sqrt(r436338);
        double r436340 = cbrt(r436339);
        double r436341 = fabs(r436340);
        double r436342 = cbrt(r436323);
        double r436343 = r436341 / r436342;
        double r436344 = r436343 / r436342;
        double r436345 = r436344 / r436333;
        double r436346 = r436333 / r436345;
        double r436347 = sqrt(r436339);
        double r436348 = r436327 / r436347;
        double r436349 = r436346 * r436348;
        double r436350 = sqrt(r436340);
        double r436351 = r436350 / r436342;
        double r436352 = r436333 / r436351;
        double r436353 = r436349 * r436352;
        double r436354 = r436332 ? r436353 : r436329;
        double r436355 = r436325 ? r436330 : r436354;
        return r436355;
}

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 +o rules:numerics
(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)))))