Average Error: 2.1 → 1.2
Time: 2.7m
Precision: 64
Internal Precision: 576
\[\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}{y}\]
↓
\[\begin{array}{l} \mathbf{if}\;\frac{\left({z}^{y} \cdot {a}^{\left(-1.0\right)}\right) \cdot {a}^{t}}{\sqrt[3]{{\left(e^{b}\right)}^{3} \cdot {\left(\frac{y}{x}\right)}^{3}}} \le +\infty:\\ \;\;\;\;\frac{\left({z}^{y} \cdot {a}^{\left(-1.0\right)}\right) \cdot {a}^{t}}{\frac{y}{x} \cdot e^{b}}\\ \mathbf{else}:\\ \;\;\;\;\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}{y}\\ \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

Derivation

  1. Split input into 2 regimes

  2. if (/ (* (* (pow z y) (pow a (- 1.0))) (pow a t)) (cbrt (* (pow (exp b) 3) (pow (/ y x) 3)))) < +inf.0

    1. Initial program 3.8

      \[\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}{y}\]

    2. Taylor expanded around inf 3.8

      \[\leadsto \frac{x \cdot \color{blue}{e^{1.0 \cdot \log \left(\frac{1}{a}\right) - \left(y \cdot \log \left(\frac{1}{z}\right) + \left(b + t \cdot \log \left(\frac{1}{a}\right)\right)\right)}}}{y}\]

    3. Applied simplify2.1

      \[\leadsto \color{blue}{\frac{\left({z}^{y} \cdot {a}^{\left(-1.0\right)}\right) \cdot {a}^{t}}{\frac{y}{x} \cdot e^{b}}}\]

    if +inf.0 < (/ (* (* (pow z y) (pow a (- 1.0))) (pow a t)) (cbrt (* (pow (exp b) 3) (pow (/ y x) 3))))

    1. Initial program 0.0

      \[\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}{y}\]
  3. Recombined 2 regimes into one program.

Runtime

Time bar (total: 2.7m)Debug logProfile

herbie shell --seed '#(1071501266 3581234924 1086666455 2685055582 1243441566 1802958749)' 
(FPCore (x y z t a b)
  :name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2"
  (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))