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

Average Error: 25.2 → 6.3

Time: 53.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 code(double x, double y, double z, double t, double a) {
	return (((x * y) * z) / sqrt(((z * z) - (t * a))));
}

↓

double code(double x, double y, double z, double t, double a) {
	double temp;
	if ((z <= -4.611098198946549e+108)) {
		temp = (-1.0 * (x * y));
	} else {
		double temp_1;
		if ((z <= 3.100709508274363e+95)) {
			temp_1 = (((cbrt(y) / (((fabs(cbrt(sqrt(((z * z) - (t * a))))) / cbrt(z)) / cbrt(z)) / cbrt(y))) * (x / sqrt(sqrt(((z * z) - (t * a)))))) * (cbrt(y) / (sqrt(cbrt(sqrt(((z * z) - (t * a))))) / cbrt(z))));
		} else {
			temp_1 = (x * y);
		}
		temp = temp_1;
	}
	return temp;
}

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)))))