\left(3 \cdot \sqrt{x}\right) \cdot \left(\left(y + \frac{1}{x \cdot 9}\right) - 1\right)\sqrt{x} \cdot \left(\left(\left(\frac{1}{x \cdot 9} + y\right) - 1\right) \cdot 3\right) + \left(\sqrt{x} \cdot \left(1 \cdot 0\right)\right) \cdot 3double f(double x, double y) {
double r382264 = 3.0;
double r382265 = x;
double r382266 = sqrt(r382265);
double r382267 = r382264 * r382266;
double r382268 = y;
double r382269 = 1.0;
double r382270 = 9.0;
double r382271 = r382265 * r382270;
double r382272 = r382269 / r382271;
double r382273 = r382268 + r382272;
double r382274 = r382273 - r382269;
double r382275 = r382267 * r382274;
return r382275;
}
double f(double x, double y) {
double r382276 = x;
double r382277 = sqrt(r382276);
double r382278 = 1.0;
double r382279 = 9.0;
double r382280 = r382276 * r382279;
double r382281 = r382278 / r382280;
double r382282 = y;
double r382283 = r382281 + r382282;
double r382284 = r382283 - r382278;
double r382285 = 3.0;
double r382286 = r382284 * r382285;
double r382287 = r382277 * r382286;
double r382288 = 0.0;
double r382289 = r382278 * r382288;
double r382290 = r382277 * r382289;
double r382291 = r382290 * r382285;
double r382292 = r382287 + r382291;
return r382292;
}




Bits error versus x




Bits error versus y
Results
| Original | 0.4 |
|---|---|
| Target | 0.4 |
| Herbie | 0.4 |
Initial program 0.4
rmApplied add-cube-cbrt0.4
Applied add-sqr-sqrt16.0
Applied prod-diff16.0
Applied distribute-lft-in16.0
Simplified0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2020045 +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)))