\frac{1}{3} \cdot \cos^{-1} \left(\frac{3 \cdot \frac{x}{y \cdot 27}}{z \cdot 2} \cdot \sqrt{t}\right)\frac{\sqrt{1}}{\sqrt[3]{3} \cdot \sqrt[3]{3}} \cdot \left(\frac{\sqrt{1}}{\sqrt[3]{3}} \cdot \cos^{-1} \left(\frac{3 \cdot \frac{x}{y \cdot 27}}{z \cdot 2} \cdot \sqrt{t}\right)\right)double f(double x, double y, double z, double t) {
double r534848 = 1.0;
double r534849 = 3.0;
double r534850 = r534848 / r534849;
double r534851 = x;
double r534852 = y;
double r534853 = 27.0;
double r534854 = r534852 * r534853;
double r534855 = r534851 / r534854;
double r534856 = r534849 * r534855;
double r534857 = z;
double r534858 = 2.0;
double r534859 = r534857 * r534858;
double r534860 = r534856 / r534859;
double r534861 = t;
double r534862 = sqrt(r534861);
double r534863 = r534860 * r534862;
double r534864 = acos(r534863);
double r534865 = r534850 * r534864;
return r534865;
}
double f(double x, double y, double z, double t) {
double r534866 = 1.0;
double r534867 = sqrt(r534866);
double r534868 = 3.0;
double r534869 = cbrt(r534868);
double r534870 = r534869 * r534869;
double r534871 = r534867 / r534870;
double r534872 = r534867 / r534869;
double r534873 = x;
double r534874 = y;
double r534875 = 27.0;
double r534876 = r534874 * r534875;
double r534877 = r534873 / r534876;
double r534878 = r534868 * r534877;
double r534879 = z;
double r534880 = 2.0;
double r534881 = r534879 * r534880;
double r534882 = r534878 / r534881;
double r534883 = t;
double r534884 = sqrt(r534883);
double r534885 = r534882 * r534884;
double r534886 = acos(r534885);
double r534887 = r534872 * r534886;
double r534888 = r534871 * r534887;
return r534888;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 1.3 |
|---|---|
| Target | 1.3 |
| Herbie | 0.4 |
Initial program 1.3
rmApplied add-cube-cbrt1.3
Applied add-sqr-sqrt1.3
Applied times-frac0.4
Applied associate-*l*0.4
Final simplification0.4
herbie shell --seed 2019209 +o rules:numerics
(FPCore (x y z t)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, D"
:precision binary64
:herbie-target
(/ (acos (* (/ (/ x 27) (* y z)) (/ (sqrt t) (/ 2 3)))) 3)
(* (/ 1 3) (acos (* (/ (* 3 (/ x (* y 27))) (* z 2)) (sqrt t)))))