x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right)\left(\left(x1 + \left(x1 \cdot \left(x1 \cdot x1\right) + \left(\left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + x2 \cdot 2\right) - x1}{1 + x1 \cdot x1} + \left(1 + x1 \cdot x1\right) \cdot \left(\left(\left(x1 \cdot x1\right) \cdot \left(\sqrt[3]{4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + x2 \cdot 2\right) - x1}{1 + x1 \cdot x1} - 6} \cdot \sqrt[3]{4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + x2 \cdot 2\right) - x1}{1 + x1 \cdot x1} - 6}\right)\right) \cdot \sqrt[3]{4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + x2 \cdot 2\right) - x1}{1 + x1 \cdot x1} - 6} + \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + x2 \cdot 2\right) - x1}{1 + x1 \cdot x1} - 3\right) \cdot \left(\left(x1 \cdot 2\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + x2 \cdot 2\right) - x1}{1 + x1 \cdot x1}\right)\right)\right)\right)\right) + \frac{\left(\left(3 \cdot x1\right) \cdot x1 - x2 \cdot 2\right) - x1}{1 + x1 \cdot x1} \cdot 3\right) + x1double f(double x1, double x2) {
double r34231228 = x1;
double r34231229 = 2.0;
double r34231230 = r34231229 * r34231228;
double r34231231 = 3.0;
double r34231232 = r34231231 * r34231228;
double r34231233 = r34231232 * r34231228;
double r34231234 = x2;
double r34231235 = r34231229 * r34231234;
double r34231236 = r34231233 + r34231235;
double r34231237 = r34231236 - r34231228;
double r34231238 = r34231228 * r34231228;
double r34231239 = 1.0;
double r34231240 = r34231238 + r34231239;
double r34231241 = r34231237 / r34231240;
double r34231242 = r34231230 * r34231241;
double r34231243 = r34231241 - r34231231;
double r34231244 = r34231242 * r34231243;
double r34231245 = 4.0;
double r34231246 = r34231245 * r34231241;
double r34231247 = 6.0;
double r34231248 = r34231246 - r34231247;
double r34231249 = r34231238 * r34231248;
double r34231250 = r34231244 + r34231249;
double r34231251 = r34231250 * r34231240;
double r34231252 = r34231233 * r34231241;
double r34231253 = r34231251 + r34231252;
double r34231254 = r34231238 * r34231228;
double r34231255 = r34231253 + r34231254;
double r34231256 = r34231255 + r34231228;
double r34231257 = r34231233 - r34231235;
double r34231258 = r34231257 - r34231228;
double r34231259 = r34231258 / r34231240;
double r34231260 = r34231231 * r34231259;
double r34231261 = r34231256 + r34231260;
double r34231262 = r34231228 + r34231261;
return r34231262;
}
double f(double x1, double x2) {
double r34231263 = x1;
double r34231264 = r34231263 * r34231263;
double r34231265 = r34231263 * r34231264;
double r34231266 = 3.0;
double r34231267 = r34231266 * r34231263;
double r34231268 = r34231267 * r34231263;
double r34231269 = x2;
double r34231270 = 2.0;
double r34231271 = r34231269 * r34231270;
double r34231272 = r34231268 + r34231271;
double r34231273 = r34231272 - r34231263;
double r34231274 = 1.0;
double r34231275 = r34231274 + r34231264;
double r34231276 = r34231273 / r34231275;
double r34231277 = r34231268 * r34231276;
double r34231278 = 4.0;
double r34231279 = r34231278 * r34231276;
double r34231280 = 6.0;
double r34231281 = r34231279 - r34231280;
double r34231282 = cbrt(r34231281);
double r34231283 = r34231282 * r34231282;
double r34231284 = r34231264 * r34231283;
double r34231285 = r34231284 * r34231282;
double r34231286 = r34231276 - r34231266;
double r34231287 = r34231263 * r34231270;
double r34231288 = r34231287 * r34231276;
double r34231289 = r34231286 * r34231288;
double r34231290 = r34231285 + r34231289;
double r34231291 = r34231275 * r34231290;
double r34231292 = r34231277 + r34231291;
double r34231293 = r34231265 + r34231292;
double r34231294 = r34231263 + r34231293;
double r34231295 = r34231268 - r34231271;
double r34231296 = r34231295 - r34231263;
double r34231297 = r34231296 / r34231275;
double r34231298 = r34231297 * r34231266;
double r34231299 = r34231294 + r34231298;
double r34231300 = r34231299 + r34231263;
return r34231300;
}



Bits error versus x1



Bits error versus x2
Results
Initial program 0.5
rmApplied add-cube-cbrt0.6
Applied associate-*r*0.6
Final simplification0.6
herbie shell --seed 2019119
(FPCore (x1 x2)
:name "Rosa's FloatVsDoubleBenchmark"
(+ x1 (+ (+ (+ (+ (* (+ (* (* (* 2 x1) (/ (- (+ (* (* 3 x1) x1) (* 2 x2)) x1) (+ (* x1 x1) 1))) (- (/ (- (+ (* (* 3 x1) x1) (* 2 x2)) x1) (+ (* x1 x1) 1)) 3)) (* (* x1 x1) (- (* 4 (/ (- (+ (* (* 3 x1) x1) (* 2 x2)) x1) (+ (* x1 x1) 1))) 6))) (+ (* x1 x1) 1)) (* (* (* 3 x1) x1) (/ (- (+ (* (* 3 x1) x1) (* 2 x2)) x1) (+ (* x1 x1) 1)))) (* (* x1 x1) x1)) x1) (* 3 (/ (- (- (* (* 3 x1) x1) (* 2 x2)) x1) (+ (* x1 x1) 1))))))