\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(\cos^{-1} \left(\left(\frac{3}{z} \cdot \frac{x}{2 \cdot \left(y \cdot 27\right)}\right) \cdot \sqrt{t}\right) \cdot \frac{1}{\sqrt[3]{3}}\right)(FPCore (x y z t) :precision binary64 (* (/ 1.0 3.0) (acos (* (/ (* 3.0 (/ x (* y 27.0))) (* z 2.0)) (sqrt t)))))
(FPCore (x y z t) :precision binary64 (* (/ 1.0 (* (cbrt 3.0) (cbrt 3.0))) (* (acos (* (* (/ 3.0 z) (/ x (* 2.0 (* y 27.0)))) (sqrt t))) (/ 1.0 (cbrt 3.0)))))
double code(double x, double y, double z, double t) {
return ((double) ((1.0 / 3.0) * ((double) acos(((double) ((((double) (3.0 * (x / ((double) (y * 27.0))))) / ((double) (z * 2.0))) * ((double) sqrt(t))))))));
}
double code(double x, double y, double z, double t) {
return ((double) ((1.0 / ((double) (((double) cbrt(3.0)) * ((double) cbrt(3.0))))) * ((double) (((double) acos(((double) (((double) ((3.0 / z) * (x / ((double) (2.0 * ((double) (y * 27.0))))))) * ((double) sqrt(t)))))) * (1.0 / ((double) cbrt(3.0)))))));
}




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 Error: 1.3 bits
rmApplied add-cube-cbrtError: 1.3 bits
Applied *-un-lft-identityError: 1.3 bits
Applied times-fracError: 0.4 bits
Applied associate-*l*Error: 0.3 bits
SimplifiedError: 0.3 bits
Final simplificationError: 0.3 bits
herbie shell --seed 2020204
(FPCore (x y z t)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, D"
:precision binary64
: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)))))