
(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 19 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}
a_m = (fabs.f64 a)
(FPCore (a_m b angle)
:precision binary64
(if (<= a_m 1.55e+163)
(*
(+ a_m b)
(*
(- b a_m)
(sin
(*
(pow PI 0.6666666666666666)
(* (cbrt PI) (* angle 0.011111111111111112))))))
(*
(+ a_m b)
(*
(- b a_m)
(*
angle
(fma
-2.2862368541380886e-7
(* (* angle angle) (* PI (* PI PI)))
(* PI 0.011111111111111112)))))))a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if (a_m <= 1.55e+163) {
tmp = (a_m + b) * ((b - a_m) * sin((pow(((double) M_PI), 0.6666666666666666) * (cbrt(((double) M_PI)) * (angle * 0.011111111111111112)))));
} else {
tmp = (a_m + b) * ((b - a_m) * (angle * fma(-2.2862368541380886e-7, ((angle * angle) * (((double) M_PI) * (((double) M_PI) * ((double) M_PI)))), (((double) M_PI) * 0.011111111111111112))));
}
return tmp;
}
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (a_m <= 1.55e+163) tmp = Float64(Float64(a_m + b) * Float64(Float64(b - a_m) * sin(Float64((pi ^ 0.6666666666666666) * Float64(cbrt(pi) * Float64(angle * 0.011111111111111112)))))); else tmp = Float64(Float64(a_m + b) * Float64(Float64(b - a_m) * Float64(angle * fma(-2.2862368541380886e-7, Float64(Float64(angle * angle) * Float64(pi * Float64(pi * pi))), Float64(pi * 0.011111111111111112))))); end return tmp end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[a$95$m, 1.55e+163], N[(N[(a$95$m + b), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * N[Sin[N[(N[Power[Pi, 0.6666666666666666], $MachinePrecision] * N[(N[Power[Pi, 1/3], $MachinePrecision] * N[(angle * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(a$95$m + b), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * N[(angle * N[(-2.2862368541380886e-7 * N[(N[(angle * angle), $MachinePrecision] * N[(Pi * N[(Pi * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;a\_m \leq 1.55 \cdot 10^{+163}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(\left(b - a\_m\right) \cdot \sin \left({\pi}^{0.6666666666666666} \cdot \left(\sqrt[3]{\pi} \cdot \left(angle \cdot 0.011111111111111112\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(\left(b - a\_m\right) \cdot \left(angle \cdot \mathsf{fma}\left(-2.2862368541380886 \cdot 10^{-7}, \left(angle \cdot angle\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), \pi \cdot 0.011111111111111112\right)\right)\right)\\
\end{array}
\end{array}
if a < 1.55000000000000014e163Initial program 55.4%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr65.3%
associate-*l*N/A
add-cube-cbrtN/A
associate-*l*N/A
*-lowering-*.f64N/A
pow1/3N/A
pow1/3N/A
pow-sqrN/A
pow-lowering-pow.f64N/A
PI-lowering-PI.f64N/A
metadata-evalN/A
*-commutativeN/A
*-lowering-*.f64N/A
cbrt-lowering-cbrt.f64N/A
PI-lowering-PI.f64N/A
*-commutativeN/A
*-lowering-*.f6464.4
Applied egg-rr64.4%
if 1.55000000000000014e163 < a Initial program 34.4%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr79.0%
Taylor expanded in angle around 0
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6488.3
Simplified88.3%
Final simplification68.4%
a_m = (fabs.f64 a)
(FPCore (a_m b angle)
:precision binary64
(let* ((t_0 (- (pow b 2.0) (pow a_m 2.0))) (t_1 (* PI (* PI PI))))
(if (<= t_0 -2e-96)
(*
(+ a_m b)
(*
(- b a_m)
(*
angle
(fma
-2.2862368541380886e-7
(* (* angle angle) t_1)
(* PI 0.011111111111111112)))))
(if (<= t_0 1e+288)
(fma (* 0.011111111111111112 (* PI angle)) (* b b) 0.0)
(*
(+ a_m b)
(*
angle
(fma
0.011111111111111112
(* (- b a_m) PI)
(*
(* -2.2862368541380886e-7 (* angle angle))
(* (- b a_m) t_1)))))))))a_m = fabs(a);
double code(double a_m, double b, double angle) {
double t_0 = pow(b, 2.0) - pow(a_m, 2.0);
double t_1 = ((double) M_PI) * (((double) M_PI) * ((double) M_PI));
double tmp;
if (t_0 <= -2e-96) {
tmp = (a_m + b) * ((b - a_m) * (angle * fma(-2.2862368541380886e-7, ((angle * angle) * t_1), (((double) M_PI) * 0.011111111111111112))));
} else if (t_0 <= 1e+288) {
tmp = fma((0.011111111111111112 * (((double) M_PI) * angle)), (b * b), 0.0);
} else {
tmp = (a_m + b) * (angle * fma(0.011111111111111112, ((b - a_m) * ((double) M_PI)), ((-2.2862368541380886e-7 * (angle * angle)) * ((b - a_m) * t_1))));
}
return tmp;
}
a_m = abs(a) function code(a_m, b, angle) t_0 = Float64((b ^ 2.0) - (a_m ^ 2.0)) t_1 = Float64(pi * Float64(pi * pi)) tmp = 0.0 if (t_0 <= -2e-96) tmp = Float64(Float64(a_m + b) * Float64(Float64(b - a_m) * Float64(angle * fma(-2.2862368541380886e-7, Float64(Float64(angle * angle) * t_1), Float64(pi * 0.011111111111111112))))); elseif (t_0 <= 1e+288) tmp = fma(Float64(0.011111111111111112 * Float64(pi * angle)), Float64(b * b), 0.0); else tmp = Float64(Float64(a_m + b) * Float64(angle * fma(0.011111111111111112, Float64(Float64(b - a_m) * pi), Float64(Float64(-2.2862368541380886e-7 * Float64(angle * angle)) * Float64(Float64(b - a_m) * t_1))))); end return tmp end
a_m = N[Abs[a], $MachinePrecision]
code[a$95$m_, b_, angle_] := Block[{t$95$0 = N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a$95$m, 2.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(Pi * N[(Pi * Pi), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-96], N[(N[(a$95$m + b), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * N[(angle * N[(-2.2862368541380886e-7 * N[(N[(angle * angle), $MachinePrecision] * t$95$1), $MachinePrecision] + N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 1e+288], N[(N[(0.011111111111111112 * N[(Pi * angle), $MachinePrecision]), $MachinePrecision] * N[(b * b), $MachinePrecision] + 0.0), $MachinePrecision], N[(N[(a$95$m + b), $MachinePrecision] * N[(angle * N[(0.011111111111111112 * N[(N[(b - a$95$m), $MachinePrecision] * Pi), $MachinePrecision] + N[(N[(-2.2862368541380886e-7 * N[(angle * angle), $MachinePrecision]), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
t_0 := {b}^{2} - {a\_m}^{2}\\
t_1 := \pi \cdot \left(\pi \cdot \pi\right)\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-96}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(\left(b - a\_m\right) \cdot \left(angle \cdot \mathsf{fma}\left(-2.2862368541380886 \cdot 10^{-7}, \left(angle \cdot angle\right) \cdot t\_1, \pi \cdot 0.011111111111111112\right)\right)\right)\\
\mathbf{elif}\;t\_0 \leq 10^{+288}:\\
\;\;\;\;\mathsf{fma}\left(0.011111111111111112 \cdot \left(\pi \cdot angle\right), b \cdot b, 0\right)\\
\mathbf{else}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(angle \cdot \mathsf{fma}\left(0.011111111111111112, \left(b - a\_m\right) \cdot \pi, \left(-2.2862368541380886 \cdot 10^{-7} \cdot \left(angle \cdot angle\right)\right) \cdot \left(\left(b - a\_m\right) \cdot t\_1\right)\right)\right)\\
\end{array}
\end{array}
if (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) < -1.9999999999999998e-96Initial program 50.7%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr62.1%
Taylor expanded in angle around 0
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6466.1
Simplified66.1%
if -1.9999999999999998e-96 < (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) < 1e288Initial program 65.4%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6461.4
Simplified61.4%
Taylor expanded in b around inf
Simplified62.2%
if 1e288 < (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) Initial program 33.9%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr78.3%
Taylor expanded in angle around 0
*-lowering-*.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
--lowering--.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
PI-lowering-PI.f6478.4
Simplified78.4%
Final simplification67.7%
a_m = (fabs.f64 a)
(FPCore (a_m b angle)
:precision binary64
(let* ((t_0 (- (pow b 2.0) (pow a_m 2.0)))
(t_1
(*
(+ a_m b)
(*
(- b a_m)
(*
angle
(fma
-2.2862368541380886e-7
(* (* angle angle) (* PI (* PI PI)))
(* PI 0.011111111111111112)))))))
(if (<= t_0 -2e-96)
t_1
(if (<= t_0 1e+288)
(fma (* 0.011111111111111112 (* PI angle)) (* b b) 0.0)
t_1))))a_m = fabs(a);
double code(double a_m, double b, double angle) {
double t_0 = pow(b, 2.0) - pow(a_m, 2.0);
double t_1 = (a_m + b) * ((b - a_m) * (angle * fma(-2.2862368541380886e-7, ((angle * angle) * (((double) M_PI) * (((double) M_PI) * ((double) M_PI)))), (((double) M_PI) * 0.011111111111111112))));
double tmp;
if (t_0 <= -2e-96) {
tmp = t_1;
} else if (t_0 <= 1e+288) {
tmp = fma((0.011111111111111112 * (((double) M_PI) * angle)), (b * b), 0.0);
} else {
tmp = t_1;
}
return tmp;
}
a_m = abs(a) function code(a_m, b, angle) t_0 = Float64((b ^ 2.0) - (a_m ^ 2.0)) t_1 = Float64(Float64(a_m + b) * Float64(Float64(b - a_m) * Float64(angle * fma(-2.2862368541380886e-7, Float64(Float64(angle * angle) * Float64(pi * Float64(pi * pi))), Float64(pi * 0.011111111111111112))))) tmp = 0.0 if (t_0 <= -2e-96) tmp = t_1; elseif (t_0 <= 1e+288) tmp = fma(Float64(0.011111111111111112 * Float64(pi * angle)), Float64(b * b), 0.0); else tmp = t_1; end return tmp end
a_m = N[Abs[a], $MachinePrecision]
code[a$95$m_, b_, angle_] := Block[{t$95$0 = N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a$95$m, 2.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(a$95$m + b), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * N[(angle * N[(-2.2862368541380886e-7 * N[(N[(angle * angle), $MachinePrecision] * N[(Pi * N[(Pi * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-96], t$95$1, If[LessEqual[t$95$0, 1e+288], N[(N[(0.011111111111111112 * N[(Pi * angle), $MachinePrecision]), $MachinePrecision] * N[(b * b), $MachinePrecision] + 0.0), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
t_0 := {b}^{2} - {a\_m}^{2}\\
t_1 := \left(a\_m + b\right) \cdot \left(\left(b - a\_m\right) \cdot \left(angle \cdot \mathsf{fma}\left(-2.2862368541380886 \cdot 10^{-7}, \left(angle \cdot angle\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), \pi \cdot 0.011111111111111112\right)\right)\right)\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-96}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 10^{+288}:\\
\;\;\;\;\mathsf{fma}\left(0.011111111111111112 \cdot \left(\pi \cdot angle\right), b \cdot b, 0\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) < -1.9999999999999998e-96 or 1e288 < (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) Initial program 44.1%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr68.5%
Taylor expanded in angle around 0
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6470.9
Simplified70.9%
if -1.9999999999999998e-96 < (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) < 1e288Initial program 65.4%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6461.4
Simplified61.4%
Taylor expanded in b around inf
Simplified62.2%
Final simplification67.7%
a_m = (fabs.f64 a) (FPCore (a_m b angle) :precision binary64 (if (<= (- (pow b 2.0) (pow a_m 2.0)) -2e-172) (* (* a_m PI) (* -0.011111111111111112 (* a_m angle))) (* (* 0.011111111111111112 (* PI angle)) (* b (- b a_m)))))
a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if ((pow(b, 2.0) - pow(a_m, 2.0)) <= -2e-172) {
tmp = (a_m * ((double) M_PI)) * (-0.011111111111111112 * (a_m * angle));
} else {
tmp = (0.011111111111111112 * (((double) M_PI) * angle)) * (b * (b - a_m));
}
return tmp;
}
a_m = Math.abs(a);
public static double code(double a_m, double b, double angle) {
double tmp;
if ((Math.pow(b, 2.0) - Math.pow(a_m, 2.0)) <= -2e-172) {
tmp = (a_m * Math.PI) * (-0.011111111111111112 * (a_m * angle));
} else {
tmp = (0.011111111111111112 * (Math.PI * angle)) * (b * (b - a_m));
}
return tmp;
}
a_m = math.fabs(a) def code(a_m, b, angle): tmp = 0 if (math.pow(b, 2.0) - math.pow(a_m, 2.0)) <= -2e-172: tmp = (a_m * math.pi) * (-0.011111111111111112 * (a_m * angle)) else: tmp = (0.011111111111111112 * (math.pi * angle)) * (b * (b - a_m)) return tmp
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (Float64((b ^ 2.0) - (a_m ^ 2.0)) <= -2e-172) tmp = Float64(Float64(a_m * pi) * Float64(-0.011111111111111112 * Float64(a_m * angle))); else tmp = Float64(Float64(0.011111111111111112 * Float64(pi * angle)) * Float64(b * Float64(b - a_m))); end return tmp end
a_m = abs(a); function tmp_2 = code(a_m, b, angle) tmp = 0.0; if (((b ^ 2.0) - (a_m ^ 2.0)) <= -2e-172) tmp = (a_m * pi) * (-0.011111111111111112 * (a_m * angle)); else tmp = (0.011111111111111112 * (pi * angle)) * (b * (b - a_m)); end tmp_2 = tmp; end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a$95$m, 2.0], $MachinePrecision]), $MachinePrecision], -2e-172], N[(N[(a$95$m * Pi), $MachinePrecision] * N[(-0.011111111111111112 * N[(a$95$m * angle), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.011111111111111112 * N[(Pi * angle), $MachinePrecision]), $MachinePrecision] * N[(b * N[(b - a$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;{b}^{2} - {a\_m}^{2} \leq -2 \cdot 10^{-172}:\\
\;\;\;\;\left(a\_m \cdot \pi\right) \cdot \left(-0.011111111111111112 \cdot \left(a\_m \cdot angle\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.011111111111111112 \cdot \left(\pi \cdot angle\right)\right) \cdot \left(b \cdot \left(b - a\_m\right)\right)\\
\end{array}
\end{array}
if (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) < -2.0000000000000001e-172Initial program 50.2%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6439.1
Simplified39.1%
Taylor expanded in b around 0
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6439.1
Simplified39.1%
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f6450.3
Applied egg-rr50.3%
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6451.2
Applied egg-rr51.2%
if -2.0000000000000001e-172 < (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) Initial program 53.0%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6459.6
Simplified59.6%
Taylor expanded in b around inf
Simplified59.1%
Final simplification55.8%
a_m = (fabs.f64 a) (FPCore (a_m b angle) :precision binary64 (if (<= (- (pow b 2.0) (pow a_m 2.0)) -2e-96) (* (* a_m PI) (* -0.011111111111111112 (* a_m angle))) (fma (* 0.011111111111111112 (* PI angle)) (* b b) 0.0)))
a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if ((pow(b, 2.0) - pow(a_m, 2.0)) <= -2e-96) {
tmp = (a_m * ((double) M_PI)) * (-0.011111111111111112 * (a_m * angle));
} else {
tmp = fma((0.011111111111111112 * (((double) M_PI) * angle)), (b * b), 0.0);
}
return tmp;
}
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (Float64((b ^ 2.0) - (a_m ^ 2.0)) <= -2e-96) tmp = Float64(Float64(a_m * pi) * Float64(-0.011111111111111112 * Float64(a_m * angle))); else tmp = fma(Float64(0.011111111111111112 * Float64(pi * angle)), Float64(b * b), 0.0); end return tmp end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a$95$m, 2.0], $MachinePrecision]), $MachinePrecision], -2e-96], N[(N[(a$95$m * Pi), $MachinePrecision] * N[(-0.011111111111111112 * N[(a$95$m * angle), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.011111111111111112 * N[(Pi * angle), $MachinePrecision]), $MachinePrecision] * N[(b * b), $MachinePrecision] + 0.0), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;{b}^{2} - {a\_m}^{2} \leq -2 \cdot 10^{-96}:\\
\;\;\;\;\left(a\_m \cdot \pi\right) \cdot \left(-0.011111111111111112 \cdot \left(a\_m \cdot angle\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(0.011111111111111112 \cdot \left(\pi \cdot angle\right), b \cdot b, 0\right)\\
\end{array}
\end{array}
if (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) < -1.9999999999999998e-96Initial program 50.7%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6439.4
Simplified39.4%
Taylor expanded in b around 0
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6439.4
Simplified39.4%
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f6451.6
Applied egg-rr51.6%
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6452.5
Applied egg-rr52.5%
if -1.9999999999999998e-96 < (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) Initial program 52.6%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6458.4
Simplified58.4%
Taylor expanded in b around inf
Simplified57.5%
Final simplification55.6%
a_m = (fabs.f64 a) (FPCore (a_m b angle) :precision binary64 (if (<= (- (pow b 2.0) (pow a_m 2.0)) -2e-96) (* (* a_m PI) (* -0.011111111111111112 (* a_m angle))) (* 0.011111111111111112 (* angle (* PI (* b b))))))
a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if ((pow(b, 2.0) - pow(a_m, 2.0)) <= -2e-96) {
tmp = (a_m * ((double) M_PI)) * (-0.011111111111111112 * (a_m * angle));
} else {
tmp = 0.011111111111111112 * (angle * (((double) M_PI) * (b * b)));
}
return tmp;
}
a_m = Math.abs(a);
public static double code(double a_m, double b, double angle) {
double tmp;
if ((Math.pow(b, 2.0) - Math.pow(a_m, 2.0)) <= -2e-96) {
tmp = (a_m * Math.PI) * (-0.011111111111111112 * (a_m * angle));
} else {
tmp = 0.011111111111111112 * (angle * (Math.PI * (b * b)));
}
return tmp;
}
a_m = math.fabs(a) def code(a_m, b, angle): tmp = 0 if (math.pow(b, 2.0) - math.pow(a_m, 2.0)) <= -2e-96: tmp = (a_m * math.pi) * (-0.011111111111111112 * (a_m * angle)) else: tmp = 0.011111111111111112 * (angle * (math.pi * (b * b))) return tmp
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (Float64((b ^ 2.0) - (a_m ^ 2.0)) <= -2e-96) tmp = Float64(Float64(a_m * pi) * Float64(-0.011111111111111112 * Float64(a_m * angle))); else tmp = Float64(0.011111111111111112 * Float64(angle * Float64(pi * Float64(b * b)))); end return tmp end
a_m = abs(a); function tmp_2 = code(a_m, b, angle) tmp = 0.0; if (((b ^ 2.0) - (a_m ^ 2.0)) <= -2e-96) tmp = (a_m * pi) * (-0.011111111111111112 * (a_m * angle)); else tmp = 0.011111111111111112 * (angle * (pi * (b * b))); end tmp_2 = tmp; end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a$95$m, 2.0], $MachinePrecision]), $MachinePrecision], -2e-96], N[(N[(a$95$m * Pi), $MachinePrecision] * N[(-0.011111111111111112 * N[(a$95$m * angle), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.011111111111111112 * N[(angle * N[(Pi * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;{b}^{2} - {a\_m}^{2} \leq -2 \cdot 10^{-96}:\\
\;\;\;\;\left(a\_m \cdot \pi\right) \cdot \left(-0.011111111111111112 \cdot \left(a\_m \cdot angle\right)\right)\\
\mathbf{else}:\\
\;\;\;\;0.011111111111111112 \cdot \left(angle \cdot \left(\pi \cdot \left(b \cdot b\right)\right)\right)\\
\end{array}
\end{array}
if (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) < -1.9999999999999998e-96Initial program 50.7%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6439.4
Simplified39.4%
Taylor expanded in b around 0
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6439.4
Simplified39.4%
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f6451.6
Applied egg-rr51.6%
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6452.5
Applied egg-rr52.5%
if -1.9999999999999998e-96 < (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) Initial program 52.6%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6458.4
Simplified58.4%
Taylor expanded in b around inf
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f6457.5
Simplified57.5%
Final simplification55.5%
a_m = (fabs.f64 a) (FPCore (a_m b angle) :precision binary64 (if (<= (- (pow b 2.0) (pow a_m 2.0)) -2e-96) (* -0.011111111111111112 (* a_m (* angle (* a_m PI)))) (* 0.011111111111111112 (* angle (* PI (* b b))))))
a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if ((pow(b, 2.0) - pow(a_m, 2.0)) <= -2e-96) {
tmp = -0.011111111111111112 * (a_m * (angle * (a_m * ((double) M_PI))));
} else {
tmp = 0.011111111111111112 * (angle * (((double) M_PI) * (b * b)));
}
return tmp;
}
a_m = Math.abs(a);
public static double code(double a_m, double b, double angle) {
double tmp;
if ((Math.pow(b, 2.0) - Math.pow(a_m, 2.0)) <= -2e-96) {
tmp = -0.011111111111111112 * (a_m * (angle * (a_m * Math.PI)));
} else {
tmp = 0.011111111111111112 * (angle * (Math.PI * (b * b)));
}
return tmp;
}
a_m = math.fabs(a) def code(a_m, b, angle): tmp = 0 if (math.pow(b, 2.0) - math.pow(a_m, 2.0)) <= -2e-96: tmp = -0.011111111111111112 * (a_m * (angle * (a_m * math.pi))) else: tmp = 0.011111111111111112 * (angle * (math.pi * (b * b))) return tmp
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (Float64((b ^ 2.0) - (a_m ^ 2.0)) <= -2e-96) tmp = Float64(-0.011111111111111112 * Float64(a_m * Float64(angle * Float64(a_m * pi)))); else tmp = Float64(0.011111111111111112 * Float64(angle * Float64(pi * Float64(b * b)))); end return tmp end
a_m = abs(a); function tmp_2 = code(a_m, b, angle) tmp = 0.0; if (((b ^ 2.0) - (a_m ^ 2.0)) <= -2e-96) tmp = -0.011111111111111112 * (a_m * (angle * (a_m * pi))); else tmp = 0.011111111111111112 * (angle * (pi * (b * b))); end tmp_2 = tmp; end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a$95$m, 2.0], $MachinePrecision]), $MachinePrecision], -2e-96], N[(-0.011111111111111112 * N[(a$95$m * N[(angle * N[(a$95$m * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.011111111111111112 * N[(angle * N[(Pi * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;{b}^{2} - {a\_m}^{2} \leq -2 \cdot 10^{-96}:\\
\;\;\;\;-0.011111111111111112 \cdot \left(a\_m \cdot \left(angle \cdot \left(a\_m \cdot \pi\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;0.011111111111111112 \cdot \left(angle \cdot \left(\pi \cdot \left(b \cdot b\right)\right)\right)\\
\end{array}
\end{array}
if (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) < -1.9999999999999998e-96Initial program 50.7%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6439.4
Simplified39.4%
Taylor expanded in b around 0
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6439.4
Simplified39.4%
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f6451.6
Applied egg-rr51.6%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6451.7
Applied egg-rr51.7%
if -1.9999999999999998e-96 < (-.f64 (pow.f64 b #s(literal 2 binary64)) (pow.f64 a #s(literal 2 binary64))) Initial program 52.6%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6458.4
Simplified58.4%
Taylor expanded in b around inf
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f6457.5
Simplified57.5%
Final simplification55.2%
a_m = (fabs.f64 a)
(FPCore (a_m b angle)
:precision binary64
(if (<= (pow b 2.0) 5e-134)
(* a_m (* (- b a_m) (sin (* angle (* PI 0.011111111111111112)))))
(*
(+ a_m b)
(*
(- b a_m)
(*
angle
(fma
-2.2862368541380886e-7
(* (* angle angle) (* PI (* PI PI)))
(* PI 0.011111111111111112)))))))a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if (pow(b, 2.0) <= 5e-134) {
tmp = a_m * ((b - a_m) * sin((angle * (((double) M_PI) * 0.011111111111111112))));
} else {
tmp = (a_m + b) * ((b - a_m) * (angle * fma(-2.2862368541380886e-7, ((angle * angle) * (((double) M_PI) * (((double) M_PI) * ((double) M_PI)))), (((double) M_PI) * 0.011111111111111112))));
}
return tmp;
}
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if ((b ^ 2.0) <= 5e-134) tmp = Float64(a_m * Float64(Float64(b - a_m) * sin(Float64(angle * Float64(pi * 0.011111111111111112))))); else tmp = Float64(Float64(a_m + b) * Float64(Float64(b - a_m) * Float64(angle * fma(-2.2862368541380886e-7, Float64(Float64(angle * angle) * Float64(pi * Float64(pi * pi))), Float64(pi * 0.011111111111111112))))); end return tmp end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[N[Power[b, 2.0], $MachinePrecision], 5e-134], N[(a$95$m * N[(N[(b - a$95$m), $MachinePrecision] * N[Sin[N[(angle * N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(a$95$m + b), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * N[(angle * N[(-2.2862368541380886e-7 * N[(N[(angle * angle), $MachinePrecision] * N[(Pi * N[(Pi * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;{b}^{2} \leq 5 \cdot 10^{-134}:\\
\;\;\;\;a\_m \cdot \left(\left(b - a\_m\right) \cdot \sin \left(angle \cdot \left(\pi \cdot 0.011111111111111112\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(\left(b - a\_m\right) \cdot \left(angle \cdot \mathsf{fma}\left(-2.2862368541380886 \cdot 10^{-7}, \left(angle \cdot angle\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), \pi \cdot 0.011111111111111112\right)\right)\right)\\
\end{array}
\end{array}
if (pow.f64 b #s(literal 2 binary64)) < 5.0000000000000003e-134Initial program 58.2%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr66.3%
*-commutativeN/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6467.4
Applied egg-rr67.4%
Taylor expanded in b around 0
Simplified66.9%
if 5.0000000000000003e-134 < (pow.f64 b #s(literal 2 binary64)) Initial program 47.1%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr68.6%
Taylor expanded in angle around 0
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6468.5
Simplified68.5%
Final simplification67.8%
a_m = (fabs.f64 a)
(FPCore (a_m b angle)
:precision binary64
(if (<= (pow a_m 2.0) 2e-104)
(* (sin (* 0.011111111111111112 (* PI angle))) (* b b))
(*
(+ a_m b)
(*
(- b a_m)
(*
angle
(fma
-2.2862368541380886e-7
(* (* angle angle) (* PI (* PI PI)))
(* PI 0.011111111111111112)))))))a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if (pow(a_m, 2.0) <= 2e-104) {
tmp = sin((0.011111111111111112 * (((double) M_PI) * angle))) * (b * b);
} else {
tmp = (a_m + b) * ((b - a_m) * (angle * fma(-2.2862368541380886e-7, ((angle * angle) * (((double) M_PI) * (((double) M_PI) * ((double) M_PI)))), (((double) M_PI) * 0.011111111111111112))));
}
return tmp;
}
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if ((a_m ^ 2.0) <= 2e-104) tmp = Float64(sin(Float64(0.011111111111111112 * Float64(pi * angle))) * Float64(b * b)); else tmp = Float64(Float64(a_m + b) * Float64(Float64(b - a_m) * Float64(angle * fma(-2.2862368541380886e-7, Float64(Float64(angle * angle) * Float64(pi * Float64(pi * pi))), Float64(pi * 0.011111111111111112))))); end return tmp end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[N[Power[a$95$m, 2.0], $MachinePrecision], 2e-104], N[(N[Sin[N[(0.011111111111111112 * N[(Pi * angle), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(b * b), $MachinePrecision]), $MachinePrecision], N[(N[(a$95$m + b), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * N[(angle * N[(-2.2862368541380886e-7 * N[(N[(angle * angle), $MachinePrecision] * N[(Pi * N[(Pi * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;{a\_m}^{2} \leq 2 \cdot 10^{-104}:\\
\;\;\;\;\sin \left(0.011111111111111112 \cdot \left(\pi \cdot angle\right)\right) \cdot \left(b \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(\left(b - a\_m\right) \cdot \left(angle \cdot \mathsf{fma}\left(-2.2862368541380886 \cdot 10^{-7}, \left(angle \cdot angle\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), \pi \cdot 0.011111111111111112\right)\right)\right)\\
\end{array}
\end{array}
if (pow.f64 a #s(literal 2 binary64)) < 1.99999999999999985e-104Initial program 62.6%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr67.7%
Taylor expanded in b around inf
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f6464.3
Simplified64.3%
if 1.99999999999999985e-104 < (pow.f64 a #s(literal 2 binary64)) Initial program 44.4%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr67.5%
Taylor expanded in angle around 0
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6469.9
Simplified69.9%
Final simplification67.6%
a_m = (fabs.f64 a)
(FPCore (a_m b angle)
:precision binary64
(if (<= (/ angle 180.0) 2e+78)
(*
(+ a_m b)
(*
(- b a_m)
(*
angle
(fma
-2.2862368541380886e-7
(* (* angle angle) (* PI (* PI PI)))
(* PI 0.011111111111111112)))))
(* (sin (* 0.011111111111111112 (* PI angle))) (* (+ a_m b) (- b a_m)))))a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if ((angle / 180.0) <= 2e+78) {
tmp = (a_m + b) * ((b - a_m) * (angle * fma(-2.2862368541380886e-7, ((angle * angle) * (((double) M_PI) * (((double) M_PI) * ((double) M_PI)))), (((double) M_PI) * 0.011111111111111112))));
} else {
tmp = sin((0.011111111111111112 * (((double) M_PI) * angle))) * ((a_m + b) * (b - a_m));
}
return tmp;
}
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (Float64(angle / 180.0) <= 2e+78) tmp = Float64(Float64(a_m + b) * Float64(Float64(b - a_m) * Float64(angle * fma(-2.2862368541380886e-7, Float64(Float64(angle * angle) * Float64(pi * Float64(pi * pi))), Float64(pi * 0.011111111111111112))))); else tmp = Float64(sin(Float64(0.011111111111111112 * Float64(pi * angle))) * Float64(Float64(a_m + b) * Float64(b - a_m))); end return tmp end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[N[(angle / 180.0), $MachinePrecision], 2e+78], N[(N[(a$95$m + b), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * N[(angle * N[(-2.2862368541380886e-7 * N[(N[(angle * angle), $MachinePrecision] * N[(Pi * N[(Pi * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Sin[N[(0.011111111111111112 * N[(Pi * angle), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(N[(a$95$m + b), $MachinePrecision] * N[(b - a$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;\frac{angle}{180} \leq 2 \cdot 10^{+78}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(\left(b - a\_m\right) \cdot \left(angle \cdot \mathsf{fma}\left(-2.2862368541380886 \cdot 10^{-7}, \left(angle \cdot angle\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), \pi \cdot 0.011111111111111112\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin \left(0.011111111111111112 \cdot \left(\pi \cdot angle\right)\right) \cdot \left(\left(a\_m + b\right) \cdot \left(b - a\_m\right)\right)\\
\end{array}
\end{array}
if (/.f64 angle #s(literal 180 binary64)) < 2.00000000000000002e78Initial program 57.5%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr76.5%
Taylor expanded in angle around 0
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6475.3
Simplified75.3%
if 2.00000000000000002e78 < (/.f64 angle #s(literal 180 binary64)) Initial program 30.3%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
sin-lowering-sin.f64N/A
associate-*r/N/A
div-invN/A
associate-*r/N/A
div-invN/A
Applied egg-rr33.6%
Final simplification66.7%
a_m = (fabs.f64 a)
(FPCore (a_m b angle)
:precision binary64
(if (<= a_m 1.35e+165)
(* (+ a_m b) (* (- b a_m) (sin (* 0.011111111111111112 (* PI angle)))))
(*
(+ a_m b)
(*
(- b a_m)
(*
angle
(fma
-2.2862368541380886e-7
(* (* angle angle) (* PI (* PI PI)))
(* PI 0.011111111111111112)))))))a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if (a_m <= 1.35e+165) {
tmp = (a_m + b) * ((b - a_m) * sin((0.011111111111111112 * (((double) M_PI) * angle))));
} else {
tmp = (a_m + b) * ((b - a_m) * (angle * fma(-2.2862368541380886e-7, ((angle * angle) * (((double) M_PI) * (((double) M_PI) * ((double) M_PI)))), (((double) M_PI) * 0.011111111111111112))));
}
return tmp;
}
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (a_m <= 1.35e+165) tmp = Float64(Float64(a_m + b) * Float64(Float64(b - a_m) * sin(Float64(0.011111111111111112 * Float64(pi * angle))))); else tmp = Float64(Float64(a_m + b) * Float64(Float64(b - a_m) * Float64(angle * fma(-2.2862368541380886e-7, Float64(Float64(angle * angle) * Float64(pi * Float64(pi * pi))), Float64(pi * 0.011111111111111112))))); end return tmp end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[a$95$m, 1.35e+165], N[(N[(a$95$m + b), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * N[Sin[N[(0.011111111111111112 * N[(Pi * angle), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(a$95$m + b), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * N[(angle * N[(-2.2862368541380886e-7 * N[(N[(angle * angle), $MachinePrecision] * N[(Pi * N[(Pi * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;a\_m \leq 1.35 \cdot 10^{+165}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(\left(b - a\_m\right) \cdot \sin \left(0.011111111111111112 \cdot \left(\pi \cdot angle\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(\left(b - a\_m\right) \cdot \left(angle \cdot \mathsf{fma}\left(-2.2862368541380886 \cdot 10^{-7}, \left(angle \cdot angle\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), \pi \cdot 0.011111111111111112\right)\right)\right)\\
\end{array}
\end{array}
if a < 1.35e165Initial program 55.4%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr65.3%
if 1.35e165 < a Initial program 34.4%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr79.0%
Taylor expanded in angle around 0
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6488.3
Simplified88.3%
Final simplification69.2%
a_m = (fabs.f64 a)
(FPCore (a_m b angle)
:precision binary64
(if (<= a_m 5e+164)
(* (+ a_m b) (* (- b a_m) (sin (* angle (* PI 0.011111111111111112)))))
(*
(+ a_m b)
(*
(- b a_m)
(*
angle
(fma
-2.2862368541380886e-7
(* (* angle angle) (* PI (* PI PI)))
(* PI 0.011111111111111112)))))))a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if (a_m <= 5e+164) {
tmp = (a_m + b) * ((b - a_m) * sin((angle * (((double) M_PI) * 0.011111111111111112))));
} else {
tmp = (a_m + b) * ((b - a_m) * (angle * fma(-2.2862368541380886e-7, ((angle * angle) * (((double) M_PI) * (((double) M_PI) * ((double) M_PI)))), (((double) M_PI) * 0.011111111111111112))));
}
return tmp;
}
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (a_m <= 5e+164) tmp = Float64(Float64(a_m + b) * Float64(Float64(b - a_m) * sin(Float64(angle * Float64(pi * 0.011111111111111112))))); else tmp = Float64(Float64(a_m + b) * Float64(Float64(b - a_m) * Float64(angle * fma(-2.2862368541380886e-7, Float64(Float64(angle * angle) * Float64(pi * Float64(pi * pi))), Float64(pi * 0.011111111111111112))))); end return tmp end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[a$95$m, 5e+164], N[(N[(a$95$m + b), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * N[Sin[N[(angle * N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(a$95$m + b), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * N[(angle * N[(-2.2862368541380886e-7 * N[(N[(angle * angle), $MachinePrecision] * N[(Pi * N[(Pi * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;a\_m \leq 5 \cdot 10^{+164}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(\left(b - a\_m\right) \cdot \sin \left(angle \cdot \left(\pi \cdot 0.011111111111111112\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(\left(b - a\_m\right) \cdot \left(angle \cdot \mathsf{fma}\left(-2.2862368541380886 \cdot 10^{-7}, \left(angle \cdot angle\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), \pi \cdot 0.011111111111111112\right)\right)\right)\\
\end{array}
\end{array}
if a < 4.9999999999999995e164Initial program 55.4%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr65.3%
*-commutativeN/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6465.8
Applied egg-rr65.8%
if 4.9999999999999995e164 < a Initial program 34.4%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr79.0%
Taylor expanded in angle around 0
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6488.3
Simplified88.3%
Final simplification69.6%
a_m = (fabs.f64 a)
(FPCore (a_m b angle)
:precision binary64
(if (<= a_m 1.05e-86)
(* (* b 2.0) (* b (sin (* PI (* angle 0.005555555555555556)))))
(if (<= a_m 3.1e+163)
(* (* (- b a_m) (* PI angle)) (* (+ a_m b) 0.011111111111111112))
(*
(+ a_m b)
(*
(- b a_m)
(*
angle
(fma
-2.2862368541380886e-7
(* (* angle angle) (* PI (* PI PI)))
(* PI 0.011111111111111112))))))))a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if (a_m <= 1.05e-86) {
tmp = (b * 2.0) * (b * sin((((double) M_PI) * (angle * 0.005555555555555556))));
} else if (a_m <= 3.1e+163) {
tmp = ((b - a_m) * (((double) M_PI) * angle)) * ((a_m + b) * 0.011111111111111112);
} else {
tmp = (a_m + b) * ((b - a_m) * (angle * fma(-2.2862368541380886e-7, ((angle * angle) * (((double) M_PI) * (((double) M_PI) * ((double) M_PI)))), (((double) M_PI) * 0.011111111111111112))));
}
return tmp;
}
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (a_m <= 1.05e-86) tmp = Float64(Float64(b * 2.0) * Float64(b * sin(Float64(pi * Float64(angle * 0.005555555555555556))))); elseif (a_m <= 3.1e+163) tmp = Float64(Float64(Float64(b - a_m) * Float64(pi * angle)) * Float64(Float64(a_m + b) * 0.011111111111111112)); else tmp = Float64(Float64(a_m + b) * Float64(Float64(b - a_m) * Float64(angle * fma(-2.2862368541380886e-7, Float64(Float64(angle * angle) * Float64(pi * Float64(pi * pi))), Float64(pi * 0.011111111111111112))))); end return tmp end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[a$95$m, 1.05e-86], N[(N[(b * 2.0), $MachinePrecision] * N[(b * N[Sin[N[(Pi * N[(angle * 0.005555555555555556), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a$95$m, 3.1e+163], N[(N[(N[(b - a$95$m), $MachinePrecision] * N[(Pi * angle), $MachinePrecision]), $MachinePrecision] * N[(N[(a$95$m + b), $MachinePrecision] * 0.011111111111111112), $MachinePrecision]), $MachinePrecision], N[(N[(a$95$m + b), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * N[(angle * N[(-2.2862368541380886e-7 * N[(N[(angle * angle), $MachinePrecision] * N[(Pi * N[(Pi * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;a\_m \leq 1.05 \cdot 10^{-86}:\\
\;\;\;\;\left(b \cdot 2\right) \cdot \left(b \cdot \sin \left(\pi \cdot \left(angle \cdot 0.005555555555555556\right)\right)\right)\\
\mathbf{elif}\;a\_m \leq 3.1 \cdot 10^{+163}:\\
\;\;\;\;\left(\left(b - a\_m\right) \cdot \left(\pi \cdot angle\right)\right) \cdot \left(\left(a\_m + b\right) \cdot 0.011111111111111112\right)\\
\mathbf{else}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(\left(b - a\_m\right) \cdot \left(angle \cdot \mathsf{fma}\left(-2.2862368541380886 \cdot 10^{-7}, \left(angle \cdot angle\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), \pi \cdot 0.011111111111111112\right)\right)\right)\\
\end{array}
\end{array}
if a < 1.05e-86Initial program 55.2%
Taylor expanded in b around inf
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
cos-lowering-cos.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
Simplified46.1%
Taylor expanded in angle around 0
Simplified49.1%
*-lft-identityN/A
associate-*l*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
div-invN/A
*-commutativeN/A
associate-*r/N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
div-invN/A
metadata-evalN/A
*-lowering-*.f6453.3
Applied egg-rr53.3%
if 1.05e-86 < a < 3.10000000000000029e163Initial program 55.8%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6451.3
Simplified51.3%
remove-double-divN/A
un-div-invN/A
*-commutativeN/A
*-commutativeN/A
associate-/r*N/A
div-invN/A
times-fracN/A
un-div-invN/A
clear-numN/A
/-rgt-identityN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
--lowering--.f64N/A
un-div-invN/A
remove-double-divN/A
*-lowering-*.f64N/A
+-lowering-+.f6458.0
Applied egg-rr58.0%
if 3.10000000000000029e163 < a Initial program 34.4%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr79.0%
Taylor expanded in angle around 0
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6488.3
Simplified88.3%
Final simplification60.2%
a_m = (fabs.f64 a)
(FPCore (a_m b angle)
:precision binary64
(if (<= a_m 1.06e-52)
(* (+ a_m b) (* b (sin (* angle (* PI 0.011111111111111112)))))
(*
(+ a_m b)
(*
(- b a_m)
(*
angle
(fma
-2.2862368541380886e-7
(* (* angle angle) (* PI (* PI PI)))
(* PI 0.011111111111111112)))))))a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if (a_m <= 1.06e-52) {
tmp = (a_m + b) * (b * sin((angle * (((double) M_PI) * 0.011111111111111112))));
} else {
tmp = (a_m + b) * ((b - a_m) * (angle * fma(-2.2862368541380886e-7, ((angle * angle) * (((double) M_PI) * (((double) M_PI) * ((double) M_PI)))), (((double) M_PI) * 0.011111111111111112))));
}
return tmp;
}
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (a_m <= 1.06e-52) tmp = Float64(Float64(a_m + b) * Float64(b * sin(Float64(angle * Float64(pi * 0.011111111111111112))))); else tmp = Float64(Float64(a_m + b) * Float64(Float64(b - a_m) * Float64(angle * fma(-2.2862368541380886e-7, Float64(Float64(angle * angle) * Float64(pi * Float64(pi * pi))), Float64(pi * 0.011111111111111112))))); end return tmp end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[a$95$m, 1.06e-52], N[(N[(a$95$m + b), $MachinePrecision] * N[(b * N[Sin[N[(angle * N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(a$95$m + b), $MachinePrecision] * N[(N[(b - a$95$m), $MachinePrecision] * N[(angle * N[(-2.2862368541380886e-7 * N[(N[(angle * angle), $MachinePrecision] * N[(Pi * N[(Pi * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(Pi * 0.011111111111111112), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;a\_m \leq 1.06 \cdot 10^{-52}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(b \cdot \sin \left(angle \cdot \left(\pi \cdot 0.011111111111111112\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(a\_m + b\right) \cdot \left(\left(b - a\_m\right) \cdot \left(angle \cdot \mathsf{fma}\left(-2.2862368541380886 \cdot 10^{-7}, \left(angle \cdot angle\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), \pi \cdot 0.011111111111111112\right)\right)\right)\\
\end{array}
\end{array}
if a < 1.06e-52Initial program 56.0%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr65.6%
*-commutativeN/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6466.8
Applied egg-rr66.8%
Taylor expanded in b around inf
Simplified51.9%
if 1.06e-52 < a Initial program 43.9%
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
associate-*l*N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
2-sinN/A
count-2N/A
Applied egg-rr71.5%
Taylor expanded in angle around 0
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6474.8
Simplified74.8%
Final simplification59.7%
a_m = (fabs.f64 a)
(FPCore (a_m b angle)
:precision binary64
(if (<= a_m 4e-184)
(fma (* 0.011111111111111112 (* PI angle)) (* b b) 0.0)
(if (<= a_m 5.05e+222)
(* (* (- b a_m) (* PI angle)) (* (+ a_m b) 0.011111111111111112))
(*
-2.0
(*
(* a_m a_m)
(*
angle
(fma
(* (* angle angle) -2.8577960676726107e-8)
(* PI (* PI PI))
(* PI 0.005555555555555556))))))))a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if (a_m <= 4e-184) {
tmp = fma((0.011111111111111112 * (((double) M_PI) * angle)), (b * b), 0.0);
} else if (a_m <= 5.05e+222) {
tmp = ((b - a_m) * (((double) M_PI) * angle)) * ((a_m + b) * 0.011111111111111112);
} else {
tmp = -2.0 * ((a_m * a_m) * (angle * fma(((angle * angle) * -2.8577960676726107e-8), (((double) M_PI) * (((double) M_PI) * ((double) M_PI))), (((double) M_PI) * 0.005555555555555556))));
}
return tmp;
}
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (a_m <= 4e-184) tmp = fma(Float64(0.011111111111111112 * Float64(pi * angle)), Float64(b * b), 0.0); elseif (a_m <= 5.05e+222) tmp = Float64(Float64(Float64(b - a_m) * Float64(pi * angle)) * Float64(Float64(a_m + b) * 0.011111111111111112)); else tmp = Float64(-2.0 * Float64(Float64(a_m * a_m) * Float64(angle * fma(Float64(Float64(angle * angle) * -2.8577960676726107e-8), Float64(pi * Float64(pi * pi)), Float64(pi * 0.005555555555555556))))); end return tmp end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[a$95$m, 4e-184], N[(N[(0.011111111111111112 * N[(Pi * angle), $MachinePrecision]), $MachinePrecision] * N[(b * b), $MachinePrecision] + 0.0), $MachinePrecision], If[LessEqual[a$95$m, 5.05e+222], N[(N[(N[(b - a$95$m), $MachinePrecision] * N[(Pi * angle), $MachinePrecision]), $MachinePrecision] * N[(N[(a$95$m + b), $MachinePrecision] * 0.011111111111111112), $MachinePrecision]), $MachinePrecision], N[(-2.0 * N[(N[(a$95$m * a$95$m), $MachinePrecision] * N[(angle * N[(N[(N[(angle * angle), $MachinePrecision] * -2.8577960676726107e-8), $MachinePrecision] * N[(Pi * N[(Pi * Pi), $MachinePrecision]), $MachinePrecision] + N[(Pi * 0.005555555555555556), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;a\_m \leq 4 \cdot 10^{-184}:\\
\;\;\;\;\mathsf{fma}\left(0.011111111111111112 \cdot \left(\pi \cdot angle\right), b \cdot b, 0\right)\\
\mathbf{elif}\;a\_m \leq 5.05 \cdot 10^{+222}:\\
\;\;\;\;\left(\left(b - a\_m\right) \cdot \left(\pi \cdot angle\right)\right) \cdot \left(\left(a\_m + b\right) \cdot 0.011111111111111112\right)\\
\mathbf{else}:\\
\;\;\;\;-2 \cdot \left(\left(a\_m \cdot a\_m\right) \cdot \left(angle \cdot \mathsf{fma}\left(\left(angle \cdot angle\right) \cdot -2.8577960676726107 \cdot 10^{-8}, \pi \cdot \left(\pi \cdot \pi\right), \pi \cdot 0.005555555555555556\right)\right)\right)\\
\end{array}
\end{array}
if a < 4.0000000000000002e-184Initial program 56.8%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6454.7
Simplified54.7%
Taylor expanded in b around inf
Simplified49.6%
if 4.0000000000000002e-184 < a < 5.0499999999999998e222Initial program 47.5%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6447.2
Simplified47.2%
remove-double-divN/A
un-div-invN/A
*-commutativeN/A
*-commutativeN/A
associate-/r*N/A
div-invN/A
times-fracN/A
un-div-invN/A
clear-numN/A
/-rgt-identityN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
--lowering--.f64N/A
un-div-invN/A
remove-double-divN/A
*-lowering-*.f64N/A
+-lowering-+.f6459.4
Applied egg-rr59.4%
if 5.0499999999999998e222 < a Initial program 44.3%
flip--N/A
clear-numN/A
/-lowering-/.f64N/A
clear-numN/A
flip--N/A
/-lowering-/.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6469.3
Applied egg-rr69.3%
Taylor expanded in b around 0
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
sin-lowering-sin.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6469.3
Simplified69.3%
Taylor expanded in angle around 0
unpow2N/A
*-lowering-*.f6458.6
Simplified58.6%
Taylor expanded in angle around 0
*-lowering-*.f64N/A
associate-*r*N/A
accelerator-lowering-fma.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
PI-lowering-PI.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
PI-lowering-PI.f6472.9
Simplified72.9%
Final simplification56.0%
a_m = (fabs.f64 a) (FPCore (a_m b angle) :precision binary64 (if (<= (/ angle 180.0) 100.0) (* (* (- b a_m) (* PI angle)) (* (+ a_m b) 0.011111111111111112)) (* (* 0.011111111111111112 (* PI angle)) (* b (+ a_m b)))))
a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if ((angle / 180.0) <= 100.0) {
tmp = ((b - a_m) * (((double) M_PI) * angle)) * ((a_m + b) * 0.011111111111111112);
} else {
tmp = (0.011111111111111112 * (((double) M_PI) * angle)) * (b * (a_m + b));
}
return tmp;
}
a_m = Math.abs(a);
public static double code(double a_m, double b, double angle) {
double tmp;
if ((angle / 180.0) <= 100.0) {
tmp = ((b - a_m) * (Math.PI * angle)) * ((a_m + b) * 0.011111111111111112);
} else {
tmp = (0.011111111111111112 * (Math.PI * angle)) * (b * (a_m + b));
}
return tmp;
}
a_m = math.fabs(a) def code(a_m, b, angle): tmp = 0 if (angle / 180.0) <= 100.0: tmp = ((b - a_m) * (math.pi * angle)) * ((a_m + b) * 0.011111111111111112) else: tmp = (0.011111111111111112 * (math.pi * angle)) * (b * (a_m + b)) return tmp
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (Float64(angle / 180.0) <= 100.0) tmp = Float64(Float64(Float64(b - a_m) * Float64(pi * angle)) * Float64(Float64(a_m + b) * 0.011111111111111112)); else tmp = Float64(Float64(0.011111111111111112 * Float64(pi * angle)) * Float64(b * Float64(a_m + b))); end return tmp end
a_m = abs(a); function tmp_2 = code(a_m, b, angle) tmp = 0.0; if ((angle / 180.0) <= 100.0) tmp = ((b - a_m) * (pi * angle)) * ((a_m + b) * 0.011111111111111112); else tmp = (0.011111111111111112 * (pi * angle)) * (b * (a_m + b)); end tmp_2 = tmp; end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[N[(angle / 180.0), $MachinePrecision], 100.0], N[(N[(N[(b - a$95$m), $MachinePrecision] * N[(Pi * angle), $MachinePrecision]), $MachinePrecision] * N[(N[(a$95$m + b), $MachinePrecision] * 0.011111111111111112), $MachinePrecision]), $MachinePrecision], N[(N[(0.011111111111111112 * N[(Pi * angle), $MachinePrecision]), $MachinePrecision] * N[(b * N[(a$95$m + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;\frac{angle}{180} \leq 100:\\
\;\;\;\;\left(\left(b - a\_m\right) \cdot \left(\pi \cdot angle\right)\right) \cdot \left(\left(a\_m + b\right) \cdot 0.011111111111111112\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.011111111111111112 \cdot \left(\pi \cdot angle\right)\right) \cdot \left(b \cdot \left(a\_m + b\right)\right)\\
\end{array}
\end{array}
if (/.f64 angle #s(literal 180 binary64)) < 100Initial program 58.8%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6461.5
Simplified61.5%
remove-double-divN/A
un-div-invN/A
*-commutativeN/A
*-commutativeN/A
associate-/r*N/A
div-invN/A
times-fracN/A
un-div-invN/A
clear-numN/A
/-rgt-identityN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
--lowering--.f64N/A
un-div-invN/A
remove-double-divN/A
*-lowering-*.f64N/A
+-lowering-+.f6470.3
Applied egg-rr70.3%
if 100 < (/.f64 angle #s(literal 180 binary64)) Initial program 32.1%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6421.4
Simplified21.4%
Taylor expanded in b around inf
Simplified20.4%
Final simplification57.2%
a_m = (fabs.f64 a) (FPCore (a_m b angle) :precision binary64 (if (<= a_m 9e-181) (* -0.011111111111111112 (* PI (* angle (* a_m a_m)))) (* -0.011111111111111112 (* a_m (* angle (* a_m PI))))))
a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if (a_m <= 9e-181) {
tmp = -0.011111111111111112 * (((double) M_PI) * (angle * (a_m * a_m)));
} else {
tmp = -0.011111111111111112 * (a_m * (angle * (a_m * ((double) M_PI))));
}
return tmp;
}
a_m = Math.abs(a);
public static double code(double a_m, double b, double angle) {
double tmp;
if (a_m <= 9e-181) {
tmp = -0.011111111111111112 * (Math.PI * (angle * (a_m * a_m)));
} else {
tmp = -0.011111111111111112 * (a_m * (angle * (a_m * Math.PI)));
}
return tmp;
}
a_m = math.fabs(a) def code(a_m, b, angle): tmp = 0 if a_m <= 9e-181: tmp = -0.011111111111111112 * (math.pi * (angle * (a_m * a_m))) else: tmp = -0.011111111111111112 * (a_m * (angle * (a_m * math.pi))) return tmp
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (a_m <= 9e-181) tmp = Float64(-0.011111111111111112 * Float64(pi * Float64(angle * Float64(a_m * a_m)))); else tmp = Float64(-0.011111111111111112 * Float64(a_m * Float64(angle * Float64(a_m * pi)))); end return tmp end
a_m = abs(a); function tmp_2 = code(a_m, b, angle) tmp = 0.0; if (a_m <= 9e-181) tmp = -0.011111111111111112 * (pi * (angle * (a_m * a_m))); else tmp = -0.011111111111111112 * (a_m * (angle * (a_m * pi))); end tmp_2 = tmp; end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[a$95$m, 9e-181], N[(-0.011111111111111112 * N[(Pi * N[(angle * N[(a$95$m * a$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-0.011111111111111112 * N[(a$95$m * N[(angle * N[(a$95$m * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;a\_m \leq 9 \cdot 10^{-181}:\\
\;\;\;\;-0.011111111111111112 \cdot \left(\pi \cdot \left(angle \cdot \left(a\_m \cdot a\_m\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;-0.011111111111111112 \cdot \left(a\_m \cdot \left(angle \cdot \left(a\_m \cdot \pi\right)\right)\right)\\
\end{array}
\end{array}
if a < 8.9999999999999998e-181Initial program 56.8%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6454.7
Simplified54.7%
Taylor expanded in b around 0
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6435.0
Simplified35.0%
if 8.9999999999999998e-181 < a Initial program 46.8%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6447.3
Simplified47.3%
Taylor expanded in b around 0
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6430.3
Simplified30.3%
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f6436.8
Applied egg-rr36.8%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6436.8
Applied egg-rr36.8%
Final simplification35.9%
a_m = (fabs.f64 a) (FPCore (a_m b angle) :precision binary64 (if (<= a_m 2e+95) (* -0.011111111111111112 (* PI (* angle (* a_m a_m)))) (* -0.011111111111111112 (* (* a_m PI) (* a_m angle)))))
a_m = fabs(a);
double code(double a_m, double b, double angle) {
double tmp;
if (a_m <= 2e+95) {
tmp = -0.011111111111111112 * (((double) M_PI) * (angle * (a_m * a_m)));
} else {
tmp = -0.011111111111111112 * ((a_m * ((double) M_PI)) * (a_m * angle));
}
return tmp;
}
a_m = Math.abs(a);
public static double code(double a_m, double b, double angle) {
double tmp;
if (a_m <= 2e+95) {
tmp = -0.011111111111111112 * (Math.PI * (angle * (a_m * a_m)));
} else {
tmp = -0.011111111111111112 * ((a_m * Math.PI) * (a_m * angle));
}
return tmp;
}
a_m = math.fabs(a) def code(a_m, b, angle): tmp = 0 if a_m <= 2e+95: tmp = -0.011111111111111112 * (math.pi * (angle * (a_m * a_m))) else: tmp = -0.011111111111111112 * ((a_m * math.pi) * (a_m * angle)) return tmp
a_m = abs(a) function code(a_m, b, angle) tmp = 0.0 if (a_m <= 2e+95) tmp = Float64(-0.011111111111111112 * Float64(pi * Float64(angle * Float64(a_m * a_m)))); else tmp = Float64(-0.011111111111111112 * Float64(Float64(a_m * pi) * Float64(a_m * angle))); end return tmp end
a_m = abs(a); function tmp_2 = code(a_m, b, angle) tmp = 0.0; if (a_m <= 2e+95) tmp = -0.011111111111111112 * (pi * (angle * (a_m * a_m))); else tmp = -0.011111111111111112 * ((a_m * pi) * (a_m * angle)); end tmp_2 = tmp; end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := If[LessEqual[a$95$m, 2e+95], N[(-0.011111111111111112 * N[(Pi * N[(angle * N[(a$95$m * a$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-0.011111111111111112 * N[(N[(a$95$m * Pi), $MachinePrecision] * N[(a$95$m * angle), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
a_m = \left|a\right|
\\
\begin{array}{l}
\mathbf{if}\;a\_m \leq 2 \cdot 10^{+95}:\\
\;\;\;\;-0.011111111111111112 \cdot \left(\pi \cdot \left(angle \cdot \left(a\_m \cdot a\_m\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;-0.011111111111111112 \cdot \left(\left(a\_m \cdot \pi\right) \cdot \left(a\_m \cdot angle\right)\right)\\
\end{array}
\end{array}
if a < 2.00000000000000004e95Initial program 55.3%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6451.7
Simplified51.7%
Taylor expanded in b around 0
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6431.0
Simplified31.0%
if 2.00000000000000004e95 < a Initial program 40.0%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6448.7
Simplified48.7%
Taylor expanded in b around 0
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6438.3
Simplified38.3%
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f6452.7
Applied egg-rr52.7%
Final simplification35.9%
a_m = (fabs.f64 a) (FPCore (a_m b angle) :precision binary64 (* -0.011111111111111112 (* (* a_m PI) (* a_m angle))))
a_m = fabs(a);
double code(double a_m, double b, double angle) {
return -0.011111111111111112 * ((a_m * ((double) M_PI)) * (a_m * angle));
}
a_m = Math.abs(a);
public static double code(double a_m, double b, double angle) {
return -0.011111111111111112 * ((a_m * Math.PI) * (a_m * angle));
}
a_m = math.fabs(a) def code(a_m, b, angle): return -0.011111111111111112 * ((a_m * math.pi) * (a_m * angle))
a_m = abs(a) function code(a_m, b, angle) return Float64(-0.011111111111111112 * Float64(Float64(a_m * pi) * Float64(a_m * angle))) end
a_m = abs(a); function tmp = code(a_m, b, angle) tmp = -0.011111111111111112 * ((a_m * pi) * (a_m * angle)); end
a_m = N[Abs[a], $MachinePrecision] code[a$95$m_, b_, angle_] := N[(-0.011111111111111112 * N[(N[(a$95$m * Pi), $MachinePrecision] * N[(a$95$m * angle), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
a_m = \left|a\right|
\\
-0.011111111111111112 \cdot \left(\left(a\_m \cdot \pi\right) \cdot \left(a\_m \cdot angle\right)\right)
\end{array}
Initial program 51.9%
Taylor expanded in angle around 0
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
unpow2N/A
unpow2N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
--lowering--.f6451.0
Simplified51.0%
Taylor expanded in b around 0
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
PI-lowering-PI.f6432.7
Simplified32.7%
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f6435.7
Applied egg-rr35.7%
Final simplification35.7%
herbie shell --seed 2024199
(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)))))