\frac{1}{3} \cdot \cos^{-1} \left(\frac{3 \cdot \frac{x}{y \cdot 27}}{z \cdot 2} \cdot \sqrt{t}\right)\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{\sqrt[3]{3} \cdot \sqrt[3]{3}} \cdot \left(\frac{\sqrt[3]{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 r741503 = 1.0;
double r741504 = 3.0;
double r741505 = r741503 / r741504;
double r741506 = x;
double r741507 = y;
double r741508 = 27.0;
double r741509 = r741507 * r741508;
double r741510 = r741506 / r741509;
double r741511 = r741504 * r741510;
double r741512 = z;
double r741513 = 2.0;
double r741514 = r741512 * r741513;
double r741515 = r741511 / r741514;
double r741516 = t;
double r741517 = sqrt(r741516);
double r741518 = r741515 * r741517;
double r741519 = acos(r741518);
double r741520 = r741505 * r741519;
return r741520;
}
double f(double x, double y, double z, double t) {
double r741521 = 1.0;
double r741522 = cbrt(r741521);
double r741523 = r741522 * r741522;
double r741524 = 3.0;
double r741525 = cbrt(r741524);
double r741526 = r741525 * r741525;
double r741527 = r741523 / r741526;
double r741528 = r741522 / r741525;
double r741529 = x;
double r741530 = y;
double r741531 = 27.0;
double r741532 = r741530 * r741531;
double r741533 = r741529 / r741532;
double r741534 = r741524 * r741533;
double r741535 = z;
double r741536 = 2.0;
double r741537 = r741535 * r741536;
double r741538 = r741534 / r741537;
double r741539 = t;
double r741540 = sqrt(r741539);
double r741541 = r741538 * r741540;
double r741542 = acos(r741541);
double r741543 = r741528 * r741542;
double r741544 = r741527 * r741543;
return r741544;
}




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