Average Error: 13.0 → 0.4
Time: 5.3m
Precision: 64
Internal Precision: 384
\[\left(\left(3 + \frac{2}{r \cdot r}\right) - \frac{\left(0.125 \cdot \left(3 - 2 \cdot v\right)\right) \cdot \left(\left(\left(w \cdot w\right) \cdot r\right) \cdot r\right)}{1 - v}\right) - 4.5\]
↓
\[\begin{array}{l} \mathbf{if}\;v \le -1.8061603950641245 \cdot 10^{+32}:\\ \;\;\;\;\left(\frac{2}{r \cdot r} + \left(3 - 4.5\right)\right) - \left(\left(\left(w \cdot r\right) \cdot \left(w \cdot r\right)\right) \cdot 0.25 - \frac{\left(w \cdot r\right) \cdot \left(w \cdot r\right)}{v} \cdot \left(0.125 + \frac{0.125}{v}\right)\right)\\ \mathbf{if}\;v \le 4.754640524038839 \cdot 10^{-19}:\\ \;\;\;\;\left(\left(3 + \frac{\frac{2}{r}}{r}\right) - \frac{0.125 \cdot \left(\left(r \cdot w\right) \cdot \left(r \cdot w\right)\right)}{\frac{1 - v \cdot v}{3 - \left(v + v\right)}} \cdot \left(1 + v\right)\right) - 4.5\\ \mathbf{else}:\\ \;\;\;\;\left(\frac{2}{r \cdot r} + \left(3 - 4.5\right)\right) - \left(\left(\left(w \cdot r\right) \cdot \left(w \cdot r\right)\right) \cdot 0.25 - \frac{\left(w \cdot r\right) \cdot \left(w \cdot r\right)}{v} \cdot \left(0.125 + \frac{0.125}{v}\right)\right)\\ \end{array}\]

Error

Bits error versus v

Bits error versus w

Bits error versus r

Derivation

  1. Split input into 2 regimes

  2. if v < -1.8061603950641245e+32 or 4.754640524038839e-19 < v

    1. Initial program 17.4

      \[\left(\left(3 + \frac{2}{r \cdot r}\right) - \frac{\left(0.125 \cdot \left(3 - 2 \cdot v\right)\right) \cdot \left(\left(\left(w \cdot w\right) \cdot r\right) \cdot r\right)}{1 - v}\right) - 4.5\]

    2. Taylor expanded around inf 17.2

      \[\leadsto \left(\left(3 + \frac{2}{r \cdot r}\right) - \color{blue}{\left(0.25 \cdot \left({w}^{2} \cdot {r}^{2}\right) - \left(0.125 \cdot \frac{{w}^{2} \cdot {r}^{2}}{v} + 0.125 \cdot \frac{{w}^{2} \cdot {r}^{2}}{{v}^{2}}\right)\right)}\right) - 4.5\]

    3. Applied simplify0.3

      \[\leadsto \color{blue}{\left(\frac{2}{r \cdot r} + \left(3 - 4.5\right)\right) - \left(\left(\left(w \cdot r\right) \cdot \left(w \cdot r\right)\right) \cdot 0.25 - \frac{\left(w \cdot r\right) \cdot \left(w \cdot r\right)}{v} \cdot \left(0.125 + \frac{0.125}{v}\right)\right)}\]

    if -1.8061603950641245e+32 < v < 4.754640524038839e-19

    1. Initial program 8.6

      \[\left(\left(3 + \frac{2}{r \cdot r}\right) - \frac{\left(0.125 \cdot \left(3 - 2 \cdot v\right)\right) \cdot \left(\left(\left(w \cdot w\right) \cdot r\right) \cdot r\right)}{1 - v}\right) - 4.5\]
    2. Using strategy rm

    3. Applied flip--8.6

      \[\leadsto \left(\left(3 + \frac{2}{r \cdot r}\right) - \frac{\left(0.125 \cdot \left(3 - 2 \cdot v\right)\right) \cdot \left(\left(\left(w \cdot w\right) \cdot r\right) \cdot r\right)}{\color{blue}{\frac{1 \cdot 1 - v \cdot v}{1 + v}}}\right) - 4.5\]

    4. Applied associate-/r/8.6

      \[\leadsto \left(\left(3 + \frac{2}{r \cdot r}\right) - \color{blue}{\frac{\left(0.125 \cdot \left(3 - 2 \cdot v\right)\right) \cdot \left(\left(\left(w \cdot w\right) \cdot r\right) \cdot r\right)}{1 \cdot 1 - v \cdot v} \cdot \left(1 + v\right)}\right) - 4.5\]

    5. Applied simplify0.4

      \[\leadsto \left(\left(3 + \frac{2}{r \cdot r}\right) - \color{blue}{\frac{0.125 \cdot \left(\left(r \cdot w\right) \cdot \left(r \cdot w\right)\right)}{\frac{1 - v \cdot v}{3 - \left(v + v\right)}}} \cdot \left(1 + v\right)\right) - 4.5\]
    6. Using strategy rm

    7. Applied associate-/r*0.4

      \[\leadsto \left(\left(3 + \color{blue}{\frac{\frac{2}{r}}{r}}\right) - \frac{0.125 \cdot \left(\left(r \cdot w\right) \cdot \left(r \cdot w\right)\right)}{\frac{1 - v \cdot v}{3 - \left(v + v\right)}} \cdot \left(1 + v\right)\right) - 4.5\]
  3. Recombined 2 regimes into one program.
  4. Removed slow pow expressions.

Runtime

Time bar (total: 5.3m)Debug log

herbie shell --seed '#(1567391828 2030694642 2833800258 828025724 3004380912 3532991858)' +o setup:early-exit +o reduce:binary-search
(FPCore (v w r)
  :name "Rosa's TurbineBenchmark"
  (- (- (+ 3 (/ 2 (* r r))) (/ (* (* 0.125 (- 3 (* 2 v))) (* (* (* w w) r) r)) (- 1 v))) 4.5))