Average Error: 36.5 → 0.6
Time: 15.4s
Precision: 64
Internal Precision: 128
\[\sin \left(x + \varepsilon\right) - \sin x\]
\[\begin{array}{l} \mathbf{if}\;\varepsilon \le -25.400968012366302:\\ \;\;\;\;\left(\cos x \cdot \sin \varepsilon - \sin x\right) + \sin x \cdot \cos \varepsilon\\ \mathbf{elif}\;\varepsilon \le 9.406785545346194 \cdot 10^{-09}:\\ \;\;\;\;2 \cdot \left(\sin \left(\frac{\varepsilon}{2}\right) \cdot \cos \left(\frac{\left(x + x\right) + \varepsilon}{2}\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\sin x \cdot \cos \varepsilon + \cos x \cdot \sin \varepsilon\right) - \sin 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

Target

Original36.5
Target15.1
Herbie0.6
\[2 \cdot \left(\cos \left(x + \frac{\varepsilon}{2}\right) \cdot \sin \left(\frac{\varepsilon}{2}\right)\right)\]

Derivation

  1. Split input into 3 regimes
  2. if eps < -25.400968012366302

    1. Initial program 29.8

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

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

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

    if -25.400968012366302 < eps < 9.406785545346194e-09

    1. Initial program 44.2

      \[\sin \left(x + \varepsilon\right) - \sin x\]
    2. Using strategy rm
    3. Applied diff-sin44.2

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

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

    if 9.406785545346194e-09 < eps

    1. Initial program 28.7

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

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

    \[\leadsto \begin{array}{l} \mathbf{if}\;\varepsilon \le -25.400968012366302:\\ \;\;\;\;\left(\cos x \cdot \sin \varepsilon - \sin x\right) + \sin x \cdot \cos \varepsilon\\ \mathbf{elif}\;\varepsilon \le 9.406785545346194 \cdot 10^{-09}:\\ \;\;\;\;2 \cdot \left(\sin \left(\frac{\varepsilon}{2}\right) \cdot \cos \left(\frac{\left(x + x\right) + \varepsilon}{2}\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\sin x \cdot \cos \varepsilon + \cos x \cdot \sin \varepsilon\right) - \sin x\\ \end{array}\]

Reproduce

herbie shell --seed 2019005 +o rules:numerics
(FPCore (x eps)
  :name "2sin (example 3.3)"

  :herbie-target
  (* 2 (* (cos (+ x (/ eps 2))) (sin (/ eps 2))))

  (- (sin (+ x eps)) (sin x)))

Details

Time bar (total: 14.4s)Debug log

sample208.0ms

Algorithm
intervals

simplify4.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

4.0ms
(- (sin (+ x eps)) (sin x))

prune41.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 36.5b

localize16.0ms

Local error

Found 2 expressions with local error:

3.8b
(sin (+ x eps))
1.8b
(- (sin (+ x eps)) (sin x))

rewrite9.0ms

Algorithm
rewrite-expression-head
Counts
2 → 40
Calls

2 calls. Slowest were:

6.0ms
(- (sin (+ x eps)) (sin x))
2.0ms
(sin (+ x eps))

series60.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

44.0ms
(- (sin (+ x eps)) (sin x))
16.0ms
(sin (+ x eps))

simplify914.0ms

Counts
33 → 46
Calls

33 calls. Slowest were:

407.0ms
(- eps (+ (* 1/2 (* x (pow eps 2))) (* 1/2 (* (pow x 2) eps))))
75.0ms
(- (+ x eps) (* 1/6 (pow x 3)))
40.0ms
(* (sin (/ (- (+ x eps) x) 2)) (cos (/ (+ (+ x eps) x) 2)))

prune398.0ms

Pruning

4 alts after pruning (4 fresh and 0 done)

Merged error: 0.3b

localize29.0ms

Local error

Found 2 expressions with local error:

3.9b
(cos (/ (+ (+ x x) eps) 2))
0.2b
(* (cos (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))

rewrite21.0ms

Algorithm
rewrite-expression-head
Counts
2 → 28
Calls

2 calls. Slowest were:

15.0ms
(* (cos (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))
5.0ms
(cos (/ (+ (+ x x) eps) 2))

series76.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

55.0ms
(* (cos (/ (+ (+ x x) eps) 2)) (sin (/ eps 2)))
20.0ms
(cos (/ (+ (+ x x) eps) 2))

simplify550.0ms

Counts
14 → 34
Calls

14 calls. Slowest were:

288.0ms
(- 1 (+ (* 1/2 (* x eps)) (* 1/2 (pow x 2))))
86.0ms
(- (* 1/2 eps) (+ (* 1/4 (* x (pow eps 2))) (* 1/4 (* (pow x 2) eps))))
77.0ms
(* (* (* (cos (/ (+ (+ x x) eps) 2)) (cos (/ (+ (+ x x) eps) 2))) (cos (/ (+ (+ x x) eps) 2))) (* (* (sin (/ eps 2)) (sin (/ eps 2))) (sin (/ eps 2))))

prune381.0ms

Pruning

8 alts after pruning (7 fresh and 1 done)

Merged error: 0.3b

localize13.0ms

Local error

Found 4 expressions with local error:

3.9b
(cos (/ (+ (+ x x) eps) 2))
3.9b
(cos (/ (+ (+ x x) eps) 2))
3.9b
(cos (/ (+ (+ x x) eps) 2))
0.6b
(cbrt (* (* (cos (/ (+ (+ x x) eps) 2)) (cos (/ (+ (+ x x) eps) 2))) (cos (/ (+ (+ x x) eps) 2))))

rewrite21.0ms

Algorithm
rewrite-expression-head
Counts
4 → 40
Calls

4 calls. Slowest were:

7.0ms
(cbrt (* (* (cos (/ (+ (+ x x) eps) 2)) (cos (/ (+ (+ x x) eps) 2))) (cos (/ (+ (+ x x) eps) 2))))
5.0ms
(cos (/ (+ (+ x x) eps) 2))
4.0ms
(cos (/ (+ (+ x x) eps) 2))

series79.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

23.0ms
(cos (/ (+ (+ x x) eps) 2))
21.0ms
(cos (/ (+ (+ x x) eps) 2))
17.0ms
(cos (/ (+ (+ x x) eps) 2))
17.0ms
(cbrt (* (* (cos (/ (+ (+ x x) eps) 2)) (cos (/ (+ (+ x x) eps) 2))) (cos (/ (+ (+ x x) eps) 2))))

simplify1.7s

Counts
15 → 52
Calls

15 calls. Slowest were:

390.0ms
(- 1 (+ (* 1/2 (* x eps)) (* 1/2 (pow x 2))))
370.0ms
(- 1 (+ (* 1/2 (* x eps)) (* 1/2 (pow x 2))))
345.0ms
(- 1 (+ (* 1/2 (* x eps)) (* 1/2 (pow x 2))))

prune567.0ms

Pruning

8 alts after pruning (6 fresh and 2 done)

Merged error: 0.3b

localize32.0ms

Local error

Found 4 expressions with local error:

3.9b
(cos (/ (+ (+ x x) eps) 2))
3.9b
(cos (/ (+ (+ x x) eps) 2))
3.9b
(cos (/ (+ (+ x x) eps) 2))
0.5b
(cbrt (cos (/ (+ (+ x x) eps) 2)))

rewrite20.0ms

Algorithm
rewrite-expression-head
Counts
4 → 40
Calls

4 calls. Slowest were:

10.0ms
(cos (/ (+ (+ x x) eps) 2))
5.0ms
(cos (/ (+ (+ x x) eps) 2))
4.0ms
(cos (/ (+ (+ x x) eps) 2))

series142.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

75.0ms
(cbrt (cos (/ (+ (+ x x) eps) 2)))
27.0ms
(cos (/ (+ (+ x x) eps) 2))
24.0ms
(cos (/ (+ (+ x x) eps) 2))
17.0ms
(cos (/ (+ (+ x x) eps) 2))

simplify1.4s

Counts
15 → 52
Calls

15 calls. Slowest were:

304.0ms
(- 1 (+ (* 1/6 (* x eps)) (* 1/6 (pow x 2))))
300.0ms
(- 1 (+ (* 1/2 (* x eps)) (* 1/2 (pow x 2))))
300.0ms
(- 1 (+ (* 1/2 (* x eps)) (* 1/2 (pow x 2))))

prune598.0ms

Pruning

8 alts after pruning (5 fresh and 3 done)

Merged error: 0.3b

regimes249.0ms

Accuracy

98.5% (0.2b remaining)

Error of 0.6b against oracle of 0.3b and baseline of 15.1b

bsearch210.0ms

end0.0ms

sample6.7s

Algorithm
intervals