\frac{1}{3} \cdot \cos^{-1} \left(\frac{3 \cdot \frac{x}{y \cdot 27}}{z \cdot 2} \cdot \sqrt{t}\right)\left(\frac{1}{\sqrt[3]{3}} \cdot \cos^{-1} \left(\frac{\frac{x}{27 \cdot y} \cdot 3}{2 \cdot z} \cdot \sqrt{t}\right)\right) \cdot \frac{1}{\sqrt[3]{3} \cdot \sqrt[3]{3}}double f(double x, double y, double z, double t) {
double r27694748 = 1.0;
double r27694749 = 3.0;
double r27694750 = r27694748 / r27694749;
double r27694751 = x;
double r27694752 = y;
double r27694753 = 27.0;
double r27694754 = r27694752 * r27694753;
double r27694755 = r27694751 / r27694754;
double r27694756 = r27694749 * r27694755;
double r27694757 = z;
double r27694758 = 2.0;
double r27694759 = r27694757 * r27694758;
double r27694760 = r27694756 / r27694759;
double r27694761 = t;
double r27694762 = sqrt(r27694761);
double r27694763 = r27694760 * r27694762;
double r27694764 = acos(r27694763);
double r27694765 = r27694750 * r27694764;
return r27694765;
}
double f(double x, double y, double z, double t) {
double r27694766 = 1.0;
double r27694767 = 3.0;
double r27694768 = cbrt(r27694767);
double r27694769 = r27694766 / r27694768;
double r27694770 = x;
double r27694771 = 27.0;
double r27694772 = y;
double r27694773 = r27694771 * r27694772;
double r27694774 = r27694770 / r27694773;
double r27694775 = r27694774 * r27694767;
double r27694776 = 2.0;
double r27694777 = z;
double r27694778 = r27694776 * r27694777;
double r27694779 = r27694775 / r27694778;
double r27694780 = t;
double r27694781 = sqrt(r27694780);
double r27694782 = r27694779 * r27694781;
double r27694783 = acos(r27694782);
double r27694784 = r27694769 * r27694783;
double r27694785 = 1.0;
double r27694786 = r27694768 * r27694768;
double r27694787 = r27694785 / r27694786;
double r27694788 = r27694784 * r27694787;
return r27694788;
}




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.3
Applied associate-*l*0.3
Final simplification0.3
herbie shell --seed 2019172 +o rules:numerics
(FPCore (x y z t)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, D"
:herbie-target
(/ (acos (* (/ (/ x 27.0) (* y z)) (/ (sqrt t) (/ 2.0 3.0)))) 3.0)
(* (/ 1.0 3.0) (acos (* (/ (* 3.0 (/ x (* y 27.0))) (* z 2.0)) (sqrt t)))))