\frac{1}{3} \cdot \cos^{-1} \left(\frac{3 \cdot \frac{x}{y \cdot 27}}{z \cdot 2} \cdot \sqrt{t}\right)\frac{1}{\sqrt[3]{3} \cdot \sqrt[3]{3}} \cdot \left(1 \cdot \frac{\cos^{-1} \left(0.055555555555555552 \cdot \left(\sqrt{t} \cdot \frac{x}{z \cdot y}\right)\right)}{\sqrt[3]{3}}\right)double f(double x, double y, double z, double t) {
double r701866 = 1.0;
double r701867 = 3.0;
double r701868 = r701866 / r701867;
double r701869 = x;
double r701870 = y;
double r701871 = 27.0;
double r701872 = r701870 * r701871;
double r701873 = r701869 / r701872;
double r701874 = r701867 * r701873;
double r701875 = z;
double r701876 = 2.0;
double r701877 = r701875 * r701876;
double r701878 = r701874 / r701877;
double r701879 = t;
double r701880 = sqrt(r701879);
double r701881 = r701878 * r701880;
double r701882 = acos(r701881);
double r701883 = r701868 * r701882;
return r701883;
}
double f(double x, double y, double z, double t) {
double r701884 = 1.0;
double r701885 = 3.0;
double r701886 = cbrt(r701885);
double r701887 = r701886 * r701886;
double r701888 = r701884 / r701887;
double r701889 = 1.0;
double r701890 = 0.05555555555555555;
double r701891 = t;
double r701892 = sqrt(r701891);
double r701893 = x;
double r701894 = z;
double r701895 = y;
double r701896 = r701894 * r701895;
double r701897 = r701893 / r701896;
double r701898 = r701892 * r701897;
double r701899 = r701890 * r701898;
double r701900 = acos(r701899);
double r701901 = r701900 / r701886;
double r701902 = r701889 * r701901;
double r701903 = r701888 * r701902;
return r701903;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 1.3 |
|---|---|
| Target | 1.1 |
| Herbie | 0.2 |
Initial program 1.3
rmApplied add-cube-cbrt1.3
Applied *-un-lft-identity1.3
Applied times-frac0.4
Applied associate-*l*0.4
Taylor expanded around 0 0.2
Final simplification0.2
herbie shell --seed 2020056
(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)))))