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 r6548743 = x1;
double r6548744 = 2.0;
double r6548745 = r6548744 * r6548743;
double r6548746 = 3.0;
double r6548747 = r6548746 * r6548743;
double r6548748 = r6548747 * r6548743;
double r6548749 = x2;
double r6548750 = r6548744 * r6548749;
double r6548751 = r6548748 + r6548750;
double r6548752 = r6548751 - r6548743;
double r6548753 = r6548743 * r6548743;
double r6548754 = 1.0;
double r6548755 = r6548753 + r6548754;
double r6548756 = r6548752 / r6548755;
double r6548757 = r6548745 * r6548756;
double r6548758 = r6548756 - r6548746;
double r6548759 = r6548757 * r6548758;
double r6548760 = 4.0;
double r6548761 = r6548760 * r6548756;
double r6548762 = 6.0;
double r6548763 = r6548761 - r6548762;
double r6548764 = r6548753 * r6548763;
double r6548765 = r6548759 + r6548764;
double r6548766 = r6548765 * r6548755;
double r6548767 = r6548748 * r6548756;
double r6548768 = r6548766 + r6548767;
double r6548769 = r6548753 * r6548743;
double r6548770 = r6548768 + r6548769;
double r6548771 = r6548770 + r6548743;
double r6548772 = r6548748 - r6548750;
double r6548773 = r6548772 - r6548743;
double r6548774 = r6548773 / r6548755;
double r6548775 = r6548746 * r6548774;
double r6548776 = r6548771 + r6548775;
double r6548777 = r6548743 + r6548776;
return r6548777;
}
double f(double x1, double x2) {
double r6548778 = x1;
double r6548779 = r6548778 * r6548778;
double r6548780 = r6548778 * r6548779;
double r6548781 = 3.0;
double r6548782 = r6548781 * r6548778;
double r6548783 = r6548782 * r6548778;
double r6548784 = x2;
double r6548785 = 2.0;
double r6548786 = r6548784 * r6548785;
double r6548787 = r6548783 + r6548786;
double r6548788 = r6548787 - r6548778;
double r6548789 = 1.0;
double r6548790 = r6548789 + r6548779;
double r6548791 = r6548788 / r6548790;
double r6548792 = r6548783 * r6548791;
double r6548793 = 4.0;
double r6548794 = r6548793 * r6548791;
double r6548795 = 6.0;
double r6548796 = r6548794 - r6548795;
double r6548797 = cbrt(r6548796);
double r6548798 = r6548797 * r6548797;
double r6548799 = r6548779 * r6548798;
double r6548800 = r6548799 * r6548797;
double r6548801 = r6548791 - r6548781;
double r6548802 = r6548778 * r6548785;
double r6548803 = r6548802 * r6548791;
double r6548804 = r6548801 * r6548803;
double r6548805 = r6548800 + r6548804;
double r6548806 = r6548790 * r6548805;
double r6548807 = r6548792 + r6548806;
double r6548808 = r6548780 + r6548807;
double r6548809 = r6548778 + r6548808;
double r6548810 = r6548783 - r6548786;
double r6548811 = r6548810 - r6548778;
double r6548812 = r6548811 / r6548790;
double r6548813 = r6548812 * r6548781;
double r6548814 = r6548809 + r6548813;
double r6548815 = r6548814 + r6548778;
return r6548815;
}



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 2019158
(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))))))