Average Error: 45.2 → 0.0
Time: 1.2m
Precision: 64
Internal Precision: 384
\[\frac{\frac{\left(i \cdot i\right) \cdot \left(i \cdot i\right)}{\left(2 \cdot i\right) \cdot \left(2 \cdot i\right)}}{\left(2 \cdot i\right) \cdot \left(2 \cdot i\right) - 1.0}\]
\[\begin{array}{l} \mathbf{if}\;i \le 1429.9899038984013:\\ \;\;\;\;\frac{\frac{i}{2} \cdot \frac{i}{2}}{\left(i + i\right) \cdot \left(i + i\right) - 1.0}\\ \mathbf{else}:\\ \;\;\;\;\left(\frac{0.00390625}{{i}^{4}} + \frac{\frac{0.015625}{i}}{i}\right) + \frac{1}{16}\\ \end{array}\]

Error

Bits error versus i

Derivation

  1. Split input into 2 regimes
  2. if i < 1429.9899038984013

    1. Initial program 43.8

      \[\frac{\frac{\left(i \cdot i\right) \cdot \left(i \cdot i\right)}{\left(2 \cdot i\right) \cdot \left(2 \cdot i\right)}}{\left(2 \cdot i\right) \cdot \left(2 \cdot i\right) - 1.0}\]
    2. Applied simplify0.0

      \[\leadsto \color{blue}{\frac{\frac{i}{2} \cdot \frac{i}{2}}{\left(i + i\right) \cdot \left(i + i\right) - 1.0}}\]

    if 1429.9899038984013 < i

    1. Initial program 46.6

      \[\frac{\frac{\left(i \cdot i\right) \cdot \left(i \cdot i\right)}{\left(2 \cdot i\right) \cdot \left(2 \cdot i\right)}}{\left(2 \cdot i\right) \cdot \left(2 \cdot i\right) - 1.0}\]
    2. Applied simplify30.8

      \[\leadsto \color{blue}{\frac{\frac{i}{2} \cdot \frac{i}{2}}{\left(i + i\right) \cdot \left(i + i\right) - 1.0}}\]
    3. Taylor expanded around inf 0.0

      \[\leadsto \color{blue}{0.00390625 \cdot \frac{1}{{i}^{4}} + \left(\frac{1}{16} + 0.015625 \cdot \frac{1}{{i}^{2}}\right)}\]
    4. Applied simplify0

      \[\leadsto \color{blue}{\left(\frac{0.00390625}{{i}^{4}} + \frac{\frac{0.015625}{i}}{i}\right) + \frac{1}{16}}\]
  3. Recombined 2 regimes into one program.
  4. Removed slow pow expressions.

Runtime

Time bar (total: 1.2m)Debug logProfile

herbie shell --seed '#(1063027428 1192549564 1443466578 604016274 3637110559 1698629644)' 
(FPCore (i)
  :name "Octave 3.8, jcobi/4, as called"
  :pre (and (> i 0))
  (/ (/ (* (* i i) (* i i)) (* (* 2 i) (* 2 i))) (- (* (* 2 i) (* 2 i)) 1.0)))