Average Error: 37.5 → 15.9
Time: 49.6s
Precision: 64
Internal Precision: 128
\[\tan \left(x + \varepsilon\right) - \tan x\]
\[\begin{array}{l} \mathbf{if}\;\varepsilon \le -2.388080272206088 \cdot 10^{-71} \lor \neg \left(\varepsilon \le 4.5316032501739245 \cdot 10^{-45}\right):\\ \;\;\;\;(\left(\frac{\tan \varepsilon + \tan x}{1 - \left(\tan \varepsilon \cdot \left(\tan \varepsilon \cdot \tan \varepsilon\right)\right) \cdot {\left(\tan x\right)}^{3}}\right) \cdot \left(\left(\tan x \cdot \tan \varepsilon + \left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right)\right) + 1\right) + \left(-\tan x\right))_*\\ \mathbf{else}:\\ \;\;\;\;(\left(x \cdot \varepsilon\right) \cdot \left(x + \varepsilon\right) + \varepsilon)_*\\ \end{array}\]

Error

Bits error versus x

Bits error versus eps

Target

Original37.5
Target15.1
Herbie15.9
\[\frac{\sin \varepsilon}{\cos x \cdot \cos \left(x + \varepsilon\right)}\]

Derivation

  1. Split input into 2 regimes
  2. if eps < -2.388080272206088e-71 or 4.5316032501739245e-45 < eps

    1. Initial program 30.7

      \[\tan \left(x + \varepsilon\right) - \tan x\]
    2. Using strategy rm
    3. Applied tan-sum4.7

      \[\leadsto \color{blue}{\frac{\tan x + \tan \varepsilon}{1 - \tan x \cdot \tan \varepsilon}} - \tan x\]
    4. Using strategy rm
    5. Applied flip3--4.7

      \[\leadsto \frac{\tan x + \tan \varepsilon}{\color{blue}{\frac{{1}^{3} - {\left(\tan x \cdot \tan \varepsilon\right)}^{3}}{1 \cdot 1 + \left(\left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right) + 1 \cdot \left(\tan x \cdot \tan \varepsilon\right)\right)}}} - \tan x\]
    6. Applied associate-/r/4.7

      \[\leadsto \color{blue}{\frac{\tan x + \tan \varepsilon}{{1}^{3} - {\left(\tan x \cdot \tan \varepsilon\right)}^{3}} \cdot \left(1 \cdot 1 + \left(\left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right) + 1 \cdot \left(\tan x \cdot \tan \varepsilon\right)\right)\right)} - \tan x\]
    7. Applied fma-neg4.7

      \[\leadsto \color{blue}{(\left(\frac{\tan x + \tan \varepsilon}{{1}^{3} - {\left(\tan x \cdot \tan \varepsilon\right)}^{3}}\right) \cdot \left(1 \cdot 1 + \left(\left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right) + 1 \cdot \left(\tan x \cdot \tan \varepsilon\right)\right)\right) + \left(-\tan x\right))_*}\]
    8. Using strategy rm
    9. Applied add-cbrt-cube4.8

      \[\leadsto (\left(\frac{\tan x + \tan \varepsilon}{{1}^{3} - {\left(\tan x \cdot \color{blue}{\sqrt[3]{\left(\tan \varepsilon \cdot \tan \varepsilon\right) \cdot \tan \varepsilon}}\right)}^{3}}\right) \cdot \left(1 \cdot 1 + \left(\left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right) + 1 \cdot \left(\tan x \cdot \tan \varepsilon\right)\right)\right) + \left(-\tan x\right))_*\]
    10. Applied add-cbrt-cube4.8

      \[\leadsto (\left(\frac{\tan x + \tan \varepsilon}{{1}^{3} - {\left(\color{blue}{\sqrt[3]{\left(\tan x \cdot \tan x\right) \cdot \tan x}} \cdot \sqrt[3]{\left(\tan \varepsilon \cdot \tan \varepsilon\right) \cdot \tan \varepsilon}\right)}^{3}}\right) \cdot \left(1 \cdot 1 + \left(\left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right) + 1 \cdot \left(\tan x \cdot \tan \varepsilon\right)\right)\right) + \left(-\tan x\right))_*\]
    11. Applied cbrt-unprod4.8

      \[\leadsto (\left(\frac{\tan x + \tan \varepsilon}{{1}^{3} - {\color{blue}{\left(\sqrt[3]{\left(\left(\tan x \cdot \tan x\right) \cdot \tan x\right) \cdot \left(\left(\tan \varepsilon \cdot \tan \varepsilon\right) \cdot \tan \varepsilon\right)}\right)}}^{3}}\right) \cdot \left(1 \cdot 1 + \left(\left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right) + 1 \cdot \left(\tan x \cdot \tan \varepsilon\right)\right)\right) + \left(-\tan x\right))_*\]
    12. Applied rem-cube-cbrt4.7

      \[\leadsto (\left(\frac{\tan x + \tan \varepsilon}{{1}^{3} - \color{blue}{\left(\left(\tan x \cdot \tan x\right) \cdot \tan x\right) \cdot \left(\left(\tan \varepsilon \cdot \tan \varepsilon\right) \cdot \tan \varepsilon\right)}}\right) \cdot \left(1 \cdot 1 + \left(\left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right) + 1 \cdot \left(\tan x \cdot \tan \varepsilon\right)\right)\right) + \left(-\tan x\right))_*\]
    13. Simplified4.7

      \[\leadsto (\left(\frac{\tan x + \tan \varepsilon}{{1}^{3} - \color{blue}{{\left(\tan x\right)}^{3}} \cdot \left(\left(\tan \varepsilon \cdot \tan \varepsilon\right) \cdot \tan \varepsilon\right)}\right) \cdot \left(1 \cdot 1 + \left(\left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right) + 1 \cdot \left(\tan x \cdot \tan \varepsilon\right)\right)\right) + \left(-\tan x\right))_*\]

    if -2.388080272206088e-71 < eps < 4.5316032501739245e-45

    1. Initial program 47.0

      \[\tan \left(x + \varepsilon\right) - \tan x\]
    2. Using strategy rm
    3. Applied tan-sum47.0

      \[\leadsto \color{blue}{\frac{\tan x + \tan \varepsilon}{1 - \tan x \cdot \tan \varepsilon}} - \tan x\]
    4. Taylor expanded around 0 31.6

      \[\leadsto \color{blue}{x \cdot {\varepsilon}^{2} + \left(\varepsilon + {x}^{2} \cdot \varepsilon\right)}\]
    5. Simplified31.6

      \[\leadsto \color{blue}{(\left(x \cdot \varepsilon\right) \cdot \left(\varepsilon + x\right) + \varepsilon)_*}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification15.9

    \[\leadsto \begin{array}{l} \mathbf{if}\;\varepsilon \le -2.388080272206088 \cdot 10^{-71} \lor \neg \left(\varepsilon \le 4.5316032501739245 \cdot 10^{-45}\right):\\ \;\;\;\;(\left(\frac{\tan \varepsilon + \tan x}{1 - \left(\tan \varepsilon \cdot \left(\tan \varepsilon \cdot \tan \varepsilon\right)\right) \cdot {\left(\tan x\right)}^{3}}\right) \cdot \left(\left(\tan x \cdot \tan \varepsilon + \left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right)\right) + 1\right) + \left(-\tan x\right))_*\\ \mathbf{else}:\\ \;\;\;\;(\left(x \cdot \varepsilon\right) \cdot \left(x + \varepsilon\right) + \varepsilon)_*\\ \end{array}\]

Reproduce

herbie shell --seed 2019004 +o rules:numerics
(FPCore (x eps)
  :name "2tan (problem 3.3.2)"

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

  (- (tan (+ x eps)) (tan x)))

Details

Time bar (total: 48.2s)Debug log

sample446.0ms

Algorithm
intervals

simplify8.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

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

prune10.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 38.0b

localize14.0ms

Local error

Found 2 expressions with local error:

3.6b
(tan (+ x eps))
0.8b
(- (tan (+ x eps)) (tan x))

rewrite8.0ms

Algorithm
rewrite-expression-head
Counts
2 → 41
Calls

2 calls. Slowest were:

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

series160.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

88.0ms
(- (tan (+ x eps)) (tan x))
72.0ms
(tan (+ x eps))

simplify1.1s

Counts
35 → 47
Calls

35 calls. Slowest were:

240.0ms
(- (* (+ (tan x) (tan eps)) (cos x)) (* (- 1 (* (tan x) (tan eps))) (sin x)))
196.0ms
(+ (* x (pow eps 2)) (+ eps (* (pow x 2) eps)))
155.0ms
(* (- 1 (* (tan x) (tan eps))) (cos x))

prune414.0ms

Pruning

9 alts after pruning (9 fresh and 0 done)

Merged error: 14.4b

localize17.0ms

Local error

Found 4 expressions with local error:

1.8b
(- (/ (+ (tan x) (tan eps)) (- 1 (* (tan x) (tan eps)))) (tan x))
0.3b
(* (tan x) (tan eps))
0.2b
(+ (tan x) (tan eps))
0.1b
(/ (+ (tan x) (tan eps)) (- 1 (* (tan x) (tan eps))))

rewrite58.0ms

Algorithm
rewrite-expression-head
Counts
4 → 168
Calls

4 calls. Slowest were:

35.0ms
(- (/ (+ (tan x) (tan eps)) (- 1 (* (tan x) (tan eps)))) (tan x))
12.0ms
(/ (+ (tan x) (tan eps)) (- 1 (* (tan x) (tan eps))))
4.0ms
(* (tan x) (tan eps))

series928.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

409.0ms
(- (/ (+ (tan x) (tan eps)) (- 1 (* (tan x) (tan eps)))) (tan x))
268.0ms
(/ (+ (tan x) (tan eps)) (- 1 (* (tan x) (tan eps))))
129.0ms
(* (tan x) (tan eps))
121.0ms
(+ (tan x) (tan eps))

simplify14.8s

Counts
199 → 180
Calls

199 calls. Slowest were:

560.0ms
(fma (/ 1 (* (cbrt (- 1 (* (tan x) (tan eps)))) (cbrt (- 1 (* (tan x) (tan eps)))))) (/ (+ (tan x) (tan eps)) (cbrt (- 1 (* (tan x) (tan eps))))) (- (* (sqrt (tan x)) (sqrt (tan x)))))
515.0ms
(fma (/ (+ (tan x) (tan eps)) (- (pow 1 3) (pow (* (tan x) (tan eps)) 3))) (+ (* 1 1) (+ (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* 1 (* (tan x) (tan eps))))) (- (* (cbrt (tan x)) (* (cbrt (tan x)) (cbrt (tan x))))))
503.0ms
(fma (/ 1 (* (cbrt (- 1 (* (tan x) (tan eps)))) (cbrt (- 1 (* (tan x) (tan eps)))))) (/ (+ (tan x) (tan eps)) (cbrt (- 1 (* (tan x) (tan eps))))) (- (* (sqrt (tan x)) (sqrt (tan x)))))

prune1.9s

Pruning

18 alts after pruning (18 fresh and 0 done)

Merged error: 14.2b

localize32.0ms

Local error

Found 4 expressions with local error:

1.9b
(fma (/ (+ (tan x) (tan eps)) (- (pow 1 3) (pow (* (tan x) (tan eps)) 3))) (+ (* 1 1) (+ (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* 1 (* (tan x) (tan eps))))) (- (tan x)))
0.3b
(pow (* (tan x) (tan eps)) 3)
0.3b
(* (* (tan x) (tan eps)) (* (tan x) (tan eps)))
0.3b
(* (tan x) (tan eps))

rewrite27.0ms

Algorithm
rewrite-expression-head
Counts
4 → 101
Calls

4 calls. Slowest were:

17.0ms
(* (* (tan x) (tan eps)) (* (tan x) (tan eps)))
4.0ms
(* (tan x) (tan eps))
4.0ms
(pow (* (tan x) (tan eps)) 3)

series1.3s

Counts
4 → 12
Calls

4 calls. Slowest were:

947.0ms
(fma (/ (+ (tan x) (tan eps)) (- (pow 1 3) (pow (* (tan x) (tan eps)) 3))) (+ (* 1 1) (+ (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* 1 (* (tan x) (tan eps))))) (- (tan x)))
153.0ms
(pow (* (tan x) (tan eps)) 3)
142.0ms
(* (* (tan x) (tan eps)) (* (tan x) (tan eps)))
101.0ms
(* (tan x) (tan eps))

simplify5.3s

Counts
83 → 113
Calls

83 calls. Slowest were:

683.0ms
(* (* (* (* (tan x) (tan x)) (tan x)) (* (* (tan eps) (tan eps)) (tan eps))) (* (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* (tan x) (tan eps))))
620.0ms
(* (* (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* (tan x) (tan eps))) (* (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* (tan x) (tan eps))))
567.0ms
(* (* (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* (tan x) (tan eps))) (* (* (* (tan x) (tan x)) (tan x)) (* (* (tan eps) (tan eps)) (tan eps))))

prune2.6s

Pruning

17 alts after pruning (17 fresh and 0 done)

Merged error: 14.2b

localize21.0ms

Local error

Found 4 expressions with local error:

1.9b
(fma (/ (+ (tan x) (tan eps)) (- (pow 1 3) (* (pow (tan x) 3) (* (* (tan eps) (tan eps)) (tan eps))))) (+ (* 1 1) (+ (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* 1 (* (tan x) (tan eps))))) (- (tan x)))
0.4b
(pow (tan x) 3)
0.3b
(* (* (tan x) (tan eps)) (* (tan x) (tan eps)))
0.3b
(* (* (tan eps) (tan eps)) (tan eps))

rewrite38.0ms

Algorithm
rewrite-expression-head
Counts
4 → 112
Calls

4 calls. Slowest were:

22.0ms
(* (* (tan x) (tan eps)) (* (tan x) (tan eps)))
11.0ms
(* (* (tan eps) (tan eps)) (tan eps))
1.0ms
(pow (tan x) 3)

series1.2s

Counts
4 → 12
Calls

4 calls. Slowest were:

972.0ms
(fma (/ (+ (tan x) (tan eps)) (- (pow 1 3) (* (pow (tan x) 3) (* (* (tan eps) (tan eps)) (tan eps))))) (+ (* 1 1) (+ (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* 1 (* (tan x) (tan eps))))) (- (tan x)))
120.0ms
(* (* (tan x) (tan eps)) (* (tan x) (tan eps)))
55.0ms
(* (* (tan eps) (tan eps)) (tan eps))
55.0ms
(pow (tan x) 3)

simplify4.2s

Counts
100 → 124
Calls

100 calls. Slowest were:

693.0ms
(* (* (* (* (tan x) (tan x)) (tan x)) (* (* (tan eps) (tan eps)) (tan eps))) (* (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* (tan x) (tan eps))))
645.0ms
(* (* (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* (tan x) (tan eps))) (* (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* (tan x) (tan eps))))
481.0ms
(* (* (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* (tan x) (tan eps))) (* (* (* (tan x) (tan x)) (tan x)) (* (* (tan eps) (tan eps)) (tan eps))))

prune2.0s

Pruning

17 alts after pruning (16 fresh and 1 done)

Merged error: 14.2b

regimes232.0ms

Accuracy

88.3% (0.9b remaining)

Error of 15.9b against oracle of 15.1b and baseline of 22.4b

bsearch92.0ms

end0.0ms

sample11.1s

Algorithm
intervals