Average Error: 42.6 → 9.7
Time: 1.1m
Precision: 64
Internal Precision: 1344
\[\frac{\sqrt{2} \cdot t}{\sqrt{\frac{x + 1}{x - 1} \cdot \left(\ell \cdot \ell + 2 \cdot \left(t \cdot t\right)\right) - \ell \cdot \ell}}\]
↓
\[\begin{array}{l} \mathbf{if}\;t \le -1.6349727824419413 \cdot 10^{+115}:\\ \;\;\;\;\frac{\sqrt{2} \cdot t}{-2 \cdot \frac{t}{x \cdot \sqrt{2}} + \sqrt{2} \cdot \left(-t\right)}\\ \mathbf{elif}\;t \le -2.501180277482156 \cdot 10^{-138}:\\ \;\;\;\;\frac{\sqrt{2} \cdot t}{\sqrt{(2 \cdot \left((\left(\frac{\ell}{x}\right) \cdot \ell + \left(t \cdot t\right))_*\right) + \left(\left(t \cdot 4\right) \cdot \frac{t}{x}\right))_*}}\\ \mathbf{elif}\;t \le -7.751097215014646 \cdot 10^{-226}:\\ \;\;\;\;\frac{\sqrt{2} \cdot t}{-2 \cdot \frac{t}{x \cdot \sqrt{2}} + \sqrt{2} \cdot \left(-t\right)}\\ \mathbf{elif}\;t \le 2.2486994535592485 \cdot 10^{+50}:\\ \;\;\;\;\frac{\sqrt{2} \cdot t}{\sqrt{(2 \cdot \left((\left(\frac{\ell}{x}\right) \cdot \ell + \left(t \cdot t\right))_*\right) + \left(\left(t \cdot 4\right) \cdot \frac{t}{x}\right))_*}}\\ \mathbf{else}:\\ \;\;\;\;\frac{\sqrt{2} \cdot t}{(\left(\frac{\frac{2}{x}}{\sqrt{2}}\right) \cdot t + \left(\sqrt{2} \cdot t\right))_*}\\ \end{array}\]

Error

Bits error versus x

Bits error versus l

Bits error versus t

Derivation

  1. Split input into 3 regimes

  2. if t < -1.6349727824419413e+115 or -2.501180277482156e-138 < t < -7.751097215014646e-226

    1. Initial program 53.8

      \[\frac{\sqrt{2} \cdot t}{\sqrt{\frac{x + 1}{x - 1} \cdot \left(\ell \cdot \ell + 2 \cdot \left(t \cdot t\right)\right) - \ell \cdot \ell}}\]

    2. Initial simplification53.8

      \[\leadsto \frac{t \cdot \sqrt{2}}{\sqrt{(\left((\left(2 \cdot t\right) \cdot t + \left(\ell \cdot \ell\right))_*\right) \cdot \left(\frac{1 + x}{-1 + x}\right) + \left(\ell \cdot \left(-\ell\right)\right))_*}}\]

    3. Taylor expanded around -inf 9.2

      \[\leadsto \frac{t \cdot \sqrt{2}}{\color{blue}{-\left(t \cdot \sqrt{2} + 2 \cdot \frac{t}{\sqrt{2} \cdot x}\right)}}\]

    if -1.6349727824419413e+115 < t < -2.501180277482156e-138 or -7.751097215014646e-226 < t < 2.2486994535592485e+50

    1. Initial program 36.1

      \[\frac{\sqrt{2} \cdot t}{\sqrt{\frac{x + 1}{x - 1} \cdot \left(\ell \cdot \ell + 2 \cdot \left(t \cdot t\right)\right) - \ell \cdot \ell}}\]

    2. Initial simplification36.0

      \[\leadsto \frac{t \cdot \sqrt{2}}{\sqrt{(\left((\left(2 \cdot t\right) \cdot t + \left(\ell \cdot \ell\right))_*\right) \cdot \left(\frac{1 + x}{-1 + x}\right) + \left(\ell \cdot \left(-\ell\right)\right))_*}}\]

    3. Taylor expanded around inf 16.7

      \[\leadsto \frac{t \cdot \sqrt{2}}{\sqrt{\color{blue}{2 \cdot {t}^{2} + \left(4 \cdot \frac{{t}^{2}}{x} + 2 \cdot \frac{{\ell}^{2}}{x}\right)}}}\]

    4. Simplified12.4

      \[\leadsto \frac{t \cdot \sqrt{2}}{\sqrt{\color{blue}{(2 \cdot \left((\left(\frac{\ell}{x}\right) \cdot \ell + \left(t \cdot t\right))_*\right) + \left(\left(t \cdot 4\right) \cdot \frac{t}{x}\right))_*}}}\]

    if 2.2486994535592485e+50 < t

    1. Initial program 44.8

      \[\frac{\sqrt{2} \cdot t}{\sqrt{\frac{x + 1}{x - 1} \cdot \left(\ell \cdot \ell + 2 \cdot \left(t \cdot t\right)\right) - \ell \cdot \ell}}\]

    2. Initial simplification44.8

      \[\leadsto \frac{t \cdot \sqrt{2}}{\sqrt{(\left((\left(2 \cdot t\right) \cdot t + \left(\ell \cdot \ell\right))_*\right) \cdot \left(\frac{1 + x}{-1 + x}\right) + \left(\ell \cdot \left(-\ell\right)\right))_*}}\]

    3. Taylor expanded around inf 4.6

      \[\leadsto \frac{t \cdot \sqrt{2}}{\color{blue}{t \cdot \sqrt{2} + 2 \cdot \frac{t}{\sqrt{2} \cdot x}}}\]

    4. Simplified4.6

      \[\leadsto \frac{t \cdot \sqrt{2}}{\color{blue}{(\left(\frac{\frac{2}{x}}{\sqrt{2}}\right) \cdot t + \left(t \cdot \sqrt{2}\right))_*}}\]
  3. Recombined 3 regimes into one program.

  4. Final simplification9.7

    \[\leadsto \begin{array}{l} \mathbf{if}\;t \le -1.6349727824419413 \cdot 10^{+115}:\\ \;\;\;\;\frac{\sqrt{2} \cdot t}{-2 \cdot \frac{t}{x \cdot \sqrt{2}} + \sqrt{2} \cdot \left(-t\right)}\\ \mathbf{elif}\;t \le -2.501180277482156 \cdot 10^{-138}:\\ \;\;\;\;\frac{\sqrt{2} \cdot t}{\sqrt{(2 \cdot \left((\left(\frac{\ell}{x}\right) \cdot \ell + \left(t \cdot t\right))_*\right) + \left(\left(t \cdot 4\right) \cdot \frac{t}{x}\right))_*}}\\ \mathbf{elif}\;t \le -7.751097215014646 \cdot 10^{-226}:\\ \;\;\;\;\frac{\sqrt{2} \cdot t}{-2 \cdot \frac{t}{x \cdot \sqrt{2}} + \sqrt{2} \cdot \left(-t\right)}\\ \mathbf{elif}\;t \le 2.2486994535592485 \cdot 10^{+50}:\\ \;\;\;\;\frac{\sqrt{2} \cdot t}{\sqrt{(2 \cdot \left((\left(\frac{\ell}{x}\right) \cdot \ell + \left(t \cdot t\right))_*\right) + \left(\left(t \cdot 4\right) \cdot \frac{t}{x}\right))_*}}\\ \mathbf{else}:\\ \;\;\;\;\frac{\sqrt{2} \cdot t}{(\left(\frac{\frac{2}{x}}{\sqrt{2}}\right) \cdot t + \left(\sqrt{2} \cdot t\right))_*}\\ \end{array}\]

Runtime

Time bar (total: 1.1m)Debug logProfile

herbie shell --seed 2018252 +o rules:numerics
(FPCore (x l t)
  :name "Toniolo and Linder, Equation (7)"
  (/ (* (sqrt 2) t) (sqrt (- (* (/ (+ x 1) (- x 1)) (+ (* l l) (* 2 (* t t)))) (* l l)))))