\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 r764882 = x;
double r764883 = 0.5;
double r764884 = r764882 * r764883;
double r764885 = y;
double r764886 = r764884 - r764885;
double r764887 = z;
double r764888 = 2.0;
double r764889 = r764887 * r764888;
double r764890 = sqrt(r764889);
double r764891 = r764886 * r764890;
double r764892 = t;
double r764893 = r764892 * r764892;
double r764894 = r764893 / r764888;
double r764895 = exp(r764894);
double r764896 = r764891 * r764895;
return r764896;
}
double f(double x, double y, double z, double t) {
double r764897 = x;
double r764898 = 0.5;
double r764899 = r764897 * r764898;
double r764900 = y;
double r764901 = r764899 - r764900;
double r764902 = z;
double r764903 = 2.0;
double r764904 = r764902 * r764903;
double r764905 = sqrt(r764904);
double r764906 = r764901 * r764905;
double r764907 = t;
double r764908 = exp(r764907);
double r764909 = cbrt(r764908);
double r764910 = r764909 * r764909;
double r764911 = r764907 / r764903;
double r764912 = pow(r764910, r764911);
double r764913 = r764906 * r764912;
double r764914 = pow(r764909, r764911);
double r764915 = r764913 * r764914;
return r764915;
}




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))))