x + \frac{e^{y \cdot \log \left(\frac{y}{z + y}\right)}}{y}\begin{array}{l}
\mathbf{if}\;y \le -9.2285838323510326 \cdot 10^{106}:\\
\;\;\;\;x + \frac{e^{-1 \cdot z}}{y}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{e^{\sqrt[3]{{\left(2 \cdot \log \left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)\right)}^{3}} \cdot y}}{\frac{y}{\frac{{\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{\left(\frac{y}{2}\right)} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{\left(\frac{y}{2}\right)}}{1}}}\\
\end{array}double f(double x, double y, double z) {
double r440170 = x;
double r440171 = y;
double r440172 = z;
double r440173 = r440172 + r440171;
double r440174 = r440171 / r440173;
double r440175 = log(r440174);
double r440176 = r440171 * r440175;
double r440177 = exp(r440176);
double r440178 = r440177 / r440171;
double r440179 = r440170 + r440178;
return r440179;
}
double f(double x, double y, double z) {
double r440180 = y;
double r440181 = -9.228583832351033e+106;
bool r440182 = r440180 <= r440181;
double r440183 = x;
double r440184 = -1.0;
double r440185 = z;
double r440186 = r440184 * r440185;
double r440187 = exp(r440186);
double r440188 = r440187 / r440180;
double r440189 = r440183 + r440188;
double r440190 = 2.0;
double r440191 = cbrt(r440180);
double r440192 = r440185 + r440180;
double r440193 = cbrt(r440192);
double r440194 = r440191 / r440193;
double r440195 = log(r440194);
double r440196 = r440190 * r440195;
double r440197 = 3.0;
double r440198 = pow(r440196, r440197);
double r440199 = cbrt(r440198);
double r440200 = r440199 * r440180;
double r440201 = exp(r440200);
double r440202 = r440180 / r440190;
double r440203 = pow(r440194, r440202);
double r440204 = r440203 * r440203;
double r440205 = 1.0;
double r440206 = r440204 / r440205;
double r440207 = r440180 / r440206;
double r440208 = r440201 / r440207;
double r440209 = r440183 + r440208;
double r440210 = r440182 ? r440189 : r440209;
return r440210;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.2 |
|---|---|
| Target | 1.1 |
| Herbie | 0.7 |
if y < -9.228583832351033e+106Initial program 2.4
Taylor expanded around inf 0.0
if -9.228583832351033e+106 < y Initial program 7.0
rmApplied add-cube-cbrt16.8
Applied add-cube-cbrt7.0
Applied times-frac7.0
Applied log-prod2.0
Applied distribute-rgt-in2.0
Applied exp-sum2.0
Applied associate-/l*2.0
Simplified2.0
rmApplied add-cbrt-cube2.0
Simplified0.8
rmApplied sqr-pow0.8
Final simplification0.7
herbie shell --seed 2020049
(FPCore (x y z)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, G"
:precision binary64
:herbie-target
(if (< (/ y (+ z y)) 7.1154157597908e-315) (+ x (/ (exp (/ -1 z)) y)) (+ x (/ (exp (log (pow (/ y (+ y z)) y))) y)))
(+ x (/ (exp (* y (log (/ y (+ z y))))) y)))