Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, G

Average Error: 6.0 → 0.6

Time: 21.8s

Precision: 64

Math

\[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{e^{y \cdot \left(2 \cdot \log \left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)\right)} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\\ \end{array}\]

C

double f(double x, double y, double z) {
        double r261271 = x;
        double r261272 = y;
        double r261273 = z;
        double r261274 = r261273 + r261272;
        double r261275 = r261272 / r261274;
        double r261276 = log(r261275);
        double r261277 = r261272 * r261276;
        double r261278 = exp(r261277);
        double r261279 = r261278 / r261272;
        double r261280 = r261271 + r261279;
        return r261280;
}

↓

double f(double x, double y, double z) {
        double r261281 = y;
        double r261282 = -5.528608899369224e+60;
        bool r261283 = r261281 <= r261282;
        double r261284 = x;
        double r261285 = z;
        double r261286 = -r261285;
        double r261287 = exp(r261286);
        double r261288 = r261287 / r261281;
        double r261289 = r261284 + r261288;
        double r261290 = 2.0;
        double r261291 = cbrt(r261281);
        double r261292 = r261285 + r261281;
        double r261293 = cbrt(r261292);
        double r261294 = r261291 / r261293;
        double r261295 = log(r261294);
        double r261296 = r261290 * r261295;
        double r261297 = r261281 * r261296;
        double r261298 = exp(r261297);
        double r261299 = pow(r261294, r261281);
        double r261300 = r261298 * r261299;
        double r261301 = r261300 / r261281;
        double r261302 = r261284 + r261301;
        double r261303 = r261283 ? r261289 : r261302;
        return r261303;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Target

Original	6.0
Target	1.0
Herbie	0.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. Simplified 2.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. Simplified 6.9

      \[\leadsto \color{blue}{x + \frac{{\left(\frac{y}{z + y}\right)}^{y}}{y}}\]
   3. Using strategy rm

   4. Applied add-cube-cbrt 15.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-cbrt 6.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-frac 6.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-down 1.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 add-exp-log 33.0

      \[\leadsto x + \frac{{\left(\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{z + y} \cdot \color{blue}{e^{\log \left(\sqrt[3]{z + y}\right)}}}\right)}^{y} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\]

   10. Applied add-exp-log 33.0

       \[\leadsto x + \frac{{\left(\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\color{blue}{e^{\log \left(\sqrt[3]{z + y}\right)}} \cdot e^{\log \left(\sqrt[3]{z + y}\right)}}\right)}^{y} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\]

   11. Applied prod-exp 33.1

       \[\leadsto x + \frac{{\left(\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\color{blue}{e^{\log \left(\sqrt[3]{z + y}\right) + \log \left(\sqrt[3]{z + y}\right)}}}\right)}^{y} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\]

   12. Applied add-exp-log 33.0

       \[\leadsto x + \frac{{\left(\frac{\sqrt[3]{y} \cdot \color{blue}{e^{\log \left(\sqrt[3]{y}\right)}}}{e^{\log \left(\sqrt[3]{z + y}\right) + \log \left(\sqrt[3]{z + y}\right)}}\right)}^{y} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\]

   13. Applied add-exp-log 27.9

       \[\leadsto x + \frac{{\left(\frac{\color{blue}{e^{\log \left(\sqrt[3]{y}\right)}} \cdot e^{\log \left(\sqrt[3]{y}\right)}}{e^{\log \left(\sqrt[3]{z + y}\right) + \log \left(\sqrt[3]{z + y}\right)}}\right)}^{y} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\]

   14. Applied prod-exp 22.8

       \[\leadsto x + \frac{{\left(\frac{\color{blue}{e^{\log \left(\sqrt[3]{y}\right) + \log \left(\sqrt[3]{y}\right)}}}{e^{\log \left(\sqrt[3]{z + y}\right) + \log \left(\sqrt[3]{z + y}\right)}}\right)}^{y} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\]

   15. Applied div-exp 22.8

       \[\leadsto x + \frac{{\color{blue}{\left(e^{\left(\log \left(\sqrt[3]{y}\right) + \log \left(\sqrt[3]{y}\right)\right) - \left(\log \left(\sqrt[3]{z + y}\right) + \log \left(\sqrt[3]{z + y}\right)\right)}\right)}}^{y} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\]

   16. Applied pow-exp 22.2

       \[\leadsto x + \frac{\color{blue}{e^{\left(\left(\log \left(\sqrt[3]{y}\right) + \log \left(\sqrt[3]{y}\right)\right) - \left(\log \left(\sqrt[3]{z + y}\right) + \log \left(\sqrt[3]{z + y}\right)\right)\right) \cdot y}} \cdot {\left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)}^{y}}{y}\]

   17. Simplified 0.7

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

4. Final simplification 0.6

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

Reproduce

herbie shell --seed 2019326 
(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)))