\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({\left(e^{\sqrt[3]{t} \cdot \sqrt[3]{t}}\right)}^{\left(\sqrt[3]{t} \cdot \frac{t}{2}\right)} \cdot \sqrt{z \cdot 2}\right)double f(double x, double y, double z, double t) {
double r744374 = x;
double r744375 = 0.5;
double r744376 = r744374 * r744375;
double r744377 = y;
double r744378 = r744376 - r744377;
double r744379 = z;
double r744380 = 2.0;
double r744381 = r744379 * r744380;
double r744382 = sqrt(r744381);
double r744383 = r744378 * r744382;
double r744384 = t;
double r744385 = r744384 * r744384;
double r744386 = r744385 / r744380;
double r744387 = exp(r744386);
double r744388 = r744383 * r744387;
return r744388;
}
double f(double x, double y, double z, double t) {
double r744389 = x;
double r744390 = 0.5;
double r744391 = r744389 * r744390;
double r744392 = y;
double r744393 = r744391 - r744392;
double r744394 = t;
double r744395 = cbrt(r744394);
double r744396 = r744395 * r744395;
double r744397 = exp(r744396);
double r744398 = 2.0;
double r744399 = r744394 / r744398;
double r744400 = r744395 * r744399;
double r744401 = pow(r744397, r744400);
double r744402 = z;
double r744403 = r744402 * r744398;
double r744404 = sqrt(r744403);
double r744405 = r744401 * r744404;
double r744406 = r744393 * r744405;
return r744406;
}




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
Simplified0.3
rmApplied add-cube-cbrt0.3
Applied exp-prod0.3
Applied pow-pow0.3
Final simplification0.3
herbie shell --seed 2019304 +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))))