\left(3 \cdot \sqrt{x}\right) \cdot \left(\left(y + \frac{1}{x \cdot 9}\right) - 1\right)\left(3 \cdot \sqrt{x}\right) \cdot \left(y + \frac{\frac{1}{x}}{9}\right) + \left(\left(\sqrt[3]{3} \cdot \sqrt[3]{3}\right) \cdot \left(\left(\sqrt[3]{\sqrt[3]{3}} \cdot \sqrt[3]{\sqrt[3]{3}}\right) \cdot \left(\sqrt[3]{\sqrt[3]{3}} \cdot \sqrt{x}\right)\right)\right) \cdot \left(-1\right)double f(double x, double y) {
double r394822 = 3.0;
double r394823 = x;
double r394824 = sqrt(r394823);
double r394825 = r394822 * r394824;
double r394826 = y;
double r394827 = 1.0;
double r394828 = 9.0;
double r394829 = r394823 * r394828;
double r394830 = r394827 / r394829;
double r394831 = r394826 + r394830;
double r394832 = r394831 - r394827;
double r394833 = r394825 * r394832;
return r394833;
}
double f(double x, double y) {
double r394834 = 3.0;
double r394835 = x;
double r394836 = sqrt(r394835);
double r394837 = r394834 * r394836;
double r394838 = y;
double r394839 = 1.0;
double r394840 = r394839 / r394835;
double r394841 = 9.0;
double r394842 = r394840 / r394841;
double r394843 = r394838 + r394842;
double r394844 = r394837 * r394843;
double r394845 = cbrt(r394834);
double r394846 = r394845 * r394845;
double r394847 = cbrt(r394845);
double r394848 = r394847 * r394847;
double r394849 = r394847 * r394836;
double r394850 = r394848 * r394849;
double r394851 = r394846 * r394850;
double r394852 = -r394839;
double r394853 = r394851 * r394852;
double r394854 = r394844 + r394853;
return r394854;
}




Bits error versus x




Bits error versus y
Results
| Original | 0.4 |
|---|---|
| Target | 0.4 |
| Herbie | 0.5 |
Initial program 0.4
rmApplied associate-/r*0.4
rmApplied sub-neg0.4
Applied distribute-lft-in0.4
rmApplied add-cube-cbrt0.4
Applied associate-*l*0.5
rmApplied add-cube-cbrt0.5
Applied associate-*l*0.5
Final simplification0.5
herbie shell --seed 2020083 +o rules:numerics
(FPCore (x y)
:name "Numeric.SpecFunctions:incompleteGamma from math-functions-0.1.5.2, B"
:precision binary64
:herbie-target
(* 3 (+ (* y (sqrt x)) (* (- (/ 1 (* x 9)) 1) (sqrt x))))
(* (* 3 (sqrt x)) (- (+ y (/ 1 (* x 9))) 1)))