Average Error: 15.1 → 0.3
Time: 27.4s
Precision: 64
Internal Precision: 128
\[\frac{r \cdot \sin b}{\cos \left(a + b\right)}\]
\[r \cdot \frac{\sin b}{\cos a \cdot \cos b - \sin a \cdot \sin b}\]

Error

Bits error versus r

Bits error versus a

Bits error versus b

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 15.1

    \[\frac{r \cdot \sin b}{\cos \left(a + b\right)}\]
  2. Using strategy rm
  3. Applied cos-sum0.3

    \[\leadsto \frac{r \cdot \sin b}{\color{blue}{\cos a \cdot \cos b - \sin a \cdot \sin b}}\]
  4. Using strategy rm
  5. Applied *-un-lft-identity0.3

    \[\leadsto \frac{r \cdot \sin b}{\color{blue}{1 \cdot \left(\cos a \cdot \cos b - \sin a \cdot \sin b\right)}}\]
  6. Applied times-frac0.3

    \[\leadsto \color{blue}{\frac{r}{1} \cdot \frac{\sin b}{\cos a \cdot \cos b - \sin a \cdot \sin b}}\]
  7. Simplified0.3

    \[\leadsto \color{blue}{r} \cdot \frac{\sin b}{\cos a \cdot \cos b - \sin a \cdot \sin b}\]
  8. Taylor expanded around inf 0.3

    \[\leadsto r \cdot \frac{\sin b}{\color{blue}{\cos a \cdot \cos b} - \sin a \cdot \sin b}\]
  9. Final simplification0.3

    \[\leadsto r \cdot \frac{\sin b}{\cos a \cdot \cos b - \sin a \cdot \sin b}\]

Reproduce

herbie shell --seed 2018362 +o rules:numerics
(FPCore (r a b)
  :name "r*sin(b)/cos(a+b), A"
  (/ (* r (sin b)) (cos (+ a b))))

Details

Time bar (total: 15.7s)Debug log

start674.0ms

Algorithm
intervals

setup600.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 12.4b

localize42.0ms

Local error

Found 3 expressions with local error:

3.0b
(cos (+ a b))
0.2b
(/ (* r (sin b)) (cos (+ a b)))
0.2b
(* r (sin b))

rewrite20.0ms

Algorithm
rewrite-expression-head
Counts
3 → 44
Calls

3 calls. Slowest were:

9.0ms
(/ (* r (sin b)) (cos (+ a b)))
6.0ms
(cos (+ a b))
2.0ms
(* r (sin b))

series252.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

157.0ms
(/ (* r (sin b)) (cos (+ a b)))
59.0ms
(* r (sin b))
35.0ms
(cos (+ a b))

simplify1.0s

Counts
24 → 53
Calls

24 calls. Slowest were:

543.0ms
(/ (* (* (* r (sin b)) (* r (sin b))) (* r (sin b))) (* (* (cos (+ a b)) (cos (+ a b))) (cos (+ a b))))
233.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))
162.0ms
(- 1 (+ (* 1/2 (pow a 2)) (* a b)))

prune599.0ms

Pruning

8 alts after pruning (8 fresh and 0 done)

Merged error: 0.1b

localize30.0ms

Local error

Found 4 expressions with local error:

0.2b
(/ (* r (sin b)) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
0.2b
(* (sin a) (sin b))
0.2b
(* r (sin b))
0.1b
(* (cos a) (cos b))

rewrite22.0ms

Algorithm
rewrite-expression-head
Counts
4 → 78
Calls

4 calls. Slowest were:

10.0ms
(/ (* r (sin b)) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
7.0ms
(* (sin a) (sin b))
3.0ms
(* (cos a) (cos b))

series441.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

257.0ms
(/ (* r (sin b)) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
73.0ms
(* (sin a) (sin b))
57.0ms
(* (cos a) (cos b))
54.0ms
(* r (sin b))

simplify1.5s

Counts
47 → 90
Calls

47 calls. Slowest were:

270.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))
219.0ms
(+ (* (* (cos a) (cos b)) (* (cos a) (cos b))) (+ (* (* (sin a) (sin b)) (* (sin a) (sin b))) (* (* (cos a) (cos b)) (* (sin a) (sin b)))))
165.0ms
(- 1 (+ (* 1/2 (pow a 2)) (* 1/2 (pow b 2))))

prune1.3s

Pruning

11 alts after pruning (11 fresh and 0 done)

Merged error: 0.0b

localize6.0ms

Local error

Found 4 expressions with local error:

0.2b
(* (sin a) (sin b))
0.2b
(/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
0.1b
(* r (/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b)))))
0.1b
(* (cos a) (cos b))

rewrite35.0ms

Algorithm
rewrite-expression-head
Counts
4 → 101
Calls

4 calls. Slowest were:

16.0ms
(/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
11.0ms
(* r (/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b)))))
4.0ms
(* (sin a) (sin b))

series420.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

220.0ms
(* r (/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b)))))
100.0ms
(/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
51.0ms
(* (cos a) (cos b))
49.0ms
(* (sin a) (sin b))

simplify2.3s

Counts
75 → 113
Calls

75 calls. Slowest were:

253.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))
231.0ms
(+ b (+ (* 1/3 (pow b 3)) (* a (pow b 2))))
216.0ms
(/ (* (* (sin b) (sin b)) (sin b)) (* (* (- (* (cos a) (cos b)) (* (sin a) (sin b))) (- (* (cos a) (cos b)) (* (sin a) (sin b)))) (- (* (cos a) (cos b)) (* (sin a) (sin b)))))

prune1.8s

Pruning

10 alts after pruning (10 fresh and 0 done)

Merged error: 0.0b

localize5.0ms

Local error

Found 4 expressions with local error:

0.2b
(* (sin a) (sin b))
0.2b
(/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
0.1b
(* r (/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b)))))
0.1b
(* (cos a) (cos b))

rewrite29.0ms

Algorithm
rewrite-expression-head
Counts
4 → 101
Calls

4 calls. Slowest were:

11.0ms
(* r (/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b)))))
10.0ms
(/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
4.0ms
(* (sin a) (sin b))

series465.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

264.0ms
(* r (/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b)))))
95.0ms
(/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
55.0ms
(* (sin a) (sin b))
51.0ms
(* (cos a) (cos b))

simplify2.2s

Counts
75 → 113
Calls

75 calls. Slowest were:

311.0ms
(+ (* (* (cos a) (cos b)) (* (cos a) (cos b))) (+ (* (* (sin a) (sin b)) (* (sin a) (sin b))) (* (* (cos a) (cos b)) (* (sin a) (sin b)))))
266.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))
229.0ms
(- 1 (+ (* 1/2 (pow a 2)) (* 1/2 (pow b 2))))

prune1.5s

Pruning

10 alts after pruning (9 fresh and 1 done)

Merged error: 0.0b

regimes351.0ms

Accuracy

0% (0.3b remaining)

Error of 0.3b against oracle of 0.1b and baseline of 0.3b

bsearch2.0ms