Average Error: 42.1 → 10.2
Time: 4.0m
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 -883.0321444790968:\\ \;\;\;\;\frac{t \cdot \sqrt{2}}{\frac{\frac{t}{\sqrt{2}}}{x \cdot x} \cdot \left(1 - 2\right) - t \cdot \left(\sqrt{2} + \frac{\frac{2}{x}}{\sqrt{2}}\right)}\\ \mathbf{if}\;t \le -2.464172910231899 \cdot 10^{-173}:\\ \;\;\;\;\frac{t \cdot \sqrt{2}}{\sqrt{\sqrt{\left(\frac{4}{x} + 2\right) \cdot \left(t \cdot t\right) + \frac{2 \cdot \ell}{\frac{x}{\ell}}}} \cdot \sqrt{\sqrt{\left(\frac{4}{x} + 2\right) \cdot \left(t \cdot t\right) + \frac{2 \cdot \ell}{\frac{x}{\ell}}}}}\\ \mathbf{if}\;t \le -1.636775027050718 \cdot 10^{-207}:\\ \;\;\;\;\frac{t \cdot \sqrt{2}}{\frac{\frac{t}{\sqrt{2}}}{x \cdot x} \cdot \left(1 - 2\right) - t \cdot \left(\sqrt{2} + \frac{\frac{2}{x}}{\sqrt{2}}\right)}\\ \mathbf{if}\;t \le 1.7573613059582121 \cdot 10^{-302}:\\ \;\;\;\;\frac{t \cdot \sqrt{2}}{\sqrt{\sqrt{\left(\frac{4}{x} + 2\right) \cdot \left(t \cdot t\right) + \frac{2 \cdot \ell}{\frac{x}{\ell}}}} \cdot \sqrt{\sqrt{\left(\frac{4}{x} + 2\right) \cdot \left(t \cdot t\right) + \frac{2 \cdot \ell}{\frac{x}{\ell}}}}}\\ \mathbf{if}\;t \le 4.3345259186715104 \cdot 10^{-157}:\\ \;\;\;\;\frac{t \cdot \sqrt{2}}{t \cdot \left(\sqrt{2} + \frac{\frac{2}{x}}{\sqrt{2}}\right) + \frac{t}{\left(x \cdot x\right) \cdot \sqrt{2}} \cdot \left(2 - 1\right)}\\ \mathbf{if}\;t \le 4.910927045495491 \cdot 10^{+82}:\\ \;\;\;\;\frac{t \cdot \sqrt{2}}{\sqrt{\sqrt{\left(\frac{4}{x} + 2\right) \cdot \left(t \cdot t\right) + \frac{2 \cdot \ell}{\frac{x}{\ell}}}} \cdot \sqrt{\sqrt{\left(\frac{4}{x} + 2\right) \cdot \left(t \cdot t\right) + \frac{2 \cdot \ell}{\frac{x}{\ell}}}}}\\ \mathbf{else}:\\ \;\;\;\;\frac{t \cdot \sqrt{2}}{t \cdot \left(\sqrt{2} + \frac{\frac{2}{x}}{\sqrt{2}}\right) + \frac{t}{\left(x \cdot x\right) \cdot \sqrt{2}} \cdot \left(2 - 1\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 < -883.0321444790968 or -2.464172910231899e-173 < t < -1.636775027050718e-207

    1. Initial program 42.7

      \[\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. Taylor expanded around -inf 7.8

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

    3. Applied simplify7.8

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

    if -883.0321444790968 < t < -2.464172910231899e-173 or -1.636775027050718e-207 < t < 1.7573613059582121e-302 or 4.3345259186715104e-157 < t < 4.910927045495491e+82

    1. Initial program 33.9

      \[\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. Taylor expanded around inf 14.1

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

    3. Applied simplify10.1

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

    5. Applied add-sqr-sqrt10.1

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

    6. Applied sqrt-prod10.2

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

    if 1.7573613059582121e-302 < t < 4.3345259186715104e-157 or 4.910927045495491e+82 < t

    1. Initial program 51.6

      \[\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. Taylor expanded around inf 12.6

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

    3. Applied simplify12.6

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

Runtime

Time bar (total: 4.0m)Debug logProfile

herbie shell --seed '#(1071852389 864846987 1238109217 3425890003 4124793586 650694553)' 
(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)))))