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

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

    \[r \cdot \frac{\sin b}{\cos \left(a + b\right)}\]
  2. Initial simplification15.1

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

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

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

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

    \[\leadsto \color{blue}{r} \cdot \frac{\sin b}{\cos b \cdot \cos a - \sin b \cdot \sin a}\]
  9. Using strategy rm
  10. Applied expm1-log1p-u0.3

    \[\leadsto r \cdot \frac{\sin b}{\cos b \cdot \cos a - \color{blue}{(e^{\log_* (1 + \sin b \cdot \sin a)} - 1)^*}}\]
  11. Using strategy rm
  12. Applied div-inv0.4

    \[\leadsto r \cdot \color{blue}{\left(\sin b \cdot \frac{1}{\cos b \cdot \cos a - (e^{\log_* (1 + \sin b \cdot \sin a)} - 1)^*}\right)}\]
  13. Applied associate-*r*0.4

    \[\leadsto \color{blue}{\left(r \cdot \sin b\right) \cdot \frac{1}{\cos b \cdot \cos a - (e^{\log_* (1 + \sin b \cdot \sin a)} - 1)^*}}\]
  14. Simplified0.4

    \[\leadsto \left(r \cdot \sin b\right) \cdot \color{blue}{\frac{1}{\cos a \cdot \cos b - \sin a \cdot \sin b}}\]
  15. Final simplification0.4

    \[\leadsto \frac{1}{\cos b \cdot \cos a - \sin b \cdot \sin a} \cdot \left(r \cdot \sin b\right)\]

Reproduce

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

Details

Time bar (total: 17.1s)Debug log

start407.0ms

Algorithm
intervals

setup63.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 15.6b

localize22.0ms

Local error

Found 3 expressions with local error:

2.3b
(cos (+ b a))
0.2b
(/ (* r (sin b)) (cos (+ b a)))
0.1b
(* r (sin b))

rewrite19.0ms

Algorithm
rewrite-expression-head
Counts
3 → 44
Calls

3 calls. Slowest were:

8.0ms
(/ (* r (sin b)) (cos (+ b a)))
5.0ms
(cos (+ b a))
2.0ms
(* r (sin b))

series523.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

437.0ms
(/ (* r (sin b)) (cos (+ b a)))
51.0ms
(* r (sin b))
35.0ms
(cos (+ b a))

simplify1.1s

Counts
24 → 53
Calls

24 calls. Slowest were:

580.0ms
(/ (* (* (* r (sin b)) (* r (sin b))) (* r (sin b))) (* (* (cos (+ b a)) (cos (+ b a))) (cos (+ b a))))
283.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))
131.0ms
(- 1 (+ (* 1/2 (pow b 2)) (* a b)))

prune731.0ms

Pruning

9 alts after pruning (9 fresh and 0 done)

Merged error: 0.2b

localize21.0ms

Local error

Found 4 expressions with local error:

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

rewrite25.0ms

Algorithm
rewrite-expression-head
Counts
4 → 83
Calls

4 calls. Slowest were:

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

series484.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

261.0ms
(/ (* r (sin b)) (- (* (cos b) (cos a)) (* (sin b) (sin a))))
111.0ms
(- (* (cos b) (cos a)) (* (sin b) (sin a)))
75.0ms
(* (sin b) (sin a))
36.0ms
(* (cos b) (cos a))

simplify2.1s

Counts
51 → 95
Calls

51 calls. Slowest were:

406.0ms
(- (* (+ (cos (+ b a)) (cos (- b a))) 2) (* 2 (- (cos (- b a)) (cos (+ b a)))))
321.0ms
(- (+ (cos (+ b a)) (cos (- b a))) (- (cos (- b a)) (cos (+ b a))))
256.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))

prune1.4s

Pruning

12 alts after pruning (12 fresh and 0 done)

Merged error: 0.0b

localize12.0ms

Local error

Found 4 expressions with local error:

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

rewrite41.0ms

Algorithm
rewrite-expression-head
Counts
4 → 99
Calls

4 calls. Slowest were:

18.0ms
(/ (sin b) (- (* (cos b) (cos a)) (* (sin b) (sin a))))
11.0ms
(* r (/ (sin b) (- (* (cos b) (cos a)) (* (sin b) (sin a)))))
7.0ms
(- (* (cos b) (cos a)) (* (sin b) (sin a)))

series561.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

237.0ms
(* r (/ (sin b) (- (* (cos b) (cos a)) (* (sin b) (sin a)))))
166.0ms
(/ (sin b) (- (* (cos b) (cos a)) (* (sin b) (sin a))))
91.0ms
(- (* (cos b) (cos a)) (* (sin b) (sin a)))
66.0ms
(* (sin b) (sin a))

simplify3.0s

Counts
73 → 111
Calls

73 calls. Slowest were:

426.0ms
(- (* (+ (cos (+ b a)) (cos (- b a))) 2) (* 2 (- (cos (- b a)) (cos (+ b a)))))
377.0ms
(- (+ (cos (+ b a)) (cos (- b a))) (- (cos (- b a)) (cos (+ b a))))
249.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))

prune1.6s

Pruning

10 alts after pruning (10 fresh and 0 done)

Merged error: 0.0b

localize13.0ms

Local error

Found 4 expressions with local error:

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

rewrite44.0ms

Algorithm
rewrite-expression-head
Counts
4 → 96
Calls

4 calls. Slowest were:

17.0ms
(/ (sin b) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a))))))
14.0ms
(* r (/ (sin b) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a)))))))
7.0ms
(* (sin b) (sin a))

series652.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

292.0ms
(* r (/ (sin b) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a)))))))
151.0ms
(/ (sin b) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a))))))
149.0ms
(- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a)))))
59.0ms
(* (sin b) (sin a))

simplify2.3s

Counts
72 → 108
Calls

72 calls. Slowest were:

263.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))
179.0ms
(+ (* (* (cos b) (cos a)) (* (cos b) (cos a))) (+ (* (expm1 (log1p (* (sin b) (sin a)))) (expm1 (log1p (* (sin b) (sin a))))) (* (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a)))))))
167.0ms
(/ (* (* (sin b) (sin b)) (sin b)) (* (* (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a))))) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a)))))) (- (* (cos b) (cos a)) (expm1 (log1p (* (sin b) (sin a)))))))

prune1.6s

Pruning

9 alts after pruning (9 fresh and 0 done)

Merged error: 0.0b

regimes346.0ms

Accuracy

0% (0.3b remaining)

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

bsearch2.0ms