\[\frac{\varepsilon \cdot \left(e^{\left(a + b\right) \cdot \varepsilon} - 1\right)}{\left(e^{a \cdot \varepsilon} - 1\right) \cdot \left(e^{b \cdot \varepsilon} - 1\right)}\]
Test:
NMSE problem 3.4.2
Bits:
128 bits
Bits error versus a
Bits error versus b
Bits error versus eps
Time: 56.4 s
Input Error: 62.0
Output Error: 0.0
Log: ⚲
Profile: 🕒
\(\begin{cases} \frac{1}{b} + \frac{1}{a} & \text{when } \frac{\varepsilon \cdot \left(e^{\left(a + b\right) \cdot \varepsilon} - 1\right)}{\left(e^{a \cdot \varepsilon} - 1\right) \cdot \left(e^{b \cdot \varepsilon} - 1\right)} \le -7.102832082435936 \cdot 10^{-280} \\ \frac{1}{b} + \frac{1}{a} & \text{when } \frac{\varepsilon \cdot \left(e^{\left(a + b\right) \cdot \varepsilon} - 1\right)}{\left(e^{a \cdot \varepsilon} - 1\right) \cdot \left(e^{b \cdot \varepsilon} - 1\right)} \le 1.4168119176958962 \cdot 10^{-79} \\ \frac{1}{b} + \frac{1}{a} & \text{otherwise} \end{cases}\)

    if (/ (* eps (- (exp (* (+ a b) eps)) 1)) (* (- (exp (* a eps)) 1) (- (exp (* b eps)) 1))) < -7.102832082435936e-280 or 1.4168119176958962e-79 < (/ (* eps (- (exp (* (+ a b) eps)) 1)) (* (- (exp (* a eps)) 1) (- (exp (* b eps)) 1)))

    1. Started with
      \[\frac{\varepsilon \cdot \left(e^{\left(a + b\right) \cdot \varepsilon} - 1\right)}{\left(e^{a \cdot \varepsilon} - 1\right) \cdot \left(e^{b \cdot \varepsilon} - 1\right)}\]
      61.6
    2. Applied taylor to get
      \[\frac{\varepsilon \cdot \left(e^{\left(a + b\right) \cdot \varepsilon} - 1\right)}{\left(e^{a \cdot \varepsilon} - 1\right) \cdot \left(e^{b \cdot \varepsilon} - 1\right)} \leadsto \frac{1}{b} + \frac{1}{a}\]
      0.0
    3. Taylor expanded around 0 to get
      \[\color{red}{\frac{1}{b} + \frac{1}{a}} \leadsto \color{blue}{\frac{1}{b} + \frac{1}{a}}\]
      0.0

    if -7.102832082435936e-280 < (/ (* eps (- (exp (* (+ a b) eps)) 1)) (* (- (exp (* a eps)) 1) (- (exp (* b eps)) 1))) < 1.4168119176958962e-79

    1. Started with
      \[\frac{\varepsilon \cdot \left(e^{\left(a + b\right) \cdot \varepsilon} - 1\right)}{\left(e^{a \cdot \varepsilon} - 1\right) \cdot \left(e^{b \cdot \varepsilon} - 1\right)}\]
      62.1
    2. Applied taylor to get
      \[\frac{\varepsilon \cdot \left(e^{\left(a + b\right) \cdot \varepsilon} - 1\right)}{\left(e^{a \cdot \varepsilon} - 1\right) \cdot \left(e^{b \cdot \varepsilon} - 1\right)} \leadsto \frac{1}{b} + \frac{1}{a}\]
      0.0
    3. Taylor expanded around 0 to get
      \[\color{red}{\frac{1}{b} + \frac{1}{a}} \leadsto \color{blue}{\frac{1}{b} + \frac{1}{a}}\]
      0.0
  1. Removed slow pow expressions

Original test:


(lambda ((a default) (b default) (eps default))
  #:name "NMSE problem 3.4.2"
  (/ (* eps (- (exp (* (+ a b) eps)) 1)) (* (- (exp (* a eps)) 1) (- (exp (* b eps)) 1)))
  #:target
  (/ (+ a b) (* a b)))