Average Error: 15.4 → 0.3
Time: 25.9s
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.4

    \[\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. Final simplification0.3

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

Reproduce

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

Details

Time bar (total: 24.5s)Debug log

sample292.0ms

Algorithm
intervals
Results
256×(pre true 80)
98×(body real 1280)
71×(body real 80)
51×(body real 640)
20×(body real 320)
16×(body real 160)

simplify25.0ms

Counts
1 → 1
Calls
1 calls:
Slowest
25.0ms
(/ (* r (sin b)) (cos (+ a b)))

prune14.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 15.7b

localize37.0ms

Local error

Found 3 expressions with local error:

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

rewrite453.0ms

Algorithm
rewrite-expression-head
Rules
add-cube-cbrt
*-un-lft-identity
add-sqr-sqrt
add-exp-log
add-cbrt-cube
add-log-exp
associate-*r*
times-frac
associate-/r*
pow1
associate-/l*
div-inv
*-commutative
div-exp
frac-2neg
clear-num
cos-sum
cbrt-undiv
Counts
3 → 38
Calls
3 calls:
Slowest
10.0ms
(/ (* r (sin b)) (cos (+ a b)))
6.0ms
(cos (+ a b))
2.0ms
(* r (sin b))

series176.0ms

Counts
3 → 9
Calls
3 calls:
Slowest
113.0ms
(/ (* r (sin b)) (cos (+ a b)))
38.0ms
(* r (sin b))
24.0ms
(cos (+ a b))

simplify945.0ms

Counts
24 → 47
Calls
24 calls:
Slowest
550.0ms
(/ (* (* (* r (sin b)) (* r (sin b))) (* r (sin b))) (* (* (cos (+ a b)) (cos (+ a b))) (cos (+ a b))))
206.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))
82.0ms
(- 1 (+ (* 1/2 (pow a 2)) (* a b)))
21.0ms
(- (log (* r (sin b))) (log (cos (+ a b))))
20.0ms
(/ (* (sin b) r) (cos (+ a b)))

prune543.0ms

Pruning

8 alts after pruning (8 fresh and 0 done)

Merged error: 0.2b

localize27.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.1b
(- (* (cos a) (cos b)) (* (sin a) (sin b)))
0.1b
(* (cos a) (cos b))

rewrite38.0ms

Algorithm
rewrite-expression-head
Rules
10×add-cube-cbrt
10×add-exp-log
10×add-cbrt-cube
10×*-un-lft-identity
10×add-sqr-sqrt
pow1
add-log-exp
associate-*l*
associate-*r*
sin-mult
cos-mult
associate-/r/
times-frac
associate-/r*
sub-div
flip--
cbrt-unprod
frac-sub
*-commutative
prod-exp
pow-prod-down
flip3--
associate-/l*
div-inv
div-exp
diff-log
frac-2neg
sub-neg
clear-num
cbrt-undiv
Counts
4 → 73
Calls
4 calls:
Slowest
16.0ms
(- (* (cos a) (cos b)) (* (sin a) (sin b)))
10.0ms
(/ (* r (sin b)) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
7.0ms
(* (cos a) (cos b))
3.0ms
(* (sin a) (sin b))

series548.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
326.0ms
(/ (* r (sin b)) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
112.0ms
(- (* (cos a) (cos b)) (* (sin a) (sin b)))
62.0ms
(* (cos a) (cos b))
49.0ms
(* (sin a) (sin b))

simplify2.0s

Counts
48 → 85
Calls
48 calls:
Slowest
426.0ms
(- (+ (cos (+ a b)) (cos (- a b))) (- (cos (- a b)) (cos (+ a b))))
300.0ms
(- (* (+ (cos (+ a b)) (cos (- a b))) 2) (* 2 (- (cos (- a b)) (cos (+ a b)))))
196.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))
191.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)))))
126.0ms
(- 1 (+ (* 1/2 (pow a 2)) (* a b)))

prune1.1s

Pruning

10 alts after pruning (10 fresh and 0 done)

Merged error: 0.0b

localize10.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)))

rewrite33.0ms

Algorithm
rewrite-expression-head
Rules
21×add-cube-cbrt
21×*-un-lft-identity
21×add-sqr-sqrt
20×associate-*r*
18×times-frac
associate-/r/
add-exp-log
add-cbrt-cube
sin-mult
add-log-exp
cos-mult
pow1
sub-div
associate-/l*
flip--
frac-sub
associate-*l*
flip3--
associate-/r*
div-inv
*-commutative
cbrt-unprod
associate-*r/
prod-exp
pow-prod-down
div-exp
diff-log
frac-2neg
sub-neg
clear-num
cbrt-undiv
Counts
4 → 89
Calls
4 calls:
Slowest
13.0ms
(* r (/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b)))))
9.0ms
(/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
6.0ms
(- (* (cos a) (cos b)) (* (sin a) (sin b)))
3.0ms
(* (sin a) (sin b))

series565.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
278.0ms
(* r (/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b)))))
135.0ms
(/ (sin b) (- (* (cos a) (cos b)) (* (sin a) (sin b))))
99.0ms
(- (* (cos a) (cos b)) (* (sin a) (sin b)))
52.0ms
(* (sin a) (sin b))

simplify2.4s

Counts
70 → 101
Calls
70 calls:
Slowest
386.0ms
(- (+ (cos (+ a b)) (cos (- a b))) (- (cos (- a b)) (cos (+ a b))))
284.0ms
(- (* (+ (cos (+ a b)) (cos (- a b))) 2) (* 2 (- (cos (- a b)) (cos (+ a b)))))
196.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)))))
184.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)))))
176.0ms
(+ (* a (* r (pow b 2))) (+ (* 1/3 (* r (pow b 3))) (* r b)))

prune1.4s

Pruning

10 alts after pruning (9 fresh and 1 done)

Merged error: 0.0b

localize12.0ms

Local error

Found 4 expressions with local error:

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

rewrite51.0ms

Algorithm
rewrite-expression-head
Rules
21×add-cube-cbrt
21×*-un-lft-identity
21×add-sqr-sqrt
18×times-frac
16×associate-/r*
add-exp-log
add-cbrt-cube
add-log-exp
pow1
sin-mult
associate-/l/
cos-mult
associate-/l*
div-inv
associate-*l*
associate-*r*
sub-div
flip--
frac-sub
flip3--
frac-2neg
clear-num
cbrt-unprod
*-commutative
associate-/r/
prod-exp
pow-prod-down
div-exp
div-sub
diff-log
sub-neg
cbrt-undiv
Counts
4 → 88
Calls
4 calls:
Slowest
25.0ms
(/ (- (* (cos a) (cos b)) (* (sin a) (sin b))) (sin b))
13.0ms
(- (* (cos a) (cos b)) (* (sin a) (sin b)))
7.0ms
(* (sin a) (sin b))
4.0ms
(/ r (/ (- (* (cos a) (cos b)) (* (sin a) (sin b))) (sin b)))

series668.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
307.0ms
(/ r (/ (- (* (cos a) (cos b)) (* (sin a) (sin b))) (sin b)))
169.0ms
(/ (- (* (cos a) (cos b)) (* (sin a) (sin b))) (sin b))
117.0ms
(- (* (cos a) (cos b)) (* (sin a) (sin b)))
74.0ms
(* (sin a) (sin b))

simplify3.9s

Counts
69 → 100
Calls
69 calls:
Slowest
349.0ms
(- (pow (* (cos a) (cos b)) 3) (pow (* (sin a) (sin b)) 3))
302.0ms
(- (* (+ (cos (+ a b)) (cos (- a b))) 2) (* 2 (- (cos (- a b)) (cos (+ a b)))))
301.0ms
(- (* (+ (cos (+ a b)) (cos (- a b))) 2) (* 2 (- (cos (- a b)) (cos (+ a b)))))
290.0ms
(- (+ (cos (+ a b)) (cos (- a b))) (- (cos (- a b)) (cos (+ a b))))
281.0ms
(- (+ (cos (+ a b)) (cos (- a b))) (- (cos (- a b)) (cos (+ a b))))

prune1.4s

Pruning

10 alts after pruning (8 fresh and 2 done)

Merged error: 0.0b

regimes352.0ms

Accuracy

0% (0.3b remaining)

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

bsearch3.0ms

end0.0ms

sample7.5s

Algorithm
intervals
Results
8000×(pre true 80)
2933×(body real 1280)
2098×(body real 80)
1786×(body real 640)
813×(body real 320)
370×(body real 160)