\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(\sqrt{z \cdot 2} \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)}\right)double f(double x, double y, double z, double t) {
double r514845 = x;
double r514846 = 0.5;
double r514847 = r514845 * r514846;
double r514848 = y;
double r514849 = r514847 - r514848;
double r514850 = z;
double r514851 = 2.0;
double r514852 = r514850 * r514851;
double r514853 = sqrt(r514852);
double r514854 = r514849 * r514853;
double r514855 = t;
double r514856 = r514855 * r514855;
double r514857 = r514856 / r514851;
double r514858 = exp(r514857);
double r514859 = r514854 * r514858;
return r514859;
}
double f(double x, double y, double z, double t) {
double r514860 = x;
double r514861 = 0.5;
double r514862 = r514860 * r514861;
double r514863 = y;
double r514864 = r514862 - r514863;
double r514865 = z;
double r514866 = 2.0;
double r514867 = r514865 * r514866;
double r514868 = sqrt(r514867);
double r514869 = t;
double r514870 = exp(r514869);
double r514871 = cbrt(r514870);
double r514872 = r514871 * r514871;
double r514873 = r514869 / r514866;
double r514874 = pow(r514872, r514873);
double r514875 = r514868 * r514874;
double r514876 = pow(r514871, r514873);
double r514877 = r514875 * r514876;
double r514878 = r514864 * r514877;
return r514878;
}




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
rmApplied add-cube-cbrt0.3
Applied unpow-prod-down0.3
Applied associate-*r*0.3
Final simplification0.3
herbie shell --seed 2019323
(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))))