Average Error: 39.2 → 0.7
Time: 19.6s
Precision: 64
Internal Precision: 128
\[\cos \left(x + \varepsilon\right) - \cos x\]
\[\begin{array}{l} \mathbf{if}\;\varepsilon \le -0.005208537887765567:\\ \;\;\;\;\cos x \cdot \cos \varepsilon - \left(\cos x + \sin x \cdot \sin \varepsilon\right)\\ \mathbf{elif}\;\varepsilon \le 9.788702847368253 \cdot 10^{-08}:\\ \;\;\;\;(e^{\log_* (1 + \sin \left(\frac{\varepsilon}{2}\right) \cdot \sin \left(\frac{\left(x + x\right) + \varepsilon}{2}\right))} - 1)^* \cdot -2\\ \mathbf{else}:\\ \;\;\;\;\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right) - \cos x\\ \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 3 regimes
  2. if eps < -0.005208537887765567

    1. Initial program 29.9

      \[\cos \left(x + \varepsilon\right) - \cos x\]
    2. Using strategy rm
    3. Applied cos-sum0.8

      \[\leadsto \color{blue}{\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right)} - \cos x\]
    4. Applied associate--l-0.8

      \[\leadsto \color{blue}{\cos x \cdot \cos \varepsilon - \left(\sin x \cdot \sin \varepsilon + \cos x\right)}\]

    if -0.005208537887765567 < eps < 9.788702847368253e-08

    1. Initial program 48.7

      \[\cos \left(x + \varepsilon\right) - \cos x\]
    2. Using strategy rm
    3. Applied diff-cos36.7

      \[\leadsto \color{blue}{-2 \cdot \left(\sin \left(\frac{\left(x + \varepsilon\right) - x}{2}\right) \cdot \sin \left(\frac{\left(x + \varepsilon\right) + x}{2}\right)\right)}\]
    4. Simplified0.5

      \[\leadsto -2 \cdot \color{blue}{\left(\sin \left(\frac{\left(x + x\right) + \varepsilon}{2}\right) \cdot \sin \left(\frac{\varepsilon}{2}\right)\right)}\]
    5. Using strategy rm
    6. Applied expm1-log1p-u0.5

      \[\leadsto -2 \cdot \color{blue}{(e^{\log_* (1 + \sin \left(\frac{\left(x + x\right) + \varepsilon}{2}\right) \cdot \sin \left(\frac{\varepsilon}{2}\right))} - 1)^*}\]

    if 9.788702847368253e-08 < eps

    1. Initial program 30.4

      \[\cos \left(x + \varepsilon\right) - \cos x\]
    2. Using strategy rm
    3. Applied cos-sum1.1

      \[\leadsto \color{blue}{\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right)} - \cos x\]
  3. Recombined 3 regimes into one program.
  4. Final simplification0.7

    \[\leadsto \begin{array}{l} \mathbf{if}\;\varepsilon \le -0.005208537887765567:\\ \;\;\;\;\cos x \cdot \cos \varepsilon - \left(\cos x + \sin x \cdot \sin \varepsilon\right)\\ \mathbf{elif}\;\varepsilon \le 9.788702847368253 \cdot 10^{-08}:\\ \;\;\;\;(e^{\log_* (1 + \sin \left(\frac{\varepsilon}{2}\right) \cdot \sin \left(\frac{\left(x + x\right) + \varepsilon}{2}\right))} - 1)^* \cdot -2\\ \mathbf{else}:\\ \;\;\;\;\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right) - \cos x\\ \end{array}\]

Reproduce

herbie shell --seed 2018365 +o rules:numerics
(FPCore (x eps)
  :name "2cos (problem 3.3.5)"
  (- (cos (+ x eps)) (cos x)))

Details

Time bar (total: 9.5s)Debug log

start258.0ms

Algorithm
intervals

setup23.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 38.5b

localize28.0ms

Local error

Found 2 expressions with local error:

3.2b
(cos (+ x eps))
0.9b
(- (cos (+ x eps)) (cos x))

rewrite15.0ms

Algorithm
rewrite-expression-head
Counts
2 → 40
Calls

2 calls. Slowest were:

10.0ms
(- (cos (+ x eps)) (cos x))
3.0ms
(cos (+ x eps))

series108.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

87.0ms
(- (cos (+ x eps)) (cos x))
21.0ms
(cos (+ x eps))

simplify787.0ms

Counts
33 → 46
Calls

33 calls. Slowest were:

263.0ms
(- (* 1/6 (* (pow x 3) eps)) (+ (* 1/2 (pow eps 2)) (* x eps)))
176.0ms
(- 1 (+ (* x eps) (* 1/2 (pow x 2))))
24.0ms
(fma (- (cbrt (cos x))) (* (cbrt (cos x)) (cbrt (cos x))) (* (cbrt (cos x)) (* (cbrt (cos x)) (cbrt (cos x)))))

prune497.0ms

Pruning

10 alts after pruning (10 fresh and 0 done)

Merged error: 0.4b

localize32.0ms

Local error

Found 3 expressions with local error:

3.1b
(sin (/ (+ (+ x x) eps) 2))
0.3b
(* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))
0.0b
(* -2 (* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2))))

rewrite28.0ms

Algorithm
rewrite-expression-head
Counts
3 → 41
Calls

3 calls. Slowest were:

12.0ms
(* -2 (* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2))))
10.0ms
(* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))
4.0ms
(sin (/ (+ (+ x x) eps) 2))

series198.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

86.0ms
(* -2 (* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2))))
73.0ms
(* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))
38.0ms
(sin (/ (+ (+ x x) eps) 2))

simplify975.0ms

Counts
19 → 50
Calls

19 calls. Slowest were:

378.0ms
(- (+ (* 1/4 (pow eps 2)) (* 1/2 (* x eps))) (* 1/12 (* (pow x 3) eps)))
213.0ms
(- (* 1/6 (* (pow x 3) eps)) (+ (* 1/2 (pow eps 2)) (* x eps)))
114.0ms
(- (+ x (* 1/2 eps)) (* 1/6 (pow x 3)))

prune698.0ms

Pruning

14 alts after pruning (14 fresh and 0 done)

Merged error: 0.4b

localize19.0ms

Local error

Found 4 expressions with local error:

3.1b
(sin (/ (+ (+ x x) eps) 2))
0.3b
(* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))
0.2b
(expm1 (log1p (* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))))
0.1b
(log1p (* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2))))

rewrite29.0ms

Algorithm
rewrite-expression-head
Counts
4 → 51
Calls

4 calls. Slowest were:

19.0ms
(* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))
9.0ms
(sin (/ (+ (+ x x) eps) 2))
1.0ms
(log1p (* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2))))

series510.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

217.0ms
(expm1 (log1p (* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))))
182.0ms
(log1p (* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2))))
80.0ms
(* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))
31.0ms
(sin (/ (+ (+ x x) eps) 2))

simplify1.5s

Counts
23 → 63
Calls

23 calls. Slowest were:

372.0ms
(- (+ (* 1/4 (pow eps 2)) (* 1/2 (* x eps))) (* 1/12 (* (pow x 3) eps)))
371.0ms
(- (+ (* 1/4 (pow eps 2)) (* 1/2 (* x eps))) (* 1/12 (* (pow x 3) eps)))
339.0ms
(- (+ (* 1/4 (pow eps 2)) (* 1/2 (* x eps))) (* 1/12 (* (pow x 3) eps)))

prune925.0ms

Pruning

13 alts after pruning (12 fresh and 1 done)

Merged error: 0.4b

localize11.0ms

Local error

Found 4 expressions with local error:

3.1b
(sin (/ (+ (+ x x) eps) 2))
0.4b
(log1p (sin (/ (+ (+ x x) eps) 2)))
0.3b
(* (expm1 (log1p (sin (/ (+ (+ x x) eps) 2)))) (sin (/ eps 2)))
0.2b
(expm1 (log1p (sin (/ (+ (+ x x) eps) 2))))

rewrite17.0ms

Algorithm
rewrite-expression-head
Counts
4 → 50
Calls

4 calls. Slowest were:

11.0ms
(* (expm1 (log1p (sin (/ (+ (+ x x) eps) 2)))) (sin (/ eps 2)))
4.0ms
(sin (/ (+ (+ x x) eps) 2))
0.0ms
(expm1 (log1p (sin (/ (+ (+ x x) eps) 2))))

series306.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

129.0ms
(* (expm1 (log1p (sin (/ (+ (+ x x) eps) 2)))) (sin (/ eps 2)))
80.0ms
(expm1 (log1p (sin (/ (+ (+ x x) eps) 2))))
61.0ms
(log1p (sin (/ (+ (+ x x) eps) 2)))
36.0ms
(sin (/ (+ (+ x x) eps) 2))

simplify894.0ms

Counts
22 → 62
Calls

22 calls. Slowest were:

351.0ms
(- (+ (* 1/4 (pow eps 2)) (* 1/2 (* x eps))) (* 1/12 (* (pow x 3) eps)))
100.0ms
(- (+ x (* 1/2 eps)) (* 1/6 (pow x 3)))
84.0ms
(- (+ x (* 1/2 eps)) (* 1/6 (pow x 3)))

prune958.0ms

Pruning

12 alts after pruning (11 fresh and 1 done)

Merged error: 0.4b

regimes364.0ms

Accuracy

97.7% (0.3b remaining)

Error of 0.7b against oracle of 0.4b and baseline of 15.1b

bsearch307.0ms