Average Error: 61.3 → 2.3
Time: 22.9min
Precision: binary64
\[-\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(0.083333333333333343 \cdot \left({f}^{2} \cdot \pi\right) + \left(\frac{4}{\pi} \cdot \left(\log \left(\frac{4}{\pi}\right) - \log f\right) - \left({\left({\left(\sqrt[3]{\pi} \cdot \sqrt[3]{\pi}\right)}^{\left(\sqrt{3}\right)} \cdot {\left(\sqrt[3]{\pi}\right)}^{\left(\sqrt{3}\right)}\right)}^{\left(\sqrt{3}\right)} \cdot {f}^{4}\right) \cdot \left(3.472222222222224 \cdot 10^{-4} + \frac{0.0138888888888888899}{{4}^{2}}\right)\right)\right)\]

Error

Bits error versus f

Derivation

  1. Initial program 61.3

    \[-\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.4

    \[\leadsto -\frac{1}{\frac{\pi}{4}} \cdot \log \left(\frac{e^{\frac{\pi}{4} \cdot f} + e^{-\frac{\pi}{4} \cdot f}}{\color{blue}{0.5 \cdot \left(f \cdot \pi\right) + \left(0.00520833333333333304 \cdot \left({f}^{3} \cdot {\pi}^{3}\right) + 1.62760416666666664 \cdot 10^{-5} \cdot \left({f}^{5} \cdot {\pi}^{5}\right)\right)}}\right)\]
  3. Taylor expanded around 0 2.3

    \[\leadsto -\color{blue}{\left(\left(4 \cdot \frac{\log \left(\frac{4}{\pi}\right)}{\pi} + 0.083333333333333343 \cdot \left({f}^{2} \cdot \pi\right)\right) - \left(0.0138888888888888899 \cdot \frac{{\pi}^{3} \cdot {f}^{4}}{{4}^{2}} + \left(3.472222222222224 \cdot 10^{-4} \cdot \left({\pi}^{3} \cdot {f}^{4}\right) + 4 \cdot \frac{\log f}{\pi}\right)\right)\right)}\]
  4. Simplified2.3

    \[\leadsto -\color{blue}{\left(0.083333333333333343 \cdot \left({f}^{2} \cdot \pi\right) + \left(\frac{4}{\pi} \cdot \left(\log \left(\frac{4}{\pi}\right) - \log f\right) - \left({\pi}^{3} \cdot {f}^{4}\right) \cdot \left(3.472222222222224 \cdot 10^{-4} + \frac{0.0138888888888888899}{{4}^{2}}\right)\right)\right)}\]
  5. Using strategy rm
  6. Applied add-sqr-sqrt2.3

    \[\leadsto -\left(0.083333333333333343 \cdot \left({f}^{2} \cdot \pi\right) + \left(\frac{4}{\pi} \cdot \left(\log \left(\frac{4}{\pi}\right) - \log f\right) - \left({\pi}^{\color{blue}{\left(\sqrt{3} \cdot \sqrt{3}\right)}} \cdot {f}^{4}\right) \cdot \left(3.472222222222224 \cdot 10^{-4} + \frac{0.0138888888888888899}{{4}^{2}}\right)\right)\right)\]
  7. Applied pow-unpow2.3

    \[\leadsto -\left(0.083333333333333343 \cdot \left({f}^{2} \cdot \pi\right) + \left(\frac{4}{\pi} \cdot \left(\log \left(\frac{4}{\pi}\right) - \log f\right) - \left(\color{blue}{{\left({\pi}^{\left(\sqrt{3}\right)}\right)}^{\left(\sqrt{3}\right)}} \cdot {f}^{4}\right) \cdot \left(3.472222222222224 \cdot 10^{-4} + \frac{0.0138888888888888899}{{4}^{2}}\right)\right)\right)\]
  8. Using strategy rm
  9. Applied add-cube-cbrt2.3

    \[\leadsto -\left(0.083333333333333343 \cdot \left({f}^{2} \cdot \pi\right) + \left(\frac{4}{\pi} \cdot \left(\log \left(\frac{4}{\pi}\right) - \log f\right) - \left({\left({\color{blue}{\left(\left(\sqrt[3]{\pi} \cdot \sqrt[3]{\pi}\right) \cdot \sqrt[3]{\pi}\right)}}^{\left(\sqrt{3}\right)}\right)}^{\left(\sqrt{3}\right)} \cdot {f}^{4}\right) \cdot \left(3.472222222222224 \cdot 10^{-4} + \frac{0.0138888888888888899}{{4}^{2}}\right)\right)\right)\]
  10. Applied unpow-prod-down2.3

    \[\leadsto -\left(0.083333333333333343 \cdot \left({f}^{2} \cdot \pi\right) + \left(\frac{4}{\pi} \cdot \left(\log \left(\frac{4}{\pi}\right) - \log f\right) - \left({\color{blue}{\left({\left(\sqrt[3]{\pi} \cdot \sqrt[3]{\pi}\right)}^{\left(\sqrt{3}\right)} \cdot {\left(\sqrt[3]{\pi}\right)}^{\left(\sqrt{3}\right)}\right)}}^{\left(\sqrt{3}\right)} \cdot {f}^{4}\right) \cdot \left(3.472222222222224 \cdot 10^{-4} + \frac{0.0138888888888888899}{{4}^{2}}\right)\right)\right)\]
  11. Final simplification2.3

    \[\leadsto -\left(0.083333333333333343 \cdot \left({f}^{2} \cdot \pi\right) + \left(\frac{4}{\pi} \cdot \left(\log \left(\frac{4}{\pi}\right) - \log f\right) - \left({\left({\left(\sqrt[3]{\pi} \cdot \sqrt[3]{\pi}\right)}^{\left(\sqrt{3}\right)} \cdot {\left(\sqrt[3]{\pi}\right)}^{\left(\sqrt{3}\right)}\right)}^{\left(\sqrt{3}\right)} \cdot {f}^{4}\right) \cdot \left(3.472222222222224 \cdot 10^{-4} + \frac{0.0138888888888888899}{{4}^{2}}\right)\right)\right)\]

Reproduce

herbie shell --seed 2020182 
(FPCore (f)
  :name "VandenBroeck and Keller, Equation (20)"
  :precision binary64
  (neg (* (/ 1.0 (/ PI 4.0)) (log (/ (+ (exp (* (/ PI 4.0) f)) (exp (neg (* (/ PI 4.0) f)))) (- (exp (* (/ PI 4.0) f)) (exp (neg (* (/ PI 4.0) f)))))))))