\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot e^{\frac{t \cdot t}{2}}\left(\left(\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z}\right) \cdot \left(\sqrt[3]{\sqrt{2}} \cdot \sqrt[3]{\sqrt{2}}\right)\right) \cdot \sqrt[3]{\sqrt{2}}\right) \cdot e^{\frac{t \cdot t}{2}}double f(double x, double y, double z, double t) {
double r839096 = x;
double r839097 = 0.5;
double r839098 = r839096 * r839097;
double r839099 = y;
double r839100 = r839098 - r839099;
double r839101 = z;
double r839102 = 2.0;
double r839103 = r839101 * r839102;
double r839104 = sqrt(r839103);
double r839105 = r839100 * r839104;
double r839106 = t;
double r839107 = r839106 * r839106;
double r839108 = r839107 / r839102;
double r839109 = exp(r839108);
double r839110 = r839105 * r839109;
return r839110;
}
double f(double x, double y, double z, double t) {
double r839111 = x;
double r839112 = 0.5;
double r839113 = r839111 * r839112;
double r839114 = y;
double r839115 = r839113 - r839114;
double r839116 = z;
double r839117 = sqrt(r839116);
double r839118 = r839115 * r839117;
double r839119 = 2.0;
double r839120 = sqrt(r839119);
double r839121 = cbrt(r839120);
double r839122 = r839121 * r839121;
double r839123 = r839118 * r839122;
double r839124 = r839123 * r839121;
double r839125 = t;
double r839126 = r839125 * r839125;
double r839127 = r839126 / r839119;
double r839128 = exp(r839127);
double r839129 = r839124 * r839128;
return r839129;
}




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.5 |
Initial program 0.3
rmApplied sqrt-prod0.5
Applied associate-*r*0.5
rmApplied add-cube-cbrt0.5
Applied associate-*r*0.5
Final simplification0.5
herbie shell --seed 2020002 +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))))