Average Error: 6.0 → 0.6
Time: 18.1s
Precision: 64
\[x + \frac{e^{y \cdot \log \left(\frac{y}{z + y}\right)}}{y}\]
\[\begin{array}{l} \mathbf{if}\;y \le -5.528608899369224207630348298103310615572 \cdot 10^{60}:\\ \;\;\;\;x + \frac{e^{-z}}{y}\\ \mathbf{else}:\\ \;\;\;\;x + \frac{\left({\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}\right) \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\\ \end{array}\]
x + \frac{e^{y \cdot \log \left(\frac{y}{z + y}\right)}}{y}
\begin{array}{l}
\mathbf{if}\;y \le -5.528608899369224207630348298103310615572 \cdot 10^{60}:\\
\;\;\;\;x + \frac{e^{-z}}{y}\\

\mathbf{else}:\\
\;\;\;\;x + \frac{\left({\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}\right) \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\\

\end{array}
double f(double x, double y, double z) {
        double r271963 = x;
        double r271964 = y;
        double r271965 = z;
        double r271966 = r271965 + r271964;
        double r271967 = r271964 / r271966;
        double r271968 = log(r271967);
        double r271969 = r271964 * r271968;
        double r271970 = exp(r271969);
        double r271971 = r271970 / r271964;
        double r271972 = r271963 + r271971;
        return r271972;
}

double f(double x, double y, double z) {
        double r271973 = y;
        double r271974 = -5.528608899369224e+60;
        bool r271975 = r271973 <= r271974;
        double r271976 = x;
        double r271977 = z;
        double r271978 = -r271977;
        double r271979 = exp(r271978);
        double r271980 = r271979 / r271973;
        double r271981 = r271976 + r271980;
        double r271982 = cbrt(r271973);
        double r271983 = r271977 + r271973;
        double r271984 = cbrt(r271983);
        double r271985 = r271982 / r271984;
        double r271986 = pow(r271985, r271973);
        double r271987 = r271986 * r271986;
        double r271988 = r271987 * r271986;
        double r271989 = r271988 / r271973;
        double r271990 = r271976 + r271989;
        double r271991 = r271975 ? r271981 : r271990;
        return r271991;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Target

Original6.0
Target1.0
Herbie0.6
\[\begin{array}{l} \mathbf{if}\;\frac{y}{z + y} \lt 7.115415759790762719541517221498726780517 \cdot 10^{-315}:\\ \;\;\;\;x + \frac{e^{\frac{-1}{z}}}{y}\\ \mathbf{else}:\\ \;\;\;\;x + \frac{e^{\log \left({\left(\frac{y}{y + z}\right)}^{y}\right)}}{y}\\ \end{array}\]

Derivation

  1. Split input into 2 regimes
  2. if y < -5.528608899369224e+60

    1. Initial program 2.4

      \[x + \frac{e^{y \cdot \log \left(\frac{y}{z + y}\right)}}{y}\]
    2. Simplified2.4

      \[\leadsto \color{blue}{x + \frac{{\left(\frac{y}{z + y}\right)}^{y}}{y}}\]
    3. Taylor expanded around inf 0.1

      \[\leadsto x + \frac{\color{blue}{e^{-z}}}{y}\]

    if -5.528608899369224e+60 < y

    1. Initial program 6.9

      \[x + \frac{e^{y \cdot \log \left(\frac{y}{z + y}\right)}}{y}\]
    2. Simplified6.9

      \[\leadsto \color{blue}{x + \frac{{\left(\frac{y}{z + y}\right)}^{y}}{y}}\]
    3. Using strategy rm
    4. Applied add-cube-cbrt15.0

      \[\leadsto x + \frac{{\left(\frac{y}{\color{blue}{\left(\sqrt[3]{z + y} \cdot \sqrt[3]{z + y}\right) \cdot \sqrt[3]{z + y}}}\right)}^{y}}{y}\]
    5. Applied add-cube-cbrt6.9

      \[\leadsto x + \frac{{\left(\frac{\color{blue}{\left(\sqrt[3]{y} \cdot \sqrt[3]{y}\right) \cdot \sqrt[3]{y}}}{\left(\sqrt[3]{z + y} \cdot \sqrt[3]{z + y}\right) \cdot \sqrt[3]{z + y}}\right)}^{y}}{y}\]
    6. Applied times-frac6.9

      \[\leadsto x + \frac{{\color{blue}{\left(\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{z + y} \cdot \sqrt[3]{z + y}} \cdot \frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}}^{y}}{y}\]
    7. Applied unpow-prod-down1.9

      \[\leadsto x + \frac{\color{blue}{{\left(\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{z + y} \cdot \sqrt[3]{z + y}}\right)}^{y} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}}{y}\]
    8. Using strategy rm
    9. Applied times-frac1.9

      \[\leadsto x + \frac{{\color{blue}{\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}} \cdot \frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}}^{y} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\]
    10. Applied unpow-prod-down0.7

      \[\leadsto x + \frac{\color{blue}{\left({\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}\right)} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification0.6

    \[\leadsto \begin{array}{l} \mathbf{if}\;y \le -5.528608899369224207630348298103310615572 \cdot 10^{60}:\\ \;\;\;\;x + \frac{e^{-z}}{y}\\ \mathbf{else}:\\ \;\;\;\;x + \frac{\left({\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}\right) \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\\ \end{array}\]

Reproduce

herbie shell --seed 2019326 +o rules:numerics
(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)))