Average Error: 59.5 → 2.3
Time: 1.9m
Precision: 64
Internal Precision: 1408
\[-\frac{1}{\frac{\pi}{4}} \cdot \log \left(\frac{e^{\frac{\pi}{4} \cdot f} + e^{-\frac{\pi}{4} \cdot f}}{e^{\frac{\pi}{4} \cdot f} - e^{-\frac{\pi}{4} \cdot f}}\right)\]
↓
\[-\left(\left(4 \cdot \frac{\log \left(\sqrt[3]{\frac{4}{\pi}} \cdot \sqrt[3]{\frac{4}{\pi}}\right) + \log \left(\sqrt[3]{\frac{4}{\pi}}\right)}{\pi} + \frac{1}{12} \cdot \left(\pi \cdot {f}^{2}\right)\right) - \left(\frac{7}{5760} \cdot \left({\pi}^{3} \cdot {f}^{4}\right) + 4 \cdot \frac{\log f}{\pi}\right)\right)\]

Error

Bits error versus f

Derivation

  1. Initial program 59.5

    \[-\frac{1}{\frac{\pi}{4}} \cdot \log \left(\frac{e^{\frac{\pi}{4} \cdot f} + e^{-\frac{\pi}{4} \cdot f}}{e^{\frac{\pi}{4} \cdot f} - e^{-\frac{\pi}{4} \cdot f}}\right)\]

  2. Taylor expanded around 0 2.3

    \[\leadsto -\color{blue}{\left(\left(4 \cdot \frac{\log \left(\frac{4}{\pi}\right)}{\pi} + \frac{1}{12} \cdot \left(\pi \cdot {f}^{2}\right)\right) - \left(\frac{7}{5760} \cdot \left({\pi}^{3} \cdot {f}^{4}\right) + 4 \cdot \frac{\log f}{\pi}\right)\right)}\]
  3. Using strategy rm

  4. Applied add-cube-cbrt2.3

    \[\leadsto -\left(\left(4 \cdot \frac{\log \color{blue}{\left(\left(\sqrt[3]{\frac{4}{\pi}} \cdot \sqrt[3]{\frac{4}{\pi}}\right) \cdot \sqrt[3]{\frac{4}{\pi}}\right)}}{\pi} + \frac{1}{12} \cdot \left(\pi \cdot {f}^{2}\right)\right) - \left(\frac{7}{5760} \cdot \left({\pi}^{3} \cdot {f}^{4}\right) + 4 \cdot \frac{\log f}{\pi}\right)\right)\]

  5. Applied log-prod2.3

    \[\leadsto -\left(\left(4 \cdot \frac{\color{blue}{\log \left(\sqrt[3]{\frac{4}{\pi}} \cdot \sqrt[3]{\frac{4}{\pi}}\right) + \log \left(\sqrt[3]{\frac{4}{\pi}}\right)}}{\pi} + \frac{1}{12} \cdot \left(\pi \cdot {f}^{2}\right)\right) - \left(\frac{7}{5760} \cdot \left({\pi}^{3} \cdot {f}^{4}\right) + 4 \cdot \frac{\log f}{\pi}\right)\right)\]

Runtime

Time bar (total: 1.9m)Debug logProfile

herbie shell --seed '#(1064173506 2580572819 2847706409 4129882574 1125180799 1845288547)' +o rules:numerics
(FPCore (f)
  :name "VandenBroeck and Keller, Equation (20)"
  (- (* (/ 1 (/ PI 4)) (log (/ (+ (exp (* (/ PI 4) f)) (exp (- (* (/ PI 4) f)))) (- (exp (* (/ PI 4) f)) (exp (- (* (/ PI 4) f)))))))))