Average Error: 0.5 → 0.6
Time: 19.5s
Precision: binary64
\[\]
↓
\[\]
↓
double code(double x1, double x2) {
	return ((double) (x1 + ((double) (((double) (((double) (((double) (((double) (((double) (((double) (((double) (((double) (2.0 * x1)) * ((double) (((double) (((double) (((double) (((double) (3.0 * x1)) * x1)) + ((double) (2.0 * x2)))) - x1)) / ((double) (((double) (x1 * x1)) + 1.0)))))) * ((double) (((double) (((double) (((double) (((double) (((double) (3.0 * x1)) * x1)) + ((double) (2.0 * x2)))) - x1)) / ((double) (((double) (x1 * x1)) + 1.0)))) - 3.0)))) + ((double) (((double) (x1 * x1)) * ((double) (((double) (4.0 * ((double) (((double) (((double) (((double) (((double) (3.0 * x1)) * x1)) + ((double) (2.0 * x2)))) - x1)) / ((double) (((double) (x1 * x1)) + 1.0)))))) - 6.0)))))) * ((double) (((double) (x1 * x1)) + 1.0)))) + ((double) (((double) (((double) (3.0 * x1)) * x1)) * ((double) (((double) (((double) (((double) (((double) (3.0 * x1)) * x1)) + ((double) (2.0 * x2)))) - x1)) / ((double) (((double) (x1 * x1)) + 1.0)))))))) + ((double) (((double) (x1 * x1)) * x1)))) + x1)) + ((double) (3.0 * ((double) (((double) (((double) (((double) (((double) (3.0 * x1)) * x1)) - ((double) (2.0 * x2)))) - x1)) / ((double) (((double) (x1 * x1)) + 1.0))))))))));
}

↓
double code(double x1, double x2) {
	return ((double) (x1 + ((double) (((double) (x1 + ((double) (((double) (((double) (((double) (((double) (x1 * x1)) + 1.0)) * ((double) (((double) (((double) (((double) (x1 * 2.0)) * ((double) (((double) (((double) (((double) (x1 * ((double) (x1 * 3.0)))) + ((double) (2.0 * x2)))) - x1)) / ((double) (((double) (x1 * x1)) + 1.0)))))) * ((double) (((double) cbrt(((double) (((double) (((double) (((double) (3.0 * ((double) (x1 * x1)))) + ((double) (((double) (2.0 * x2)) - x1)))) / ((double) (((double) (x1 * x1)) + 1.0)))) - 3.0)))) * ((double) (((double) cbrt(((double) (((double) (((double) (((double) (3.0 * ((double) (x1 * x1)))) + ((double) (((double) (2.0 * x2)) - x1)))) / ((double) (((double) (x1 * x1)) + 1.0)))) - 3.0)))) * ((double) cbrt(((double) (((double) (((double) (((double) (3.0 * ((double) (x1 * x1)))) + ((double) (((double) (2.0 * x2)) - x1)))) / ((double) (((double) (x1 * x1)) + 1.0)))) - 3.0)))))))))) + ((double) (((double) (x1 * x1)) * ((double) (((double) (((double) (((double) (((double) (((double) (x1 * ((double) (x1 * 3.0)))) + ((double) (2.0 * x2)))) - x1)) / ((double) (((double) (x1 * x1)) + 1.0)))) * 4.0)) - 6.0)))))))) + ((double) (((double) (x1 * ((double) (x1 * 3.0)))) * ((double) (((double) (((double) (((double) (x1 * ((double) (x1 * 3.0)))) + ((double) (2.0 * x2)))) - x1)) / ((double) (((double) (x1 * x1)) + 1.0)))))))) + ((double) (x1 * ((double) (x1 * x1)))))))) + ((double) (3.0 * ((double) (((double) (((double) (((double) (x1 * ((double) (x1 * 3.0)))) - ((double) (2.0 * x2)))) - x1)) / ((double) (((double) (x1 * x1)) + 1.0))))))))));
}

Error
Bits error versus x1
Bits error versus x2
Try it out
Your Program's Arguments
Results
Enter valid numbers for all inputs
Derivation
  1. Initial program 0.5
    \[\]
  2. Using strategy rm
  3. Applied add-cube-cbrt0.6
    \[\leadsto \]
  4. Simplified0.6
    \[\leadsto \]
  5. Simplified0.6
    \[\leadsto \]
  6. Final simplification0.6
    \[\leadsto \]
Reproduce
herbie shell --seed 2020190 
(FPCore (x1 x2)
  :name "Rosa's FloatVsDoubleBenchmark"
  :precision binary64
  (+ x1 (+ (+ (+ (+ (* (+ (* (* (* 2.0 x1) (/ (- (+ (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0))) (- (/ (- (+ (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0)) 3.0)) (* (* x1 x1) (- (* 4.0 (/ (- (+ (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0))) 6.0))) (+ (* x1 x1) 1.0)) (* (* (* 3.0 x1) x1) (/ (- (+ (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0)))) (* (* x1 x1) x1)) x1) (* 3.0 (/ (- (- (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0))))))