Average Error: 4.1 → 0.8
Time: 3.5m
Precision: 64
Internal Precision: 2624
\[\frac{{\left(\frac{1}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-s}}\right)}^{c_n}}{{\left(\frac{1}{1 + e^{-t}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}\]
\[e^{(c_n \cdot \left(\log_* (1 + \frac{-1}{1 + e^{-s}})\right) + \left((\left((\frac{1}{8} \cdot t + \frac{-1}{2})_*\right) \cdot t + \left(\log 2\right))_* \cdot c_p\right))_* - (\left((t \cdot \left((\frac{-1}{8} \cdot t + \frac{-1}{2})_*\right) + \left(\log \frac{1}{2}\right))_*\right) \cdot c_n + \left(\log_* (1 + e^{-s}) \cdot c_p\right))_*}\]

Error

Bits error versus c_p

Bits error versus c_n

Bits error versus t

Bits error versus s

Target

Original4.1
Target2.1
Herbie0.8
\[{\left(\frac{1 + e^{-t}}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(\frac{1 + e^{t}}{1 + e^{s}}\right)}^{c_n}\]

Derivation

  1. Initial program 4.1

    \[\frac{{\left(\frac{1}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-s}}\right)}^{c_n}}{{\left(\frac{1}{1 + e^{-t}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}\]
  2. Initial simplification4.1

    \[\leadsto \frac{{\left(1 - \frac{1}{e^{-s} + 1}\right)}^{c_n}}{{\left(1 - \frac{1}{e^{-t} + 1}\right)}^{c_n}} \cdot \frac{{\left(\frac{1}{e^{-s} + 1}\right)}^{c_p}}{{\left(\frac{1}{e^{-t} + 1}\right)}^{c_p}}\]
  3. Using strategy rm
  4. Applied add-exp-log4.1

    \[\leadsto \frac{{\left(1 - \frac{1}{e^{-s} + 1}\right)}^{c_n}}{{\left(1 - \frac{1}{e^{-t} + 1}\right)}^{c_n}} \cdot \color{blue}{e^{\log \left(\frac{{\left(\frac{1}{e^{-s} + 1}\right)}^{c_p}}{{\left(\frac{1}{e^{-t} + 1}\right)}^{c_p}}\right)}}\]
  5. Applied add-exp-log4.1

    \[\leadsto \frac{{\left(1 - \frac{1}{e^{-s} + 1}\right)}^{c_n}}{\color{blue}{e^{\log \left({\left(1 - \frac{1}{e^{-t} + 1}\right)}^{c_n}\right)}}} \cdot e^{\log \left(\frac{{\left(\frac{1}{e^{-s} + 1}\right)}^{c_p}}{{\left(\frac{1}{e^{-t} + 1}\right)}^{c_p}}\right)}\]
  6. Applied add-exp-log4.1

    \[\leadsto \frac{\color{blue}{e^{\log \left({\left(1 - \frac{1}{e^{-s} + 1}\right)}^{c_n}\right)}}}{e^{\log \left({\left(1 - \frac{1}{e^{-t} + 1}\right)}^{c_n}\right)}} \cdot e^{\log \left(\frac{{\left(\frac{1}{e^{-s} + 1}\right)}^{c_p}}{{\left(\frac{1}{e^{-t} + 1}\right)}^{c_p}}\right)}\]
  7. Applied div-exp4.1

    \[\leadsto \color{blue}{e^{\log \left({\left(1 - \frac{1}{e^{-s} + 1}\right)}^{c_n}\right) - \log \left({\left(1 - \frac{1}{e^{-t} + 1}\right)}^{c_n}\right)}} \cdot e^{\log \left(\frac{{\left(\frac{1}{e^{-s} + 1}\right)}^{c_p}}{{\left(\frac{1}{e^{-t} + 1}\right)}^{c_p}}\right)}\]
  8. Applied prod-exp4.1

    \[\leadsto \color{blue}{e^{\left(\log \left({\left(1 - \frac{1}{e^{-s} + 1}\right)}^{c_n}\right) - \log \left({\left(1 - \frac{1}{e^{-t} + 1}\right)}^{c_n}\right)\right) + \log \left(\frac{{\left(\frac{1}{e^{-s} + 1}\right)}^{c_p}}{{\left(\frac{1}{e^{-t} + 1}\right)}^{c_p}}\right)}}\]
  9. Simplified1.8

    \[\leadsto e^{\color{blue}{(c_n \cdot \left(\log_* (1 + \frac{-1}{e^{-s} + 1})\right) + \left(\log_* (1 + e^{-t}) \cdot c_p\right))_* - (\left(\log_* (1 + \frac{-1}{e^{-t} + 1})\right) \cdot c_n + \left(c_p \cdot \log_* (1 + e^{-s})\right))_*}}\]
  10. Taylor expanded around 0 0.7

    \[\leadsto e^{(c_n \cdot \left(\log_* (1 + \frac{-1}{e^{-s} + 1})\right) + \left(\color{blue}{\left(\left(\log 2 + \frac{1}{8} \cdot {t}^{2}\right) - \frac{1}{2} \cdot t\right)} \cdot c_p\right))_* - (\left(\log_* (1 + \frac{-1}{e^{-t} + 1})\right) \cdot c_n + \left(c_p \cdot \log_* (1 + e^{-s})\right))_*}\]
  11. Simplified0.7

    \[\leadsto e^{(c_n \cdot \left(\log_* (1 + \frac{-1}{e^{-s} + 1})\right) + \left(\color{blue}{(\left((\frac{1}{8} \cdot t + \frac{-1}{2})_*\right) \cdot t + \left(\log 2\right))_*} \cdot c_p\right))_* - (\left(\log_* (1 + \frac{-1}{e^{-t} + 1})\right) \cdot c_n + \left(c_p \cdot \log_* (1 + e^{-s})\right))_*}\]
  12. Taylor expanded around 0 0.8

    \[\leadsto e^{(c_n \cdot \left(\log_* (1 + \frac{-1}{e^{-s} + 1})\right) + \left((\left((\frac{1}{8} \cdot t + \frac{-1}{2})_*\right) \cdot t + \left(\log 2\right))_* \cdot c_p\right))_* - (\color{blue}{\left(\log \frac{1}{2} - \left(\frac{1}{2} \cdot t + \frac{1}{8} \cdot {t}^{2}\right)\right)} \cdot c_n + \left(c_p \cdot \log_* (1 + e^{-s})\right))_*}\]
  13. Simplified0.8

    \[\leadsto e^{(c_n \cdot \left(\log_* (1 + \frac{-1}{e^{-s} + 1})\right) + \left((\left((\frac{1}{8} \cdot t + \frac{-1}{2})_*\right) \cdot t + \left(\log 2\right))_* \cdot c_p\right))_* - (\color{blue}{\left((t \cdot \left((\frac{-1}{8} \cdot t + \frac{-1}{2})_*\right) + \left(\log \frac{1}{2}\right))_*\right)} \cdot c_n + \left(c_p \cdot \log_* (1 + e^{-s})\right))_*}\]
  14. Final simplification0.8

    \[\leadsto e^{(c_n \cdot \left(\log_* (1 + \frac{-1}{1 + e^{-s}})\right) + \left((\left((\frac{1}{8} \cdot t + \frac{-1}{2})_*\right) \cdot t + \left(\log 2\right))_* \cdot c_p\right))_* - (\left((t \cdot \left((\frac{-1}{8} \cdot t + \frac{-1}{2})_*\right) + \left(\log \frac{1}{2}\right))_*\right) \cdot c_n + \left(\log_* (1 + e^{-s}) \cdot c_p\right))_*}\]

Runtime

Time bar (total: 3.5m)Debug logProfile

herbie shell --seed 2018251 +o rules:numerics
(FPCore (c_p c_n t s)
  :name "Harley's example"
  :pre (and (< 0 c_p) (< 0 c_n))

  :herbie-target
  (* (pow (/ (+ 1 (exp (- t))) (+ 1 (exp (- s)))) c_p) (pow (/ (+ 1 (exp t)) (+ 1 (exp s))) c_n))

  (/ (* (pow (/ 1 (+ 1 (exp (- s)))) c_p) (pow (- 1 (/ 1 (+ 1 (exp (- s))))) c_n)) (* (pow (/ 1 (+ 1 (exp (- t)))) c_p) (pow (- 1 (/ 1 (+ 1 (exp (- t))))) c_n))))