Average Error: 39.5 → 0.7
Time: 18.5s
Precision: 64
Internal Precision: 128
\[\cos \left(x + \varepsilon\right) - \cos x\]
\[\begin{array}{l} \mathbf{if}\;\varepsilon \le -3.3020901580082745 \cdot 10^{-05} \lor \neg \left(\varepsilon \le 0.02011113661874671\right):\\ \;\;\;\;\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right) - \cos x\\ \mathbf{else}:\\ \;\;\;\;-2 \cdot \left(\sin \left(\frac{\varepsilon}{2}\right) \cdot \sin \left(\frac{\left(x + x\right) + \varepsilon}{2}\right)\right)\\ \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 eps < -3.3020901580082745e-05 or 0.02011113661874671 < eps

    1. Initial program 29.8

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

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

    if -3.3020901580082745e-05 < eps < 0.02011113661874671

    1. Initial program 49.5

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

      \[\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.6

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

    \[\leadsto \begin{array}{l} \mathbf{if}\;\varepsilon \le -3.3020901580082745 \cdot 10^{-05} \lor \neg \left(\varepsilon \le 0.02011113661874671\right):\\ \;\;\;\;\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right) - \cos x\\ \mathbf{else}:\\ \;\;\;\;-2 \cdot \left(\sin \left(\frac{\varepsilon}{2}\right) \cdot \sin \left(\frac{\left(x + x\right) + \varepsilon}{2}\right)\right)\\ \end{array}\]

Reproduce

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

Details

Time bar (total: 17.0s)Debug log

sample240.0ms

Algorithm
intervals

simplify8.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

8.0ms
(- (cos (+ x eps)) (cos x))

prune10.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 37.0b

localize28.0ms

Local error

Found 2 expressions with local error:

3.4b
(cos (+ x eps))
1.5b
(- (cos (+ x eps)) (cos x))

rewrite457.0ms

Algorithm
rewrite-expression-head
Counts
2 → 40
Calls

2 calls. Slowest were:

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

series81.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

53.0ms
(- (cos (+ x eps)) (cos x))
28.0ms
(cos (+ x eps))

simplify901.0ms

Counts
33 → 46
Calls

33 calls. Slowest were:

265.0ms
(- (* 1/6 (* (pow x 3) eps)) (+ (* 1/2 (pow eps 2)) (* x eps)))
173.0ms
(- 1 (+ (* x eps) (* 1/2 (pow x 2))))
38.0ms
(fma 1 (cos (+ x eps)) (- (* (cos x) 1)))

prune533.0ms

Pruning

8 alts after pruning (8 fresh and 0 done)

Merged error: 0.2b

localize41.0ms

Local error

Found 2 expressions with local error:

3.4b
(sin (/ (+ (+ x x) eps) 2))
0.3b
(* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))

rewrite27.0ms

Algorithm
rewrite-expression-head
Counts
2 → 29
Calls

2 calls. Slowest were:

18.0ms
(* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))
8.0ms
(sin (/ (+ (+ x x) eps) 2))

series127.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

91.0ms
(* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))
36.0ms
(sin (/ (+ (+ x x) eps) 2))

simplify681.0ms

Counts
15 → 35
Calls

15 calls. Slowest were:

344.0ms
(- (+ (* 1/4 (pow eps 2)) (* 1/2 (* x eps))) (* 1/12 (* (pow x 3) eps)))
102.0ms
(- (+ x (* 1/2 eps)) (* 1/6 (pow x 3)))
81.0ms
(* (* (* (sin (/ (+ (+ x x) eps) 2)) (sin (/ (+ (+ x x) eps) 2))) (sin (/ (+ (+ x x) eps) 2))) (* (* (sin (/ eps 2)) (sin (/ eps 2))) (sin (/ eps 2))))

prune442.0ms

Pruning

14 alts after pruning (13 fresh and 1 done)

Merged error: 0.2b

localize9.0ms

Local error

Found 4 expressions with local error:

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

rewrite15.0ms

Algorithm
rewrite-expression-head
Counts
4 → 51
Calls

4 calls. Slowest were:

10.0ms
(* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))
5.0ms
(sin (/ (+ (+ x x) eps) 2))
0.0ms
(expm1 (* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2))))

series471.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

185.0ms
(log1p (expm1 (* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))))
166.0ms
(expm1 (* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2))))
95.0ms
(* (sin (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))
25.0ms
(sin (/ (+ (+ x x) eps) 2))

simplify1.7s

Counts
23 → 63
Calls

23 calls. Slowest were:

445.0ms
(- (+ (* 1/4 (pow eps 2)) (* 1/2 (* x eps))) (* 1/12 (* (pow x 3) eps)))
404.0ms
(- (+ (* 1/4 (pow eps 2)) (* 1/2 (* x eps))) (* 1/12 (* (pow x 3) eps)))
370.0ms
(- (+ (* 1/4 (pow eps 2)) (* 1/2 (* x eps))) (* 1/12 (* (pow x 3) eps)))

prune857.0ms

Pruning

12 alts after pruning (11 fresh and 1 done)

Merged error: 0.2b

localize33.0ms

Local error

Found 4 expressions with local error:

3.4b
(sin (/ (+ (+ x x) eps) 2))
3.4b
(sin (/ (+ (+ x x) eps) 2))
3.4b
(sin (/ (+ (+ x x) eps) 2))
0.6b
(cbrt (sin (/ (+ (+ x x) eps) 2)))

rewrite18.0ms

Algorithm
rewrite-expression-head
Counts
4 → 40
Calls

4 calls. Slowest were:

5.0ms
(sin (/ (+ (+ x x) eps) 2))
5.0ms
(sin (/ (+ (+ x x) eps) 2))
5.0ms
(sin (/ (+ (+ x x) eps) 2))

series146.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

75.0ms
(cbrt (sin (/ (+ (+ x x) eps) 2)))
29.0ms
(sin (/ (+ (+ x x) eps) 2))
21.0ms
(sin (/ (+ (+ x x) eps) 2))
20.0ms
(sin (/ (+ (+ x x) eps) 2))

simplify672.0ms

Counts
15 → 52
Calls

15 calls. Slowest were:

151.0ms
(- (+ (* 1/3 (* x (exp (* 1/3 (- (log 4) (* 2 (log eps))))))) (exp (* 1/3 (+ (log 1/2) (log eps))))) (* 7/81 (* (pow x 2) (exp (* 1/3 (+ (log 1/2) (log eps)))))))
128.0ms
(- (+ x (* 1/2 eps)) (* 1/6 (pow x 3)))
117.0ms
(- (+ x (* 1/2 eps)) (* 1/6 (pow x 3)))

prune753.0ms

Pruning

12 alts after pruning (11 fresh and 1 done)

Merged error: 0.2b

regimes363.0ms

Accuracy

97.8% (0.3b remaining)

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

bsearch245.0ms

end0.0ms

sample8.1s

Algorithm
intervals