\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z}{\sqrt[3]{t} \cdot \sqrt[3]{t}} \cdot \frac{\sqrt{t + a}}{\sqrt[3]{t}} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}double f(double x, double y, double z, double t, double a, double b, double c) {
double r443265 = x;
double r443266 = y;
double r443267 = 2.0;
double r443268 = z;
double r443269 = t;
double r443270 = a;
double r443271 = r443269 + r443270;
double r443272 = sqrt(r443271);
double r443273 = r443268 * r443272;
double r443274 = r443273 / r443269;
double r443275 = b;
double r443276 = c;
double r443277 = r443275 - r443276;
double r443278 = 5.0;
double r443279 = 6.0;
double r443280 = r443278 / r443279;
double r443281 = r443270 + r443280;
double r443282 = 3.0;
double r443283 = r443269 * r443282;
double r443284 = r443267 / r443283;
double r443285 = r443281 - r443284;
double r443286 = r443277 * r443285;
double r443287 = r443274 - r443286;
double r443288 = r443267 * r443287;
double r443289 = exp(r443288);
double r443290 = r443266 * r443289;
double r443291 = r443265 + r443290;
double r443292 = r443265 / r443291;
return r443292;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r443293 = x;
double r443294 = y;
double r443295 = 2.0;
double r443296 = z;
double r443297 = t;
double r443298 = cbrt(r443297);
double r443299 = r443298 * r443298;
double r443300 = r443296 / r443299;
double r443301 = a;
double r443302 = r443297 + r443301;
double r443303 = sqrt(r443302);
double r443304 = r443303 / r443298;
double r443305 = r443300 * r443304;
double r443306 = b;
double r443307 = c;
double r443308 = r443306 - r443307;
double r443309 = 5.0;
double r443310 = 6.0;
double r443311 = r443309 / r443310;
double r443312 = r443301 + r443311;
double r443313 = 3.0;
double r443314 = r443297 * r443313;
double r443315 = r443295 / r443314;
double r443316 = r443312 - r443315;
double r443317 = r443308 * r443316;
double r443318 = r443305 - r443317;
double r443319 = r443295 * r443318;
double r443320 = exp(r443319);
double r443321 = r443294 * r443320;
double r443322 = r443293 + r443321;
double r443323 = r443293 / r443322;
return r443323;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 3.8 |
|---|---|
| Target | 2.7 |
| Herbie | 2.8 |
Initial program 3.8
rmApplied add-cube-cbrt3.8
Applied times-frac2.8
Final simplification2.8
herbie shell --seed 2020034
(FPCore (x y z t a b c)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, I"
:precision binary64
:herbie-target
(if (< t -2.118326644891581e-50) (/ x (+ x (* y (exp (* 2 (- (+ (* a c) (* 0.8333333333333334 c)) (* a b))))))) (if (< t 5.196588770651547e-123) (/ x (+ x (* y (exp (* 2 (/ (- (* (* z (sqrt (+ t a))) (* (* 3 t) (- a (/ 5 6)))) (* (- (* (+ (/ 5 6) a) (* 3 t)) 2) (* (- a (/ 5 6)) (* (- b c) t)))) (* (* (* t t) 3) (- a (/ 5 6))))))))) (/ x (+ x (* y (exp (* 2 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5 6)) (/ 2 (* t 3))))))))))))
(/ x (+ x (* y (exp (* 2 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5 6)) (/ 2 (* t 3)))))))))))