\frac{1}{3} \cdot \cos^{-1} \left(\frac{3 \cdot \frac{x}{y \cdot 27}}{z \cdot 2} \cdot \sqrt{t}\right)\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{\sqrt[3]{3} \cdot \sqrt[3]{3}} \cdot \left(\frac{\sqrt[3]{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 r690825 = 1.0;
double r690826 = 3.0;
double r690827 = r690825 / r690826;
double r690828 = x;
double r690829 = y;
double r690830 = 27.0;
double r690831 = r690829 * r690830;
double r690832 = r690828 / r690831;
double r690833 = r690826 * r690832;
double r690834 = z;
double r690835 = 2.0;
double r690836 = r690834 * r690835;
double r690837 = r690833 / r690836;
double r690838 = t;
double r690839 = sqrt(r690838);
double r690840 = r690837 * r690839;
double r690841 = acos(r690840);
double r690842 = r690827 * r690841;
return r690842;
}
double f(double x, double y, double z, double t) {
double r690843 = 1.0;
double r690844 = cbrt(r690843);
double r690845 = r690844 * r690844;
double r690846 = 3.0;
double r690847 = cbrt(r690846);
double r690848 = r690847 * r690847;
double r690849 = r690845 / r690848;
double r690850 = r690844 / r690847;
double r690851 = x;
double r690852 = y;
double r690853 = 27.0;
double r690854 = r690852 * r690853;
double r690855 = r690851 / r690854;
double r690856 = r690846 * r690855;
double r690857 = z;
double r690858 = 2.0;
double r690859 = r690857 * r690858;
double r690860 = r690856 / r690859;
double r690861 = t;
double r690862 = sqrt(r690861);
double r690863 = r690860 * r690862;
double r690864 = acos(r690863);
double r690865 = r690850 * r690864;
double r690866 = r690849 * r690865;
return r690866;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 1.3 |
|---|---|
| Target | 1.2 |
| Herbie | 0.4 |
Initial program 1.3
rmApplied add-cube-cbrt1.3
Applied add-cube-cbrt1.3
Applied times-frac0.4
Applied associate-*l*0.4
Final simplification0.4
herbie shell --seed 2020083 +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)))))