Average Error: 59.1 → 3.0
Time: 26.6s
Precision: 64
Internal Precision: 128
\[\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)}\]
\[\begin{array}{l} \mathbf{if}\;\varepsilon \le -1.7224963786636187 \cdot 10^{-95}:\\ \;\;\;\;\frac{(e^{\left(a + b\right) \cdot \varepsilon} - 1)^*}{(e^{b \cdot \varepsilon} - 1)^*} \cdot \frac{1}{\frac{(e^{a \cdot \varepsilon} - 1)^*}{\varepsilon}}\\ \mathbf{else}:\\ \;\;\;\;\frac{1}{a} + \frac{1}{b}\\ \end{array}\]

Error

Bits error versus a

Bits error versus b

Bits error versus eps

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Target

Original59.1
Target13.8
Herbie3.0
\[\frac{a + b}{a \cdot b}\]

Derivation

  1. Split input into 2 regimes
  2. if eps < -1.7224963786636187e-95

    1. Initial program 53.5

      \[\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)}\]
    2. Simplified8.9

      \[\leadsto \color{blue}{\frac{(e^{\left(a + b\right) \cdot \varepsilon} - 1)^*}{(e^{\varepsilon \cdot b} - 1)^*} \cdot \frac{\varepsilon}{(e^{\varepsilon \cdot a} - 1)^*}}\]
    3. Using strategy rm
    4. Applied clear-num8.9

      \[\leadsto \frac{(e^{\left(a + b\right) \cdot \varepsilon} - 1)^*}{(e^{\varepsilon \cdot b} - 1)^*} \cdot \color{blue}{\frac{1}{\frac{(e^{\varepsilon \cdot a} - 1)^*}{\varepsilon}}}\]

    if -1.7224963786636187e-95 < eps

    1. Initial program 59.9

      \[\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)}\]
    2. Simplified30.7

      \[\leadsto \color{blue}{\frac{(e^{\left(a + b\right) \cdot \varepsilon} - 1)^*}{(e^{\varepsilon \cdot b} - 1)^*} \cdot \frac{\varepsilon}{(e^{\varepsilon \cdot a} - 1)^*}}\]
    3. Taylor expanded around 0 2.2

      \[\leadsto \color{blue}{\frac{1}{a} + \frac{1}{b}}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification3.0

    \[\leadsto \begin{array}{l} \mathbf{if}\;\varepsilon \le -1.7224963786636187 \cdot 10^{-95}:\\ \;\;\;\;\frac{(e^{\left(a + b\right) \cdot \varepsilon} - 1)^*}{(e^{b \cdot \varepsilon} - 1)^*} \cdot \frac{1}{\frac{(e^{a \cdot \varepsilon} - 1)^*}{\varepsilon}}\\ \mathbf{else}:\\ \;\;\;\;\frac{1}{a} + \frac{1}{b}\\ \end{array}\]

Reproduce

herbie shell --seed 2019002 +o rules:numerics
(FPCore (a b eps)
  :name "expq3 (problem 3.4.2)"
  :pre (and (< -1 eps) (< eps 1))

  :herbie-target
  (/ (+ a b) (* a b))

  (/ (* eps (- (exp (* (+ a b) eps)) 1)) (* (- (exp (* a eps)) 1) (- (exp (* b eps)) 1))))

Details

Time bar (total: 26.1s)Debug log

sample545.0ms

Algorithm
intervals

simplify218.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

218.0ms
(/ (* eps (- (exp (* (+ a b) eps)) 1)) (* (- (exp (* a eps)) 1) (- (exp (* b eps)) 1)))

prune22.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 29.7b

localize37.0ms

Local error

Found 4 expressions with local error:

19.5b
(/ eps (expm1 (* eps a)))
13.6b
(/ (expm1 (* (+ a b) eps)) (expm1 (* eps b)))
4.3b
(* (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b))) (/ eps (expm1 (* eps a))))
0.0b
(* (+ a b) eps)

rewrite13.0ms

Algorithm
rewrite-expression-head
Counts
4 → 95
Calls

4 calls. Slowest were:

6.0ms
(* (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b))) (/ eps (expm1 (* eps a))))
3.0ms
(/ (expm1 (* (+ a b) eps)) (expm1 (* eps b)))
2.0ms
(* (+ a b) eps)

series588.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

447.0ms
(* (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b))) (/ eps (expm1 (* eps a))))
85.0ms
(/ (expm1 (* (+ a b) eps)) (expm1 (* eps b)))
32.0ms
(/ eps (expm1 (* eps a)))
24.0ms
(* (+ a b) eps)

simplify2.6s

Counts
70 → 107
Calls

70 calls. Slowest were:

433.0ms
(* (* (* (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b))) (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b)))) (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b)))) (* (* (/ eps (expm1 (* eps a))) (/ eps (expm1 (* eps a)))) (/ eps (expm1 (* eps a)))))
363.0ms
(+ (- (log (expm1 (* (+ a b) eps))) (log (expm1 (* eps b)))) (log (/ eps (expm1 (* eps a)))))
240.0ms
(+ (log (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b)))) (log (/ eps (expm1 (* eps a)))))

prune1.3s

Pruning

5 alts after pruning (5 fresh and 0 done)

Merged error: 0b

localize8.0ms

Local error

Found 1 expressions with local error:

0.0b
(+ (/ 1 a) (/ 1 b))

rewrite4.0ms

Algorithm
rewrite-expression-head
Counts
1 → 22
Calls

1 calls. Slowest were:

4.0ms
(+ (/ 1 a) (/ 1 b))

series7.0ms

Counts
1 → 3
Calls

1 calls. Slowest were:

7.0ms
(+ (/ 1 a) (/ 1 b))

simplify31.0ms

Counts
12 → 25
Calls

12 calls. Slowest were:

5.0ms
(* (exp (/ 1 a)) (exp (/ 1 b)))
4.0ms
(+ (* 1 b) (* a 1))
2.0ms
(+ (/ 1 a) (/ 1 b))

prune160.0ms

Pruning

5 alts after pruning (4 fresh and 1 done)

Merged error: 0b

localize7.0ms

Local error

Found 1 expressions with local error:

13.3b
(/ (+ b a) (* a b))

rewrite3.0ms

Algorithm
rewrite-expression-head
Counts
1 → 23
Calls

1 calls. Slowest were:

3.0ms
(/ (+ b a) (* a b))

series28.0ms

Counts
1 → 3
Calls

1 calls. Slowest were:

28.0ms
(/ (+ b a) (* a b))

simplify145.0ms

Counts
15 → 26
Calls

15 calls. Slowest were:

71.0ms
(/ (* (* (+ b a) (+ b a)) (+ b a)) (* (* (* a b) (* a b)) (* a b)))
29.0ms
(- (log (+ b a)) (log (* a b)))
12.0ms
(- (* b b) (* a a))

prune217.0ms

Pruning

4 alts after pruning (3 fresh and 1 done)

Merged error: 0b

localize11.0ms

Local error

Found 2 expressions with local error:

8.0b
(/ (/ (+ b a) a) b)
0.0b
(/ (+ b a) a)

rewrite14.0ms

Algorithm
rewrite-expression-head
Counts
2 → 34
Calls

2 calls. Slowest were:

9.0ms
(/ (/ (+ b a) a) b)
5.0ms
(/ (+ b a) a)

series48.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

35.0ms
(/ (/ (+ b a) a) b)
12.0ms
(/ (+ b a) a)

simplify47.0ms

Counts
14 → 40
Calls

14 calls. Slowest were:

12.0ms
(- (* b b) (* a a))
5.0ms
(+ (pow b 3) (pow a 3))
4.0ms
(+ 1 (/ b a))

prune278.0ms

Pruning

4 alts after pruning (2 fresh and 2 done)

Merged error: 0b

regimes95.0ms

Accuracy

-3.9% (3.0b remaining)

Error of 3.0b against oracle of 0.0b and baseline of 2.9b

bsearch403.0ms

end0.0ms

sample19.2s

Algorithm
intervals