\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(\frac{1}{\sqrt[3]{3}} \cdot \cos^{-1} \left(\frac{3 \cdot x}{\left(z \cdot 2\right) \cdot \left(y \cdot 27\right)} \cdot \sqrt{t}\right)\right)double f(double x, double y, double z, double t) {
double r770828 = 1.0;
double r770829 = 3.0;
double r770830 = r770828 / r770829;
double r770831 = x;
double r770832 = y;
double r770833 = 27.0;
double r770834 = r770832 * r770833;
double r770835 = r770831 / r770834;
double r770836 = r770829 * r770835;
double r770837 = z;
double r770838 = 2.0;
double r770839 = r770837 * r770838;
double r770840 = r770836 / r770839;
double r770841 = t;
double r770842 = sqrt(r770841);
double r770843 = r770840 * r770842;
double r770844 = acos(r770843);
double r770845 = r770830 * r770844;
return r770845;
}
double f(double x, double y, double z, double t) {
double r770846 = 1.0;
double r770847 = 3.0;
double r770848 = cbrt(r770847);
double r770849 = r770848 * r770848;
double r770850 = r770846 / r770849;
double r770851 = 1.0;
double r770852 = r770851 / r770848;
double r770853 = x;
double r770854 = r770847 * r770853;
double r770855 = z;
double r770856 = 2.0;
double r770857 = r770855 * r770856;
double r770858 = y;
double r770859 = 27.0;
double r770860 = r770858 * r770859;
double r770861 = r770857 * r770860;
double r770862 = r770854 / r770861;
double r770863 = t;
double r770864 = sqrt(r770863);
double r770865 = r770862 * r770864;
double r770866 = acos(r770865);
double r770867 = r770852 * r770866;
double r770868 = r770850 * r770867;
return r770868;
}




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.3 |
Initial program 1.3
rmApplied add-cube-cbrt1.3
Applied *-un-lft-identity1.3
Applied times-frac0.4
Applied associate-*l*0.4
rmApplied associate-*r/0.4
Applied associate-/l/0.3
Final simplification0.3
herbie shell --seed 2020024
(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)))))