Average Error: 31.0 → 0.1
Time: 53.7s
Precision: 64
Internal Precision: 128
\[\frac{x - \sin x}{x - \tan x}\]
\[\begin{array}{l} \mathbf{if}\;x \le -1.5592990590494797:\\ \;\;\;\;(\left(\sqrt{\frac{x}{x - \tan x}}\right) \cdot \left(\sqrt{\frac{x}{x - \tan x}}\right) + \left(\frac{-\sin x}{x - \tan x}\right))_*\\ \mathbf{elif}\;x \le 0.028176034411070774:\\ \;\;\;\;(\left(\frac{9}{40} \cdot x\right) \cdot x + \left((\frac{-27}{2800} \cdot \left({x}^{4}\right) + \frac{-1}{2})_*\right))_*\\ \mathbf{else}:\\ \;\;\;\;\frac{x}{x - \tan x} - \frac{\sin x}{x - \tan x}\\ \end{array}\]

Error

Bits error versus x

Derivation

  1. Split input into 3 regimes
  2. if x < -1.5592990590494797

    1. Initial program 0.0

      \[\frac{x - \sin x}{x - \tan x}\]
    2. Using strategy rm
    3. Applied div-sub0.0

      \[\leadsto \color{blue}{\frac{x}{x - \tan x} - \frac{\sin x}{x - \tan x}}\]
    4. Using strategy rm
    5. Applied add-sqr-sqrt0.1

      \[\leadsto \color{blue}{\sqrt{\frac{x}{x - \tan x}} \cdot \sqrt{\frac{x}{x - \tan x}}} - \frac{\sin x}{x - \tan x}\]
    6. Applied fma-neg0.1

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

    if -1.5592990590494797 < x < 0.028176034411070774

    1. Initial program 62.4

      \[\frac{x - \sin x}{x - \tan x}\]
    2. Using strategy rm
    3. Applied div-sub62.3

      \[\leadsto \color{blue}{\frac{x}{x - \tan x} - \frac{\sin x}{x - \tan x}}\]
    4. Taylor expanded around 0 0.1

      \[\leadsto \color{blue}{\frac{9}{40} \cdot {x}^{2} - \left(\frac{27}{2800} \cdot {x}^{4} + \frac{1}{2}\right)}\]
    5. Simplified0.1

      \[\leadsto \color{blue}{(\left(\frac{9}{40} \cdot x\right) \cdot x + \left((\frac{-27}{2800} \cdot \left({x}^{4}\right) + \frac{-1}{2})_*\right))_*}\]

    if 0.028176034411070774 < x

    1. Initial program 0.0

      \[\frac{x - \sin x}{x - \tan x}\]
    2. Using strategy rm
    3. Applied div-sub0.0

      \[\leadsto \color{blue}{\frac{x}{x - \tan x} - \frac{\sin x}{x - \tan x}}\]
  3. Recombined 3 regimes into one program.
  4. Final simplification0.1

    \[\leadsto \begin{array}{l} \mathbf{if}\;x \le -1.5592990590494797:\\ \;\;\;\;(\left(\sqrt{\frac{x}{x - \tan x}}\right) \cdot \left(\sqrt{\frac{x}{x - \tan x}}\right) + \left(\frac{-\sin x}{x - \tan x}\right))_*\\ \mathbf{elif}\;x \le 0.028176034411070774:\\ \;\;\;\;(\left(\frac{9}{40} \cdot x\right) \cdot x + \left((\frac{-27}{2800} \cdot \left({x}^{4}\right) + \frac{-1}{2})_*\right))_*\\ \mathbf{else}:\\ \;\;\;\;\frac{x}{x - \tan x} - \frac{\sin x}{x - \tan x}\\ \end{array}\]

Reproduce

herbie shell --seed 2019007 +o rules:numerics
(FPCore (x)
  :name "sintan (problem 3.4.5)"
  (/ (- x (sin x)) (- x (tan x))))

Details

Time bar (total: 53.1s)Debug log

sample721.0ms

Algorithm
intervals

simplify21.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

20.0ms
(/ (- x (sin x)) (- x (tan x)))

prune6.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 31.4b

localize44.0ms

Local error

Found 3 expressions with local error:

1.5b
(- x (sin x))
1.5b
(- x (tan x))
0.0b
(/ (- x (sin x)) (- x (tan x)))

rewrite16.0ms

Algorithm
rewrite-expression-head
Counts
3 → 60
Calls

3 calls. Slowest were:

8.0ms
(/ (- x (sin x)) (- x (tan x)))
3.0ms
(- x (sin x))
1.0ms
(- x (tan x))

series230.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

89.0ms
(/ (- x (sin x)) (- x (tan x)))
77.0ms
(- x (tan x))
64.0ms
(- x (sin x))

simplify1.4s

Counts
39 → 69
Calls

39 calls. Slowest were:

269.0ms
(- (+ 1 (+ (/ (pow (sin x) 2) (* (pow (cos x) 2) (pow x 2))) (/ (sin x) (* (cos x) x)))) (+ (/ (sin x) x) (/ (pow (sin x) 2) (* (cos x) (pow x 2)))))
251.0ms
(- (+ 1 (+ (/ (pow (sin x) 2) (* (pow (cos x) 2) (pow x 2))) (/ (sin x) (* (cos x) x)))) (+ (/ (sin x) x) (/ (pow (sin x) 2) (* (cos x) (pow x 2)))))
184.0ms
(- (+ (* 1/3 (pow x 3)) (+ (* 2/15 (pow x 5)) (* 17/315 (pow x 7)))))

prune666.0ms

Pruning

3 alts after pruning (3 fresh and 0 done)

Merged error: 0.0b

localize17.0ms

Local error

Found 4 expressions with local error:

1.5b
(- (/ x (- x (tan x))) (/ (sin x) (- x (tan x))))
1.5b
(- x (tan x))
1.5b
(- x (tan x))
0.1b
(/ (sin x) (- x (tan x)))

rewrite29.0ms

Algorithm
rewrite-expression-head
Counts
4 → 175
Calls

4 calls. Slowest were:

15.0ms
(- (/ x (- x (tan x))) (/ (sin x) (- x (tan x))))
7.0ms
(/ (sin x) (- x (tan x)))
1.0ms
(- x (tan x))

series527.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

281.0ms
(- (/ x (- x (tan x))) (/ (sin x) (- x (tan x))))
94.0ms
(/ (sin x) (- x (tan x)))
83.0ms
(- x (tan x))
68.0ms
(- x (tan x))

simplify29.8s

Counts
239 → 187
Calls

239 calls. Slowest were:

747.0ms
(fma (/ x (- (pow x 3) (pow (tan x) 3))) (+ (* x x) (+ (* (tan x) (tan x)) (* x (tan x)))) (- (* (cbrt (/ (sin x) (- x (tan x)))) (* (cbrt (/ (sin x) (- x (tan x)))) (cbrt (/ (sin x) (- x (tan x))))))))
605.0ms
(fma (/ x (- (* x x) (* (tan x) (tan x)))) (+ x (tan x)) (- (* (/ (sin x) (- x (tan x))) 1)))
533.0ms
(fma x (/ 1 (- x (tan x))) (- (* (/ (cbrt (sin x)) (- x (tan x))) (/ (* (cbrt (sin x)) (cbrt (sin x))) 1))))

prune1.2s

Pruning

3 alts after pruning (2 fresh and 1 done)

Merged error: 0.0b

localize15.0ms

Local error

Found 3 expressions with local error:

0.5b
(fma -27/2800 (pow x 4) -1/2)
0.2b
(* 9/40 x)
0.0b
(fma (* 9/40 x) x (fma -27/2800 (pow x 4) -1/2))

rewrite1.0ms

Algorithm
rewrite-expression-head
Counts
3 → 30
Calls

3 calls. Slowest were:

0.0ms
(* 9/40 x)
0.0ms
(fma -27/2800 (pow x 4) -1/2)
0.0ms
(fma (* 9/40 x) x (fma -27/2800 (pow x 4) -1/2))

series89.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

46.0ms
(fma -27/2800 (pow x 4) -1/2)
32.0ms
(fma (* 9/40 x) x (fma -27/2800 (pow x 4) -1/2))
11.0ms
(* 9/40 x)

simplify245.0ms

Counts
9 → 39
Calls

9 calls. Slowest were:

87.0ms
(- (* 9/40 (pow x 2)) (+ (* 27/2800 (pow x 4)) 1/2))
66.0ms
(- (* 9/40 (pow x 2)) (+ (* 27/2800 (pow x 4)) 1/2))
54.0ms
(- (* 9/40 (pow x 2)) (+ (* 27/2800 (pow x 4)) 1/2))

prune286.0ms

Pruning

4 alts after pruning (2 fresh and 2 done)

Merged error: 0.0b

localize21.0ms

Local error

Found 4 expressions with local error:

5.9b
(cbrt (* (* (fma -27/2800 (pow x 4) -1/2) (fma -27/2800 (pow x 4) -1/2)) (fma -27/2800 (pow x 4) -1/2)))
0.5b
(fma -27/2800 (pow x 4) -1/2)
0.5b
(fma -27/2800 (pow x 4) -1/2)
0.5b
(fma -27/2800 (pow x 4) -1/2)

rewrite12.0ms

Algorithm
rewrite-expression-head
Counts
4 → 42
Calls

4 calls. Slowest were:

10.0ms
(cbrt (* (* (fma -27/2800 (pow x 4) -1/2) (fma -27/2800 (pow x 4) -1/2)) (fma -27/2800 (pow x 4) -1/2)))
0.0ms
(fma -27/2800 (pow x 4) -1/2)
0.0ms
(fma -27/2800 (pow x 4) -1/2)

series213.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

57.0ms
(cbrt (* (* (fma -27/2800 (pow x 4) -1/2) (fma -27/2800 (pow x 4) -1/2)) (fma -27/2800 (pow x 4) -1/2)))
53.0ms
(fma -27/2800 (pow x 4) -1/2)
52.0ms
(fma -27/2800 (pow x 4) -1/2)
52.0ms
(fma -27/2800 (pow x 4) -1/2)

simplify248.0ms

Counts
13 → 54
Calls

13 calls. Slowest were:

58.0ms
(cbrt (* (fma -27/2800 (pow x 4) -1/2) (fma -27/2800 (pow x 4) -1/2)))
20.0ms
(- (+ (* 27/2800 (pow x 4)) 1/2))
18.0ms
(- (+ (* 27/2800 (pow x 4)) 1/2))

prune340.0ms

Pruning

4 alts after pruning (2 fresh and 2 done)

Merged error: 0.0b

regimes39.0ms

Accuracy

99.8% (0.0b remaining)

Error of 0.1b against oracle of 0.0b and baseline of 31.0b

bsearch113.0ms

end0.0ms

sample16.7s

Algorithm
intervals