
(FPCore (A B C) :precision binary64 (* 180.0 (/ (atan (* (/ 1.0 B) (- (- C A) (sqrt (+ (pow (- A C) 2.0) (pow B 2.0)))))) PI)))
double code(double A, double B, double C) {
return 180.0 * (atan(((1.0 / B) * ((C - A) - sqrt((pow((A - C), 2.0) + pow(B, 2.0)))))) / ((double) M_PI));
}
public static double code(double A, double B, double C) {
return 180.0 * (Math.atan(((1.0 / B) * ((C - A) - Math.sqrt((Math.pow((A - C), 2.0) + Math.pow(B, 2.0)))))) / Math.PI);
}
def code(A, B, C): return 180.0 * (math.atan(((1.0 / B) * ((C - A) - math.sqrt((math.pow((A - C), 2.0) + math.pow(B, 2.0)))))) / math.pi)
function code(A, B, C) return Float64(180.0 * Float64(atan(Float64(Float64(1.0 / B) * Float64(Float64(C - A) - sqrt(Float64((Float64(A - C) ^ 2.0) + (B ^ 2.0)))))) / pi)) end
function tmp = code(A, B, C) tmp = 180.0 * (atan(((1.0 / B) * ((C - A) - sqrt((((A - C) ^ 2.0) + (B ^ 2.0)))))) / pi); end
code[A_, B_, C_] := N[(180.0 * N[(N[ArcTan[N[(N[(1.0 / B), $MachinePrecision] * N[(N[(C - A), $MachinePrecision] - N[Sqrt[N[(N[Power[N[(A - C), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[B, 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]
180 \cdot \frac{\tan^{-1} \left(\frac{1}{B} \cdot \left(\left(C - A\right) - \sqrt{{\left(A - C\right)}^{2} + {B}^{2}}\right)\right)}{\pi}
Herbie found 15 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (A B C) :precision binary64 (* 180.0 (/ (atan (* (/ 1.0 B) (- (- C A) (sqrt (+ (pow (- A C) 2.0) (pow B 2.0)))))) PI)))
double code(double A, double B, double C) {
return 180.0 * (atan(((1.0 / B) * ((C - A) - sqrt((pow((A - C), 2.0) + pow(B, 2.0)))))) / ((double) M_PI));
}
public static double code(double A, double B, double C) {
return 180.0 * (Math.atan(((1.0 / B) * ((C - A) - Math.sqrt((Math.pow((A - C), 2.0) + Math.pow(B, 2.0)))))) / Math.PI);
}
def code(A, B, C): return 180.0 * (math.atan(((1.0 / B) * ((C - A) - math.sqrt((math.pow((A - C), 2.0) + math.pow(B, 2.0)))))) / math.pi)
function code(A, B, C) return Float64(180.0 * Float64(atan(Float64(Float64(1.0 / B) * Float64(Float64(C - A) - sqrt(Float64((Float64(A - C) ^ 2.0) + (B ^ 2.0)))))) / pi)) end
function tmp = code(A, B, C) tmp = 180.0 * (atan(((1.0 / B) * ((C - A) - sqrt((((A - C) ^ 2.0) + (B ^ 2.0)))))) / pi); end
code[A_, B_, C_] := N[(180.0 * N[(N[ArcTan[N[(N[(1.0 / B), $MachinePrecision] * N[(N[(C - A), $MachinePrecision] - N[Sqrt[N[(N[Power[N[(A - C), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[B, 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]
180 \cdot \frac{\tan^{-1} \left(\frac{1}{B} \cdot \left(\left(C - A\right) - \sqrt{{\left(A - C\right)}^{2} + {B}^{2}}\right)\right)}{\pi}
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<=
(*
180.0
(/
(atan
(*
(/ 1.0 (fabs B))
(- (- C A) (sqrt (+ (pow (- A C) 2.0) (pow (fabs B) 2.0))))))
PI))
-4e-62)
(/ (* (atan (/ (- (- C A) (hypot (- C A) (fabs B))) (fabs B))) 180.0) PI)
(/
(*
(atan (fma -1.0 (/ (+ A (* -1.0 A)) (fabs B)) (* -0.5 (/ (fabs B) C))))
180.0)
PI))))double code(double A, double B, double C) {
double tmp;
if ((180.0 * (atan(((1.0 / fabs(B)) * ((C - A) - sqrt((pow((A - C), 2.0) + pow(fabs(B), 2.0)))))) / ((double) M_PI))) <= -4e-62) {
tmp = (atan((((C - A) - hypot((C - A), fabs(B))) / fabs(B))) * 180.0) / ((double) M_PI);
} else {
tmp = (atan(fma(-1.0, ((A + (-1.0 * A)) / fabs(B)), (-0.5 * (fabs(B) / C)))) * 180.0) / ((double) M_PI);
}
return copysign(1.0, B) * tmp;
}
function code(A, B, C) tmp = 0.0 if (Float64(180.0 * Float64(atan(Float64(Float64(1.0 / abs(B)) * Float64(Float64(C - A) - sqrt(Float64((Float64(A - C) ^ 2.0) + (abs(B) ^ 2.0)))))) / pi)) <= -4e-62) tmp = Float64(Float64(atan(Float64(Float64(Float64(C - A) - hypot(Float64(C - A), abs(B))) / abs(B))) * 180.0) / pi); else tmp = Float64(Float64(atan(fma(-1.0, Float64(Float64(A + Float64(-1.0 * A)) / abs(B)), Float64(-0.5 * Float64(abs(B) / C)))) * 180.0) / pi); end return Float64(copysign(1.0, B) * tmp) end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[N[(180.0 * N[(N[ArcTan[N[(N[(1.0 / N[Abs[B], $MachinePrecision]), $MachinePrecision] * N[(N[(C - A), $MachinePrecision] - N[Sqrt[N[(N[Power[N[(A - C), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[N[Abs[B], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], -4e-62], N[(N[(N[ArcTan[N[(N[(N[(C - A), $MachinePrecision] - N[Sqrt[N[(C - A), $MachinePrecision] ^ 2 + N[Abs[B], $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision], N[(N[(N[ArcTan[N[(-1.0 * N[(N[(A + N[(-1.0 * A), $MachinePrecision]), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision] + N[(-0.5 * N[(N[Abs[B], $MachinePrecision] / C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;180 \cdot \frac{\tan^{-1} \left(\frac{1}{\left|B\right|} \cdot \left(\left(C - A\right) - \sqrt{{\left(A - C\right)}^{2} + {\left(\left|B\right|\right)}^{2}}\right)\right)}{\pi} \leq -4 \cdot 10^{-62}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{\left(C - A\right) - \mathsf{hypot}\left(C - A, \left|B\right|\right)}{\left|B\right|}\right) \cdot 180}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\frac{\tan^{-1} \left(\mathsf{fma}\left(-1, \frac{A + -1 \cdot A}{\left|B\right|}, -0.5 \cdot \frac{\left|B\right|}{C}\right)\right) \cdot 180}{\pi}\\
\end{array}
if (*.f64 #s(literal 180 binary64) (/.f64 (atan.f64 (*.f64 (/.f64 #s(literal 1 binary64) B) (-.f64 (-.f64 C A) (sqrt.f64 (+.f64 (pow.f64 (-.f64 A C) #s(literal 2 binary64)) (pow.f64 B #s(literal 2 binary64))))))) (PI.f64))) < -4.0000000000000002e-62Initial program 53.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites53.8%
lift-sqrt.f64N/A
lift-fma.f64N/A
sqr-neg-revN/A
lift--.f64N/A
sub-negate-revN/A
lift--.f64N/A
lift--.f64N/A
sub-negate-revN/A
lift--.f64N/A
lift-*.f64N/A
lower-hypot.f6478.0%
Applied rewrites78.0%
if -4.0000000000000002e-62 < (*.f64 #s(literal 180 binary64) (/.f64 (atan.f64 (*.f64 (/.f64 #s(literal 1 binary64) B) (-.f64 (-.f64 C A) (sqrt.f64 (+.f64 (pow.f64 (-.f64 A C) #s(literal 2 binary64)) (pow.f64 B #s(literal 2 binary64))))))) (PI.f64))) Initial program 53.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites53.8%
lift-sqrt.f64N/A
lift-fma.f64N/A
sqr-neg-revN/A
lift--.f64N/A
sub-negate-revN/A
lift--.f64N/A
lift--.f64N/A
sub-negate-revN/A
lift--.f64N/A
lift-*.f64N/A
lower-hypot.f6478.0%
Applied rewrites78.0%
Taylor expanded in C around inf
lower-fma.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6425.7%
Applied rewrites25.7%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= C 3.8e+77)
(* (/ (atan (- (/ (- C A) (fabs B)) 1.0)) PI) 180.0)
(/
(*
(atan (fma -1.0 (/ (+ A (* -1.0 A)) (fabs B)) (* -0.5 (/ (fabs B) C))))
180.0)
PI))))double code(double A, double B, double C) {
double tmp;
if (C <= 3.8e+77) {
tmp = (atan((((C - A) / fabs(B)) - 1.0)) / ((double) M_PI)) * 180.0;
} else {
tmp = (atan(fma(-1.0, ((A + (-1.0 * A)) / fabs(B)), (-0.5 * (fabs(B) / C)))) * 180.0) / ((double) M_PI);
}
return copysign(1.0, B) * tmp;
}
function code(A, B, C) tmp = 0.0 if (C <= 3.8e+77) tmp = Float64(Float64(atan(Float64(Float64(Float64(C - A) / abs(B)) - 1.0)) / pi) * 180.0); else tmp = Float64(Float64(atan(fma(-1.0, Float64(Float64(A + Float64(-1.0 * A)) / abs(B)), Float64(-0.5 * Float64(abs(B) / C)))) * 180.0) / pi); end return Float64(copysign(1.0, B) * tmp) end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[C, 3.8e+77], N[(N[(N[ArcTan[N[(N[(N[(C - A), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision] * 180.0), $MachinePrecision], N[(N[(N[ArcTan[N[(-1.0 * N[(N[(A + N[(-1.0 * A), $MachinePrecision]), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision] + N[(-0.5 * N[(N[Abs[B], $MachinePrecision] / C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;C \leq 3.8 \cdot 10^{+77}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{C - A}{\left|B\right|} - 1\right)}{\pi} \cdot 180\\
\mathbf{else}:\\
\;\;\;\;\frac{\tan^{-1} \left(\mathsf{fma}\left(-1, \frac{A + -1 \cdot A}{\left|B\right|}, -0.5 \cdot \frac{\left|B\right|}{C}\right)\right) \cdot 180}{\pi}\\
\end{array}
if C < 3.8000000000000001e77Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6449.7%
Applied rewrites50.7%
if 3.8000000000000001e77 < C Initial program 53.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites53.8%
lift-sqrt.f64N/A
lift-fma.f64N/A
sqr-neg-revN/A
lift--.f64N/A
sub-negate-revN/A
lift--.f64N/A
lift--.f64N/A
sub-negate-revN/A
lift--.f64N/A
lift-*.f64N/A
lower-hypot.f6478.0%
Applied rewrites78.0%
Taylor expanded in C around inf
lower-fma.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6425.7%
Applied rewrites25.7%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= C 3.8e+77)
(* (/ (atan (- (/ (- C A) (fabs B)) 1.0)) PI) 180.0)
(*
180.0
(/
(atan (fma -1.0 (/ (+ A (* -1.0 A)) (fabs B)) (* -0.5 (/ (fabs B) C))))
PI)))))double code(double A, double B, double C) {
double tmp;
if (C <= 3.8e+77) {
tmp = (atan((((C - A) / fabs(B)) - 1.0)) / ((double) M_PI)) * 180.0;
} else {
tmp = 180.0 * (atan(fma(-1.0, ((A + (-1.0 * A)) / fabs(B)), (-0.5 * (fabs(B) / C)))) / ((double) M_PI));
}
return copysign(1.0, B) * tmp;
}
function code(A, B, C) tmp = 0.0 if (C <= 3.8e+77) tmp = Float64(Float64(atan(Float64(Float64(Float64(C - A) / abs(B)) - 1.0)) / pi) * 180.0); else tmp = Float64(180.0 * Float64(atan(fma(-1.0, Float64(Float64(A + Float64(-1.0 * A)) / abs(B)), Float64(-0.5 * Float64(abs(B) / C)))) / pi)); end return Float64(copysign(1.0, B) * tmp) end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[C, 3.8e+77], N[(N[(N[ArcTan[N[(N[(N[(C - A), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision] * 180.0), $MachinePrecision], N[(180.0 * N[(N[ArcTan[N[(-1.0 * N[(N[(A + N[(-1.0 * A), $MachinePrecision]), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision] + N[(-0.5 * N[(N[Abs[B], $MachinePrecision] / C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;C \leq 3.8 \cdot 10^{+77}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{C - A}{\left|B\right|} - 1\right)}{\pi} \cdot 180\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\mathsf{fma}\left(-1, \frac{A + -1 \cdot A}{\left|B\right|}, -0.5 \cdot \frac{\left|B\right|}{C}\right)\right)}{\pi}\\
\end{array}
if C < 3.8000000000000001e77Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6449.7%
Applied rewrites50.7%
if 3.8000000000000001e77 < C Initial program 53.8%
Taylor expanded in C around inf
lower-fma.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6425.7%
Applied rewrites25.7%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= A -2.5e+59)
(* (atan (* 0.5 (/ (fma (/ (fabs B) A) C (fabs B)) A))) (/ 180.0 PI))
(if (<= A 2.8e-174)
(/ (* (atan (/ (- C (hypot C (fabs B))) (fabs B))) 180.0) PI)
(* (/ (atan (- (/ (- C A) (fabs B)) 1.0)) PI) 180.0)))))double code(double A, double B, double C) {
double tmp;
if (A <= -2.5e+59) {
tmp = atan((0.5 * (fma((fabs(B) / A), C, fabs(B)) / A))) * (180.0 / ((double) M_PI));
} else if (A <= 2.8e-174) {
tmp = (atan(((C - hypot(C, fabs(B))) / fabs(B))) * 180.0) / ((double) M_PI);
} else {
tmp = (atan((((C - A) / fabs(B)) - 1.0)) / ((double) M_PI)) * 180.0;
}
return copysign(1.0, B) * tmp;
}
function code(A, B, C) tmp = 0.0 if (A <= -2.5e+59) tmp = Float64(atan(Float64(0.5 * Float64(fma(Float64(abs(B) / A), C, abs(B)) / A))) * Float64(180.0 / pi)); elseif (A <= 2.8e-174) tmp = Float64(Float64(atan(Float64(Float64(C - hypot(C, abs(B))) / abs(B))) * 180.0) / pi); else tmp = Float64(Float64(atan(Float64(Float64(Float64(C - A) / abs(B)) - 1.0)) / pi) * 180.0); end return Float64(copysign(1.0, B) * tmp) end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[A, -2.5e+59], N[(N[ArcTan[N[(0.5 * N[(N[(N[(N[Abs[B], $MachinePrecision] / A), $MachinePrecision] * C + N[Abs[B], $MachinePrecision]), $MachinePrecision] / A), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(180.0 / Pi), $MachinePrecision]), $MachinePrecision], If[LessEqual[A, 2.8e-174], N[(N[(N[ArcTan[N[(N[(C - N[Sqrt[C ^ 2 + N[Abs[B], $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision], N[(N[(N[ArcTan[N[(N[(N[(C - A), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision] * 180.0), $MachinePrecision]]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;A \leq -2.5 \cdot 10^{+59}:\\
\;\;\;\;\tan^{-1} \left(0.5 \cdot \frac{\mathsf{fma}\left(\frac{\left|B\right|}{A}, C, \left|B\right|\right)}{A}\right) \cdot \frac{180}{\pi}\\
\mathbf{elif}\;A \leq 2.8 \cdot 10^{-174}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{C - \mathsf{hypot}\left(C, \left|B\right|\right)}{\left|B\right|}\right) \cdot 180}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{C - A}{\left|B\right|} - 1\right)}{\pi} \cdot 180\\
\end{array}
if A < -2.4999999999999999e59Initial program 53.8%
Taylor expanded in A around -inf
lower-*.f64N/A
lower-/.f64N/A
lower-fma.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f6432.6%
Applied rewrites32.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites33.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
lower-*.f64N/A
Applied rewrites33.5%
if -2.4999999999999999e59 < A < 2.79999999999999999e-174Initial program 53.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites53.8%
lift-sqrt.f64N/A
lift-fma.f64N/A
sqr-neg-revN/A
lift--.f64N/A
sub-negate-revN/A
lift--.f64N/A
lift--.f64N/A
sub-negate-revN/A
lift--.f64N/A
lift-*.f64N/A
lower-hypot.f6478.0%
Applied rewrites78.0%
Taylor expanded in A around 0
Applied rewrites72.6%
Taylor expanded in A around 0
Applied rewrites63.5%
if 2.79999999999999999e-174 < A Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6449.7%
Applied rewrites50.7%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= A -3.8e-9)
(* (atan (* 0.5 (/ (fma (/ (fabs B) A) C (fabs B)) A))) (/ 180.0 PI))
(* (/ (atan (- (/ (- C A) (fabs B)) 1.0)) PI) 180.0))))double code(double A, double B, double C) {
double tmp;
if (A <= -3.8e-9) {
tmp = atan((0.5 * (fma((fabs(B) / A), C, fabs(B)) / A))) * (180.0 / ((double) M_PI));
} else {
tmp = (atan((((C - A) / fabs(B)) - 1.0)) / ((double) M_PI)) * 180.0;
}
return copysign(1.0, B) * tmp;
}
function code(A, B, C) tmp = 0.0 if (A <= -3.8e-9) tmp = Float64(atan(Float64(0.5 * Float64(fma(Float64(abs(B) / A), C, abs(B)) / A))) * Float64(180.0 / pi)); else tmp = Float64(Float64(atan(Float64(Float64(Float64(C - A) / abs(B)) - 1.0)) / pi) * 180.0); end return Float64(copysign(1.0, B) * tmp) end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[A, -3.8e-9], N[(N[ArcTan[N[(0.5 * N[(N[(N[(N[Abs[B], $MachinePrecision] / A), $MachinePrecision] * C + N[Abs[B], $MachinePrecision]), $MachinePrecision] / A), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(180.0 / Pi), $MachinePrecision]), $MachinePrecision], N[(N[(N[ArcTan[N[(N[(N[(C - A), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision] * 180.0), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;A \leq -3.8 \cdot 10^{-9}:\\
\;\;\;\;\tan^{-1} \left(0.5 \cdot \frac{\mathsf{fma}\left(\frac{\left|B\right|}{A}, C, \left|B\right|\right)}{A}\right) \cdot \frac{180}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{C - A}{\left|B\right|} - 1\right)}{\pi} \cdot 180\\
\end{array}
if A < -3.80000000000000011e-9Initial program 53.8%
Taylor expanded in A around -inf
lower-*.f64N/A
lower-/.f64N/A
lower-fma.f64N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f6432.6%
Applied rewrites32.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites33.5%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
lower-*.f64N/A
Applied rewrites33.5%
if -3.80000000000000011e-9 < A Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6449.7%
Applied rewrites50.7%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= A -2.5e+59)
(* (atan (* (/ (fabs B) A) 0.5)) (* (/ 1.0 PI) 180.0))
(* (/ (atan (- (/ (- C A) (fabs B)) 1.0)) PI) 180.0))))double code(double A, double B, double C) {
double tmp;
if (A <= -2.5e+59) {
tmp = atan(((fabs(B) / A) * 0.5)) * ((1.0 / ((double) M_PI)) * 180.0);
} else {
tmp = (atan((((C - A) / fabs(B)) - 1.0)) / ((double) M_PI)) * 180.0;
}
return copysign(1.0, B) * tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (A <= -2.5e+59) {
tmp = Math.atan(((Math.abs(B) / A) * 0.5)) * ((1.0 / Math.PI) * 180.0);
} else {
tmp = (Math.atan((((C - A) / Math.abs(B)) - 1.0)) / Math.PI) * 180.0;
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if A <= -2.5e+59: tmp = math.atan(((math.fabs(B) / A) * 0.5)) * ((1.0 / math.pi) * 180.0) else: tmp = (math.atan((((C - A) / math.fabs(B)) - 1.0)) / math.pi) * 180.0 return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (A <= -2.5e+59) tmp = Float64(atan(Float64(Float64(abs(B) / A) * 0.5)) * Float64(Float64(1.0 / pi) * 180.0)); else tmp = Float64(Float64(atan(Float64(Float64(Float64(C - A) / abs(B)) - 1.0)) / pi) * 180.0); end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (A <= -2.5e+59) tmp = atan(((abs(B) / A) * 0.5)) * ((1.0 / pi) * 180.0); else tmp = (atan((((C - A) / abs(B)) - 1.0)) / pi) * 180.0; end tmp_2 = (sign(B) * abs(1.0)) * tmp; end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[A, -2.5e+59], N[(N[ArcTan[N[(N[(N[Abs[B], $MachinePrecision] / A), $MachinePrecision] * 0.5), $MachinePrecision]], $MachinePrecision] * N[(N[(1.0 / Pi), $MachinePrecision] * 180.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[ArcTan[N[(N[(N[(C - A), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision] * 180.0), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;A \leq -2.5 \cdot 10^{+59}:\\
\;\;\;\;\tan^{-1} \left(\frac{\left|B\right|}{A} \cdot 0.5\right) \cdot \left(\frac{1}{\pi} \cdot 180\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{C - A}{\left|B\right|} - 1\right)}{\pi} \cdot 180\\
\end{array}
if A < -2.4999999999999999e59Initial program 53.8%
Taylor expanded in A around -inf
lower-*.f64N/A
lower-/.f6425.8%
Applied rewrites25.8%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
mult-flipN/A
associate-*l*N/A
lower-*.f64N/A
Applied rewrites25.8%
if -2.4999999999999999e59 < A Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6449.7%
Applied rewrites50.7%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= A -2.5e+59)
(/ (* (atan (* (/ (fabs B) A) 0.5)) 180.0) PI)
(* (/ (atan (- (/ (- C A) (fabs B)) 1.0)) PI) 180.0))))double code(double A, double B, double C) {
double tmp;
if (A <= -2.5e+59) {
tmp = (atan(((fabs(B) / A) * 0.5)) * 180.0) / ((double) M_PI);
} else {
tmp = (atan((((C - A) / fabs(B)) - 1.0)) / ((double) M_PI)) * 180.0;
}
return copysign(1.0, B) * tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (A <= -2.5e+59) {
tmp = (Math.atan(((Math.abs(B) / A) * 0.5)) * 180.0) / Math.PI;
} else {
tmp = (Math.atan((((C - A) / Math.abs(B)) - 1.0)) / Math.PI) * 180.0;
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if A <= -2.5e+59: tmp = (math.atan(((math.fabs(B) / A) * 0.5)) * 180.0) / math.pi else: tmp = (math.atan((((C - A) / math.fabs(B)) - 1.0)) / math.pi) * 180.0 return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (A <= -2.5e+59) tmp = Float64(Float64(atan(Float64(Float64(abs(B) / A) * 0.5)) * 180.0) / pi); else tmp = Float64(Float64(atan(Float64(Float64(Float64(C - A) / abs(B)) - 1.0)) / pi) * 180.0); end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (A <= -2.5e+59) tmp = (atan(((abs(B) / A) * 0.5)) * 180.0) / pi; else tmp = (atan((((C - A) / abs(B)) - 1.0)) / pi) * 180.0; end tmp_2 = (sign(B) * abs(1.0)) * tmp; end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[A, -2.5e+59], N[(N[(N[ArcTan[N[(N[(N[Abs[B], $MachinePrecision] / A), $MachinePrecision] * 0.5), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision], N[(N[(N[ArcTan[N[(N[(N[(C - A), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision] * 180.0), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;A \leq -2.5 \cdot 10^{+59}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{\left|B\right|}{A} \cdot 0.5\right) \cdot 180}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{C - A}{\left|B\right|} - 1\right)}{\pi} \cdot 180\\
\end{array}
if A < -2.4999999999999999e59Initial program 53.8%
Taylor expanded in A around -inf
lower-*.f64N/A
lower-/.f6425.8%
Applied rewrites25.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites25.8%
if -2.4999999999999999e59 < A Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6449.7%
Applied rewrites50.7%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= A -2.5e+59)
(/ (* (atan (* (/ (fabs B) A) 0.5)) 180.0) PI)
(if (<= A 9.2e+172)
(* 180.0 (/ (atan (- (/ C (fabs B)) 1.0)) PI))
(/ (* (atan (/ (* -2.0 A) (fabs B))) 180.0) PI)))))double code(double A, double B, double C) {
double tmp;
if (A <= -2.5e+59) {
tmp = (atan(((fabs(B) / A) * 0.5)) * 180.0) / ((double) M_PI);
} else if (A <= 9.2e+172) {
tmp = 180.0 * (atan(((C / fabs(B)) - 1.0)) / ((double) M_PI));
} else {
tmp = (atan(((-2.0 * A) / fabs(B))) * 180.0) / ((double) M_PI);
}
return copysign(1.0, B) * tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (A <= -2.5e+59) {
tmp = (Math.atan(((Math.abs(B) / A) * 0.5)) * 180.0) / Math.PI;
} else if (A <= 9.2e+172) {
tmp = 180.0 * (Math.atan(((C / Math.abs(B)) - 1.0)) / Math.PI);
} else {
tmp = (Math.atan(((-2.0 * A) / Math.abs(B))) * 180.0) / Math.PI;
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if A <= -2.5e+59: tmp = (math.atan(((math.fabs(B) / A) * 0.5)) * 180.0) / math.pi elif A <= 9.2e+172: tmp = 180.0 * (math.atan(((C / math.fabs(B)) - 1.0)) / math.pi) else: tmp = (math.atan(((-2.0 * A) / math.fabs(B))) * 180.0) / math.pi return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (A <= -2.5e+59) tmp = Float64(Float64(atan(Float64(Float64(abs(B) / A) * 0.5)) * 180.0) / pi); elseif (A <= 9.2e+172) tmp = Float64(180.0 * Float64(atan(Float64(Float64(C / abs(B)) - 1.0)) / pi)); else tmp = Float64(Float64(atan(Float64(Float64(-2.0 * A) / abs(B))) * 180.0) / pi); end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (A <= -2.5e+59) tmp = (atan(((abs(B) / A) * 0.5)) * 180.0) / pi; elseif (A <= 9.2e+172) tmp = 180.0 * (atan(((C / abs(B)) - 1.0)) / pi); else tmp = (atan(((-2.0 * A) / abs(B))) * 180.0) / pi; end tmp_2 = (sign(B) * abs(1.0)) * tmp; end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[A, -2.5e+59], N[(N[(N[ArcTan[N[(N[(N[Abs[B], $MachinePrecision] / A), $MachinePrecision] * 0.5), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision], If[LessEqual[A, 9.2e+172], N[(180.0 * N[(N[ArcTan[N[(N[(C / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(N[(N[ArcTan[N[(N[(-2.0 * A), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision]]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;A \leq -2.5 \cdot 10^{+59}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{\left|B\right|}{A} \cdot 0.5\right) \cdot 180}{\pi}\\
\mathbf{elif}\;A \leq 9.2 \cdot 10^{+172}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{C}{\left|B\right|} - 1\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{-2 \cdot A}{\left|B\right|}\right) \cdot 180}{\pi}\\
\end{array}
if A < -2.4999999999999999e59Initial program 53.8%
Taylor expanded in A around -inf
lower-*.f64N/A
lower-/.f6425.8%
Applied rewrites25.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites25.8%
if -2.4999999999999999e59 < A < 9.2000000000000003e172Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
Taylor expanded in A around 0
Applied rewrites38.9%
if 9.2000000000000003e172 < A Initial program 53.8%
Taylor expanded in A around inf
lower-*.f6423.7%
Applied rewrites23.7%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites23.7%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= A -2.5e+59)
(/ (* (atan (* (/ (fabs B) A) 0.5)) 180.0) PI)
(if (<= A 9.2e+172)
(* 180.0 (/ (atan (- (/ C (fabs B)) 1.0)) PI))
(* (/ (atan (/ (- A) (fabs B))) PI) 180.0)))))double code(double A, double B, double C) {
double tmp;
if (A <= -2.5e+59) {
tmp = (atan(((fabs(B) / A) * 0.5)) * 180.0) / ((double) M_PI);
} else if (A <= 9.2e+172) {
tmp = 180.0 * (atan(((C / fabs(B)) - 1.0)) / ((double) M_PI));
} else {
tmp = (atan((-A / fabs(B))) / ((double) M_PI)) * 180.0;
}
return copysign(1.0, B) * tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (A <= -2.5e+59) {
tmp = (Math.atan(((Math.abs(B) / A) * 0.5)) * 180.0) / Math.PI;
} else if (A <= 9.2e+172) {
tmp = 180.0 * (Math.atan(((C / Math.abs(B)) - 1.0)) / Math.PI);
} else {
tmp = (Math.atan((-A / Math.abs(B))) / Math.PI) * 180.0;
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if A <= -2.5e+59: tmp = (math.atan(((math.fabs(B) / A) * 0.5)) * 180.0) / math.pi elif A <= 9.2e+172: tmp = 180.0 * (math.atan(((C / math.fabs(B)) - 1.0)) / math.pi) else: tmp = (math.atan((-A / math.fabs(B))) / math.pi) * 180.0 return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (A <= -2.5e+59) tmp = Float64(Float64(atan(Float64(Float64(abs(B) / A) * 0.5)) * 180.0) / pi); elseif (A <= 9.2e+172) tmp = Float64(180.0 * Float64(atan(Float64(Float64(C / abs(B)) - 1.0)) / pi)); else tmp = Float64(Float64(atan(Float64(Float64(-A) / abs(B))) / pi) * 180.0); end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (A <= -2.5e+59) tmp = (atan(((abs(B) / A) * 0.5)) * 180.0) / pi; elseif (A <= 9.2e+172) tmp = 180.0 * (atan(((C / abs(B)) - 1.0)) / pi); else tmp = (atan((-A / abs(B))) / pi) * 180.0; end tmp_2 = (sign(B) * abs(1.0)) * tmp; end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[A, -2.5e+59], N[(N[(N[ArcTan[N[(N[(N[Abs[B], $MachinePrecision] / A), $MachinePrecision] * 0.5), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision], If[LessEqual[A, 9.2e+172], N[(180.0 * N[(N[ArcTan[N[(N[(C / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(N[(N[ArcTan[N[((-A) / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision] * 180.0), $MachinePrecision]]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;A \leq -2.5 \cdot 10^{+59}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{\left|B\right|}{A} \cdot 0.5\right) \cdot 180}{\pi}\\
\mathbf{elif}\;A \leq 9.2 \cdot 10^{+172}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{C}{\left|B\right|} - 1\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{-A}{\left|B\right|}\right)}{\pi} \cdot 180\\
\end{array}
if A < -2.4999999999999999e59Initial program 53.8%
Taylor expanded in A around -inf
lower-*.f64N/A
lower-/.f6425.8%
Applied rewrites25.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites25.8%
if -2.4999999999999999e59 < A < 9.2000000000000003e172Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
Taylor expanded in A around 0
Applied rewrites38.9%
if 9.2000000000000003e172 < A Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
Taylor expanded in A around inf
lower-*.f64N/A
lower-/.f6423.6%
Applied rewrites23.6%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6423.6%
Applied rewrites23.6%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= A -2.5e+59)
(* (/ (atan (* (/ (fabs B) A) 0.5)) PI) 180.0)
(if (<= A 9.2e+172)
(* 180.0 (/ (atan (- (/ C (fabs B)) 1.0)) PI))
(* (/ (atan (/ (- A) (fabs B))) PI) 180.0)))))double code(double A, double B, double C) {
double tmp;
if (A <= -2.5e+59) {
tmp = (atan(((fabs(B) / A) * 0.5)) / ((double) M_PI)) * 180.0;
} else if (A <= 9.2e+172) {
tmp = 180.0 * (atan(((C / fabs(B)) - 1.0)) / ((double) M_PI));
} else {
tmp = (atan((-A / fabs(B))) / ((double) M_PI)) * 180.0;
}
return copysign(1.0, B) * tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (A <= -2.5e+59) {
tmp = (Math.atan(((Math.abs(B) / A) * 0.5)) / Math.PI) * 180.0;
} else if (A <= 9.2e+172) {
tmp = 180.0 * (Math.atan(((C / Math.abs(B)) - 1.0)) / Math.PI);
} else {
tmp = (Math.atan((-A / Math.abs(B))) / Math.PI) * 180.0;
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if A <= -2.5e+59: tmp = (math.atan(((math.fabs(B) / A) * 0.5)) / math.pi) * 180.0 elif A <= 9.2e+172: tmp = 180.0 * (math.atan(((C / math.fabs(B)) - 1.0)) / math.pi) else: tmp = (math.atan((-A / math.fabs(B))) / math.pi) * 180.0 return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (A <= -2.5e+59) tmp = Float64(Float64(atan(Float64(Float64(abs(B) / A) * 0.5)) / pi) * 180.0); elseif (A <= 9.2e+172) tmp = Float64(180.0 * Float64(atan(Float64(Float64(C / abs(B)) - 1.0)) / pi)); else tmp = Float64(Float64(atan(Float64(Float64(-A) / abs(B))) / pi) * 180.0); end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (A <= -2.5e+59) tmp = (atan(((abs(B) / A) * 0.5)) / pi) * 180.0; elseif (A <= 9.2e+172) tmp = 180.0 * (atan(((C / abs(B)) - 1.0)) / pi); else tmp = (atan((-A / abs(B))) / pi) * 180.0; end tmp_2 = (sign(B) * abs(1.0)) * tmp; end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[A, -2.5e+59], N[(N[(N[ArcTan[N[(N[(N[Abs[B], $MachinePrecision] / A), $MachinePrecision] * 0.5), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision] * 180.0), $MachinePrecision], If[LessEqual[A, 9.2e+172], N[(180.0 * N[(N[ArcTan[N[(N[(C / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(N[(N[ArcTan[N[((-A) / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision] * 180.0), $MachinePrecision]]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;A \leq -2.5 \cdot 10^{+59}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{\left|B\right|}{A} \cdot 0.5\right)}{\pi} \cdot 180\\
\mathbf{elif}\;A \leq 9.2 \cdot 10^{+172}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{C}{\left|B\right|} - 1\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{-A}{\left|B\right|}\right)}{\pi} \cdot 180\\
\end{array}
if A < -2.4999999999999999e59Initial program 53.8%
Taylor expanded in A around -inf
lower-*.f64N/A
lower-/.f6425.8%
Applied rewrites25.8%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6425.8%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6425.8%
Applied rewrites25.8%
if -2.4999999999999999e59 < A < 9.2000000000000003e172Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
Taylor expanded in A around 0
Applied rewrites38.9%
if 9.2000000000000003e172 < A Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
Taylor expanded in A around inf
lower-*.f64N/A
lower-/.f6423.6%
Applied rewrites23.6%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6423.6%
Applied rewrites23.6%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= A -1.4e+176)
(* 180.0 (/ (atan 0.0) PI))
(if (<= A 9.2e+172)
(* 180.0 (/ (atan (- (/ C (fabs B)) 1.0)) PI))
(* (/ (atan (/ (- A) (fabs B))) PI) 180.0)))))double code(double A, double B, double C) {
double tmp;
if (A <= -1.4e+176) {
tmp = 180.0 * (atan(0.0) / ((double) M_PI));
} else if (A <= 9.2e+172) {
tmp = 180.0 * (atan(((C / fabs(B)) - 1.0)) / ((double) M_PI));
} else {
tmp = (atan((-A / fabs(B))) / ((double) M_PI)) * 180.0;
}
return copysign(1.0, B) * tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (A <= -1.4e+176) {
tmp = 180.0 * (Math.atan(0.0) / Math.PI);
} else if (A <= 9.2e+172) {
tmp = 180.0 * (Math.atan(((C / Math.abs(B)) - 1.0)) / Math.PI);
} else {
tmp = (Math.atan((-A / Math.abs(B))) / Math.PI) * 180.0;
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if A <= -1.4e+176: tmp = 180.0 * (math.atan(0.0) / math.pi) elif A <= 9.2e+172: tmp = 180.0 * (math.atan(((C / math.fabs(B)) - 1.0)) / math.pi) else: tmp = (math.atan((-A / math.fabs(B))) / math.pi) * 180.0 return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (A <= -1.4e+176) tmp = Float64(180.0 * Float64(atan(0.0) / pi)); elseif (A <= 9.2e+172) tmp = Float64(180.0 * Float64(atan(Float64(Float64(C / abs(B)) - 1.0)) / pi)); else tmp = Float64(Float64(atan(Float64(Float64(-A) / abs(B))) / pi) * 180.0); end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (A <= -1.4e+176) tmp = 180.0 * (atan(0.0) / pi); elseif (A <= 9.2e+172) tmp = 180.0 * (atan(((C / abs(B)) - 1.0)) / pi); else tmp = (atan((-A / abs(B))) / pi) * 180.0; end tmp_2 = (sign(B) * abs(1.0)) * tmp; end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[A, -1.4e+176], N[(180.0 * N[(N[ArcTan[0.0], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], If[LessEqual[A, 9.2e+172], N[(180.0 * N[(N[ArcTan[N[(N[(C / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(N[(N[ArcTan[N[((-A) / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision] * 180.0), $MachinePrecision]]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;A \leq -1.4 \cdot 10^{+176}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} 0}{\pi}\\
\mathbf{elif}\;A \leq 9.2 \cdot 10^{+172}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{C}{\left|B\right|} - 1\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{-A}{\left|B\right|}\right)}{\pi} \cdot 180\\
\end{array}
if A < -1.4000000000000001e176Initial program 53.8%
Taylor expanded in C around inf
lower-*.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f6412.6%
Applied rewrites12.6%
Taylor expanded in A around 0
Applied rewrites12.6%
if -1.4000000000000001e176 < A < 9.2000000000000003e172Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
Taylor expanded in A around 0
Applied rewrites38.9%
if 9.2000000000000003e172 < A Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
Taylor expanded in A around inf
lower-*.f64N/A
lower-/.f6423.6%
Applied rewrites23.6%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6423.6%
Applied rewrites23.6%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= (fabs B) 1.26e+94)
(* 180.0 (/ (atan (/ (- C A) (fabs B))) PI))
(* 180.0 (/ (atan -1.0) PI)))))double code(double A, double B, double C) {
double tmp;
if (fabs(B) <= 1.26e+94) {
tmp = 180.0 * (atan(((C - A) / fabs(B))) / ((double) M_PI));
} else {
tmp = 180.0 * (atan(-1.0) / ((double) M_PI));
}
return copysign(1.0, B) * tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (Math.abs(B) <= 1.26e+94) {
tmp = 180.0 * (Math.atan(((C - A) / Math.abs(B))) / Math.PI);
} else {
tmp = 180.0 * (Math.atan(-1.0) / Math.PI);
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if math.fabs(B) <= 1.26e+94: tmp = 180.0 * (math.atan(((C - A) / math.fabs(B))) / math.pi) else: tmp = 180.0 * (math.atan(-1.0) / math.pi) return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (abs(B) <= 1.26e+94) tmp = Float64(180.0 * Float64(atan(Float64(Float64(C - A) / abs(B))) / pi)); else tmp = Float64(180.0 * Float64(atan(-1.0) / pi)); end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (abs(B) <= 1.26e+94) tmp = 180.0 * (atan(((C - A) / abs(B))) / pi); else tmp = 180.0 * (atan(-1.0) / pi); end tmp_2 = (sign(B) * abs(1.0)) * tmp; end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[N[Abs[B], $MachinePrecision], 1.26e+94], N[(180.0 * N[(N[ArcTan[N[(N[(C - A), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(180.0 * N[(N[ArcTan[-1.0], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;\left|B\right| \leq 1.26 \cdot 10^{+94}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{C - A}{\left|B\right|}\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} -1}{\pi}\\
\end{array}
if B < 1.25999999999999997e94Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
Taylor expanded in B around 0
lower-/.f64N/A
lower--.f6435.3%
Applied rewrites35.3%
if 1.25999999999999997e94 < B Initial program 53.8%
Taylor expanded in B around inf
Applied rewrites21.2%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= C -6e+47)
(* 180.0 (/ (atan (/ C (fabs B))) PI))
(if (<= C 1.85e+108)
(* 180.0 (/ (atan -1.0) PI))
(* 180.0 (/ (atan 0.0) PI))))))double code(double A, double B, double C) {
double tmp;
if (C <= -6e+47) {
tmp = 180.0 * (atan((C / fabs(B))) / ((double) M_PI));
} else if (C <= 1.85e+108) {
tmp = 180.0 * (atan(-1.0) / ((double) M_PI));
} else {
tmp = 180.0 * (atan(0.0) / ((double) M_PI));
}
return copysign(1.0, B) * tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (C <= -6e+47) {
tmp = 180.0 * (Math.atan((C / Math.abs(B))) / Math.PI);
} else if (C <= 1.85e+108) {
tmp = 180.0 * (Math.atan(-1.0) / Math.PI);
} else {
tmp = 180.0 * (Math.atan(0.0) / Math.PI);
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if C <= -6e+47: tmp = 180.0 * (math.atan((C / math.fabs(B))) / math.pi) elif C <= 1.85e+108: tmp = 180.0 * (math.atan(-1.0) / math.pi) else: tmp = 180.0 * (math.atan(0.0) / math.pi) return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (C <= -6e+47) tmp = Float64(180.0 * Float64(atan(Float64(C / abs(B))) / pi)); elseif (C <= 1.85e+108) tmp = Float64(180.0 * Float64(atan(-1.0) / pi)); else tmp = Float64(180.0 * Float64(atan(0.0) / pi)); end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (C <= -6e+47) tmp = 180.0 * (atan((C / abs(B))) / pi); elseif (C <= 1.85e+108) tmp = 180.0 * (atan(-1.0) / pi); else tmp = 180.0 * (atan(0.0) / pi); end tmp_2 = (sign(B) * abs(1.0)) * tmp; end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[C, -6e+47], N[(180.0 * N[(N[ArcTan[N[(C / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], If[LessEqual[C, 1.85e+108], N[(180.0 * N[(N[ArcTan[-1.0], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(180.0 * N[(N[ArcTan[0.0], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;C \leq -6 \cdot 10^{+47}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{C}{\left|B\right|}\right)}{\pi}\\
\mathbf{elif}\;C \leq 1.85 \cdot 10^{+108}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} -1}{\pi}\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} 0}{\pi}\\
\end{array}
if C < -6.0000000000000003e47Initial program 53.8%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6449.7%
Applied rewrites49.7%
Taylor expanded in A around inf
lower-*.f64N/A
lower-/.f6423.6%
Applied rewrites23.6%
Taylor expanded in C around inf
lower-/.f6423.1%
Applied rewrites23.1%
if -6.0000000000000003e47 < C < 1.8499999999999999e108Initial program 53.8%
Taylor expanded in B around inf
Applied rewrites21.2%
if 1.8499999999999999e108 < C Initial program 53.8%
Taylor expanded in C around inf
lower-*.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f6412.6%
Applied rewrites12.6%
Taylor expanded in A around 0
Applied rewrites12.6%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= (fabs B) 2e-127)
(* 180.0 (/ (atan 0.0) PI))
(* 180.0 (/ (atan -1.0) PI)))))double code(double A, double B, double C) {
double tmp;
if (fabs(B) <= 2e-127) {
tmp = 180.0 * (atan(0.0) / ((double) M_PI));
} else {
tmp = 180.0 * (atan(-1.0) / ((double) M_PI));
}
return copysign(1.0, B) * tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (Math.abs(B) <= 2e-127) {
tmp = 180.0 * (Math.atan(0.0) / Math.PI);
} else {
tmp = 180.0 * (Math.atan(-1.0) / Math.PI);
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if math.fabs(B) <= 2e-127: tmp = 180.0 * (math.atan(0.0) / math.pi) else: tmp = 180.0 * (math.atan(-1.0) / math.pi) return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (abs(B) <= 2e-127) tmp = Float64(180.0 * Float64(atan(0.0) / pi)); else tmp = Float64(180.0 * Float64(atan(-1.0) / pi)); end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (abs(B) <= 2e-127) tmp = 180.0 * (atan(0.0) / pi); else tmp = 180.0 * (atan(-1.0) / pi); end tmp_2 = (sign(B) * abs(1.0)) * tmp; end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * If[LessEqual[N[Abs[B], $MachinePrecision], 2e-127], N[(180.0 * N[(N[ArcTan[0.0], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(180.0 * N[(N[ArcTan[-1.0], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;\left|B\right| \leq 2 \cdot 10^{-127}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} 0}{\pi}\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} -1}{\pi}\\
\end{array}
if B < 2.0000000000000001e-127Initial program 53.8%
Taylor expanded in C around inf
lower-*.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f6412.6%
Applied rewrites12.6%
Taylor expanded in A around 0
Applied rewrites12.6%
if 2.0000000000000001e-127 < B Initial program 53.8%
Taylor expanded in B around inf
Applied rewrites21.2%
(FPCore (A B C) :precision binary64 (* (copysign 1.0 B) (* 180.0 (/ (atan -1.0) PI))))
double code(double A, double B, double C) {
return copysign(1.0, B) * (180.0 * (atan(-1.0) / ((double) M_PI)));
}
public static double code(double A, double B, double C) {
return Math.copySign(1.0, B) * (180.0 * (Math.atan(-1.0) / Math.PI));
}
def code(A, B, C): return math.copysign(1.0, B) * (180.0 * (math.atan(-1.0) / math.pi))
function code(A, B, C) return Float64(copysign(1.0, B) * Float64(180.0 * Float64(atan(-1.0) / pi))) end
function tmp = code(A, B, C) tmp = (sign(B) * abs(1.0)) * (180.0 * (atan(-1.0) / pi)); end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * N[(180.0 * N[(N[ArcTan[-1.0], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \left(180 \cdot \frac{\tan^{-1} -1}{\pi}\right)
Initial program 53.8%
Taylor expanded in B around inf
Applied rewrites21.2%
herbie shell --seed 2025183
(FPCore (A B C)
:name "ABCF->ab-angle angle"
:precision binary64
(* 180.0 (/ (atan (* (/ 1.0 B) (- (- C A) (sqrt (+ (pow (- A C) 2.0) (pow B 2.0)))))) PI)))