Average Error: 29.4 → 0.1
Time: 24.3s
Precision: 64
Internal Precision: 1408
\[\frac{x}{x + 1} - \frac{x + 1}{x - 1}\]
↓
\[\begin{array}{l} \mathbf{if}\;\left(\sqrt[3]{(\left(1 + \frac{3}{x}\right) \cdot \left(\frac{-1}{x \cdot x}\right) + \left(\frac{-3}{x}\right))_*} \cdot \sqrt[3]{(\left(1 + \frac{3}{x}\right) \cdot \left(\frac{-1}{x \cdot x}\right) + \left(\frac{-3}{x}\right))_*}\right) \cdot \sqrt[3]{(\left(1 + \frac{3}{x}\right) \cdot \left(\frac{-1}{x \cdot x}\right) + \left(\frac{-3}{x}\right))_*} \le -0.010637352795128258:\\ \;\;\;\;(x \cdot \left(\frac{1}{x + 1}\right) + \left(-\frac{x + 1}{x - 1}\right))_*\\ \mathbf{if}\;\left(\sqrt[3]{(\left(1 + \frac{3}{x}\right) \cdot \left(\frac{-1}{x \cdot x}\right) + \left(\frac{-3}{x}\right))_*} \cdot \sqrt[3]{(\left(1 + \frac{3}{x}\right) \cdot \left(\frac{-1}{x \cdot x}\right) + \left(\frac{-3}{x}\right))_*}\right) \cdot \sqrt[3]{(\left(1 + \frac{3}{x}\right) \cdot \left(\frac{-1}{x \cdot x}\right) + \left(\frac{-3}{x}\right))_*} \le 0.0002862999886700314:\\ \;\;\;\;(\left(1 + \frac{3}{x}\right) \cdot \left(\frac{-1}{x \cdot x}\right) + \left(\frac{-3}{x}\right))_*\\ \mathbf{if}\;\left(\sqrt[3]{(\left(1 + \frac{3}{x}\right) \cdot \left(\frac{-1}{x \cdot x}\right) + \left(\frac{-3}{x}\right))_*} \cdot \sqrt[3]{(\left(1 + \frac{3}{x}\right) \cdot \left(\frac{-1}{x \cdot x}\right) + \left(\frac{-3}{x}\right))_*}\right) \cdot \sqrt[3]{(\left(1 + \frac{3}{x}\right) \cdot \left(\frac{-1}{x \cdot x}\right) + \left(\frac{-3}{x}\right))_*} \le +\infty:\\ \;\;\;\;\sqrt{\frac{x}{x + 1} - \frac{x + 1}{x - 1}} \cdot \sqrt{\frac{x}{x + 1} - \frac{x + 1}{x - 1}}\\ \mathbf{else}:\\ \;\;\;\;\sqrt{\frac{x}{x + 1} - \frac{x + 1}{x - 1}} \cdot \sqrt{\frac{x}{x + 1} - \frac{x + 1}{x - 1}}\\ \end{array}\]

Error

Bits error versus x

Derivation

  1. Split input into 3 regimes

  2. if (* (* (cbrt (fma (+ 1 (/ 3 x)) (/ (- 1) (* x x)) (/ (- 3) x))) (cbrt (fma (+ 1 (/ 3 x)) (/ (- 1) (* x x)) (/ (- 3) x)))) (cbrt (fma (+ 1 (/ 3 x)) (/ (- 1) (* x x)) (/ (- 3) x)))) < -0.010637352795128258

    1. Initial program 0.1

      \[\frac{x}{x + 1} - \frac{x + 1}{x - 1}\]
    2. Using strategy rm

    3. Applied div-inv0.0

      \[\leadsto \color{blue}{x \cdot \frac{1}{x + 1}} - \frac{x + 1}{x - 1}\]

    4. Applied fma-neg0.0

      \[\leadsto \color{blue}{(x \cdot \left(\frac{1}{x + 1}\right) + \left(-\frac{x + 1}{x - 1}\right))_*}\]

    if -0.010637352795128258 < (* (* (cbrt (fma (+ 1 (/ 3 x)) (/ (- 1) (* x x)) (/ (- 3) x))) (cbrt (fma (+ 1 (/ 3 x)) (/ (- 1) (* x x)) (/ (- 3) x)))) (cbrt (fma (+ 1 (/ 3 x)) (/ (- 1) (* x x)) (/ (- 3) x)))) < 0.0002862999886700314

    1. Initial program 59.0

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

    2. Taylor expanded around inf 0.4

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

    3. Applied simplify0.1

      \[\leadsto \color{blue}{(\left(1 + \frac{3}{x}\right) \cdot \left(\frac{-1}{x \cdot x}\right) + \left(\frac{-3}{x}\right))_*}\]

    if 0.0002862999886700314 < (* (* (cbrt (fma (+ 1 (/ 3 x)) (/ (- 1) (* x x)) (/ (- 3) x))) (cbrt (fma (+ 1 (/ 3 x)) (/ (- 1) (* x x)) (/ (- 3) x)))) (cbrt (fma (+ 1 (/ 3 x)) (/ (- 1) (* x x)) (/ (- 3) x))))

    1. Initial program 0.1

      \[\frac{x}{x + 1} - \frac{x + 1}{x - 1}\]
    2. Using strategy rm

    3. Applied add-sqr-sqrt0.2

      \[\leadsto \color{blue}{\sqrt{\frac{x}{x + 1} - \frac{x + 1}{x - 1}} \cdot \sqrt{\frac{x}{x + 1} - \frac{x + 1}{x - 1}}}\]
  3. Recombined 3 regimes into one program.

Runtime

Time bar (total: 24.3s)Debug logProfile

herbie shell --seed '#(1070131407 1246090267 3027482374 2150728003 2026520792 2347815650)' +o rules:numerics
(FPCore (x)
  :name "Asymptote C"
  (- (/ x (+ x 1)) (/ (+ x 1) (- x 1))))