Average Error: 15.0 → 0.0
Time: 26.2s
Precision: 64
Internal Precision: 384
\[\frac{x}{x \cdot x + 1}\]
↓
\[\begin{array}{l} \mathbf{if}\;x \le -1.0689697925879056 \cdot 10^{+65}:\\ \;\;\;\;\left(\frac{1}{{x}^{5}} + \frac{1}{x}\right) - \frac{1}{{x}^{3}}\\ \mathbf{if}\;x \le 2458290.8697445355:\\ \;\;\;\;\frac{x}{x \cdot x + 1}\\ \mathbf{else}:\\ \;\;\;\;\left(\frac{1}{{x}^{5}} + \frac{1}{x}\right) - \frac{1}{{x}^{3}}\\ \end{array}\]

Error

Bits error versus x

Target

Original15.0
Target0.1
Herbie0.0
\[\frac{1}{x + \frac{1}{x}}\]

Derivation

  1. Split input into 2 regimes

  2. if x < -1.0689697925879056e+65 or 2458290.8697445355 < x

    1. Initial program 34.6

      \[\frac{x}{x \cdot x + 1}\]

    2. Taylor expanded around inf 0.0

      \[\leadsto \color{blue}{\left(\frac{1}{{x}^{5}} + \frac{1}{x}\right) - \frac{1}{{x}^{3}}}\]

    if -1.0689697925879056e+65 < x < 2458290.8697445355

    1. Initial program 0.0

      \[\frac{x}{x \cdot x + 1}\]
  3. Recombined 2 regimes into one program.

Runtime

Time bar (total: 26.2s)Debug logProfile

herbie shell --seed '#(1070706311 3771791028 4128836681 4194990999 2341756049 504035650)' +o rules:numerics
(FPCore (x)
  :name "x / (x^2 + 1)"

  :herbie-target
  (/ 1 (+ x (/ 1 x)))

  (/ x (+ (* x x) 1)))