Average Error: 3.9 → 2.3
Time: 40.8s
Precision: 64
Internal Precision: 128
\[\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 -2.913213379994578 \cdot 10^{-175} \lor \neg \left(t \le 5.391777885674001 \cdot 10^{-268}\right):\\ \;\;\;\;\frac{x}{(y \cdot \left(e^{(\left(\left(\frac{5.0}{6.0} + a\right) - \frac{\frac{2.0}{t}}{3.0}\right) \cdot \left(2.0 \cdot \left(c - b\right)\right) + \left(\sqrt{t + a} \cdot \frac{2.0 \cdot z}{t}\right))_*}\right) + x)_*}\\ \mathbf{else}:\\ \;\;\;\;\frac{x}{(y \cdot \left(e^{(\left(\frac{1.3333333333333333}{t}\right) \cdot \left(b - c\right) + \left(1.6666666666666667 \cdot c\right))_*}\right) + x)_*}\\ \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

Derivation

  1. Split input into 2 regimes

  2. if t < -2.913213379994578e-175 or 5.391777885674001e-268 < t

    1. Initial program 3.0

      \[\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. Initial simplification1.1

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

    if -2.913213379994578e-175 < t < 5.391777885674001e-268

    1. Initial program 9.7

      \[\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. Initial simplification6.6

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

    3. Taylor expanded around 0 15.6

      \[\leadsto \frac{x}{(y \cdot \left(e^{\color{blue}{\left(1.3333333333333333 \cdot \frac{b}{t} + 1.6666666666666667 \cdot c\right) - 1.3333333333333333 \cdot \frac{c}{t}}}\right) + x)_*}\]

    4. Simplified10.2

      \[\leadsto \frac{x}{(y \cdot \left(e^{\color{blue}{(\left(\frac{1.3333333333333333}{t}\right) \cdot \left(b - c\right) + \left(1.6666666666666667 \cdot c\right))_*}}\right) + x)_*}\]
  3. Recombined 2 regimes into one program.

  4. Final simplification2.3

    \[\leadsto \begin{array}{l} \mathbf{if}\;t \le -2.913213379994578 \cdot 10^{-175} \lor \neg \left(t \le 5.391777885674001 \cdot 10^{-268}\right):\\ \;\;\;\;\frac{x}{(y \cdot \left(e^{(\left(\left(\frac{5.0}{6.0} + a\right) - \frac{\frac{2.0}{t}}{3.0}\right) \cdot \left(2.0 \cdot \left(c - b\right)\right) + \left(\sqrt{t + a} \cdot \frac{2.0 \cdot z}{t}\right))_*}\right) + x)_*}\\ \mathbf{else}:\\ \;\;\;\;\frac{x}{(y \cdot \left(e^{(\left(\frac{1.3333333333333333}{t}\right) \cdot \left(b - c\right) + \left(1.6666666666666667 \cdot c\right))_*}\right) + x)_*}\\ \end{array}\]

Runtime

Time bar (total: 40.8s)Debug logProfile

herbie shell --seed 2018274 +o rules:numerics
(FPCore (x y z t a b c)
  :name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2"
  (/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))