\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot e^{\frac{t \cdot t}{2}}\left(\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot {\left(\sqrt[3]{e^{t}} \cdot \sqrt[3]{e^{t}}\right)}^{\left(\frac{t}{2}\right)}\right) \cdot {\left(\sqrt[3]{e^{t}}\right)}^{\left(\frac{t}{2}\right)}double f(double x, double y, double z, double t) {
double r757227 = x;
double r757228 = 0.5;
double r757229 = r757227 * r757228;
double r757230 = y;
double r757231 = r757229 - r757230;
double r757232 = z;
double r757233 = 2.0;
double r757234 = r757232 * r757233;
double r757235 = sqrt(r757234);
double r757236 = r757231 * r757235;
double r757237 = t;
double r757238 = r757237 * r757237;
double r757239 = r757238 / r757233;
double r757240 = exp(r757239);
double r757241 = r757236 * r757240;
return r757241;
}
double f(double x, double y, double z, double t) {
double r757242 = x;
double r757243 = 0.5;
double r757244 = r757242 * r757243;
double r757245 = y;
double r757246 = r757244 - r757245;
double r757247 = z;
double r757248 = 2.0;
double r757249 = r757247 * r757248;
double r757250 = sqrt(r757249);
double r757251 = r757246 * r757250;
double r757252 = t;
double r757253 = exp(r757252);
double r757254 = cbrt(r757253);
double r757255 = r757254 * r757254;
double r757256 = r757252 / r757248;
double r757257 = pow(r757255, r757256);
double r757258 = r757251 * r757257;
double r757259 = pow(r757254, r757256);
double r757260 = r757258 * r757259;
return r757260;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 0.3 |
|---|---|
| Target | 0.3 |
| Herbie | 0.3 |
Initial program 0.3
rmApplied *-un-lft-identity0.3
Applied times-frac0.3
Applied exp-prod0.3
Simplified0.3
rmApplied add-cube-cbrt0.3
Applied unpow-prod-down0.3
Applied associate-*r*0.3
Final simplification0.3
herbie shell --seed 2020001 +o rules:numerics
(FPCore (x y z t)
:name "Data.Number.Erf:$cinvnormcdf from erf-2.0.0.0, A"
:precision binary64
:herbie-target
(* (* (- (* x 0.5) y) (sqrt (* z 2))) (pow (exp 1) (/ (* t t) 2)))
(* (* (- (* x 0.5) y) (sqrt (* z 2))) (exp (/ (* t t) 2))))