\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 r416160 = 3.0;
double r416161 = x;
double r416162 = sqrt(r416161);
double r416163 = r416160 * r416162;
double r416164 = y;
double r416165 = 1.0;
double r416166 = 9.0;
double r416167 = r416161 * r416166;
double r416168 = r416165 / r416167;
double r416169 = r416164 + r416168;
double r416170 = r416169 - r416165;
double r416171 = r416163 * r416170;
return r416171;
}
double f(double x, double y) {
double r416172 = 3.0;
double r416173 = x;
double r416174 = sqrt(r416173);
double r416175 = r416172 * r416174;
double r416176 = y;
double r416177 = 1.0;
double r416178 = r416177 / r416173;
double r416179 = 9.0;
double r416180 = r416178 / r416179;
double r416181 = r416176 + r416180;
double r416182 = r416175 * r416181;
double r416183 = cbrt(r416172);
double r416184 = r416183 * r416183;
double r416185 = cbrt(r416183);
double r416186 = r416185 * r416185;
double r416187 = r416185 * r416174;
double r416188 = r416186 * r416187;
double r416189 = r416184 * r416188;
double r416190 = -r416177;
double r416191 = r416189 * r416190;
double r416192 = r416182 + r416191;
return r416192;
}




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