1 - \sqrt{0.5 \cdot \left(1 + \frac{1}{\mathsf{hypot}\left(1, x\right)}\right)}\frac{0.5 - \frac{0.5}{\mathsf{hypot}\left(1, x\right)}}{1 + \sqrt[3]{\sqrt[3]{{\left({\left(\sqrt{0.5 \cdot \left(1 + \frac{1}{\mathsf{hypot}\left(1, x\right)}\right)}\right)}^{3}\right)}^{3}}}}double f(double x) {
double r173748 = 1.0;
double r173749 = 0.5;
double r173750 = x;
double r173751 = hypot(r173748, r173750);
double r173752 = r173748 / r173751;
double r173753 = r173748 + r173752;
double r173754 = r173749 * r173753;
double r173755 = sqrt(r173754);
double r173756 = r173748 - r173755;
return r173756;
}
double f(double x) {
double r173757 = 0.5;
double r173758 = 1.0;
double r173759 = x;
double r173760 = hypot(r173758, r173759);
double r173761 = r173757 / r173760;
double r173762 = r173757 - r173761;
double r173763 = r173758 / r173760;
double r173764 = r173758 + r173763;
double r173765 = r173757 * r173764;
double r173766 = sqrt(r173765);
double r173767 = 3.0;
double r173768 = pow(r173766, r173767);
double r173769 = pow(r173768, r173767);
double r173770 = cbrt(r173769);
double r173771 = cbrt(r173770);
double r173772 = r173758 + r173771;
double r173773 = r173762 / r173772;
return r173773;
}



Bits error versus x
Results
Initial program 15.0
rmApplied flip--15.0
Simplified14.5
rmApplied add-cbrt-cube14.5
Simplified14.5
Taylor expanded around 0 14.5
Simplified14.5
rmApplied add-cbrt-cube14.5
Simplified14.5
Final simplification14.5
herbie shell --seed 2020045
(FPCore (x)
:name "Given's Rotation SVD example, simplified"
:precision binary64
(- 1 (sqrt (* 0.5 (+ 1 (/ 1 (hypot 1 x)))))))