Average Error: 29.7 → 1.0
Time: 1.3m
Precision: 64
Internal Precision: 128
\[\frac{\left(1 + \frac{1}{\varepsilon}\right) \cdot e^{-\left(1 - \varepsilon\right) \cdot x} - \left(\frac{1}{\varepsilon} - 1\right) \cdot e^{-\left(1 + \varepsilon\right) \cdot x}}{2}\]
\[\begin{array}{l} \mathbf{if}\;x \le 324.30491579076437:\\ \;\;\;\;\frac{\left(2 + \frac{2}{3} \cdot {x}^{3}\right) - {x}^{2}}{2}\\ \mathbf{else}:\\ \;\;\;\;\frac{\left(\frac{1}{\varepsilon} + 1\right) \cdot e^{\left(1 - \varepsilon\right) \cdot \left(-x\right)} - \left(\frac{1}{\varepsilon} - 1\right) \cdot {\left(e^{-1}\right)}^{\left(x \cdot \left(1 + \varepsilon\right)\right)}}{2}\\ \end{array}\]

Error

Bits error versus x

Bits error versus eps

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 2 regimes
  2. if x < 324.30491579076437

    1. Initial program 39.3

      \[\frac{\left(1 + \frac{1}{\varepsilon}\right) \cdot e^{-\left(1 - \varepsilon\right) \cdot x} - \left(\frac{1}{\varepsilon} - 1\right) \cdot e^{-\left(1 + \varepsilon\right) \cdot x}}{2}\]
    2. Taylor expanded around 0 1.4

      \[\leadsto \frac{\color{blue}{\left(\frac{2}{3} \cdot {x}^{3} + 2\right) - {x}^{2}}}{2}\]
    3. Using strategy rm
    4. Applied add-exp-log1.4

      \[\leadsto \frac{\color{blue}{e^{\log \left(\left(\frac{2}{3} \cdot {x}^{3} + 2\right) - {x}^{2}\right)}}}{2}\]
    5. Taylor expanded around -inf 1.4

      \[\leadsto \frac{\color{blue}{\left(\frac{2}{3} \cdot {x}^{3} + 2\right) - {x}^{2}}}{2}\]

    if 324.30491579076437 < x

    1. Initial program 0

      \[\frac{\left(1 + \frac{1}{\varepsilon}\right) \cdot e^{-\left(1 - \varepsilon\right) \cdot x} - \left(\frac{1}{\varepsilon} - 1\right) \cdot e^{-\left(1 + \varepsilon\right) \cdot x}}{2}\]
    2. Using strategy rm
    3. Applied neg-mul-10

      \[\leadsto \frac{\left(1 + \frac{1}{\varepsilon}\right) \cdot e^{-\left(1 - \varepsilon\right) \cdot x} - \left(\frac{1}{\varepsilon} - 1\right) \cdot e^{\color{blue}{-1 \cdot \left(\left(1 + \varepsilon\right) \cdot x\right)}}}{2}\]
    4. Applied exp-prod0

      \[\leadsto \frac{\left(1 + \frac{1}{\varepsilon}\right) \cdot e^{-\left(1 - \varepsilon\right) \cdot x} - \left(\frac{1}{\varepsilon} - 1\right) \cdot \color{blue}{{\left(e^{-1}\right)}^{\left(\left(1 + \varepsilon\right) \cdot x\right)}}}{2}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification1.0

    \[\leadsto \begin{array}{l} \mathbf{if}\;x \le 324.30491579076437:\\ \;\;\;\;\frac{\left(2 + \frac{2}{3} \cdot {x}^{3}\right) - {x}^{2}}{2}\\ \mathbf{else}:\\ \;\;\;\;\frac{\left(\frac{1}{\varepsilon} + 1\right) \cdot e^{\left(1 - \varepsilon\right) \cdot \left(-x\right)} - \left(\frac{1}{\varepsilon} - 1\right) \cdot {\left(e^{-1}\right)}^{\left(x \cdot \left(1 + \varepsilon\right)\right)}}{2}\\ \end{array}\]

Reproduce

herbie shell --seed 2018360 
(FPCore (x eps)
  :name "NMSE Section 6.1 mentioned, A"
  (/ (- (* (+ 1 (/ 1 eps)) (exp (- (* (- 1 eps) x)))) (* (- (/ 1 eps) 1) (exp (- (* (+ 1 eps) x))))) 2))

Details

Time bar (total: 1.0m)Debug log

start519.0ms

Algorithm
intervals

setup373.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 29.3b

localize38.0ms

Local error

Found 4 expressions with local error:

1.9b
(- (* (+ 1 (/ 1 eps)) (exp (- (* (- 1 eps) x)))) (* (- (/ 1 eps) 1) (exp (- (* (+ 1 eps) x)))))
0.0b
(exp (- (* (+ 1 eps) x)))
0.0b
(exp (- (* (- 1 eps) x)))
0.0b
(* (+ 1 (/ 1 eps)) (exp (- (* (- 1 eps) x))))

rewrite47.0ms

Algorithm
rewrite-expression-head
Counts
4 → 153
Calls

4 calls. Slowest were:

24.0ms
(- (* (+ 1 (/ 1 eps)) (exp (- (* (- 1 eps) x)))) (* (- (/ 1 eps) 1) (exp (- (* (+ 1 eps) x)))))
11.0ms
(* (+ 1 (/ 1 eps)) (exp (- (* (- 1 eps) x))))
5.0ms
(exp (- (* (+ 1 eps) x)))

series187.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

109.0ms
(- (* (+ 1 (/ 1 eps)) (exp (- (* (- 1 eps) x)))) (* (- (/ 1 eps) 1) (exp (- (* (+ 1 eps) x)))))
30.0ms
(* (+ 1 (/ 1 eps)) (exp (- (* (- 1 eps) x))))
26.0ms
(exp (- (* (- 1 eps) x)))
23.0ms
(exp (- (* (+ 1 eps) x)))

simplify49.3s

Counts
212 → 165
Calls

212 calls. Slowest were:

902.0ms
(- (* (* (- (* 1 1) (* (/ 1 eps) (/ 1 eps))) (exp (- (* (- 1 eps) x)))) (* (+ (/ 1 eps) 1) (exp (* (+ 1 eps) x)))) (* (- 1 (/ 1 eps)) (* (- (* (/ 1 eps) (/ 1 eps)) (* 1 1)) (exp 0))))
815.0ms
(- (* (* (+ (pow 1 3) (pow (/ 1 eps) 3)) 1) (* (+ (/ 1 eps) 1) (exp (* (+ 1 eps) x)))) (* (* (+ (* 1 1) (- (* (/ 1 eps) (/ 1 eps)) (* 1 (/ 1 eps)))) (exp (* (- 1 eps) x))) (* (- (* (/ 1 eps) (/ 1 eps)) (* 1 1)) (exp 0))))
784.0ms
(- (* (* (+ (pow 1 3) (pow (/ 1 eps) 3)) (exp (- (* (- 1 eps) x)))) (* (+ (/ 1 eps) 1) (exp (* (+ 1 eps) x)))) (* (+ (* 1 1) (- (* (/ 1 eps) (/ 1 eps)) (* 1 (/ 1 eps)))) (* (- (* (/ 1 eps) (/ 1 eps)) (* 1 1)) (exp 0))))

prune2.9s

Pruning

4 alts after pruning (4 fresh and 0 done)

Merged error: 0.3b

localize27.0ms

Local error

Found 2 expressions with local error:

7.0b
(- (+ (* 2/3 (pow x 3)) 2) (pow x 2))
0.1b
(* 2/3 (pow x 3))

rewrite13.0ms

Algorithm
rewrite-expression-head
Counts
2 → 29
Calls

2 calls. Slowest were:

11.0ms
(- (+ (* 2/3 (pow x 3)) 2) (pow x 2))
1.0ms
(* 2/3 (pow x 3))

series42.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

26.0ms
(* 2/3 (pow x 3))
16.0ms
(- (+ (* 2/3 (pow x 3)) 2) (pow x 2))

simplify412.0ms

Counts
17 → 35
Calls

17 calls. Slowest were:

88.0ms
(/ (exp (+ (* 2/3 (pow x 3)) 2)) (exp (pow x 2)))
59.0ms
(- (+ (* 2/3 (pow x 3)) 2) (pow x 2))
57.0ms
(- (+ (* 2/3 (pow x 3)) 2) (pow x 2))

prune542.0ms

Pruning

5 alts after pruning (5 fresh and 0 done)

Merged error: 0.3b

localize20.0ms

Local error

Found 4 expressions with local error:

10.7b
(log (- (+ (* 2/3 (pow x 3)) 2) (pow x 2)))
7.0b
(- (+ (* 2/3 (pow x 3)) 2) (pow x 2))
0.5b
(exp (log (- (+ (* 2/3 (pow x 3)) 2) (pow x 2))))
0.1b
(* 2/3 (pow x 3))

rewrite58.0ms

Algorithm
rewrite-expression-head
Counts
4 → 66
Calls

4 calls. Slowest were:

22.0ms
(- (+ (* 2/3 (pow x 3)) 2) (pow x 2))
16.0ms
(exp (log (- (+ (* 2/3 (pow x 3)) 2) (pow x 2))))
15.0ms
(log (- (+ (* 2/3 (pow x 3)) 2) (pow x 2)))

series119.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

47.0ms
(log (- (+ (* 2/3 (pow x 3)) 2) (pow x 2)))
27.0ms
(* 2/3 (pow x 3))
25.0ms
(exp (log (- (+ (* 2/3 (pow x 3)) 2) (pow x 2))))
20.0ms
(- (+ (* 2/3 (pow x 3)) 2) (pow x 2))

simplify2.7s

Counts
45 → 78
Calls

45 calls. Slowest were:

323.0ms
(- (log 2/3) (+ (* 3 (log (/ 1 x))) (+ (* 9/8 (/ 1 (pow x 2))) (* 3/2 (/ 1 x)))))
245.0ms
(exp (log (- (pow (+ (* 2/3 (pow x 3)) 2) 3) (pow (pow x 2) 3))))
206.0ms
(exp (log (- (* (+ (* 2/3 (pow x 3)) 2) (+ (* 2/3 (pow x 3)) 2)) (* (pow x 2) (pow x 2)))))

prune1.3s

Pruning

6 alts after pruning (6 fresh and 0 done)

Merged error: 0.3b

localize1.0ms

Local error

Found 2 expressions with local error:

7.0b
(- (+ (* 2/3 (pow x 3)) 2) (pow x 2))
0.1b
(* 2/3 (pow x 3))

rewrite26.0ms

Algorithm
rewrite-expression-head
Counts
2 → 29
Calls

2 calls. Slowest were:

23.0ms
(- (+ (* 2/3 (pow x 3)) 2) (pow x 2))
2.0ms
(* 2/3 (pow x 3))

series56.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

29.0ms
(- (+ (* 2/3 (pow x 3)) 2) (pow x 2))
27.0ms
(* 2/3 (pow x 3))

simplify499.0ms

Counts
17 → 35
Calls

17 calls. Slowest were:

93.0ms
(/ (exp (+ (* 2/3 (pow x 3)) 2)) (exp (pow x 2)))
84.0ms
(- (+ (* 2/3 (pow x 3)) 2) (pow x 2))
76.0ms
(- (+ (* 2/3 (pow x 3)) 2) (pow x 2))

prune559.0ms

Pruning

6 alts after pruning (5 fresh and 1 done)

Merged error: 0.3b

regimes140.0ms

Accuracy

96.2% (0.6b remaining)

Error of 1.0b against oracle of 0.4b and baseline of 16.3b

bsearch376.0ms