\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot e^{\frac{t \cdot t}{2}}\left(\sqrt{{\left(e^{t}\right)}^{\left(\frac{t}{2}\right)}} \cdot \left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right)\right) \cdot \sqrt{{\left(e^{t}\right)}^{\left(\frac{t}{2}\right)}}double f(double x, double y, double z, double t) {
double r26887457 = x;
double r26887458 = 0.5;
double r26887459 = r26887457 * r26887458;
double r26887460 = y;
double r26887461 = r26887459 - r26887460;
double r26887462 = z;
double r26887463 = 2.0;
double r26887464 = r26887462 * r26887463;
double r26887465 = sqrt(r26887464);
double r26887466 = r26887461 * r26887465;
double r26887467 = t;
double r26887468 = r26887467 * r26887467;
double r26887469 = r26887468 / r26887463;
double r26887470 = exp(r26887469);
double r26887471 = r26887466 * r26887470;
return r26887471;
}
double f(double x, double y, double z, double t) {
double r26887472 = t;
double r26887473 = exp(r26887472);
double r26887474 = 2.0;
double r26887475 = r26887472 / r26887474;
double r26887476 = pow(r26887473, r26887475);
double r26887477 = sqrt(r26887476);
double r26887478 = x;
double r26887479 = 0.5;
double r26887480 = r26887478 * r26887479;
double r26887481 = y;
double r26887482 = r26887480 - r26887481;
double r26887483 = z;
double r26887484 = r26887483 * r26887474;
double r26887485 = sqrt(r26887484);
double r26887486 = r26887482 * r26887485;
double r26887487 = r26887477 * r26887486;
double r26887488 = r26887487 * r26887477;
return r26887488;
}




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-sqr-sqrt0.3
Applied associate-*r*0.3
Final simplification0.3
herbie shell --seed 2019171
(FPCore (x y z t)
:name "Data.Number.Erf:$cinvnormcdf from erf-2.0.0.0, A"
:herbie-target
(* (* (- (* x 0.5) y) (sqrt (* z 2.0))) (pow (exp 1.0) (/ (* t t) 2.0)))
(* (* (- (* x 0.5) y) (sqrt (* z 2.0))) (exp (/ (* t t) 2.0))))