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

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. 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. Taylor expanded around -inf 0.3

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

    \[\leadsto \frac{\sin b \cdot r}{\cos a \cdot \cos b - \color{blue}{(e^{\log_* (1 + \sin a \cdot \sin b)} - 1)^*}}\]
  7. Final simplification0.3

    \[\leadsto \frac{r \cdot \sin b}{\cos a \cdot \cos b - (e^{\log_* (1 + \sin a \cdot \sin b)} - 1)^*}\]

Reproduce

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

sample575.0ms

Algorithm
intervals

simplify77.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

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

prune14.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 14.3b

localize35.0ms

Local error

Found 3 expressions with local error:

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

rewrite20.0ms

Algorithm
rewrite-expression-head
Counts
3 → 63
Calls

3 calls. Slowest were:

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

series192.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

93.0ms
(* r (/ (sin b) (cos (+ a b))))
67.0ms
(/ (sin b) (cos (+ a b)))
31.0ms
(cos (+ a b))

simplify1.1s

Counts
49 → 72
Calls

49 calls. Slowest were:

349.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))
215.0ms
(+ b (+ (* 1/3 (pow b 3)) (* a (pow b 2))))
190.0ms
(- 1 (+ (* 1/2 (pow a 2)) (* a b)))

prune845.0ms

Pruning

10 alts after pruning (10 fresh and 0 done)

Merged error: 0.2b

localize50.0ms

Local error

Found 4 expressions with local error:

0.3b
(* (sin a) (sin b))
0.2b
(/ (sin b) (- (* (cos a) (cos b)) (* (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)))

rewrite41.0ms

Algorithm
rewrite-expression-head
Counts
4 → 99
Calls

4 calls. Slowest were:

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

series556.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

274.0ms
(* r (/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b)))))
127.0ms
(/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
86.0ms
(- (* (cos a) (cos b)) (* (sin a) (sin b)))
69.0ms
(* (sin a) (sin b))

simplify3.0s

Counts
73 → 111
Calls

73 calls. Slowest were:

394.0ms
(- (+ (cos (+ a b)) (cos (- a b))) (- (cos (- a b)) (cos (+ a b))))
364.0ms
(- (* (+ (cos (+ a b)) (cos (- a b))) 2) (* 2 (- (cos (- a b)) (cos (+ a b)))))
297.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))

prune1.6s

Pruning

14 alts after pruning (14 fresh and 0 done)

Merged error: 0.0b

localize6.0ms

Local error

Found 4 expressions with local error:

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

rewrite47.0ms

Algorithm
rewrite-expression-head
Counts
4 → 76
Calls

4 calls. Slowest were:

21.0ms
(/ (* (sin b) r) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
14.0ms
(- (* (cos a) (cos b)) (* (sin a) (sin b)))
6.0ms
(* (sin a) (sin b))

series475.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

257.0ms
(/ (* (sin b) r) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
102.0ms
(- (* (cos a) (cos b)) (* (sin a) (sin b)))
78.0ms
(* (sin a) (sin b))
37.0ms
(* (sin b) r)

simplify2.0s

Counts
45 → 88
Calls

45 calls. Slowest were:

373.0ms
(- (* (+ (cos (+ a b)) (cos (- a b))) 2) (* 2 (- (cos (- a b)) (cos (+ a b)))))
359.0ms
(- (+ (cos (+ a b)) (cos (- a b))) (- (cos (- a b)) (cos (+ a b))))
281.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))

prune1.2s

Pruning

13 alts after pruning (13 fresh and 0 done)

Merged error: 0.0b

localize9.0ms

Local error

Found 4 expressions with local error:

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

rewrite11.0ms

Algorithm
rewrite-expression-head
Counts
4 → 64
Calls

4 calls. Slowest were:

4.0ms
(/ (* (sin b) r) (fma (cos a) (cos b) (- (* (sin a) (sin b)))))
4.0ms
(* (sin a) (sin b))
2.0ms
(* (sin b) r)

series520.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

298.0ms
(/ (* (sin b) r) (fma (cos a) (cos b) (- (* (sin a) (sin b)))))
99.0ms
(fma (cos a) (cos b) (- (* (sin a) (sin b))))
75.0ms
(* (sin a) (sin b))
48.0ms
(* (sin b) r)

simplify956.0ms

Counts
35 → 76
Calls

35 calls. Slowest were:

229.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))
181.0ms
(- 1 (+ (* 1/2 (pow a 2)) (* a b)))
99.0ms
(* (* (* (sin a) (sin a)) (sin a)) (* (* (sin b) (sin b)) (sin b)))

prune1.1s

Pruning

11 alts after pruning (11 fresh and 0 done)

Merged error: 0.0b

regimes315.0ms

Accuracy

0% (0.3b remaining)

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

bsearch2.0ms

end0.0ms

sample7.3s

Algorithm
intervals