Average Error: 37.1 → 15.5
Time: 57.4s
Precision: 64
Internal Precision: 128
\[\tan \left(x + \varepsilon\right) - \tan x\]
\[\begin{array}{l} \mathbf{if}\;\varepsilon \le -1.6718721488586065 \cdot 10^{-18}:\\ \;\;\;\;(\left(\tan \varepsilon + \tan x\right) \cdot \left(\sqrt[3]{\frac{1}{1 - \tan x \cdot \tan \varepsilon} \cdot \left(\frac{1}{1 - \tan x \cdot \tan \varepsilon} \cdot \frac{1}{1 - \tan x \cdot \tan \varepsilon}\right)}\right) + \left(-\tan x\right))_*\\ \mathbf{elif}\;\varepsilon \le 1.0288762116509696 \cdot 10^{-29}:\\ \;\;\;\;(\left(x \cdot \varepsilon\right) \cdot \left(x + \varepsilon\right) + \varepsilon)_*\\ \mathbf{else}:\\ \;\;\;\;(\left(\frac{\tan \varepsilon + \tan x}{1 - \left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right)}\right) \cdot \left(1 + \tan x \cdot \tan \varepsilon\right) + \left(-\tan x\right))_*\\ \end{array}\]

Error

Bits error versus x

Bits error versus eps

Target

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

Derivation

  1. Split input into 3 regimes
  2. if eps < -1.6718721488586065e-18

    1. Initial program 29.4

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

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

      \[\leadsto \color{blue}{\left(\tan x + \tan \varepsilon\right) \cdot \frac{1}{1 - \tan x \cdot \tan \varepsilon}} - \tan x\]
    6. Applied fma-neg0.9

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

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

    if -1.6718721488586065e-18 < eps < 1.0288762116509696e-29

    1. Initial program 45.7

      \[\tan \left(x + \varepsilon\right) - \tan x\]
    2. Taylor expanded around 0 31.4

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

      \[\leadsto \color{blue}{(\left(x \cdot \varepsilon\right) \cdot \left(\varepsilon + x\right) + \varepsilon)_*}\]

    if 1.0288762116509696e-29 < eps

    1. Initial program 29.5

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

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

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

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

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

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

    \[\leadsto \begin{array}{l} \mathbf{if}\;\varepsilon \le -1.6718721488586065 \cdot 10^{-18}:\\ \;\;\;\;(\left(\tan \varepsilon + \tan x\right) \cdot \left(\sqrt[3]{\frac{1}{1 - \tan x \cdot \tan \varepsilon} \cdot \left(\frac{1}{1 - \tan x \cdot \tan \varepsilon} \cdot \frac{1}{1 - \tan x \cdot \tan \varepsilon}\right)}\right) + \left(-\tan x\right))_*\\ \mathbf{elif}\;\varepsilon \le 1.0288762116509696 \cdot 10^{-29}:\\ \;\;\;\;(\left(x \cdot \varepsilon\right) \cdot \left(x + \varepsilon\right) + \varepsilon)_*\\ \mathbf{else}:\\ \;\;\;\;(\left(\frac{\tan \varepsilon + \tan x}{1 - \left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right)}\right) \cdot \left(1 + \tan x \cdot \tan \varepsilon\right) + \left(-\tan x\right))_*\\ \end{array}\]

Reproduce

herbie shell --seed 2018360 +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: 43.1s)Debug log

start560.0ms

Algorithm
intervals

setup21.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 33.6b

localize23.0ms

Local error

Found 2 expressions with local error:

3.4b
(tan (+ x eps))
1.2b
(- (tan (+ x eps)) (tan x))

rewrite14.0ms

Algorithm
rewrite-expression-head
Counts
2 → 41
Calls

2 calls. Slowest were:

11.0ms
(- (tan (+ x eps)) (tan x))
3.0ms
(tan (+ x eps))

series173.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

100.0ms
(- (tan (+ x eps)) (tan x))
74.0ms
(tan (+ x eps))

simplify1.2s

Counts
35 → 47
Calls

35 calls. Slowest were:

499.0ms
(- (* (+ (tan x) (tan eps)) (cos x)) (* (- 1 (* (tan x) (tan eps))) (sin x)))
118.0ms
(+ (* x (pow eps 2)) (+ eps (* (pow x 2) eps)))
89.0ms
(* (- 1 (* (tan x) (tan eps))) (cos x))

prune471.0ms

Pruning

8 alts after pruning (8 fresh and 0 done)

Merged error: 11.5b

localize32.0ms

Local error

Found 4 expressions with local error:

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

rewrite87.0ms

Algorithm
rewrite-expression-head
Counts
4 → 168
Calls

4 calls. Slowest were:

64.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))

series868.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

415.0ms
(- (/ (+ (tan x) (tan eps)) (- 1 (* (tan x) (tan eps)))) (tan x))
234.0ms
(/ (+ (tan x) (tan eps)) (- 1 (* (tan x) (tan eps))))
110.0ms
(+ (tan x) (tan eps))
108.0ms
(* (tan x) (tan eps))

simplify10.6s

Counts
199 → 180
Calls

199 calls. Slowest were:

359.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))))) (- (* (cbrt (tan x)) (* (cbrt (tan x)) (cbrt (tan x))))))
320.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)))))
312.0ms
(- (* (+ (tan x) (tan eps)) (cos x)) (* (- 1 (* (tan x) (tan eps))) (sin x)))

prune1.5s

Pruning

12 alts after pruning (12 fresh and 0 done)

Merged error: 11.4b

localize8.0ms

Local error

Found 4 expressions with local error:

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

rewrite10.0ms

Algorithm
rewrite-expression-head
Counts
4 → 70
Calls

4 calls. Slowest were:

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

series545.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

244.0ms
(fma (+ (tan x) (tan eps)) (/ 1 (- 1 (* (tan x) (tan eps)))) (- (tan x)))
142.0ms
(/ 1 (- 1 (* (tan x) (tan eps))))
79.0ms
(* (tan x) (tan eps))
79.0ms
(+ (tan x) (tan eps))

simplify1.4s

Counts
36 → 82
Calls

36 calls. Slowest were:

244.0ms
(+ (* 1/3 (* (pow x 3) eps)) (+ (* 1/3 (* x (pow eps 3))) (* x eps)))
197.0ms
(+ (* 1/3 (* (pow x 3) eps)) (+ (* x eps) 1))
154.0ms
(/ 1 (- 1 (/ (* (sin x) (sin eps)) (* (cos x) (cos eps)))))

prune1.5s

Pruning

13 alts after pruning (13 fresh and 0 done)

Merged error: 11.4b

localize35.0ms

Local error

Found 4 expressions with local error:

28.0b
(cbrt (* (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* (tan x) (tan eps))))
1.8b
(fma (+ (tan x) (tan eps)) (/ 1 (- 1 (cbrt (* (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* (tan x) (tan eps)))))) (- (tan x)))
0.3b
(* (* (tan x) (tan eps)) (* (tan x) (tan eps)))
0.2b
(* (tan x) (tan eps))

rewrite123.0ms

Algorithm
rewrite-expression-head
Counts
4 → 153
Calls

4 calls. Slowest were:

73.0ms
(cbrt (* (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* (tan x) (tan eps))))
35.0ms
(* (* (tan x) (tan eps)) (* (tan x) (tan eps)))
8.0ms
(* (tan x) (tan eps))

series643.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

324.0ms
(fma (+ (tan x) (tan eps)) (/ 1 (- 1 (cbrt (* (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* (tan x) (tan eps)))))) (- (tan x)))
115.0ms
(cbrt (* (* (* (tan x) (tan eps)) (* (tan x) (tan eps))) (* (tan x) (tan eps))))
104.0ms
(* (tan x) (tan eps))
101.0ms
(* (* (tan x) (tan eps)) (* (tan x) (tan eps)))

simplify20.8s

Counts
137 → 165
Calls

137 calls. Slowest were:

551.0ms
(cbrt (* (* (* (tan x) (tan eps)) (* (sin x) (sin eps))) (* (sin x) (tan eps))))
541.0ms
(cbrt (* (* (* (tan x) (tan eps)) (* (sin x) (tan eps))) (* (sin x) (sin eps))))
523.0ms
(cbrt (* (* (* (tan x) (tan eps)) (* (tan x) (sin eps))) (* (sin x) (tan eps))))

prune1.8s

Pruning

14 alts after pruning (14 fresh and 0 done)

Merged error: 11.4b

regimes203.0ms

Accuracy

90.6% (0.7b remaining)

Error of 15.5b against oracle of 14.8b and baseline of 22.2b

bsearch342.0ms