\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot e^{\frac{t \cdot t}{2}}\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot \left(\left(\sqrt[3]{{\left(e^{t}\right)}^{\left(\frac{t}{2}\right)}} \cdot \sqrt[3]{{\left(e^{t}\right)}^{\left(\frac{t}{2}\right)}}\right) \cdot \sqrt[3]{{\left(e^{t}\right)}^{\left(\frac{t}{2}\right)}}\right)double f(double x, double y, double z, double t) {
double r845064 = x;
double r845065 = 0.5;
double r845066 = r845064 * r845065;
double r845067 = y;
double r845068 = r845066 - r845067;
double r845069 = z;
double r845070 = 2.0;
double r845071 = r845069 * r845070;
double r845072 = sqrt(r845071);
double r845073 = r845068 * r845072;
double r845074 = t;
double r845075 = r845074 * r845074;
double r845076 = r845075 / r845070;
double r845077 = exp(r845076);
double r845078 = r845073 * r845077;
return r845078;
}
double f(double x, double y, double z, double t) {
double r845079 = x;
double r845080 = 0.5;
double r845081 = r845079 * r845080;
double r845082 = y;
double r845083 = r845081 - r845082;
double r845084 = z;
double r845085 = 2.0;
double r845086 = r845084 * r845085;
double r845087 = sqrt(r845086);
double r845088 = r845083 * r845087;
double r845089 = t;
double r845090 = exp(r845089);
double r845091 = r845089 / r845085;
double r845092 = pow(r845090, r845091);
double r845093 = cbrt(r845092);
double r845094 = r845093 * r845093;
double r845095 = r845094 * r845093;
double r845096 = r845088 * r845095;
return r845096;
}




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
Final simplification0.3
herbie shell --seed 2020062
(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))))