\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 \frac{x}{y \cdot 27}}{z \cdot 2} \cdot \sqrt{t}\right)\right)double f(double x, double y, double z, double t) {
double r756660 = 1.0;
double r756661 = 3.0;
double r756662 = r756660 / r756661;
double r756663 = x;
double r756664 = y;
double r756665 = 27.0;
double r756666 = r756664 * r756665;
double r756667 = r756663 / r756666;
double r756668 = r756661 * r756667;
double r756669 = z;
double r756670 = 2.0;
double r756671 = r756669 * r756670;
double r756672 = r756668 / r756671;
double r756673 = t;
double r756674 = sqrt(r756673);
double r756675 = r756672 * r756674;
double r756676 = acos(r756675);
double r756677 = r756662 * r756676;
return r756677;
}
double f(double x, double y, double z, double t) {
double r756678 = 1.0;
double r756679 = 3.0;
double r756680 = cbrt(r756679);
double r756681 = r756680 * r756680;
double r756682 = r756678 / r756681;
double r756683 = 1.0;
double r756684 = r756683 / r756680;
double r756685 = x;
double r756686 = y;
double r756687 = 27.0;
double r756688 = r756686 * r756687;
double r756689 = r756685 / r756688;
double r756690 = r756679 * r756689;
double r756691 = z;
double r756692 = 2.0;
double r756693 = r756691 * r756692;
double r756694 = r756690 / r756693;
double r756695 = t;
double r756696 = sqrt(r756695);
double r756697 = r756694 * r756696;
double r756698 = acos(r756697);
double r756699 = r756684 * r756698;
double r756700 = r756682 * r756699;
return r756700;
}




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 2020003
(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)))))