| Alternative 1 | |
|---|---|
| Accuracy | 99.5% |
| Cost | 39040 |
\[r \cdot \frac{\sin b}{\mathsf{fma}\left(\cos b, \cos a, \sin b \cdot \left(-\sin a\right)\right)}
\]
(FPCore (r a b) :precision binary64 (* r (/ (sin b) (cos (+ a b)))))
(FPCore (r a b) :precision binary64 (* (/ (sin b) (fma (sin b) (- (sin a)) (* (cos a) (cos b)))) r))
double code(double r, double a, double b) {
return r * (sin(b) / cos((a + b)));
}
double code(double r, double a, double b) {
return (sin(b) / fma(sin(b), -sin(a), (cos(a) * cos(b)))) * r;
}
function code(r, a, b) return Float64(r * Float64(sin(b) / cos(Float64(a + b)))) end
function code(r, a, b) return Float64(Float64(sin(b) / fma(sin(b), Float64(-sin(a)), Float64(cos(a) * cos(b)))) * r) end
code[r_, a_, b_] := N[(r * N[(N[Sin[b], $MachinePrecision] / N[Cos[N[(a + b), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[r_, a_, b_] := N[(N[(N[Sin[b], $MachinePrecision] / N[(N[Sin[b], $MachinePrecision] * (-N[Sin[a], $MachinePrecision]) + N[(N[Cos[a], $MachinePrecision] * N[Cos[b], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * r), $MachinePrecision]
r \cdot \frac{\sin b}{\cos \left(a + b\right)}
\frac{\sin b}{\mathsf{fma}\left(\sin b, -\sin a, \cos a \cdot \cos b\right)} \cdot r
Initial program 74.6%
Simplified74.6%
[Start]74.6 | \[ r \cdot \frac{\sin b}{\cos \left(a + b\right)}
\] |
|---|---|
+-commutative [=>]74.6 | \[ r \cdot \frac{\sin b}{\cos \color{blue}{\left(b + a\right)}}
\] |
Applied egg-rr99.5%
[Start]74.6 | \[ r \cdot \frac{\sin b}{\cos \left(b + a\right)}
\] |
|---|---|
cos-sum [=>]99.5 | \[ r \cdot \frac{\sin b}{\color{blue}{\cos b \cdot \cos a - \sin b \cdot \sin a}}
\] |
cancel-sign-sub-inv [=>]99.5 | \[ r \cdot \frac{\sin b}{\color{blue}{\cos b \cdot \cos a + \left(-\sin b\right) \cdot \sin a}}
\] |
fma-def [=>]99.5 | \[ r \cdot \frac{\sin b}{\color{blue}{\mathsf{fma}\left(\cos b, \cos a, \left(-\sin b\right) \cdot \sin a\right)}}
\] |
Taylor expanded in r around 0 99.4%
Simplified99.5%
[Start]99.4 | \[ \frac{\sin b \cdot r}{-1 \cdot \left(\sin a \cdot \sin b\right) + \cos a \cdot \cos b}
\] |
|---|---|
associate-/l* [=>]99.4 | \[ \color{blue}{\frac{\sin b}{\frac{-1 \cdot \left(\sin a \cdot \sin b\right) + \cos a \cdot \cos b}{r}}}
\] |
associate-/r/ [=>]99.5 | \[ \color{blue}{\frac{\sin b}{-1 \cdot \left(\sin a \cdot \sin b\right) + \cos a \cdot \cos b} \cdot r}
\] |
mul-1-neg [=>]99.5 | \[ \frac{\sin b}{\color{blue}{\left(-\sin a \cdot \sin b\right)} + \cos a \cdot \cos b} \cdot r
\] |
distribute-lft-neg-out [<=]99.5 | \[ \frac{\sin b}{\color{blue}{\left(-\sin a\right) \cdot \sin b} + \cos a \cdot \cos b} \cdot r
\] |
+-commutative [=>]99.5 | \[ \frac{\sin b}{\color{blue}{\cos a \cdot \cos b + \left(-\sin a\right) \cdot \sin b}} \cdot r
\] |
*-commutative [=>]99.5 | \[ \frac{\sin b}{\color{blue}{\cos b \cdot \cos a} + \left(-\sin a\right) \cdot \sin b} \cdot r
\] |
cancel-sign-sub-inv [<=]99.5 | \[ \frac{\sin b}{\color{blue}{\cos b \cdot \cos a - \sin a \cdot \sin b}} \cdot r
\] |
*-commutative [<=]99.5 | \[ \frac{\sin b}{\color{blue}{\cos a \cdot \cos b} - \sin a \cdot \sin b} \cdot r
\] |
Applied egg-rr99.5%
[Start]99.5 | \[ \frac{\sin b}{\cos a \cdot \cos b - \sin a \cdot \sin b} \cdot r
\] |
|---|---|
sub-neg [=>]99.5 | \[ \frac{\sin b}{\color{blue}{\cos a \cdot \cos b + \left(-\sin a \cdot \sin b\right)}} \cdot r
\] |
+-commutative [=>]99.5 | \[ \frac{\sin b}{\color{blue}{\left(-\sin a \cdot \sin b\right) + \cos a \cdot \cos b}} \cdot r
\] |
*-commutative [=>]99.5 | \[ \frac{\sin b}{\left(-\color{blue}{\sin b \cdot \sin a}\right) + \cos a \cdot \cos b} \cdot r
\] |
distribute-rgt-neg-in [=>]99.5 | \[ \frac{\sin b}{\color{blue}{\sin b \cdot \left(-\sin a\right)} + \cos a \cdot \cos b} \cdot r
\] |
fma-def [=>]99.5 | \[ \frac{\sin b}{\color{blue}{\mathsf{fma}\left(\sin b, -\sin a, \cos a \cdot \cos b\right)}} \cdot r
\] |
Final simplification99.5%
| Alternative 1 | |
|---|---|
| Accuracy | 99.5% |
| Cost | 39040 |
| Alternative 2 | |
|---|---|
| Accuracy | 99.5% |
| Cost | 32704 |
| Alternative 3 | |
|---|---|
| Accuracy | 98.7% |
| Cost | 26432 |
| Alternative 4 | |
|---|---|
| Accuracy | 98.8% |
| Cost | 26432 |
| Alternative 5 | |
|---|---|
| Accuracy | 77.5% |
| Cost | 26048 |
| Alternative 6 | |
|---|---|
| Accuracy | 76.4% |
| Cost | 13385 |
| Alternative 7 | |
|---|---|
| Accuracy | 76.4% |
| Cost | 13385 |
| Alternative 8 | |
|---|---|
| Accuracy | 76.4% |
| Cost | 13384 |
| Alternative 9 | |
|---|---|
| Accuracy | 76.4% |
| Cost | 13384 |
| Alternative 10 | |
|---|---|
| Accuracy | 76.4% |
| Cost | 13384 |
| Alternative 11 | |
|---|---|
| Accuracy | 76.6% |
| Cost | 13248 |
| Alternative 12 | |
|---|---|
| Accuracy | 76.6% |
| Cost | 13248 |
| Alternative 13 | |
|---|---|
| Accuracy | 54.9% |
| Cost | 13120 |
| Alternative 14 | |
|---|---|
| Accuracy | 55.1% |
| Cost | 6985 |
| Alternative 15 | |
|---|---|
| Accuracy | 55.1% |
| Cost | 6985 |
| Alternative 16 | |
|---|---|
| Accuracy | 55.1% |
| Cost | 6985 |
| Alternative 17 | |
|---|---|
| Accuracy | 38.7% |
| Cost | 6592 |
| Alternative 18 | |
|---|---|
| Accuracy | 34.3% |
| Cost | 192 |
herbie shell --seed 2023157
(FPCore (r a b)
:name "rsin B (should all be same)"
:precision binary64
(* r (/ (sin b) (cos (+ a b)))))