double code(double x, double y) {
return ((double) (((double) (((double) (((double) (333.75 * ((double) pow(y, 6.0)))) + ((double) (((double) (x * x)) * ((double) (((double) (((double) (((double) (((double) (((double) (((double) (11.0 * x)) * x)) * y)) * y)) - ((double) pow(y, 6.0)))) - ((double) (121.0 * ((double) pow(y, 4.0)))))) - 2.0)))))) + ((double) (5.5 * ((double) pow(y, 8.0)))))) + ((double) (x / ((double) (2.0 * y))))));
}
double code(double x, double y) {
return ((double) (((double) (0.5 * ((double) (x / y)))) - ((double) (x * ((double) (x * 2.0))))));
}
Results
Initial program 58.1
Simplified58.1
Taylor expanded around 0 57.1
Simplified57.1
Final simplification57.1
herbie shell --seed 2020191
(FPCore (x y)
:name "Rump's expression from Stadtherr's award speech"
:precision binary64
:pre (and (== x 77617.0) (== y 33096.0))
(+ (+ (+ (* 333.75 (pow y 6.0)) (* (* x x) (- (- (- (* (* (* (* 11.0 x) x) y) y) (pow y 6.0)) (* 121.0 (pow y 4.0))) 2.0))) (* 5.5 (pow y 8.0))) (/ x (* 2.0 y))))