\left(1 - \frac{1}{x \cdot 9}\right) - \frac{y}{3 \cdot \sqrt{x}}\left(1 + \frac{\sqrt[3]{-1} \cdot \sqrt[3]{-1}}{\frac{9}{\sqrt[3]{-1}} \cdot x}\right) - \frac{\frac{1}{\sqrt[3]{3} \cdot \sqrt[3]{3}}}{\sqrt{x} \cdot \sqrt[3]{3}} \cdot ydouble f(double x, double y) {
double r528279 = 1.0;
double r528280 = x;
double r528281 = 9.0;
double r528282 = r528280 * r528281;
double r528283 = r528279 / r528282;
double r528284 = r528279 - r528283;
double r528285 = y;
double r528286 = 3.0;
double r528287 = sqrt(r528280);
double r528288 = r528286 * r528287;
double r528289 = r528285 / r528288;
double r528290 = r528284 - r528289;
return r528290;
}
double f(double x, double y) {
double r528291 = 1.0;
double r528292 = -r528291;
double r528293 = cbrt(r528292);
double r528294 = r528293 * r528293;
double r528295 = 9.0;
double r528296 = r528295 / r528293;
double r528297 = x;
double r528298 = r528296 * r528297;
double r528299 = r528294 / r528298;
double r528300 = r528291 + r528299;
double r528301 = 1.0;
double r528302 = 3.0;
double r528303 = cbrt(r528302);
double r528304 = r528303 * r528303;
double r528305 = r528301 / r528304;
double r528306 = sqrt(r528297);
double r528307 = r528306 * r528303;
double r528308 = r528305 / r528307;
double r528309 = y;
double r528310 = r528308 * r528309;
double r528311 = r528300 - r528310;
return r528311;
}




Bits error versus x




Bits error versus y
Results
| Original | 0.2 |
|---|---|
| Target | 0.3 |
| Herbie | 0.2 |
Initial program 0.2
Simplified0.2
rmApplied sub-neg0.2
Simplified0.3
rmApplied div-inv0.3
Simplified0.3
rmApplied add-cube-cbrt0.3
Applied add-sqr-sqrt0.3
Applied times-frac0.4
Applied associate-/l*0.3
Simplified0.3
rmApplied *-un-lft-identity0.3
Applied add-cube-cbrt0.3
Applied times-frac0.3
Applied associate-/l*0.3
Simplified0.2
Final simplification0.2
herbie shell --seed 2019174
(FPCore (x y)
:name "Numeric.SpecFunctions:invIncompleteGamma from math-functions-0.1.5.2, D"
:herbie-target
(- (- 1.0 (/ (/ 1.0 x) 9.0)) (/ y (* 3.0 (sqrt x))))
(- (- 1.0 (/ 1.0 (* x 9.0))) (/ y (* 3.0 (sqrt x)))))