\left(1 - \frac{1}{x \cdot 9}\right) - \frac{y}{3 \cdot \sqrt{x}}\left(1 - \frac{\frac{1}{x}}{9}\right) - \frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{\sqrt[3]{3} \cdot \sqrt[3]{3}} \cdot \frac{\frac{y}{\sqrt{x}}}{\sqrt[3]{3}}double f(double x, double y) {
double r422473 = 1.0;
double r422474 = x;
double r422475 = 9.0;
double r422476 = r422474 * r422475;
double r422477 = r422473 / r422476;
double r422478 = r422473 - r422477;
double r422479 = y;
double r422480 = 3.0;
double r422481 = sqrt(r422474);
double r422482 = r422480 * r422481;
double r422483 = r422479 / r422482;
double r422484 = r422478 - r422483;
return r422484;
}
double f(double x, double y) {
double r422485 = 1.0;
double r422486 = x;
double r422487 = r422485 / r422486;
double r422488 = 9.0;
double r422489 = r422487 / r422488;
double r422490 = r422485 - r422489;
double r422491 = 1.0;
double r422492 = cbrt(r422491);
double r422493 = r422492 * r422492;
double r422494 = 3.0;
double r422495 = cbrt(r422494);
double r422496 = r422495 * r422495;
double r422497 = r422493 / r422496;
double r422498 = y;
double r422499 = sqrt(r422486);
double r422500 = r422498 / r422499;
double r422501 = r422500 / r422495;
double r422502 = r422497 * r422501;
double r422503 = r422490 - r422502;
return r422503;
}




Bits error versus x




Bits error versus y
Results
| Original | 0.2 |
|---|---|
| Target | 0.2 |
| Herbie | 0.3 |
Initial program 0.2
rmApplied associate-/r*0.2
rmApplied *-un-lft-identity0.2
Applied times-frac0.2
rmApplied add-cube-cbrt0.2
Applied add-cube-cbrt0.2
Applied times-frac0.4
Applied associate-*l*0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2020064
(FPCore (x y)
:name "Numeric.SpecFunctions:invIncompleteGamma from math-functions-0.1.5.2, D"
:precision binary64
:herbie-target
(- (- 1 (/ (/ 1 x) 9)) (/ y (* 3 (sqrt x))))
(- (- 1 (/ 1 (* x 9))) (/ y (* 3 (sqrt x)))))