\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot e^{\frac{t \cdot t}{2}}\left(x \cdot 0.5 - y\right) \cdot \left(\sqrt{z \cdot 2} \cdot {\left(e^{t}\right)}^{\left(\frac{t}{2}\right)}\right)double f(double x, double y, double z, double t) {
double r861507 = x;
double r861508 = 0.5;
double r861509 = r861507 * r861508;
double r861510 = y;
double r861511 = r861509 - r861510;
double r861512 = z;
double r861513 = 2.0;
double r861514 = r861512 * r861513;
double r861515 = sqrt(r861514);
double r861516 = r861511 * r861515;
double r861517 = t;
double r861518 = r861517 * r861517;
double r861519 = r861518 / r861513;
double r861520 = exp(r861519);
double r861521 = r861516 * r861520;
return r861521;
}
double f(double x, double y, double z, double t) {
double r861522 = x;
double r861523 = 0.5;
double r861524 = r861522 * r861523;
double r861525 = y;
double r861526 = r861524 - r861525;
double r861527 = z;
double r861528 = 2.0;
double r861529 = r861527 * r861528;
double r861530 = sqrt(r861529);
double r861531 = t;
double r861532 = exp(r861531);
double r861533 = r861531 / r861528;
double r861534 = pow(r861532, r861533);
double r861535 = r861530 * r861534;
double r861536 = r861526 * r861535;
return r861536;
}




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 associate-*l*0.3
Final simplification0.3
herbie shell --seed 2020047 +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))))