
(FPCore (a b angle) :precision binary64 (let* ((t_0 (* PI (/ angle 180.0)))) (* (* (* 2.0 (- (pow b 2.0) (pow a 2.0))) (sin t_0)) (cos t_0))))
double code(double a, double b, double angle) {
double t_0 = ((double) M_PI) * (angle / 180.0);
return ((2.0 * (pow(b, 2.0) - pow(a, 2.0))) * sin(t_0)) * cos(t_0);
}
public static double code(double a, double b, double angle) {
double t_0 = Math.PI * (angle / 180.0);
return ((2.0 * (Math.pow(b, 2.0) - Math.pow(a, 2.0))) * Math.sin(t_0)) * Math.cos(t_0);
}
def code(a, b, angle): t_0 = math.pi * (angle / 180.0) return ((2.0 * (math.pow(b, 2.0) - math.pow(a, 2.0))) * math.sin(t_0)) * math.cos(t_0)
function code(a, b, angle) t_0 = Float64(pi * Float64(angle / 180.0)) return Float64(Float64(Float64(2.0 * Float64((b ^ 2.0) - (a ^ 2.0))) * sin(t_0)) * cos(t_0)) end
function tmp = code(a, b, angle) t_0 = pi * (angle / 180.0); tmp = ((2.0 * ((b ^ 2.0) - (a ^ 2.0))) * sin(t_0)) * cos(t_0); end
code[a_, b_, angle_] := Block[{t$95$0 = N[(Pi * N[(angle / 180.0), $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(2.0 * N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sin[t$95$0], $MachinePrecision]), $MachinePrecision] * N[Cos[t$95$0], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \pi \cdot \frac{angle}{180}\\
\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin t\_0\right) \cdot \cos t\_0
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 15 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a b angle) :precision binary64 (let* ((t_0 (* PI (/ angle 180.0)))) (* (* (* 2.0 (- (pow b 2.0) (pow a 2.0))) (sin t_0)) (cos t_0))))
double code(double a, double b, double angle) {
double t_0 = ((double) M_PI) * (angle / 180.0);
return ((2.0 * (pow(b, 2.0) - pow(a, 2.0))) * sin(t_0)) * cos(t_0);
}
public static double code(double a, double b, double angle) {
double t_0 = Math.PI * (angle / 180.0);
return ((2.0 * (Math.pow(b, 2.0) - Math.pow(a, 2.0))) * Math.sin(t_0)) * Math.cos(t_0);
}
def code(a, b, angle): t_0 = math.pi * (angle / 180.0) return ((2.0 * (math.pow(b, 2.0) - math.pow(a, 2.0))) * math.sin(t_0)) * math.cos(t_0)
function code(a, b, angle) t_0 = Float64(pi * Float64(angle / 180.0)) return Float64(Float64(Float64(2.0 * Float64((b ^ 2.0) - (a ^ 2.0))) * sin(t_0)) * cos(t_0)) end
function tmp = code(a, b, angle) t_0 = pi * (angle / 180.0); tmp = ((2.0 * ((b ^ 2.0) - (a ^ 2.0))) * sin(t_0)) * cos(t_0); end
code[a_, b_, angle_] := Block[{t$95$0 = N[(Pi * N[(angle / 180.0), $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(2.0 * N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sin[t$95$0], $MachinePrecision]), $MachinePrecision] * N[Cos[t$95$0], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \pi \cdot \frac{angle}{180}\\
\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin t\_0\right) \cdot \cos t\_0
\end{array}
\end{array}
angle\_m = (fabs.f64 angle)
angle\_s = (copysign.f64 #s(literal 1 binary64) angle)
(FPCore (angle_s a b angle_m)
:precision binary64
(let* ((t_0 (* (+ b a) (- b a))))
(*
angle_s
(if (<= (/ angle_m 180.0) 1e+42)
(* (+ b a) (* (- b a) (sin (/ (* (* angle_m PI) 360.0) 32400.0))))
(if (<= (/ angle_m 180.0) 1e+154)
(* t_0 (sin (pow (cbrt (* (* angle_m PI) 0.011111111111111112)) 3.0)))
(*
t_0
(sin
(* PI (* (exp (log (* angle_m 0.005555555555555556))) 2.0)))))))))angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
double t_0 = (b + a) * (b - a);
double tmp;
if ((angle_m / 180.0) <= 1e+42) {
tmp = (b + a) * ((b - a) * sin((((angle_m * ((double) M_PI)) * 360.0) / 32400.0)));
} else if ((angle_m / 180.0) <= 1e+154) {
tmp = t_0 * sin(pow(cbrt(((angle_m * ((double) M_PI)) * 0.011111111111111112)), 3.0));
} else {
tmp = t_0 * sin((((double) M_PI) * (exp(log((angle_m * 0.005555555555555556))) * 2.0)));
}
return angle_s * tmp;
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
double t_0 = (b + a) * (b - a);
double tmp;
if ((angle_m / 180.0) <= 1e+42) {
tmp = (b + a) * ((b - a) * Math.sin((((angle_m * Math.PI) * 360.0) / 32400.0)));
} else if ((angle_m / 180.0) <= 1e+154) {
tmp = t_0 * Math.sin(Math.pow(Math.cbrt(((angle_m * Math.PI) * 0.011111111111111112)), 3.0));
} else {
tmp = t_0 * Math.sin((Math.PI * (Math.exp(Math.log((angle_m * 0.005555555555555556))) * 2.0)));
}
return angle_s * tmp;
}
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) t_0 = Float64(Float64(b + a) * Float64(b - a)) tmp = 0.0 if (Float64(angle_m / 180.0) <= 1e+42) tmp = Float64(Float64(b + a) * Float64(Float64(b - a) * sin(Float64(Float64(Float64(angle_m * pi) * 360.0) / 32400.0)))); elseif (Float64(angle_m / 180.0) <= 1e+154) tmp = Float64(t_0 * sin((cbrt(Float64(Float64(angle_m * pi) * 0.011111111111111112)) ^ 3.0))); else tmp = Float64(t_0 * sin(Float64(pi * Float64(exp(log(Float64(angle_m * 0.005555555555555556))) * 2.0)))); end return Float64(angle_s * tmp) end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := Block[{t$95$0 = N[(N[(b + a), $MachinePrecision] * N[(b - a), $MachinePrecision]), $MachinePrecision]}, N[(angle$95$s * If[LessEqual[N[(angle$95$m / 180.0), $MachinePrecision], 1e+42], N[(N[(b + a), $MachinePrecision] * N[(N[(b - a), $MachinePrecision] * N[Sin[N[(N[(N[(angle$95$m * Pi), $MachinePrecision] * 360.0), $MachinePrecision] / 32400.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(angle$95$m / 180.0), $MachinePrecision], 1e+154], N[(t$95$0 * N[Sin[N[Power[N[Power[N[(N[(angle$95$m * Pi), $MachinePrecision] * 0.011111111111111112), $MachinePrecision], 1/3], $MachinePrecision], 3.0], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(t$95$0 * N[Sin[N[(Pi * N[(N[Exp[N[Log[N[(angle$95$m * 0.005555555555555556), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
\begin{array}{l}
t_0 := \left(b + a\right) \cdot \left(b - a\right)\\
angle\_s \cdot \begin{array}{l}
\mathbf{if}\;\frac{angle\_m}{180} \leq 10^{+42}:\\
\;\;\;\;\left(b + a\right) \cdot \left(\left(b - a\right) \cdot \sin \left(\frac{\left(angle\_m \cdot \pi\right) \cdot 360}{32400}\right)\right)\\
\mathbf{elif}\;\frac{angle\_m}{180} \leq 10^{+154}:\\
\;\;\;\;t\_0 \cdot \sin \left({\left(\sqrt[3]{\left(angle\_m \cdot \pi\right) \cdot 0.011111111111111112}\right)}^{3}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot \sin \left(\pi \cdot \left(e^{\log \left(angle\_m \cdot 0.005555555555555556\right)} \cdot 2\right)\right)\\
\end{array}
\end{array}
\end{array}
if (/.f64 angle #s(literal 180 binary64)) < 1.00000000000000004e42Initial program 56.8%
associate-*l*56.8%
*-commutative56.8%
associate-*l*56.8%
Simplified56.8%
associate-*r*56.8%
*-commutative56.8%
associate-*l*56.8%
add-cbrt-cube39.6%
pow1/325.2%
Applied egg-rr25.3%
unpow1/340.7%
rem-cbrt-cube57.7%
unpow257.7%
unpow257.7%
difference-of-squares64.4%
*-commutative64.4%
metadata-eval64.4%
div-inv64.0%
2-sin64.0%
associate-*l*73.9%
2-sin73.9%
div-inv74.2%
metadata-eval74.2%
Applied egg-rr74.2%
add-exp-log30.4%
Applied egg-rr30.4%
add-log-exp5.4%
associate-*r*5.4%
rem-exp-log11.4%
*-commutative11.4%
exp-lft-sqr11.4%
sum-log11.4%
add-log-exp28.1%
add-sqr-sqrt28.1%
unpow228.1%
*-commutative28.1%
metadata-eval28.1%
div-inv28.1%
add-log-exp73.2%
*-commutative73.2%
metadata-eval73.2%
div-inv73.0%
Applied egg-rr74.8%
*-commutative74.8%
*-commutative74.8%
*-commutative74.8%
distribute-rgt-out74.8%
metadata-eval74.8%
Simplified74.8%
if 1.00000000000000004e42 < (/.f64 angle #s(literal 180 binary64)) < 1.00000000000000004e154Initial program 32.2%
associate-*l*32.2%
*-commutative32.2%
associate-*l*32.2%
Simplified32.2%
unpow232.2%
unpow232.2%
difference-of-squares35.8%
Applied egg-rr35.8%
2-sin35.8%
div-inv29.4%
metadata-eval29.4%
*-commutative29.4%
metadata-eval29.4%
div-inv35.8%
associate-*l*35.8%
div-inv29.4%
metadata-eval29.4%
Applied egg-rr29.4%
add-cube-cbrt38.9%
pow349.3%
associate-*l*49.3%
associate-*r*52.9%
metadata-eval52.9%
Applied egg-rr52.9%
if 1.00000000000000004e154 < (/.f64 angle #s(literal 180 binary64)) Initial program 14.5%
associate-*l*14.5%
*-commutative14.5%
associate-*l*14.5%
Simplified14.5%
unpow214.5%
unpow214.5%
difference-of-squares14.5%
Applied egg-rr14.5%
2-sin14.5%
div-inv15.1%
metadata-eval15.1%
*-commutative15.1%
metadata-eval15.1%
div-inv14.5%
associate-*l*14.5%
div-inv15.1%
metadata-eval15.1%
Applied egg-rr15.1%
add-exp-log29.9%
Applied egg-rr29.9%
Final simplification66.8%
angle\_m = (fabs.f64 angle)
angle\_s = (copysign.f64 #s(literal 1 binary64) angle)
(FPCore (angle_s a b angle_m)
:precision binary64
(let* ((t_0 (* (+ b a) (- b a))))
(*
angle_s
(if (<= (/ angle_m 180.0) 1e+42)
(* (+ b a) (* (- b a) (sin (/ (* (* angle_m PI) 360.0) 32400.0))))
(if (<= (/ angle_m 180.0) 1e+155)
(* t_0 (sin (pow (cbrt (* (* angle_m PI) 0.011111111111111112)) 3.0)))
(*
t_0
(sin (* PI (expm1 (log1p (* angle_m 0.011111111111111112)))))))))))angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
double t_0 = (b + a) * (b - a);
double tmp;
if ((angle_m / 180.0) <= 1e+42) {
tmp = (b + a) * ((b - a) * sin((((angle_m * ((double) M_PI)) * 360.0) / 32400.0)));
} else if ((angle_m / 180.0) <= 1e+155) {
tmp = t_0 * sin(pow(cbrt(((angle_m * ((double) M_PI)) * 0.011111111111111112)), 3.0));
} else {
tmp = t_0 * sin((((double) M_PI) * expm1(log1p((angle_m * 0.011111111111111112)))));
}
return angle_s * tmp;
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
double t_0 = (b + a) * (b - a);
double tmp;
if ((angle_m / 180.0) <= 1e+42) {
tmp = (b + a) * ((b - a) * Math.sin((((angle_m * Math.PI) * 360.0) / 32400.0)));
} else if ((angle_m / 180.0) <= 1e+155) {
tmp = t_0 * Math.sin(Math.pow(Math.cbrt(((angle_m * Math.PI) * 0.011111111111111112)), 3.0));
} else {
tmp = t_0 * Math.sin((Math.PI * Math.expm1(Math.log1p((angle_m * 0.011111111111111112)))));
}
return angle_s * tmp;
}
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) t_0 = Float64(Float64(b + a) * Float64(b - a)) tmp = 0.0 if (Float64(angle_m / 180.0) <= 1e+42) tmp = Float64(Float64(b + a) * Float64(Float64(b - a) * sin(Float64(Float64(Float64(angle_m * pi) * 360.0) / 32400.0)))); elseif (Float64(angle_m / 180.0) <= 1e+155) tmp = Float64(t_0 * sin((cbrt(Float64(Float64(angle_m * pi) * 0.011111111111111112)) ^ 3.0))); else tmp = Float64(t_0 * sin(Float64(pi * expm1(log1p(Float64(angle_m * 0.011111111111111112)))))); end return Float64(angle_s * tmp) end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := Block[{t$95$0 = N[(N[(b + a), $MachinePrecision] * N[(b - a), $MachinePrecision]), $MachinePrecision]}, N[(angle$95$s * If[LessEqual[N[(angle$95$m / 180.0), $MachinePrecision], 1e+42], N[(N[(b + a), $MachinePrecision] * N[(N[(b - a), $MachinePrecision] * N[Sin[N[(N[(N[(angle$95$m * Pi), $MachinePrecision] * 360.0), $MachinePrecision] / 32400.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(angle$95$m / 180.0), $MachinePrecision], 1e+155], N[(t$95$0 * N[Sin[N[Power[N[Power[N[(N[(angle$95$m * Pi), $MachinePrecision] * 0.011111111111111112), $MachinePrecision], 1/3], $MachinePrecision], 3.0], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(t$95$0 * N[Sin[N[(Pi * N[(Exp[N[Log[1 + N[(angle$95$m * 0.011111111111111112), $MachinePrecision]], $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
\begin{array}{l}
t_0 := \left(b + a\right) \cdot \left(b - a\right)\\
angle\_s \cdot \begin{array}{l}
\mathbf{if}\;\frac{angle\_m}{180} \leq 10^{+42}:\\
\;\;\;\;\left(b + a\right) \cdot \left(\left(b - a\right) \cdot \sin \left(\frac{\left(angle\_m \cdot \pi\right) \cdot 360}{32400}\right)\right)\\
\mathbf{elif}\;\frac{angle\_m}{180} \leq 10^{+155}:\\
\;\;\;\;t\_0 \cdot \sin \left({\left(\sqrt[3]{\left(angle\_m \cdot \pi\right) \cdot 0.011111111111111112}\right)}^{3}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot \sin \left(\pi \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(angle\_m \cdot 0.011111111111111112\right)\right)\right)\\
\end{array}
\end{array}
\end{array}
if (/.f64 angle #s(literal 180 binary64)) < 1.00000000000000004e42Initial program 56.8%
associate-*l*56.8%
*-commutative56.8%
associate-*l*56.8%
Simplified56.8%
associate-*r*56.8%
*-commutative56.8%
associate-*l*56.8%
add-cbrt-cube39.6%
pow1/325.2%
Applied egg-rr25.3%
unpow1/340.7%
rem-cbrt-cube57.7%
unpow257.7%
unpow257.7%
difference-of-squares64.4%
*-commutative64.4%
metadata-eval64.4%
div-inv64.0%
2-sin64.0%
associate-*l*73.9%
2-sin73.9%
div-inv74.2%
metadata-eval74.2%
Applied egg-rr74.2%
add-exp-log30.4%
Applied egg-rr30.4%
add-log-exp5.4%
associate-*r*5.4%
rem-exp-log11.4%
*-commutative11.4%
exp-lft-sqr11.4%
sum-log11.4%
add-log-exp28.1%
add-sqr-sqrt28.1%
unpow228.1%
*-commutative28.1%
metadata-eval28.1%
div-inv28.1%
add-log-exp73.2%
*-commutative73.2%
metadata-eval73.2%
div-inv73.0%
Applied egg-rr74.8%
*-commutative74.8%
*-commutative74.8%
*-commutative74.8%
distribute-rgt-out74.8%
metadata-eval74.8%
Simplified74.8%
if 1.00000000000000004e42 < (/.f64 angle #s(literal 180 binary64)) < 1.00000000000000001e155Initial program 32.2%
associate-*l*32.2%
*-commutative32.2%
associate-*l*32.2%
Simplified32.2%
unpow232.2%
unpow232.2%
difference-of-squares35.8%
Applied egg-rr35.8%
2-sin35.8%
div-inv29.4%
metadata-eval29.4%
*-commutative29.4%
metadata-eval29.4%
div-inv35.8%
associate-*l*35.8%
div-inv29.4%
metadata-eval29.4%
Applied egg-rr29.4%
add-cube-cbrt38.9%
pow349.3%
associate-*l*49.3%
associate-*r*52.9%
metadata-eval52.9%
Applied egg-rr52.9%
if 1.00000000000000001e155 < (/.f64 angle #s(literal 180 binary64)) Initial program 14.5%
associate-*l*14.5%
*-commutative14.5%
associate-*l*14.5%
Simplified14.5%
unpow214.5%
unpow214.5%
difference-of-squares14.5%
Applied egg-rr14.5%
2-sin14.5%
div-inv15.1%
metadata-eval15.1%
*-commutative15.1%
metadata-eval15.1%
div-inv14.5%
associate-*l*14.5%
div-inv15.1%
metadata-eval15.1%
Applied egg-rr15.1%
expm1-log1p-u25.9%
associate-*l*25.9%
metadata-eval25.9%
Applied egg-rr25.9%
Final simplification66.3%
angle\_m = (fabs.f64 angle)
angle\_s = (copysign.f64 #s(literal 1 binary64) angle)
(FPCore (angle_s a b angle_m)
:precision binary64
(*
angle_s
(if (<= b 0.00072)
(* (+ b a) (* (- b a) (sin (/ (* (* angle_m PI) 2.0) 180.0))))
(if (<= b 1.92e+174)
(* (+ b a) (* (- b a) (fabs (sin (* 2.0 (* (/ angle_m 180.0) PI))))))
(*
(+ b a)
(*
(- b a)
(*
angle_m
(+
(* -2.2862368541380886e-7 (* (pow angle_m 2.0) (pow PI 3.0)))
(* PI 0.011111111111111112)))))))))angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if (b <= 0.00072) {
tmp = (b + a) * ((b - a) * sin((((angle_m * ((double) M_PI)) * 2.0) / 180.0)));
} else if (b <= 1.92e+174) {
tmp = (b + a) * ((b - a) * fabs(sin((2.0 * ((angle_m / 180.0) * ((double) M_PI))))));
} else {
tmp = (b + a) * ((b - a) * (angle_m * ((-2.2862368541380886e-7 * (pow(angle_m, 2.0) * pow(((double) M_PI), 3.0))) + (((double) M_PI) * 0.011111111111111112))));
}
return angle_s * tmp;
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if (b <= 0.00072) {
tmp = (b + a) * ((b - a) * Math.sin((((angle_m * Math.PI) * 2.0) / 180.0)));
} else if (b <= 1.92e+174) {
tmp = (b + a) * ((b - a) * Math.abs(Math.sin((2.0 * ((angle_m / 180.0) * Math.PI)))));
} else {
tmp = (b + a) * ((b - a) * (angle_m * ((-2.2862368541380886e-7 * (Math.pow(angle_m, 2.0) * Math.pow(Math.PI, 3.0))) + (Math.PI * 0.011111111111111112))));
}
return angle_s * tmp;
}
angle\_m = math.fabs(angle) angle\_s = math.copysign(1.0, angle) def code(angle_s, a, b, angle_m): tmp = 0 if b <= 0.00072: tmp = (b + a) * ((b - a) * math.sin((((angle_m * math.pi) * 2.0) / 180.0))) elif b <= 1.92e+174: tmp = (b + a) * ((b - a) * math.fabs(math.sin((2.0 * ((angle_m / 180.0) * math.pi))))) else: tmp = (b + a) * ((b - a) * (angle_m * ((-2.2862368541380886e-7 * (math.pow(angle_m, 2.0) * math.pow(math.pi, 3.0))) + (math.pi * 0.011111111111111112)))) return angle_s * tmp
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) tmp = 0.0 if (b <= 0.00072) tmp = Float64(Float64(b + a) * Float64(Float64(b - a) * sin(Float64(Float64(Float64(angle_m * pi) * 2.0) / 180.0)))); elseif (b <= 1.92e+174) tmp = Float64(Float64(b + a) * Float64(Float64(b - a) * abs(sin(Float64(2.0 * Float64(Float64(angle_m / 180.0) * pi)))))); else tmp = Float64(Float64(b + a) * Float64(Float64(b - a) * Float64(angle_m * Float64(Float64(-2.2862368541380886e-7 * Float64((angle_m ^ 2.0) * (pi ^ 3.0))) + Float64(pi * 0.011111111111111112))))); end return Float64(angle_s * tmp) end
angle\_m = abs(angle); angle\_s = sign(angle) * abs(1.0); function tmp_2 = code(angle_s, a, b, angle_m) tmp = 0.0; if (b <= 0.00072) tmp = (b + a) * ((b - a) * sin((((angle_m * pi) * 2.0) / 180.0))); elseif (b <= 1.92e+174) tmp = (b + a) * ((b - a) * abs(sin((2.0 * ((angle_m / 180.0) * pi))))); else tmp = (b + a) * ((b - a) * (angle_m * ((-2.2862368541380886e-7 * ((angle_m ^ 2.0) * (pi ^ 3.0))) + (pi * 0.011111111111111112)))); end tmp_2 = angle_s * tmp; end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := N[(angle$95$s * If[LessEqual[b, 0.00072], N[(N[(b + a), $MachinePrecision] * N[(N[(b - a), $MachinePrecision] * N[Sin[N[(N[(N[(angle$95$m * Pi), $MachinePrecision] * 2.0), $MachinePrecision] / 180.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.92e+174], N[(N[(b + a), $MachinePrecision] * N[(N[(b - a), $MachinePrecision] * N[Abs[N[Sin[N[(2.0 * N[(N[(angle$95$m / 180.0), $MachinePrecision] * Pi), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(b + a), $MachinePrecision] * N[(N[(b - a), $MachinePrecision] * N[(angle$95$m * N[(N[(-2.2862368541380886e-7 * N[(N[Power[angle$95$m, 2.0], $MachinePrecision] * N[Power[Pi, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
angle\_s \cdot \begin{array}{l}
\mathbf{if}\;b \leq 0.00072:\\
\;\;\;\;\left(b + a\right) \cdot \left(\left(b - a\right) \cdot \sin \left(\frac{\left(angle\_m \cdot \pi\right) \cdot 2}{180}\right)\right)\\
\mathbf{elif}\;b \leq 1.92 \cdot 10^{+174}:\\
\;\;\;\;\left(b + a\right) \cdot \left(\left(b - a\right) \cdot \left|\sin \left(2 \cdot \left(\frac{angle\_m}{180} \cdot \pi\right)\right)\right|\right)\\
\mathbf{else}:\\
\;\;\;\;\left(b + a\right) \cdot \left(\left(b - a\right) \cdot \left(angle\_m \cdot \left(-2.2862368541380886 \cdot 10^{-7} \cdot \left({angle\_m}^{2} \cdot {\pi}^{3}\right) + \pi \cdot 0.011111111111111112\right)\right)\right)\\
\end{array}
\end{array}
if b < 7.20000000000000045e-4Initial program 52.1%
associate-*l*52.1%
*-commutative52.1%
associate-*l*52.1%
Simplified52.1%
associate-*r*52.1%
*-commutative52.1%
associate-*l*52.1%
add-cbrt-cube37.3%
pow1/326.4%
Applied egg-rr26.5%
unpow1/338.0%
rem-cbrt-cube52.7%
unpow252.7%
unpow252.7%
difference-of-squares58.4%
*-commutative58.4%
metadata-eval58.4%
div-inv58.4%
2-sin58.4%
associate-*l*61.7%
2-sin61.7%
div-inv61.7%
metadata-eval61.7%
Applied egg-rr61.7%
associate-*r*61.7%
metadata-eval61.7%
div-inv61.7%
associate-*r/62.6%
associate-*l/62.6%
Applied egg-rr62.6%
if 7.20000000000000045e-4 < b < 1.9199999999999999e174Initial program 45.7%
associate-*l*45.7%
*-commutative45.7%
associate-*l*45.7%
Simplified45.7%
associate-*r*45.7%
*-commutative45.7%
associate-*l*45.7%
add-cbrt-cube30.9%
pow1/318.4%
Applied egg-rr18.4%
unpow1/331.0%
rem-cbrt-cube45.7%
unpow245.7%
unpow245.7%
difference-of-squares48.0%
*-commutative48.0%
metadata-eval48.0%
div-inv48.0%
2-sin48.0%
associate-*l*60.6%
2-sin60.6%
div-inv60.7%
metadata-eval60.7%
Applied egg-rr60.7%
add-exp-log31.1%
Applied egg-rr31.1%
add-sqr-sqrt26.9%
sqrt-unprod16.0%
pow216.0%
associate-*r*16.0%
rem-exp-log19.6%
metadata-eval19.6%
div-inv19.5%
associate-*r/19.0%
Applied egg-rr19.0%
unpow219.0%
rem-sqrt-square34.8%
*-commutative34.8%
associate-/l*35.4%
Simplified35.4%
if 1.9199999999999999e174 < b Initial program 25.2%
associate-*l*25.2%
*-commutative25.2%
associate-*l*25.2%
Simplified25.2%
associate-*r*25.2%
*-commutative25.2%
associate-*l*25.2%
add-cbrt-cube25.2%
pow1/320.0%
Applied egg-rr20.0%
unpow1/320.5%
rem-cbrt-cube20.5%
unpow220.5%
unpow220.5%
difference-of-squares30.9%
*-commutative30.9%
metadata-eval30.9%
div-inv35.7%
2-sin35.7%
associate-*l*71.1%
2-sin71.1%
div-inv66.3%
metadata-eval66.3%
Applied egg-rr66.3%
Taylor expanded in angle around 0 90.3%
Final simplification60.3%
angle\_m = (fabs.f64 angle)
angle\_s = (copysign.f64 #s(literal 1 binary64) angle)
(FPCore (angle_s a b angle_m)
:precision binary64
(*
angle_s
(if (<= (pow b 2.0) 1e-226)
(* (+ b a) (* (sin (* (* angle_m PI) 0.011111111111111112)) (- a)))
(* (+ b a) (* 0.011111111111111112 (* angle_m (* (- b a) PI)))))))angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if (pow(b, 2.0) <= 1e-226) {
tmp = (b + a) * (sin(((angle_m * ((double) M_PI)) * 0.011111111111111112)) * -a);
} else {
tmp = (b + a) * (0.011111111111111112 * (angle_m * ((b - a) * ((double) M_PI))));
}
return angle_s * tmp;
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if (Math.pow(b, 2.0) <= 1e-226) {
tmp = (b + a) * (Math.sin(((angle_m * Math.PI) * 0.011111111111111112)) * -a);
} else {
tmp = (b + a) * (0.011111111111111112 * (angle_m * ((b - a) * Math.PI)));
}
return angle_s * tmp;
}
angle\_m = math.fabs(angle) angle\_s = math.copysign(1.0, angle) def code(angle_s, a, b, angle_m): tmp = 0 if math.pow(b, 2.0) <= 1e-226: tmp = (b + a) * (math.sin(((angle_m * math.pi) * 0.011111111111111112)) * -a) else: tmp = (b + a) * (0.011111111111111112 * (angle_m * ((b - a) * math.pi))) return angle_s * tmp
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) tmp = 0.0 if ((b ^ 2.0) <= 1e-226) tmp = Float64(Float64(b + a) * Float64(sin(Float64(Float64(angle_m * pi) * 0.011111111111111112)) * Float64(-a))); else tmp = Float64(Float64(b + a) * Float64(0.011111111111111112 * Float64(angle_m * Float64(Float64(b - a) * pi)))); end return Float64(angle_s * tmp) end
angle\_m = abs(angle); angle\_s = sign(angle) * abs(1.0); function tmp_2 = code(angle_s, a, b, angle_m) tmp = 0.0; if ((b ^ 2.0) <= 1e-226) tmp = (b + a) * (sin(((angle_m * pi) * 0.011111111111111112)) * -a); else tmp = (b + a) * (0.011111111111111112 * (angle_m * ((b - a) * pi))); end tmp_2 = angle_s * tmp; end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := N[(angle$95$s * If[LessEqual[N[Power[b, 2.0], $MachinePrecision], 1e-226], N[(N[(b + a), $MachinePrecision] * N[(N[Sin[N[(N[(angle$95$m * Pi), $MachinePrecision] * 0.011111111111111112), $MachinePrecision]], $MachinePrecision] * (-a)), $MachinePrecision]), $MachinePrecision], N[(N[(b + a), $MachinePrecision] * N[(0.011111111111111112 * N[(angle$95$m * N[(N[(b - a), $MachinePrecision] * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
angle\_s \cdot \begin{array}{l}
\mathbf{if}\;{b}^{2} \leq 10^{-226}:\\
\;\;\;\;\left(b + a\right) \cdot \left(\sin \left(\left(angle\_m \cdot \pi\right) \cdot 0.011111111111111112\right) \cdot \left(-a\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(b + a\right) \cdot \left(0.011111111111111112 \cdot \left(angle\_m \cdot \left(\left(b - a\right) \cdot \pi\right)\right)\right)\\
\end{array}
\end{array}
if (pow.f64 b #s(literal 2 binary64)) < 9.99999999999999921e-227Initial program 57.4%
associate-*l*57.4%
*-commutative57.4%
associate-*l*57.4%
Simplified57.4%
associate-*r*57.4%
*-commutative57.4%
associate-*l*57.4%
add-cbrt-cube44.6%
pow1/335.8%
Applied egg-rr36.0%
unpow1/344.8%
rem-cbrt-cube57.4%
unpow257.4%
unpow257.4%
difference-of-squares57.4%
*-commutative57.4%
metadata-eval57.4%
div-inv57.4%
2-sin57.4%
associate-*l*61.3%
2-sin61.3%
div-inv61.4%
metadata-eval61.4%
Applied egg-rr61.4%
add-exp-log35.1%
Applied egg-rr35.1%
Taylor expanded in b around 0 62.3%
mul-1-neg62.3%
*-commutative62.3%
distribute-rgt-neg-in62.3%
Simplified62.3%
if 9.99999999999999921e-227 < (pow.f64 b #s(literal 2 binary64)) Initial program 44.0%
associate-*l*44.0%
*-commutative44.0%
associate-*l*44.0%
Simplified44.0%
associate-*r*44.0%
*-commutative44.0%
associate-*l*44.0%
add-cbrt-cube30.0%
pow1/318.2%
Applied egg-rr18.2%
unpow1/330.1%
rem-cbrt-cube44.1%
unpow244.1%
unpow244.1%
difference-of-squares52.7%
*-commutative52.7%
metadata-eval52.7%
div-inv53.3%
2-sin53.3%
associate-*l*62.8%
2-sin62.8%
div-inv62.2%
metadata-eval62.2%
Applied egg-rr62.2%
Taylor expanded in angle around 0 61.6%
Final simplification61.9%
angle\_m = (fabs.f64 angle)
angle\_s = (copysign.f64 #s(literal 1 binary64) angle)
(FPCore (angle_s a b angle_m)
:precision binary64
(*
angle_s
(if (<= (/ angle_m 180.0) 1e+79)
(* (+ b a) (* (- b a) (sin (* PI (/ (* angle_m 2.0) 180.0)))))
(* (* (+ b a) (- b a)) (sin (* PI (/ 2.0 (/ 180.0 angle_m))))))))angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if ((angle_m / 180.0) <= 1e+79) {
tmp = (b + a) * ((b - a) * sin((((double) M_PI) * ((angle_m * 2.0) / 180.0))));
} else {
tmp = ((b + a) * (b - a)) * sin((((double) M_PI) * (2.0 / (180.0 / angle_m))));
}
return angle_s * tmp;
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if ((angle_m / 180.0) <= 1e+79) {
tmp = (b + a) * ((b - a) * Math.sin((Math.PI * ((angle_m * 2.0) / 180.0))));
} else {
tmp = ((b + a) * (b - a)) * Math.sin((Math.PI * (2.0 / (180.0 / angle_m))));
}
return angle_s * tmp;
}
angle\_m = math.fabs(angle) angle\_s = math.copysign(1.0, angle) def code(angle_s, a, b, angle_m): tmp = 0 if (angle_m / 180.0) <= 1e+79: tmp = (b + a) * ((b - a) * math.sin((math.pi * ((angle_m * 2.0) / 180.0)))) else: tmp = ((b + a) * (b - a)) * math.sin((math.pi * (2.0 / (180.0 / angle_m)))) return angle_s * tmp
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) tmp = 0.0 if (Float64(angle_m / 180.0) <= 1e+79) tmp = Float64(Float64(b + a) * Float64(Float64(b - a) * sin(Float64(pi * Float64(Float64(angle_m * 2.0) / 180.0))))); else tmp = Float64(Float64(Float64(b + a) * Float64(b - a)) * sin(Float64(pi * Float64(2.0 / Float64(180.0 / angle_m))))); end return Float64(angle_s * tmp) end
angle\_m = abs(angle); angle\_s = sign(angle) * abs(1.0); function tmp_2 = code(angle_s, a, b, angle_m) tmp = 0.0; if ((angle_m / 180.0) <= 1e+79) tmp = (b + a) * ((b - a) * sin((pi * ((angle_m * 2.0) / 180.0)))); else tmp = ((b + a) * (b - a)) * sin((pi * (2.0 / (180.0 / angle_m)))); end tmp_2 = angle_s * tmp; end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := N[(angle$95$s * If[LessEqual[N[(angle$95$m / 180.0), $MachinePrecision], 1e+79], N[(N[(b + a), $MachinePrecision] * N[(N[(b - a), $MachinePrecision] * N[Sin[N[(Pi * N[(N[(angle$95$m * 2.0), $MachinePrecision] / 180.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(b + a), $MachinePrecision] * N[(b - a), $MachinePrecision]), $MachinePrecision] * N[Sin[N[(Pi * N[(2.0 / N[(180.0 / angle$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
angle\_s \cdot \begin{array}{l}
\mathbf{if}\;\frac{angle\_m}{180} \leq 10^{+79}:\\
\;\;\;\;\left(b + a\right) \cdot \left(\left(b - a\right) \cdot \sin \left(\pi \cdot \frac{angle\_m \cdot 2}{180}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(b + a\right) \cdot \left(b - a\right)\right) \cdot \sin \left(\pi \cdot \frac{2}{\frac{180}{angle\_m}}\right)\\
\end{array}
\end{array}
if (/.f64 angle #s(literal 180 binary64)) < 9.99999999999999967e78Initial program 55.2%
associate-*l*55.2%
*-commutative55.2%
associate-*l*55.2%
Simplified55.2%
associate-*r*55.2%
*-commutative55.2%
associate-*l*55.2%
add-cbrt-cube38.8%
pow1/325.0%
Applied egg-rr25.1%
unpow1/339.4%
rem-cbrt-cube55.6%
unpow255.6%
unpow255.6%
difference-of-squares62.5%
*-commutative62.5%
metadata-eval62.5%
div-inv62.6%
2-sin62.6%
associate-*l*72.0%
2-sin72.0%
div-inv71.8%
metadata-eval71.8%
Applied egg-rr71.8%
metadata-eval71.8%
div-inv72.0%
associate-*l/72.0%
Applied egg-rr72.0%
if 9.99999999999999967e78 < (/.f64 angle #s(literal 180 binary64)) Initial program 22.5%
associate-*l*22.5%
*-commutative22.5%
associate-*l*22.5%
Simplified22.5%
unpow222.5%
unpow222.5%
difference-of-squares22.5%
Applied egg-rr22.5%
2-sin22.5%
div-inv21.3%
metadata-eval21.3%
*-commutative21.3%
metadata-eval21.3%
div-inv22.5%
associate-*l*22.5%
div-inv21.3%
metadata-eval21.3%
Applied egg-rr21.3%
*-commutative21.3%
metadata-eval21.3%
div-inv22.5%
clear-num28.8%
un-div-inv28.8%
Applied egg-rr28.8%
angle\_m = (fabs.f64 angle)
angle\_s = (copysign.f64 #s(literal 1 binary64) angle)
(FPCore (angle_s a b angle_m)
:precision binary64
(*
angle_s
(if (<= (/ angle_m 180.0) 5e+39)
(* (+ b a) (* (- b a) (sin (* PI (* angle_m 0.011111111111111112)))))
(* (* (+ b a) (- b a)) (sin (* PI (/ 2.0 (/ 180.0 angle_m))))))))angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if ((angle_m / 180.0) <= 5e+39) {
tmp = (b + a) * ((b - a) * sin((((double) M_PI) * (angle_m * 0.011111111111111112))));
} else {
tmp = ((b + a) * (b - a)) * sin((((double) M_PI) * (2.0 / (180.0 / angle_m))));
}
return angle_s * tmp;
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if ((angle_m / 180.0) <= 5e+39) {
tmp = (b + a) * ((b - a) * Math.sin((Math.PI * (angle_m * 0.011111111111111112))));
} else {
tmp = ((b + a) * (b - a)) * Math.sin((Math.PI * (2.0 / (180.0 / angle_m))));
}
return angle_s * tmp;
}
angle\_m = math.fabs(angle) angle\_s = math.copysign(1.0, angle) def code(angle_s, a, b, angle_m): tmp = 0 if (angle_m / 180.0) <= 5e+39: tmp = (b + a) * ((b - a) * math.sin((math.pi * (angle_m * 0.011111111111111112)))) else: tmp = ((b + a) * (b - a)) * math.sin((math.pi * (2.0 / (180.0 / angle_m)))) return angle_s * tmp
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) tmp = 0.0 if (Float64(angle_m / 180.0) <= 5e+39) tmp = Float64(Float64(b + a) * Float64(Float64(b - a) * sin(Float64(pi * Float64(angle_m * 0.011111111111111112))))); else tmp = Float64(Float64(Float64(b + a) * Float64(b - a)) * sin(Float64(pi * Float64(2.0 / Float64(180.0 / angle_m))))); end return Float64(angle_s * tmp) end
angle\_m = abs(angle); angle\_s = sign(angle) * abs(1.0); function tmp_2 = code(angle_s, a, b, angle_m) tmp = 0.0; if ((angle_m / 180.0) <= 5e+39) tmp = (b + a) * ((b - a) * sin((pi * (angle_m * 0.011111111111111112)))); else tmp = ((b + a) * (b - a)) * sin((pi * (2.0 / (180.0 / angle_m)))); end tmp_2 = angle_s * tmp; end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := N[(angle$95$s * If[LessEqual[N[(angle$95$m / 180.0), $MachinePrecision], 5e+39], N[(N[(b + a), $MachinePrecision] * N[(N[(b - a), $MachinePrecision] * N[Sin[N[(Pi * N[(angle$95$m * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(b + a), $MachinePrecision] * N[(b - a), $MachinePrecision]), $MachinePrecision] * N[Sin[N[(Pi * N[(2.0 / N[(180.0 / angle$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
angle\_s \cdot \begin{array}{l}
\mathbf{if}\;\frac{angle\_m}{180} \leq 5 \cdot 10^{+39}:\\
\;\;\;\;\left(b + a\right) \cdot \left(\left(b - a\right) \cdot \sin \left(\pi \cdot \left(angle\_m \cdot 0.011111111111111112\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(b + a\right) \cdot \left(b - a\right)\right) \cdot \sin \left(\pi \cdot \frac{2}{\frac{180}{angle\_m}}\right)\\
\end{array}
\end{array}
if (/.f64 angle #s(literal 180 binary64)) < 5.00000000000000015e39Initial program 57.3%
associate-*l*57.3%
*-commutative57.3%
associate-*l*57.3%
Simplified57.3%
associate-*r*57.3%
*-commutative57.3%
associate-*l*57.3%
add-cbrt-cube39.9%
pow1/325.5%
Applied egg-rr25.5%
unpow1/341.1%
rem-cbrt-cube58.2%
unpow258.2%
unpow258.2%
difference-of-squares65.0%
*-commutative65.0%
metadata-eval65.0%
div-inv64.6%
2-sin64.6%
associate-*l*74.5%
2-sin74.5%
div-inv74.9%
metadata-eval74.9%
Applied egg-rr74.9%
Taylor expanded in angle around 0 74.9%
if 5.00000000000000015e39 < (/.f64 angle #s(literal 180 binary64)) Initial program 22.4%
associate-*l*22.4%
*-commutative22.4%
associate-*l*22.4%
Simplified22.4%
unpow222.4%
unpow222.4%
difference-of-squares24.0%
Applied egg-rr24.0%
2-sin24.0%
div-inv21.4%
metadata-eval21.4%
*-commutative21.4%
metadata-eval21.4%
div-inv24.0%
associate-*l*24.0%
div-inv21.4%
metadata-eval21.4%
Applied egg-rr21.4%
*-commutative21.4%
metadata-eval21.4%
div-inv24.0%
clear-num28.5%
un-div-inv28.5%
Applied egg-rr28.5%
Final simplification63.7%
angle\_m = (fabs.f64 angle)
angle\_s = (copysign.f64 #s(literal 1 binary64) angle)
(FPCore (angle_s a b angle_m)
:precision binary64
(*
angle_s
(if (<= (/ angle_m 180.0) 5.5e+31)
(* (+ b a) (* (- b a) (sin (* PI (* angle_m 0.011111111111111112)))))
(* (* (+ b a) (- b a)) (sin (* (* angle_m PI) 0.011111111111111112))))))angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if ((angle_m / 180.0) <= 5.5e+31) {
tmp = (b + a) * ((b - a) * sin((((double) M_PI) * (angle_m * 0.011111111111111112))));
} else {
tmp = ((b + a) * (b - a)) * sin(((angle_m * ((double) M_PI)) * 0.011111111111111112));
}
return angle_s * tmp;
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if ((angle_m / 180.0) <= 5.5e+31) {
tmp = (b + a) * ((b - a) * Math.sin((Math.PI * (angle_m * 0.011111111111111112))));
} else {
tmp = ((b + a) * (b - a)) * Math.sin(((angle_m * Math.PI) * 0.011111111111111112));
}
return angle_s * tmp;
}
angle\_m = math.fabs(angle) angle\_s = math.copysign(1.0, angle) def code(angle_s, a, b, angle_m): tmp = 0 if (angle_m / 180.0) <= 5.5e+31: tmp = (b + a) * ((b - a) * math.sin((math.pi * (angle_m * 0.011111111111111112)))) else: tmp = ((b + a) * (b - a)) * math.sin(((angle_m * math.pi) * 0.011111111111111112)) return angle_s * tmp
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) tmp = 0.0 if (Float64(angle_m / 180.0) <= 5.5e+31) tmp = Float64(Float64(b + a) * Float64(Float64(b - a) * sin(Float64(pi * Float64(angle_m * 0.011111111111111112))))); else tmp = Float64(Float64(Float64(b + a) * Float64(b - a)) * sin(Float64(Float64(angle_m * pi) * 0.011111111111111112))); end return Float64(angle_s * tmp) end
angle\_m = abs(angle); angle\_s = sign(angle) * abs(1.0); function tmp_2 = code(angle_s, a, b, angle_m) tmp = 0.0; if ((angle_m / 180.0) <= 5.5e+31) tmp = (b + a) * ((b - a) * sin((pi * (angle_m * 0.011111111111111112)))); else tmp = ((b + a) * (b - a)) * sin(((angle_m * pi) * 0.011111111111111112)); end tmp_2 = angle_s * tmp; end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := N[(angle$95$s * If[LessEqual[N[(angle$95$m / 180.0), $MachinePrecision], 5.5e+31], N[(N[(b + a), $MachinePrecision] * N[(N[(b - a), $MachinePrecision] * N[Sin[N[(Pi * N[(angle$95$m * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(b + a), $MachinePrecision] * N[(b - a), $MachinePrecision]), $MachinePrecision] * N[Sin[N[(N[(angle$95$m * Pi), $MachinePrecision] * 0.011111111111111112), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
angle\_s \cdot \begin{array}{l}
\mathbf{if}\;\frac{angle\_m}{180} \leq 5.5 \cdot 10^{+31}:\\
\;\;\;\;\left(b + a\right) \cdot \left(\left(b - a\right) \cdot \sin \left(\pi \cdot \left(angle\_m \cdot 0.011111111111111112\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(b + a\right) \cdot \left(b - a\right)\right) \cdot \sin \left(\left(angle\_m \cdot \pi\right) \cdot 0.011111111111111112\right)\\
\end{array}
\end{array}
if (/.f64 angle #s(literal 180 binary64)) < 5.50000000000000002e31Initial program 58.2%
associate-*l*58.2%
*-commutative58.2%
associate-*l*58.2%
Simplified58.2%
associate-*r*58.2%
*-commutative58.2%
associate-*l*58.2%
add-cbrt-cube40.6%
pow1/325.8%
Applied egg-rr25.9%
unpow1/341.2%
rem-cbrt-cube58.6%
unpow258.6%
unpow258.6%
difference-of-squares65.5%
*-commutative65.5%
metadata-eval65.5%
div-inv65.6%
2-sin65.6%
associate-*l*75.7%
2-sin75.7%
div-inv75.5%
metadata-eval75.5%
Applied egg-rr75.5%
Taylor expanded in angle around 0 75.6%
if 5.50000000000000002e31 < (/.f64 angle #s(literal 180 binary64)) Initial program 21.3%
associate-*l*21.3%
*-commutative21.3%
associate-*l*21.3%
Simplified21.3%
unpow221.3%
unpow221.3%
difference-of-squares22.9%
Applied egg-rr22.9%
2-sin22.9%
div-inv21.9%
metadata-eval21.9%
*-commutative21.9%
metadata-eval21.9%
div-inv22.9%
associate-*l*22.9%
div-inv21.9%
metadata-eval21.9%
Applied egg-rr21.9%
Taylor expanded in angle around inf 30.2%
Final simplification64.0%
angle\_m = (fabs.f64 angle)
angle\_s = (copysign.f64 #s(literal 1 binary64) angle)
(FPCore (angle_s a b angle_m)
:precision binary64
(*
angle_s
(if (<= (/ angle_m 180.0) 1e-11)
(* (+ b a) (* 0.011111111111111112 (* angle_m (* (- b a) PI))))
(* (* (+ b a) (- b a)) (sin (* (* angle_m PI) 0.011111111111111112))))))angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if ((angle_m / 180.0) <= 1e-11) {
tmp = (b + a) * (0.011111111111111112 * (angle_m * ((b - a) * ((double) M_PI))));
} else {
tmp = ((b + a) * (b - a)) * sin(((angle_m * ((double) M_PI)) * 0.011111111111111112));
}
return angle_s * tmp;
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if ((angle_m / 180.0) <= 1e-11) {
tmp = (b + a) * (0.011111111111111112 * (angle_m * ((b - a) * Math.PI)));
} else {
tmp = ((b + a) * (b - a)) * Math.sin(((angle_m * Math.PI) * 0.011111111111111112));
}
return angle_s * tmp;
}
angle\_m = math.fabs(angle) angle\_s = math.copysign(1.0, angle) def code(angle_s, a, b, angle_m): tmp = 0 if (angle_m / 180.0) <= 1e-11: tmp = (b + a) * (0.011111111111111112 * (angle_m * ((b - a) * math.pi))) else: tmp = ((b + a) * (b - a)) * math.sin(((angle_m * math.pi) * 0.011111111111111112)) return angle_s * tmp
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) tmp = 0.0 if (Float64(angle_m / 180.0) <= 1e-11) tmp = Float64(Float64(b + a) * Float64(0.011111111111111112 * Float64(angle_m * Float64(Float64(b - a) * pi)))); else tmp = Float64(Float64(Float64(b + a) * Float64(b - a)) * sin(Float64(Float64(angle_m * pi) * 0.011111111111111112))); end return Float64(angle_s * tmp) end
angle\_m = abs(angle); angle\_s = sign(angle) * abs(1.0); function tmp_2 = code(angle_s, a, b, angle_m) tmp = 0.0; if ((angle_m / 180.0) <= 1e-11) tmp = (b + a) * (0.011111111111111112 * (angle_m * ((b - a) * pi))); else tmp = ((b + a) * (b - a)) * sin(((angle_m * pi) * 0.011111111111111112)); end tmp_2 = angle_s * tmp; end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := N[(angle$95$s * If[LessEqual[N[(angle$95$m / 180.0), $MachinePrecision], 1e-11], N[(N[(b + a), $MachinePrecision] * N[(0.011111111111111112 * N[(angle$95$m * N[(N[(b - a), $MachinePrecision] * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(b + a), $MachinePrecision] * N[(b - a), $MachinePrecision]), $MachinePrecision] * N[Sin[N[(N[(angle$95$m * Pi), $MachinePrecision] * 0.011111111111111112), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
angle\_s \cdot \begin{array}{l}
\mathbf{if}\;\frac{angle\_m}{180} \leq 10^{-11}:\\
\;\;\;\;\left(b + a\right) \cdot \left(0.011111111111111112 \cdot \left(angle\_m \cdot \left(\left(b - a\right) \cdot \pi\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(b + a\right) \cdot \left(b - a\right)\right) \cdot \sin \left(\left(angle\_m \cdot \pi\right) \cdot 0.011111111111111112\right)\\
\end{array}
\end{array}
if (/.f64 angle #s(literal 180 binary64)) < 9.99999999999999939e-12Initial program 58.8%
associate-*l*58.8%
*-commutative58.8%
associate-*l*58.8%
Simplified58.8%
associate-*r*58.8%
*-commutative58.8%
associate-*l*58.8%
add-cbrt-cube40.5%
pow1/326.4%
Applied egg-rr26.4%
unpow1/341.2%
rem-cbrt-cube59.2%
unpow259.2%
unpow259.2%
difference-of-squares65.8%
*-commutative65.8%
metadata-eval65.8%
div-inv66.0%
2-sin66.0%
associate-*l*76.5%
2-sin76.5%
div-inv76.3%
metadata-eval76.3%
Applied egg-rr76.3%
Taylor expanded in angle around 0 72.6%
if 9.99999999999999939e-12 < (/.f64 angle #s(literal 180 binary64)) Initial program 23.9%
associate-*l*23.9%
*-commutative23.9%
associate-*l*23.9%
Simplified23.9%
unpow223.9%
unpow223.9%
difference-of-squares26.6%
Applied egg-rr26.6%
2-sin26.6%
div-inv25.8%
metadata-eval25.8%
*-commutative25.8%
metadata-eval25.8%
div-inv26.6%
associate-*l*26.6%
div-inv25.8%
metadata-eval25.8%
Applied egg-rr25.8%
Taylor expanded in angle around inf 31.8%
Final simplification61.0%
angle\_m = (fabs.f64 angle) angle\_s = (copysign.f64 #s(literal 1 binary64) angle) (FPCore (angle_s a b angle_m) :precision binary64 (* angle_s (* (+ b a) (* (- b a) (sin (/ (* (* angle_m PI) 360.0) 32400.0))))))
angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
return angle_s * ((b + a) * ((b - a) * sin((((angle_m * ((double) M_PI)) * 360.0) / 32400.0))));
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
return angle_s * ((b + a) * ((b - a) * Math.sin((((angle_m * Math.PI) * 360.0) / 32400.0))));
}
angle\_m = math.fabs(angle) angle\_s = math.copysign(1.0, angle) def code(angle_s, a, b, angle_m): return angle_s * ((b + a) * ((b - a) * math.sin((((angle_m * math.pi) * 360.0) / 32400.0))))
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) return Float64(angle_s * Float64(Float64(b + a) * Float64(Float64(b - a) * sin(Float64(Float64(Float64(angle_m * pi) * 360.0) / 32400.0))))) end
angle\_m = abs(angle); angle\_s = sign(angle) * abs(1.0); function tmp = code(angle_s, a, b, angle_m) tmp = angle_s * ((b + a) * ((b - a) * sin((((angle_m * pi) * 360.0) / 32400.0)))); end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := N[(angle$95$s * N[(N[(b + a), $MachinePrecision] * N[(N[(b - a), $MachinePrecision] * N[Sin[N[(N[(N[(angle$95$m * Pi), $MachinePrecision] * 360.0), $MachinePrecision] / 32400.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
angle\_s \cdot \left(\left(b + a\right) \cdot \left(\left(b - a\right) \cdot \sin \left(\frac{\left(angle\_m \cdot \pi\right) \cdot 360}{32400}\right)\right)\right)
\end{array}
Initial program 48.8%
associate-*l*48.8%
*-commutative48.8%
associate-*l*48.8%
Simplified48.8%
associate-*r*48.8%
*-commutative48.8%
associate-*l*48.8%
add-cbrt-cube35.2%
pow1/324.6%
Applied egg-rr24.6%
unpow1/335.4%
rem-cbrt-cube48.9%
unpow248.9%
unpow248.9%
difference-of-squares54.4%
*-commutative54.4%
metadata-eval54.4%
div-inv54.8%
2-sin54.8%
associate-*l*62.3%
2-sin62.3%
div-inv61.9%
metadata-eval61.9%
Applied egg-rr61.9%
add-exp-log32.2%
Applied egg-rr32.2%
add-log-exp4.1%
associate-*r*4.1%
rem-exp-log8.7%
*-commutative8.7%
exp-lft-sqr8.7%
sum-log8.7%
add-log-exp21.5%
add-sqr-sqrt21.5%
unpow221.5%
*-commutative21.5%
metadata-eval21.5%
div-inv21.5%
add-log-exp62.3%
*-commutative62.3%
metadata-eval62.3%
div-inv61.3%
Applied egg-rr64.4%
*-commutative64.4%
*-commutative64.4%
*-commutative64.4%
distribute-rgt-out64.4%
metadata-eval64.4%
Simplified64.4%
angle\_m = (fabs.f64 angle)
angle\_s = (copysign.f64 #s(literal 1 binary64) angle)
(FPCore (angle_s a b angle_m)
:precision binary64
(*
angle_s
(if (<= angle_m 2.7e+58)
(* (+ b a) (* 0.011111111111111112 (* angle_m (* (- b a) PI))))
(* 0.011111111111111112 (* angle_m (* PI (* a (- (- a) b))))))))angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if (angle_m <= 2.7e+58) {
tmp = (b + a) * (0.011111111111111112 * (angle_m * ((b - a) * ((double) M_PI))));
} else {
tmp = 0.011111111111111112 * (angle_m * (((double) M_PI) * (a * (-a - b))));
}
return angle_s * tmp;
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if (angle_m <= 2.7e+58) {
tmp = (b + a) * (0.011111111111111112 * (angle_m * ((b - a) * Math.PI)));
} else {
tmp = 0.011111111111111112 * (angle_m * (Math.PI * (a * (-a - b))));
}
return angle_s * tmp;
}
angle\_m = math.fabs(angle) angle\_s = math.copysign(1.0, angle) def code(angle_s, a, b, angle_m): tmp = 0 if angle_m <= 2.7e+58: tmp = (b + a) * (0.011111111111111112 * (angle_m * ((b - a) * math.pi))) else: tmp = 0.011111111111111112 * (angle_m * (math.pi * (a * (-a - b)))) return angle_s * tmp
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) tmp = 0.0 if (angle_m <= 2.7e+58) tmp = Float64(Float64(b + a) * Float64(0.011111111111111112 * Float64(angle_m * Float64(Float64(b - a) * pi)))); else tmp = Float64(0.011111111111111112 * Float64(angle_m * Float64(pi * Float64(a * Float64(Float64(-a) - b))))); end return Float64(angle_s * tmp) end
angle\_m = abs(angle); angle\_s = sign(angle) * abs(1.0); function tmp_2 = code(angle_s, a, b, angle_m) tmp = 0.0; if (angle_m <= 2.7e+58) tmp = (b + a) * (0.011111111111111112 * (angle_m * ((b - a) * pi))); else tmp = 0.011111111111111112 * (angle_m * (pi * (a * (-a - b)))); end tmp_2 = angle_s * tmp; end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := N[(angle$95$s * If[LessEqual[angle$95$m, 2.7e+58], N[(N[(b + a), $MachinePrecision] * N[(0.011111111111111112 * N[(angle$95$m * N[(N[(b - a), $MachinePrecision] * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.011111111111111112 * N[(angle$95$m * N[(Pi * N[(a * N[((-a) - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
angle\_s \cdot \begin{array}{l}
\mathbf{if}\;angle\_m \leq 2.7 \cdot 10^{+58}:\\
\;\;\;\;\left(b + a\right) \cdot \left(0.011111111111111112 \cdot \left(angle\_m \cdot \left(\left(b - a\right) \cdot \pi\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;0.011111111111111112 \cdot \left(angle\_m \cdot \left(\pi \cdot \left(a \cdot \left(\left(-a\right) - b\right)\right)\right)\right)\\
\end{array}
\end{array}
if angle < 2.7000000000000001e58Initial program 55.7%
associate-*l*55.7%
*-commutative55.7%
associate-*l*55.7%
Simplified55.7%
associate-*r*55.7%
*-commutative55.7%
associate-*l*55.7%
add-cbrt-cube38.8%
pow1/325.2%
Applied egg-rr25.3%
unpow1/339.9%
rem-cbrt-cube56.5%
unpow256.5%
unpow256.5%
difference-of-squares63.1%
*-commutative63.1%
metadata-eval63.1%
div-inv62.8%
2-sin62.8%
associate-*l*72.4%
2-sin72.4%
div-inv72.7%
metadata-eval72.7%
Applied egg-rr72.7%
Taylor expanded in angle around 0 69.1%
if 2.7000000000000001e58 < angle Initial program 24.3%
associate-*l*24.3%
*-commutative24.3%
associate-*l*24.3%
Simplified24.3%
unpow224.3%
unpow224.3%
difference-of-squares26.1%
Applied egg-rr26.1%
2-sin26.1%
div-inv23.3%
metadata-eval23.3%
*-commutative23.3%
metadata-eval23.3%
div-inv26.1%
associate-*l*26.1%
div-inv23.3%
metadata-eval23.3%
Applied egg-rr23.3%
Taylor expanded in angle around 0 24.9%
Taylor expanded in b around 0 23.4%
neg-mul-123.4%
Simplified23.4%
Final simplification59.1%
angle\_m = (fabs.f64 angle)
angle\_s = (copysign.f64 #s(literal 1 binary64) angle)
(FPCore (angle_s a b angle_m)
:precision binary64
(*
angle_s
(if (<= angle_m 8.8e+173)
(* (* a 0.011111111111111112) (* (- b a) (* angle_m PI)))
(* 0.011111111111111112 (* angle_m (* PI (* a (- (- a) b))))))))angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if (angle_m <= 8.8e+173) {
tmp = (a * 0.011111111111111112) * ((b - a) * (angle_m * ((double) M_PI)));
} else {
tmp = 0.011111111111111112 * (angle_m * (((double) M_PI) * (a * (-a - b))));
}
return angle_s * tmp;
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
double tmp;
if (angle_m <= 8.8e+173) {
tmp = (a * 0.011111111111111112) * ((b - a) * (angle_m * Math.PI));
} else {
tmp = 0.011111111111111112 * (angle_m * (Math.PI * (a * (-a - b))));
}
return angle_s * tmp;
}
angle\_m = math.fabs(angle) angle\_s = math.copysign(1.0, angle) def code(angle_s, a, b, angle_m): tmp = 0 if angle_m <= 8.8e+173: tmp = (a * 0.011111111111111112) * ((b - a) * (angle_m * math.pi)) else: tmp = 0.011111111111111112 * (angle_m * (math.pi * (a * (-a - b)))) return angle_s * tmp
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) tmp = 0.0 if (angle_m <= 8.8e+173) tmp = Float64(Float64(a * 0.011111111111111112) * Float64(Float64(b - a) * Float64(angle_m * pi))); else tmp = Float64(0.011111111111111112 * Float64(angle_m * Float64(pi * Float64(a * Float64(Float64(-a) - b))))); end return Float64(angle_s * tmp) end
angle\_m = abs(angle); angle\_s = sign(angle) * abs(1.0); function tmp_2 = code(angle_s, a, b, angle_m) tmp = 0.0; if (angle_m <= 8.8e+173) tmp = (a * 0.011111111111111112) * ((b - a) * (angle_m * pi)); else tmp = 0.011111111111111112 * (angle_m * (pi * (a * (-a - b)))); end tmp_2 = angle_s * tmp; end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := N[(angle$95$s * If[LessEqual[angle$95$m, 8.8e+173], N[(N[(a * 0.011111111111111112), $MachinePrecision] * N[(N[(b - a), $MachinePrecision] * N[(angle$95$m * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.011111111111111112 * N[(angle$95$m * N[(Pi * N[(a * N[((-a) - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
angle\_s \cdot \begin{array}{l}
\mathbf{if}\;angle\_m \leq 8.8 \cdot 10^{+173}:\\
\;\;\;\;\left(a \cdot 0.011111111111111112\right) \cdot \left(\left(b - a\right) \cdot \left(angle\_m \cdot \pi\right)\right)\\
\mathbf{else}:\\
\;\;\;\;0.011111111111111112 \cdot \left(angle\_m \cdot \left(\pi \cdot \left(a \cdot \left(\left(-a\right) - b\right)\right)\right)\right)\\
\end{array}
\end{array}
if angle < 8.7999999999999999e173Initial program 53.3%
associate-*l*53.3%
*-commutative53.3%
associate-*l*53.3%
Simplified53.3%
unpow253.3%
unpow253.3%
difference-of-squares60.0%
Applied egg-rr60.0%
Taylor expanded in angle around inf 60.1%
associate-*r*59.6%
*-commutative59.6%
*-commutative59.6%
*-commutative59.6%
Simplified59.6%
Taylor expanded in b around 0 42.3%
Taylor expanded in angle around 0 43.5%
associate-*r*43.5%
associate-*r*43.5%
Simplified43.5%
if 8.7999999999999999e173 < angle Initial program 15.4%
associate-*l*15.4%
*-commutative15.4%
associate-*l*15.4%
Simplified15.4%
unpow215.4%
unpow215.4%
difference-of-squares15.4%
Applied egg-rr15.4%
2-sin15.4%
div-inv16.1%
metadata-eval16.1%
*-commutative16.1%
metadata-eval16.1%
div-inv15.4%
associate-*l*15.4%
div-inv16.1%
metadata-eval16.1%
Applied egg-rr16.1%
Taylor expanded in angle around 0 15.4%
Taylor expanded in b around 0 19.2%
neg-mul-119.2%
Simplified19.2%
Final simplification40.6%
angle\_m = (fabs.f64 angle) angle\_s = (copysign.f64 #s(literal 1 binary64) angle) (FPCore (angle_s a b angle_m) :precision binary64 (* angle_s (* (* (+ b a) (- b a)) (* (* angle_m PI) 0.011111111111111112))))
angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
return angle_s * (((b + a) * (b - a)) * ((angle_m * ((double) M_PI)) * 0.011111111111111112));
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
return angle_s * (((b + a) * (b - a)) * ((angle_m * Math.PI) * 0.011111111111111112));
}
angle\_m = math.fabs(angle) angle\_s = math.copysign(1.0, angle) def code(angle_s, a, b, angle_m): return angle_s * (((b + a) * (b - a)) * ((angle_m * math.pi) * 0.011111111111111112))
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) return Float64(angle_s * Float64(Float64(Float64(b + a) * Float64(b - a)) * Float64(Float64(angle_m * pi) * 0.011111111111111112))) end
angle\_m = abs(angle); angle\_s = sign(angle) * abs(1.0); function tmp = code(angle_s, a, b, angle_m) tmp = angle_s * (((b + a) * (b - a)) * ((angle_m * pi) * 0.011111111111111112)); end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := N[(angle$95$s * N[(N[(N[(b + a), $MachinePrecision] * N[(b - a), $MachinePrecision]), $MachinePrecision] * N[(N[(angle$95$m * Pi), $MachinePrecision] * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
angle\_s \cdot \left(\left(\left(b + a\right) \cdot \left(b - a\right)\right) \cdot \left(\left(angle\_m \cdot \pi\right) \cdot 0.011111111111111112\right)\right)
\end{array}
Initial program 48.8%
associate-*l*48.8%
*-commutative48.8%
associate-*l*48.8%
Simplified48.8%
unpow248.8%
unpow248.8%
difference-of-squares54.8%
Applied egg-rr54.8%
2-sin54.8%
div-inv54.4%
metadata-eval54.4%
*-commutative54.4%
metadata-eval54.4%
div-inv54.8%
associate-*l*54.8%
div-inv54.4%
metadata-eval54.4%
Applied egg-rr54.4%
Taylor expanded in angle around 0 53.4%
Final simplification53.4%
angle\_m = (fabs.f64 angle) angle\_s = (copysign.f64 #s(literal 1 binary64) angle) (FPCore (angle_s a b angle_m) :precision binary64 (* angle_s (* 0.011111111111111112 (* angle_m (* PI (* (+ b a) (- b a)))))))
angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
return angle_s * (0.011111111111111112 * (angle_m * (((double) M_PI) * ((b + a) * (b - a)))));
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
return angle_s * (0.011111111111111112 * (angle_m * (Math.PI * ((b + a) * (b - a)))));
}
angle\_m = math.fabs(angle) angle\_s = math.copysign(1.0, angle) def code(angle_s, a, b, angle_m): return angle_s * (0.011111111111111112 * (angle_m * (math.pi * ((b + a) * (b - a)))))
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) return Float64(angle_s * Float64(0.011111111111111112 * Float64(angle_m * Float64(pi * Float64(Float64(b + a) * Float64(b - a)))))) end
angle\_m = abs(angle); angle\_s = sign(angle) * abs(1.0); function tmp = code(angle_s, a, b, angle_m) tmp = angle_s * (0.011111111111111112 * (angle_m * (pi * ((b + a) * (b - a))))); end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := N[(angle$95$s * N[(0.011111111111111112 * N[(angle$95$m * N[(Pi * N[(N[(b + a), $MachinePrecision] * N[(b - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
angle\_s \cdot \left(0.011111111111111112 \cdot \left(angle\_m \cdot \left(\pi \cdot \left(\left(b + a\right) \cdot \left(b - a\right)\right)\right)\right)\right)
\end{array}
Initial program 48.8%
associate-*l*48.8%
*-commutative48.8%
associate-*l*48.8%
Simplified48.8%
unpow248.8%
unpow248.8%
difference-of-squares54.8%
Applied egg-rr54.8%
2-sin54.8%
div-inv54.4%
metadata-eval54.4%
*-commutative54.4%
metadata-eval54.4%
div-inv54.8%
associate-*l*54.8%
div-inv54.4%
metadata-eval54.4%
Applied egg-rr54.4%
Taylor expanded in angle around 0 53.4%
Final simplification53.4%
angle\_m = (fabs.f64 angle) angle\_s = (copysign.f64 #s(literal 1 binary64) angle) (FPCore (angle_s a b angle_m) :precision binary64 (* angle_s (* (* a 0.011111111111111112) (* (- b a) (* angle_m PI)))))
angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
return angle_s * ((a * 0.011111111111111112) * ((b - a) * (angle_m * ((double) M_PI))));
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
return angle_s * ((a * 0.011111111111111112) * ((b - a) * (angle_m * Math.PI)));
}
angle\_m = math.fabs(angle) angle\_s = math.copysign(1.0, angle) def code(angle_s, a, b, angle_m): return angle_s * ((a * 0.011111111111111112) * ((b - a) * (angle_m * math.pi)))
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) return Float64(angle_s * Float64(Float64(a * 0.011111111111111112) * Float64(Float64(b - a) * Float64(angle_m * pi)))) end
angle\_m = abs(angle); angle\_s = sign(angle) * abs(1.0); function tmp = code(angle_s, a, b, angle_m) tmp = angle_s * ((a * 0.011111111111111112) * ((b - a) * (angle_m * pi))); end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := N[(angle$95$s * N[(N[(a * 0.011111111111111112), $MachinePrecision] * N[(N[(b - a), $MachinePrecision] * N[(angle$95$m * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
angle\_s \cdot \left(\left(a \cdot 0.011111111111111112\right) \cdot \left(\left(b - a\right) \cdot \left(angle\_m \cdot \pi\right)\right)\right)
\end{array}
Initial program 48.8%
associate-*l*48.8%
*-commutative48.8%
associate-*l*48.8%
Simplified48.8%
unpow248.8%
unpow248.8%
difference-of-squares54.8%
Applied egg-rr54.8%
Taylor expanded in angle around inf 55.5%
associate-*r*54.8%
*-commutative54.8%
*-commutative54.8%
*-commutative54.8%
Simplified54.8%
Taylor expanded in b around 0 38.8%
Taylor expanded in angle around 0 41.3%
associate-*r*41.4%
associate-*r*41.4%
Simplified41.4%
Final simplification41.4%
angle\_m = (fabs.f64 angle) angle\_s = (copysign.f64 #s(literal 1 binary64) angle) (FPCore (angle_s a b angle_m) :precision binary64 (* angle_s (* 0.011111111111111112 (* a (* angle_m (* (- b a) PI))))))
angle\_m = fabs(angle);
angle\_s = copysign(1.0, angle);
double code(double angle_s, double a, double b, double angle_m) {
return angle_s * (0.011111111111111112 * (a * (angle_m * ((b - a) * ((double) M_PI)))));
}
angle\_m = Math.abs(angle);
angle\_s = Math.copySign(1.0, angle);
public static double code(double angle_s, double a, double b, double angle_m) {
return angle_s * (0.011111111111111112 * (a * (angle_m * ((b - a) * Math.PI))));
}
angle\_m = math.fabs(angle) angle\_s = math.copysign(1.0, angle) def code(angle_s, a, b, angle_m): return angle_s * (0.011111111111111112 * (a * (angle_m * ((b - a) * math.pi))))
angle\_m = abs(angle) angle\_s = copysign(1.0, angle) function code(angle_s, a, b, angle_m) return Float64(angle_s * Float64(0.011111111111111112 * Float64(a * Float64(angle_m * Float64(Float64(b - a) * pi))))) end
angle\_m = abs(angle); angle\_s = sign(angle) * abs(1.0); function tmp = code(angle_s, a, b, angle_m) tmp = angle_s * (0.011111111111111112 * (a * (angle_m * ((b - a) * pi)))); end
angle\_m = N[Abs[angle], $MachinePrecision]
angle\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[angle]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[angle$95$s_, a_, b_, angle$95$m_] := N[(angle$95$s * N[(0.011111111111111112 * N[(a * N[(angle$95$m * N[(N[(b - a), $MachinePrecision] * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
angle\_m = \left|angle\right|
\\
angle\_s = \mathsf{copysign}\left(1, angle\right)
\\
angle\_s \cdot \left(0.011111111111111112 \cdot \left(a \cdot \left(angle\_m \cdot \left(\left(b - a\right) \cdot \pi\right)\right)\right)\right)
\end{array}
Initial program 48.8%
associate-*l*48.8%
*-commutative48.8%
associate-*l*48.8%
Simplified48.8%
unpow248.8%
unpow248.8%
difference-of-squares54.8%
Applied egg-rr54.8%
Taylor expanded in angle around inf 55.5%
associate-*r*54.8%
*-commutative54.8%
*-commutative54.8%
*-commutative54.8%
Simplified54.8%
Taylor expanded in b around 0 38.8%
Taylor expanded in angle around 0 41.3%
Final simplification41.3%
herbie shell --seed 2024121
(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)))))