\frac{1}{3} \cdot \cos^{-1} \left(\frac{3 \cdot \frac{x}{y \cdot 27}}{z \cdot 2} \cdot \sqrt{t}\right)\frac{\sqrt{1}}{\sqrt[3]{3} \cdot \sqrt[3]{3}} \cdot \frac{\cos^{-1} \left(0.05555555555555555247160270937456516548991 \cdot \left(\sqrt{t} \cdot \frac{x}{z \cdot y}\right)\right) \cdot \sqrt{1}}{\sqrt[3]{3}}double f(double x, double y, double z, double t) {
double r458651 = 1.0;
double r458652 = 3.0;
double r458653 = r458651 / r458652;
double r458654 = x;
double r458655 = y;
double r458656 = 27.0;
double r458657 = r458655 * r458656;
double r458658 = r458654 / r458657;
double r458659 = r458652 * r458658;
double r458660 = z;
double r458661 = 2.0;
double r458662 = r458660 * r458661;
double r458663 = r458659 / r458662;
double r458664 = t;
double r458665 = sqrt(r458664);
double r458666 = r458663 * r458665;
double r458667 = acos(r458666);
double r458668 = r458653 * r458667;
return r458668;
}
double f(double x, double y, double z, double t) {
double r458669 = 1.0;
double r458670 = sqrt(r458669);
double r458671 = 3.0;
double r458672 = cbrt(r458671);
double r458673 = r458672 * r458672;
double r458674 = r458670 / r458673;
double r458675 = 0.05555555555555555;
double r458676 = t;
double r458677 = sqrt(r458676);
double r458678 = x;
double r458679 = z;
double r458680 = y;
double r458681 = r458679 * r458680;
double r458682 = r458678 / r458681;
double r458683 = r458677 * r458682;
double r458684 = r458675 * r458683;
double r458685 = acos(r458684);
double r458686 = r458685 * r458670;
double r458687 = r458686 / r458672;
double r458688 = r458674 * r458687;
return r458688;
}




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.2 |
Initial program 1.3
rmApplied add-cube-cbrt1.3
Applied add-sqr-sqrt1.3
Applied times-frac0.4
Applied associate-*l*0.4
Taylor expanded around 0 0.2
Final simplification0.2
herbie shell --seed 2019325 +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)))))