Average Error: 15.2 → 0.4
Time: 27.5s
Precision: 64
Internal Precision: 128
\[r \cdot \frac{\sin b}{\cos \left(a + b\right)}\]
\[\frac{\sin b}{\cos a \cdot \cos b - \sqrt[3]{{\left(\sin a \cdot \sin b\right)}^{3}}} \cdot r\]

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.2

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

    \[\leadsto r \cdot \frac{\sin b}{\color{blue}{\cos a \cdot \cos b - \sin a \cdot \sin b}}\]
  4. Using strategy rm
  5. Applied add-cbrt-cube0.4

    \[\leadsto r \cdot \frac{\sin b}{\cos a \cdot \cos b - \sin a \cdot \color{blue}{\sqrt[3]{\left(\sin b \cdot \sin b\right) \cdot \sin b}}}\]
  6. Applied add-cbrt-cube0.4

    \[\leadsto r \cdot \frac{\sin b}{\cos a \cdot \cos b - \color{blue}{\sqrt[3]{\left(\sin a \cdot \sin a\right) \cdot \sin a}} \cdot \sqrt[3]{\left(\sin b \cdot \sin b\right) \cdot \sin b}}\]
  7. Applied cbrt-unprod0.4

    \[\leadsto r \cdot \frac{\sin b}{\cos a \cdot \cos b - \color{blue}{\sqrt[3]{\left(\left(\sin a \cdot \sin a\right) \cdot \sin a\right) \cdot \left(\left(\sin b \cdot \sin b\right) \cdot \sin b\right)}}}\]
  8. Simplified0.4

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

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

Reproduce

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

sample393.0ms

Algorithm
intervals

simplify74.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

73.0ms
(* r (/ (sin b) (cos (+ a b))))

prune15.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 13.2b

localize21.0ms

Local error

Found 3 expressions with local error:

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

rewrite12.0ms

Algorithm
rewrite-expression-head
Counts
3 → 63
Calls

3 calls. Slowest were:

4.0ms
(/ (sin b) (cos (+ a b)))
3.0ms
(cos (+ a b))
3.0ms
(* r (/ (sin b) (cos (+ a b))))

series201.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

124.0ms
(* r (/ (sin b) (cos (+ a b))))
56.0ms
(/ (sin b) (cos (+ a b)))
22.0ms
(cos (+ a b))

simplify1.1s

Counts
49 → 72
Calls

49 calls. Slowest were:

283.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))
226.0ms
(+ b (+ (* 1/3 (pow b 3)) (* a (pow b 2))))
199.0ms
(- 1 (+ (* 1/2 (pow a 2)) (* a b)))

prune753.0ms

Pruning

8 alts after pruning (8 fresh and 0 done)

Merged error: 0.2b

localize51.0ms

Local error

Found 4 expressions with local error:

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

rewrite61.0ms

Algorithm
rewrite-expression-head
Counts
4 → 99
Calls

4 calls. Slowest were:

22.0ms
(* r (/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b)))))
22.0ms
(/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
8.0ms
(- (* (cos a) (cos b)) (* (sin a) (sin b)))

series553.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

277.0ms
(* r (/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b)))))
132.0ms
(/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
88.0ms
(- (* (cos a) (cos b)) (* (sin a) (sin b)))
56.0ms
(* (sin a) (sin b))

simplify3.1s

Counts
73 → 111
Calls

73 calls. Slowest were:

407.0ms
(- (+ (cos (+ a b)) (cos (- a b))) (- (cos (- a b)) (cos (+ a b))))
309.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))
304.0ms
(- (* (+ (cos (+ a b)) (cos (- a b))) 2) (* 2 (- (cos (- a b)) (cos (+ a b)))))

prune1.5s

Pruning

12 alts after pruning (12 fresh and 0 done)

Merged error: 0.0b

localize17.0ms

Local error

Found 4 expressions with local error:

26.8b
(cbrt (pow (* (sin b) (sin a)) 3))
0.3b
(pow (* (sin b) (sin a)) 3)
0.2b
(/ (sin b) (- (* (cos a) (cos b)) (cbrt (pow (* (sin b) (sin a)) 3))))
0.2b
(* (sin b) (sin a))

rewrite22.0ms

Algorithm
rewrite-expression-head
Counts
4 → 92
Calls

4 calls. Slowest were:

12.0ms
(/ (sin b) (- (* (cos a) (cos b)) (cbrt (pow (* (sin b) (sin a)) 3))))
4.0ms
(* (sin b) (sin a))
3.0ms
(pow (* (sin b) (sin a)) 3)

series340.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

125.0ms
(/ (sin b) (- (* (cos a) (cos b)) (cbrt (pow (* (sin b) (sin a)) 3))))
73.0ms
(* (sin b) (sin a))
72.0ms
(pow (* (sin b) (sin a)) 3)
70.0ms
(cbrt (pow (* (sin b) (sin a)) 3))

simplify4.1s

Counts
67 → 104
Calls

67 calls. Slowest were:

350.0ms
(cbrt (pow (- (cos (- b a)) (cos (+ b a))) 3))
304.0ms
(pow (- (cos (- b a)) (cos (+ b a))) 3)
292.0ms
(/ (* (* (sin b) (sin b)) (sin b)) (* (* (- (* (cos a) (cos b)) (cbrt (pow (* (sin b) (sin a)) 3))) (- (* (cos a) (cos b)) (cbrt (pow (* (sin b) (sin a)) 3)))) (- (* (cos a) (cos b)) (cbrt (pow (* (sin b) (sin a)) 3)))))

prune1.6s

Pruning

12 alts after pruning (11 fresh and 1 done)

Merged error: 0.0b

localize13.0ms

Local error

Found 4 expressions with local error:

0.4b
(expm1 (- (* (cos a) (cos b)) (* (sin a) (sin b))))
0.3b
(log1p (expm1 (- (* (cos a) (cos b)) (* (sin a) (sin b)))))
0.2b
(/ (sin b) (log1p (expm1 (- (* (cos a) (cos b)) (* (sin a) (sin b))))))
0.2b
(* (sin a) (sin b))

rewrite8.0ms

Algorithm
rewrite-expression-head
Counts
4 → 71
Calls

4 calls. Slowest were:

4.0ms
(* (sin a) (sin b))
3.0ms
(/ (sin b) (log1p (expm1 (- (* (cos a) (cos b)) (* (sin a) (sin b))))))
0.0ms
(expm1 (- (* (cos a) (cos b)) (* (sin a) (sin b))))

series553.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

175.0ms
(log1p (expm1 (- (* (cos a) (cos b)) (* (sin a) (sin b)))))
154.0ms
(/ (sin b) (log1p (expm1 (- (* (cos a) (cos b)) (* (sin a) (sin b))))))
149.0ms
(expm1 (- (* (cos a) (cos b)) (* (sin a) (sin b))))
75.0ms
(* (sin a) (sin b))

simplify1.5s

Counts
48 → 83
Calls

48 calls. Slowest were:

263.0ms
(- E (+ 1 (+ (* 1/2 (* (pow a 2) E)) (* a (* E b)))))
239.0ms
(+ b (+ (* 1/3 (pow b 3)) (* a (pow b 2))))
150.0ms
(/ (* (* (sin b) (sin b)) (sin b)) (* (* (log1p (expm1 (- (* (cos a) (cos b)) (* (sin a) (sin b))))) (log1p (expm1 (- (* (cos a) (cos b)) (* (sin a) (sin b)))))) (log1p (expm1 (- (* (cos a) (cos b)) (* (sin a) (sin b)))))))

prune1.3s

Pruning

13 alts after pruning (12 fresh and 1 done)

Merged error: 0.0b

regimes359.0ms

Accuracy

0% (0.3b remaining)

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

bsearch4.0ms

end0.0ms

sample7.7s

Algorithm
intervals