\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 r751233 = 1.0;
double r751234 = 3.0;
double r751235 = r751233 / r751234;
double r751236 = x;
double r751237 = y;
double r751238 = 27.0;
double r751239 = r751237 * r751238;
double r751240 = r751236 / r751239;
double r751241 = r751234 * r751240;
double r751242 = z;
double r751243 = 2.0;
double r751244 = r751242 * r751243;
double r751245 = r751241 / r751244;
double r751246 = t;
double r751247 = sqrt(r751246);
double r751248 = r751245 * r751247;
double r751249 = acos(r751248);
double r751250 = r751235 * r751249;
return r751250;
}
double f(double x, double y, double z, double t) {
double r751251 = 1.0;
double r751252 = 3.0;
double r751253 = cbrt(r751252);
double r751254 = r751253 * r751253;
double r751255 = r751251 / r751254;
double r751256 = 1.0;
double r751257 = r751256 / r751253;
double r751258 = x;
double r751259 = y;
double r751260 = 27.0;
double r751261 = r751259 * r751260;
double r751262 = r751258 / r751261;
double r751263 = r751252 * r751262;
double r751264 = z;
double r751265 = 2.0;
double r751266 = r751264 * r751265;
double r751267 = r751263 / r751266;
double r751268 = t;
double r751269 = sqrt(r751268);
double r751270 = r751267 * r751269;
double r751271 = acos(r751270);
double r751272 = r751257 * r751271;
double r751273 = r751255 * r751272;
return r751273;
}




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 *-un-lft-identity1.3
Applied times-frac0.4
Applied associate-*l*0.4
Final simplification0.4
herbie shell --seed 2020062 +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)))))