Average Error: 4.6 → 2.4
Time: 8.7m
Precision: 64
Internal Precision: 384
\[\frac{x}{x + y \cdot e^{2.0 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right)\right)}}\]
↓
\[\begin{array}{l} \mathbf{if}\;t \le -3.546490195958238 \cdot 10^{-104}:\\ \;\;\;\;\frac{x}{x + y \cdot {\left(e^{2.0}\right)}^{\left(\left(\frac{0.6666666666666666}{t} \cdot \left(b - c\right) - b \cdot 0.8333333333333334\right) + \frac{\sqrt{t + a}}{\frac{t}{z}}\right)}}\\ \mathbf{if}\;t \le 3.620412988201743 \cdot 10^{-245}:\\ \;\;\;\;\frac{x}{x + y \cdot e^{2.0 \cdot \frac{\left(\left(z \cdot 3.0\right) \cdot \sqrt{t + a}\right) \cdot \left(\left(a - \frac{5.0}{6.0}\right) \cdot t\right) - \left(\left(b - c\right) \cdot \left(\left(a - \frac{5.0}{6.0}\right) \cdot t\right)\right) \cdot \left(\left(3.0 \cdot t\right) \cdot \left(\frac{5.0}{6.0} + a\right) - 2.0\right)}{t \cdot \left(\left(a - \frac{5.0}{6.0}\right) \cdot \left(t \cdot 3.0\right)\right)}}}\\ \mathbf{else}:\\ \;\;\;\;\frac{x}{x + y \cdot e^{2.0 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right)\right)}}\\ \end{array}\]

Error

Bits error versus x

Bits error versus y

Bits error versus z

Bits error versus t

Bits error versus a

Bits error versus b

Bits error versus c

Target

Original4.6
Target4.0
Herbie2.4
\[\begin{array}{l} \mathbf{if}\;t \lt -2.118326644891581 \cdot 10^{-50}:\\ \;\;\;\;\frac{x}{x + y \cdot e^{2.0 \cdot \left(\left(a \cdot c + 0.8333333333333334 \cdot c\right) - a \cdot b\right)}}\\ \mathbf{if}\;t \lt 5.196588770651547 \cdot 10^{-123}:\\ \;\;\;\;\frac{x}{x + y \cdot e^{2.0 \cdot \frac{\left(z \cdot \sqrt{t + a}\right) \cdot \left(\left(3.0 \cdot t\right) \cdot \left(a - \frac{5.0}{6.0}\right)\right) - \left(\left(\frac{5.0}{6.0} + a\right) \cdot \left(3.0 \cdot t\right) - 2.0\right) \cdot \left(\left(a - \frac{5.0}{6.0}\right) \cdot \left(\left(b - c\right) \cdot t\right)\right)}{\left(\left(t \cdot t\right) \cdot 3.0\right) \cdot \left(a - \frac{5.0}{6.0}\right)}}}\\ \mathbf{else}:\\ \;\;\;\;\frac{x}{x + y \cdot e^{2.0 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right)\right)}}\\ \end{array}\]

Derivation

  1. Split input into 3 regimes

  2. if t < -3.546490195958238e-104

    1. Initial program 10.8

      \[\frac{x}{x + y \cdot e^{2.0 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right)\right)}}\]

    2. Taylor expanded around 0 1.0

      \[\leadsto \frac{x}{x + y \cdot e^{2.0 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \color{blue}{\left(\left(0.8333333333333334 \cdot b + 0.6666666666666666 \cdot \frac{c}{t}\right) - 0.6666666666666666 \cdot \frac{b}{t}\right)}\right)}}\]

    3. Applied simplify0.4

      \[\leadsto \color{blue}{\frac{x}{x + y \cdot {\left(e^{2.0}\right)}^{\left(\left(\frac{0.6666666666666666}{t} \cdot \left(b - c\right) - b \cdot 0.8333333333333334\right) + \frac{\sqrt{t + a}}{\frac{t}{z}}\right)}}}\]

    if -3.546490195958238e-104 < t < 3.620412988201743e-245

    1. Initial program 8.4

      \[\frac{x}{x + y \cdot e^{2.0 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right)\right)}}\]
    2. Using strategy rm

    3. Applied flip-+11.7

      \[\leadsto \frac{x}{x + y \cdot e^{2.0 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\color{blue}{\frac{a \cdot a - \frac{5.0}{6.0} \cdot \frac{5.0}{6.0}}{a - \frac{5.0}{6.0}}} - \frac{2.0}{t \cdot 3.0}\right)\right)}}\]

    4. Applied frac-sub11.7

      \[\leadsto \frac{x}{x + y \cdot e^{2.0 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \color{blue}{\frac{\left(a \cdot a - \frac{5.0}{6.0} \cdot \frac{5.0}{6.0}\right) \cdot \left(t \cdot 3.0\right) - \left(a - \frac{5.0}{6.0}\right) \cdot 2.0}{\left(a - \frac{5.0}{6.0}\right) \cdot \left(t \cdot 3.0\right)}}\right)}}\]

    5. Applied associate-*r/11.7

      \[\leadsto \frac{x}{x + y \cdot e^{2.0 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \color{blue}{\frac{\left(b - c\right) \cdot \left(\left(a \cdot a - \frac{5.0}{6.0} \cdot \frac{5.0}{6.0}\right) \cdot \left(t \cdot 3.0\right) - \left(a - \frac{5.0}{6.0}\right) \cdot 2.0\right)}{\left(a - \frac{5.0}{6.0}\right) \cdot \left(t \cdot 3.0\right)}}\right)}}\]

    6. Applied frac-sub8.9

      \[\leadsto \frac{x}{x + y \cdot e^{2.0 \cdot \color{blue}{\frac{\left(z \cdot \sqrt{t + a}\right) \cdot \left(\left(a - \frac{5.0}{6.0}\right) \cdot \left(t \cdot 3.0\right)\right) - t \cdot \left(\left(b - c\right) \cdot \left(\left(a \cdot a - \frac{5.0}{6.0} \cdot \frac{5.0}{6.0}\right) \cdot \left(t \cdot 3.0\right) - \left(a - \frac{5.0}{6.0}\right) \cdot 2.0\right)\right)}{t \cdot \left(\left(a - \frac{5.0}{6.0}\right) \cdot \left(t \cdot 3.0\right)\right)}}}}\]

    7. Applied simplify3.2

      \[\leadsto \frac{x}{x + y \cdot e^{2.0 \cdot \frac{\color{blue}{\left(\left(z \cdot 3.0\right) \cdot \sqrt{t + a}\right) \cdot \left(\left(a - \frac{5.0}{6.0}\right) \cdot t\right) - \left(\left(b - c\right) \cdot \left(\left(a - \frac{5.0}{6.0}\right) \cdot t\right)\right) \cdot \left(\left(3.0 \cdot t\right) \cdot \left(\frac{5.0}{6.0} + a\right) - 2.0\right)}}{t \cdot \left(\left(a - \frac{5.0}{6.0}\right) \cdot \left(t \cdot 3.0\right)\right)}}}\]

    if 3.620412988201743e-245 < t

    1. Initial program 2.5

      \[\frac{x}{x + y \cdot e^{2.0 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5.0}{6.0}\right) - \frac{2.0}{t \cdot 3.0}\right)\right)}}\]
  3. Recombined 3 regimes into one program.
  4. Removed slow pow expressions.

Runtime

Time bar (total: 8.7m)Debug log

herbie shell --seed '#(1567391828 2030694642 2833800258 828025724 3004380912 3532991858)' +o setup:early-exit
(FPCore (x y z t a b c)
  :name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, I"

  :herbie-target
  (if (< t -2.118326644891581e-50) (/ x (+ x (* y (exp (* 2.0 (- (+ (* a c) (* 0.8333333333333334 c)) (* a b))))))) (if (< t 5.196588770651547e-123) (/ x (+ x (* y (exp (* 2.0 (/ (- (* (* z (sqrt (+ t a))) (* (* 3.0 t) (- a (/ 5.0 6.0)))) (* (- (* (+ (/ 5.0 6.0) a) (* 3.0 t)) 2.0) (* (- a (/ 5.0 6.0)) (* (- b c) t)))) (* (* (* t t) 3.0) (- a (/ 5.0 6.0))))))))) (/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))))))))))

  (/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))