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 + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 + \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(\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}\right) + \left(x1 \cdot x1\right) \cdot \left(-6\right)\right)\right) \cdot \left(x1 \cdot x1 + 1\right)\right)\right) + \left(\sqrt[3]{\left(x1 \cdot x1\right) \cdot \left(x1 + \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} \cdot 3\right)} \cdot \sqrt[3]{\left(x1 \cdot x1\right) \cdot \left(x1 + \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} \cdot 3\right)}\right) \cdot \sqrt[3]{\left(x1 \cdot x1\right) \cdot \left(x1 + \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} \cdot 3\right)}double f(double x1, double x2) {
double r44156 = x1;
double r44157 = 2.0;
double r44158 = r44157 * r44156;
double r44159 = 3.0;
double r44160 = r44159 * r44156;
double r44161 = r44160 * r44156;
double r44162 = x2;
double r44163 = r44157 * r44162;
double r44164 = r44161 + r44163;
double r44165 = r44164 - r44156;
double r44166 = r44156 * r44156;
double r44167 = 1.0;
double r44168 = r44166 + r44167;
double r44169 = r44165 / r44168;
double r44170 = r44158 * r44169;
double r44171 = r44169 - r44159;
double r44172 = r44170 * r44171;
double r44173 = 4.0;
double r44174 = r44173 * r44169;
double r44175 = 6.0;
double r44176 = r44174 - r44175;
double r44177 = r44166 * r44176;
double r44178 = r44172 + r44177;
double r44179 = r44178 * r44168;
double r44180 = r44161 * r44169;
double r44181 = r44179 + r44180;
double r44182 = r44166 * r44156;
double r44183 = r44181 + r44182;
double r44184 = r44183 + r44156;
double r44185 = r44161 - r44163;
double r44186 = r44185 - r44156;
double r44187 = r44186 / r44168;
double r44188 = r44159 * r44187;
double r44189 = r44184 + r44188;
double r44190 = r44156 + r44189;
return r44190;
}
double f(double x1, double x2) {
double r44191 = x1;
double r44192 = 3.0;
double r44193 = r44192 * r44191;
double r44194 = r44193 * r44191;
double r44195 = 2.0;
double r44196 = x2;
double r44197 = r44195 * r44196;
double r44198 = r44194 - r44197;
double r44199 = r44198 - r44191;
double r44200 = r44191 * r44191;
double r44201 = 1.0;
double r44202 = r44200 + r44201;
double r44203 = r44199 / r44202;
double r44204 = r44192 * r44203;
double r44205 = r44191 + r44204;
double r44206 = r44195 * r44191;
double r44207 = r44194 + r44197;
double r44208 = r44207 - r44191;
double r44209 = r44208 / r44202;
double r44210 = r44206 * r44209;
double r44211 = r44209 - r44192;
double r44212 = r44210 * r44211;
double r44213 = 4.0;
double r44214 = r44213 * r44209;
double r44215 = r44200 * r44214;
double r44216 = 6.0;
double r44217 = -r44216;
double r44218 = r44200 * r44217;
double r44219 = r44215 + r44218;
double r44220 = r44212 + r44219;
double r44221 = r44220 * r44202;
double r44222 = r44191 + r44221;
double r44223 = r44205 + r44222;
double r44224 = r44209 * r44192;
double r44225 = r44191 + r44224;
double r44226 = r44200 * r44225;
double r44227 = cbrt(r44226);
double r44228 = r44227 * r44227;
double r44229 = r44228 * r44227;
double r44230 = r44223 + r44229;
return r44230;
}



Bits error versus x1



Bits error versus x2
Results
Initial program 0.5
Simplified0.5
rmApplied sub-neg0.5
Applied distribute-lft-in0.5
rmApplied add-cube-cbrt0.5
Final simplification0.5
herbie shell --seed 2020083
(FPCore (x1 x2)
:name "Rosa's FloatVsDoubleBenchmark"
:precision binary64
(+ 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))))))