| Alternative 1 | |
|---|---|
| Accuracy | 60.5% |
| Cost | 26956 |
(FPCore (a b angle) :precision binary64 (* (* (* 2.0 (- (pow b 2.0) (pow a 2.0))) (sin (* PI (/ angle 180.0)))) (cos (* PI (/ angle 180.0)))))
(FPCore (a b angle)
:precision binary64
(let* ((t_0 (* PI (* 0.005555555555555556 angle))) (t_1 (sin t_0)))
(if (<= (- (pow b 2.0) (pow a 2.0)) 5e-243)
(* (cos (* 0.005555555555555556 (* PI angle))) (* -2.0 (* a (* a t_1))))
(* 2.0 (* (cos t_0) (* b (* b t_1)))))))double code(double a, double b, double angle) {
return ((2.0 * (pow(b, 2.0) - pow(a, 2.0))) * sin((((double) M_PI) * (angle / 180.0)))) * cos((((double) M_PI) * (angle / 180.0)));
}
double code(double a, double b, double angle) {
double t_0 = ((double) M_PI) * (0.005555555555555556 * angle);
double t_1 = sin(t_0);
double tmp;
if ((pow(b, 2.0) - pow(a, 2.0)) <= 5e-243) {
tmp = cos((0.005555555555555556 * (((double) M_PI) * angle))) * (-2.0 * (a * (a * t_1)));
} else {
tmp = 2.0 * (cos(t_0) * (b * (b * t_1)));
}
return tmp;
}
public static double code(double a, double b, double angle) {
return ((2.0 * (Math.pow(b, 2.0) - Math.pow(a, 2.0))) * Math.sin((Math.PI * (angle / 180.0)))) * Math.cos((Math.PI * (angle / 180.0)));
}
public static double code(double a, double b, double angle) {
double t_0 = Math.PI * (0.005555555555555556 * angle);
double t_1 = Math.sin(t_0);
double tmp;
if ((Math.pow(b, 2.0) - Math.pow(a, 2.0)) <= 5e-243) {
tmp = Math.cos((0.005555555555555556 * (Math.PI * angle))) * (-2.0 * (a * (a * t_1)));
} else {
tmp = 2.0 * (Math.cos(t_0) * (b * (b * t_1)));
}
return tmp;
}
def code(a, b, angle): return ((2.0 * (math.pow(b, 2.0) - math.pow(a, 2.0))) * math.sin((math.pi * (angle / 180.0)))) * math.cos((math.pi * (angle / 180.0)))
def code(a, b, angle): t_0 = math.pi * (0.005555555555555556 * angle) t_1 = math.sin(t_0) tmp = 0 if (math.pow(b, 2.0) - math.pow(a, 2.0)) <= 5e-243: tmp = math.cos((0.005555555555555556 * (math.pi * angle))) * (-2.0 * (a * (a * t_1))) else: tmp = 2.0 * (math.cos(t_0) * (b * (b * t_1))) return tmp
function code(a, b, angle) return Float64(Float64(Float64(2.0 * Float64((b ^ 2.0) - (a ^ 2.0))) * sin(Float64(pi * Float64(angle / 180.0)))) * cos(Float64(pi * Float64(angle / 180.0)))) end
function code(a, b, angle) t_0 = Float64(pi * Float64(0.005555555555555556 * angle)) t_1 = sin(t_0) tmp = 0.0 if (Float64((b ^ 2.0) - (a ^ 2.0)) <= 5e-243) tmp = Float64(cos(Float64(0.005555555555555556 * Float64(pi * angle))) * Float64(-2.0 * Float64(a * Float64(a * t_1)))); else tmp = Float64(2.0 * Float64(cos(t_0) * Float64(b * Float64(b * t_1)))); end return tmp end
function tmp = code(a, b, angle) tmp = ((2.0 * ((b ^ 2.0) - (a ^ 2.0))) * sin((pi * (angle / 180.0)))) * cos((pi * (angle / 180.0))); end
function tmp_2 = code(a, b, angle) t_0 = pi * (0.005555555555555556 * angle); t_1 = sin(t_0); tmp = 0.0; if (((b ^ 2.0) - (a ^ 2.0)) <= 5e-243) tmp = cos((0.005555555555555556 * (pi * angle))) * (-2.0 * (a * (a * t_1))); else tmp = 2.0 * (cos(t_0) * (b * (b * t_1))); end tmp_2 = tmp; end
code[a_, b_, angle_] := N[(N[(N[(2.0 * N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sin[N[(Pi * N[(angle / 180.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(Pi * N[(angle / 180.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
code[a_, b_, angle_] := Block[{t$95$0 = N[(Pi * N[(0.005555555555555556 * angle), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Sin[t$95$0], $MachinePrecision]}, If[LessEqual[N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a, 2.0], $MachinePrecision]), $MachinePrecision], 5e-243], N[(N[Cos[N[(0.005555555555555556 * N[(Pi * angle), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(-2.0 * N[(a * N[(a * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[Cos[t$95$0], $MachinePrecision] * N[(b * N[(b * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\begin{array}{l}
t_0 := \pi \cdot \left(0.005555555555555556 \cdot angle\right)\\
t_1 := \sin t_0\\
\mathbf{if}\;{b}^{2} - {a}^{2} \leq 5 \cdot 10^{-243}:\\
\;\;\;\;\cos \left(0.005555555555555556 \cdot \left(\pi \cdot angle\right)\right) \cdot \left(-2 \cdot \left(a \cdot \left(a \cdot t_1\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\cos t_0 \cdot \left(b \cdot \left(b \cdot t_1\right)\right)\right)\\
\end{array}
Results
if (-.f64 (pow.f64 b 2) (pow.f64 a 2)) < 5e-243Initial program 55.0%
Simplified55.0%
[Start]55.0 | \[ \left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
|---|---|
*-commutative [=>]55.0 | \[ \left(\color{blue}{\left(\left({b}^{2} - {a}^{2}\right) \cdot 2\right)} \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
associate-*l* [=>]55.0 | \[ \color{blue}{\left(\left({b}^{2} - {a}^{2}\right) \cdot \left(2 \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)\right)} \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
unpow2 [=>]55.0 | \[ \left(\left(\color{blue}{b \cdot b} - {a}^{2}\right) \cdot \left(2 \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
fma-neg [=>]55.0 | \[ \left(\color{blue}{\mathsf{fma}\left(b, b, -{a}^{2}\right)} \cdot \left(2 \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
unpow2 [=>]55.0 | \[ \left(\mathsf{fma}\left(b, b, -\color{blue}{a \cdot a}\right) \cdot \left(2 \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
distribute-rgt-neg-in [=>]55.0 | \[ \left(\mathsf{fma}\left(b, b, \color{blue}{a \cdot \left(-a\right)}\right) \cdot \left(2 \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
Taylor expanded in b around 0 54.3%
Simplified54.3%
[Start]54.3 | \[ \left(-2 \cdot \left({a}^{2} \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
|---|---|
*-commutative [=>]54.3 | \[ \color{blue}{\left(\left({a}^{2} \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right) \cdot -2\right)} \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
associate-*l* [=>]54.3 | \[ \color{blue}{\left({a}^{2} \cdot \left(\sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right) \cdot -2\right)\right)} \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
unpow2 [=>]54.3 | \[ \left(\color{blue}{\left(a \cdot a\right)} \cdot \left(\sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right) \cdot -2\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
*-commutative [=>]54.3 | \[ \left(\left(a \cdot a\right) \cdot \color{blue}{\left(-2 \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right)}\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
Taylor expanded in a around 0 54.3%
Simplified67.9%
[Start]54.3 | \[ \left(-2 \cdot \left({a}^{2} \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
|---|---|
unpow2 [=>]54.3 | \[ \left(-2 \cdot \left(\color{blue}{\left(a \cdot a\right)} \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
associate-*l* [=>]67.9 | \[ \left(-2 \cdot \color{blue}{\left(a \cdot \left(a \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right)\right)}\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
associate-*r* [=>]67.9 | \[ \left(-2 \cdot \left(a \cdot \left(a \cdot \sin \color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \pi\right)}\right)\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
*-commutative [=>]67.9 | \[ \left(-2 \cdot \left(a \cdot \left(a \cdot \sin \color{blue}{\left(\pi \cdot \left(0.005555555555555556 \cdot angle\right)\right)}\right)\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
Taylor expanded in angle around inf 67.8%
if 5e-243 < (-.f64 (pow.f64 b 2) (pow.f64 a 2)) Initial program 45.9%
Simplified45.9%
[Start]45.9 | \[ \left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
|---|---|
*-commutative [=>]45.9 | \[ \color{blue}{\left(\sin \left(\pi \cdot \frac{angle}{180}\right) \cdot \left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right)\right)} \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)
\] |
associate-*l* [=>]45.9 | \[ \color{blue}{\sin \left(\pi \cdot \frac{angle}{180}\right) \cdot \left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}
\] |
unpow2 [=>]45.9 | \[ \sin \left(\pi \cdot \frac{angle}{180}\right) \cdot \left(\left(2 \cdot \left(\color{blue}{b \cdot b} - {a}^{2}\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)
\] |
fma-neg [=>]45.9 | \[ \sin \left(\pi \cdot \frac{angle}{180}\right) \cdot \left(\left(2 \cdot \color{blue}{\mathsf{fma}\left(b, b, -{a}^{2}\right)}\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)
\] |
unpow2 [=>]45.9 | \[ \sin \left(\pi \cdot \frac{angle}{180}\right) \cdot \left(\left(2 \cdot \mathsf{fma}\left(b, b, -\color{blue}{a \cdot a}\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)
\] |
Applied egg-rr23.7%
[Start]45.9 | \[ \sin \left(\pi \cdot \frac{angle}{180}\right) \cdot \left(\left(2 \cdot \mathsf{fma}\left(b, b, -a \cdot a\right)\right) \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)
\] |
|---|---|
add-sqr-sqrt [=>]23.7 | \[ \sin \left(\pi \cdot \frac{angle}{180}\right) \cdot \left(\left(2 \cdot \mathsf{fma}\left(b, b, -a \cdot a\right)\right) \cdot \cos \color{blue}{\left(\sqrt{\pi \cdot \frac{angle}{180}} \cdot \sqrt{\pi \cdot \frac{angle}{180}}\right)}\right)
\] |
pow2 [=>]23.7 | \[ \sin \left(\pi \cdot \frac{angle}{180}\right) \cdot \left(\left(2 \cdot \mathsf{fma}\left(b, b, -a \cdot a\right)\right) \cdot \cos \color{blue}{\left({\left(\sqrt{\pi \cdot \frac{angle}{180}}\right)}^{2}\right)}\right)
\] |
div-inv [=>]23.7 | \[ \sin \left(\pi \cdot \frac{angle}{180}\right) \cdot \left(\left(2 \cdot \mathsf{fma}\left(b, b, -a \cdot a\right)\right) \cdot \cos \left({\left(\sqrt{\pi \cdot \color{blue}{\left(angle \cdot \frac{1}{180}\right)}}\right)}^{2}\right)\right)
\] |
metadata-eval [=>]23.7 | \[ \sin \left(\pi \cdot \frac{angle}{180}\right) \cdot \left(\left(2 \cdot \mathsf{fma}\left(b, b, -a \cdot a\right)\right) \cdot \cos \left({\left(\sqrt{\pi \cdot \left(angle \cdot \color{blue}{0.005555555555555556}\right)}\right)}^{2}\right)\right)
\] |
Taylor expanded in b around inf 45.5%
Simplified63.5%
[Start]45.5 | \[ 2 \cdot \left({b}^{2} \cdot \left(\cos \left(angle \cdot \left({\left(\sqrt{0.005555555555555556}\right)}^{2} \cdot \pi\right)\right) \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right)\right)
\] |
|---|---|
*-commutative [<=]45.5 | \[ 2 \cdot \left({b}^{2} \cdot \color{blue}{\left(\sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right) \cdot \cos \left(angle \cdot \left({\left(\sqrt{0.005555555555555556}\right)}^{2} \cdot \pi\right)\right)\right)}\right)
\] |
associate-*r* [=>]45.5 | \[ 2 \cdot \color{blue}{\left(\left({b}^{2} \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right) \cdot \cos \left(angle \cdot \left({\left(\sqrt{0.005555555555555556}\right)}^{2} \cdot \pi\right)\right)\right)}
\] |
*-commutative [=>]45.5 | \[ 2 \cdot \color{blue}{\left(\cos \left(angle \cdot \left({\left(\sqrt{0.005555555555555556}\right)}^{2} \cdot \pi\right)\right) \cdot \left({b}^{2} \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right)\right)}
\] |
associate-*r* [=>]45.6 | \[ 2 \cdot \left(\cos \color{blue}{\left(\left(angle \cdot {\left(\sqrt{0.005555555555555556}\right)}^{2}\right) \cdot \pi\right)} \cdot \left({b}^{2} \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right)\right)
\] |
unpow2 [=>]45.6 | \[ 2 \cdot \left(\cos \left(\left(angle \cdot \color{blue}{\left(\sqrt{0.005555555555555556} \cdot \sqrt{0.005555555555555556}\right)}\right) \cdot \pi\right) \cdot \left({b}^{2} \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right)\right)
\] |
rem-square-sqrt [=>]45.6 | \[ 2 \cdot \left(\cos \left(\left(angle \cdot \color{blue}{0.005555555555555556}\right) \cdot \pi\right) \cdot \left({b}^{2} \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right)\right)
\] |
*-commutative [<=]45.6 | \[ 2 \cdot \left(\cos \color{blue}{\left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)} \cdot \left({b}^{2} \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right)\right)
\] |
unpow2 [=>]45.6 | \[ 2 \cdot \left(\cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right) \cdot \left(\color{blue}{\left(b \cdot b\right)} \cdot \sin \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right)\right)
\] |
associate-*r* [=>]45.7 | \[ 2 \cdot \left(\cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right) \cdot \left(\left(b \cdot b\right) \cdot \sin \color{blue}{\left(\left(0.005555555555555556 \cdot angle\right) \cdot \pi\right)}\right)\right)
\] |
*-commutative [<=]45.7 | \[ 2 \cdot \left(\cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right) \cdot \left(\left(b \cdot b\right) \cdot \sin \left(\color{blue}{\left(angle \cdot 0.005555555555555556\right)} \cdot \pi\right)\right)\right)
\] |
*-commutative [<=]45.7 | \[ 2 \cdot \left(\cos \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right) \cdot \left(\left(b \cdot b\right) \cdot \sin \color{blue}{\left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)}\right)\right)
\] |
Final simplification66.0%
| Alternative 1 | |
|---|---|
| Accuracy | 60.5% |
| Cost | 26956 |
| Alternative 2 | |
|---|---|
| Accuracy | 59.8% |
| Cost | 26956 |
| Alternative 3 | |
|---|---|
| Accuracy | 59.9% |
| Cost | 14096 |
| Alternative 4 | |
|---|---|
| Accuracy | 59.9% |
| Cost | 14096 |
| Alternative 5 | |
|---|---|
| Accuracy | 55.0% |
| Cost | 13704 |
| Alternative 6 | |
|---|---|
| Accuracy | 53.5% |
| Cost | 7432 |
| Alternative 7 | |
|---|---|
| Accuracy | 53.5% |
| Cost | 7432 |
| Alternative 8 | |
|---|---|
| Accuracy | 53.5% |
| Cost | 7432 |
| Alternative 9 | |
|---|---|
| Accuracy | 39.9% |
| Cost | 7177 |
| Alternative 10 | |
|---|---|
| Accuracy | 39.9% |
| Cost | 7177 |
| Alternative 11 | |
|---|---|
| Accuracy | 48.3% |
| Cost | 7177 |
| Alternative 12 | |
|---|---|
| Accuracy | 49.0% |
| Cost | 7176 |
| Alternative 13 | |
|---|---|
| Accuracy | 48.9% |
| Cost | 7176 |
| Alternative 14 | |
|---|---|
| Accuracy | 49.0% |
| Cost | 7176 |
| Alternative 15 | |
|---|---|
| Accuracy | 32.3% |
| Cost | 6912 |
herbie shell --seed 2023147
(FPCore (a b angle)
:name "ab-angle->ABCF B"
:precision binary64
(* (* (* 2.0 (- (pow b 2.0) (pow a 2.0))) (sin (* PI (/ angle 180.0)))) (cos (* PI (/ angle 180.0)))))