Average Error: 4.2 → 0.9
Time: 2.8m
Precision: 64
Internal Precision: 2368
\[\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}}\]
\[\begin{array}{l} \mathbf{if}\;\frac{1}{1 + e^{-t}} \le 1.3219482130838553 \cdot 10^{-261}:\\ \;\;\;\;\frac{{\left(\frac{1}{e^{-s} + 1}\right)}^{c_p} \cdot {\left(1 - \frac{1}{e^{-s} + 1}\right)}^{c_n}}{e^{\log_* (1 + e^{-t}) \cdot c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}\\ \mathbf{else}:\\ \;\;\;\;e^{(\left(-\log_* (1 + e^{-s})\right) \cdot c_p + \left(\log \left(1 - \frac{1}{e^{-s} + 1}\right) \cdot c_n\right))_* - (\left(-\log_* (1 + e^{-t})\right) \cdot c_p + \left(\log \left(1 - \frac{1}{1 + e^{-t}}\right) \cdot c_n\right))_*}\\ \end{array}\]

Error

Bits error versus c_p

Bits error versus c_n

Bits error versus t

Bits error versus s

Target

Original4.2
Target2.0
Herbie0.9
\[{\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. Split input into 2 regimes
  2. if (/ 1 (+ 1 (exp (- t)))) < 1.3219482130838553e-261

    1. Initial program 62.4

      \[\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. Using strategy rm
    3. Applied add-exp-log62.4

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

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

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

      \[\leadsto \frac{{\left(\frac{1}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-s}}\right)}^{c_n}}{e^{\color{blue}{c_p \cdot \log_* (1 + e^{-t})}} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}\]

    if 1.3219482130838553e-261 < (/ 1 (+ 1 (exp (- t))))

    1. Initial program 3.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. Using strategy rm
    3. Applied add-exp-log3.1

      \[\leadsto \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 {\color{blue}{\left(e^{\log \left(1 - \frac{1}{1 + e^{-t}}\right)}\right)}}^{c_n}}\]
    4. Applied pow-exp3.1

      \[\leadsto \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 \color{blue}{e^{\log \left(1 - \frac{1}{1 + e^{-t}}\right) \cdot c_n}}}\]
    5. Applied add-exp-log3.1

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

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

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

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

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

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

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

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

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

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

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

      \[\leadsto e^{\color{blue}{(\left(-\log_* (1 + e^{-s})\right) \cdot c_p + \left(\log \left(1 - \frac{1}{1 + e^{-s}}\right) \cdot c_n\right))_* - (\left(-\log_* (1 + e^{-t})\right) \cdot c_p + \left(\log \left(1 - \frac{1}{e^{-t} + 1}\right) \cdot c_n\right))_*}}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification0.9

    \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{1}{1 + e^{-t}} \le 1.3219482130838553 \cdot 10^{-261}:\\ \;\;\;\;\frac{{\left(\frac{1}{e^{-s} + 1}\right)}^{c_p} \cdot {\left(1 - \frac{1}{e^{-s} + 1}\right)}^{c_n}}{e^{\log_* (1 + e^{-t}) \cdot c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}\\ \mathbf{else}:\\ \;\;\;\;e^{(\left(-\log_* (1 + e^{-s})\right) \cdot c_p + \left(\log \left(1 - \frac{1}{e^{-s} + 1}\right) \cdot c_n\right))_* - (\left(-\log_* (1 + e^{-t})\right) \cdot c_p + \left(\log \left(1 - \frac{1}{1 + e^{-t}}\right) \cdot c_n\right))_*}\\ \end{array}\]

Runtime

Time bar (total: 2.8m)Debug logProfile

herbie shell --seed 2018214 +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))))