
(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 8 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
(let* ((t_0 (/ 1.0 (fabs B))))
(*
(copysign 1.0 B)
(if (<=
(*
180.0
(/
(atan
(*
t_0
(-
(- C A)
(sqrt (+ (pow (- A C) 2.0) (pow (fabs B) 2.0))))))
PI))
-40.0)
(*
180.0
(/ (atan (* t_0 (- (- C A) (hypot (- C A) (fabs B))))) PI))
(*
(atan (fma (/ (fabs B) C) -0.5 (/ 0.0 (fabs B))))
(* (/ 1.0 PI) 180.0))))))double code(double A, double B, double C) {
double t_0 = 1.0 / fabs(B);
double tmp;
if ((180.0 * (atan((t_0 * ((C - A) - sqrt((pow((A - C), 2.0) + pow(fabs(B), 2.0)))))) / ((double) M_PI))) <= -40.0) {
tmp = 180.0 * (atan((t_0 * ((C - A) - hypot((C - A), fabs(B))))) / ((double) M_PI));
} else {
tmp = atan(fma((fabs(B) / C), -0.5, (0.0 / fabs(B)))) * ((1.0 / ((double) M_PI)) * 180.0);
}
return copysign(1.0, B) * tmp;
}
function code(A, B, C) t_0 = Float64(1.0 / abs(B)) tmp = 0.0 if (Float64(180.0 * Float64(atan(Float64(t_0 * Float64(Float64(C - A) - sqrt(Float64((Float64(A - C) ^ 2.0) + (abs(B) ^ 2.0)))))) / pi)) <= -40.0) tmp = Float64(180.0 * Float64(atan(Float64(t_0 * Float64(Float64(C - A) - hypot(Float64(C - A), abs(B))))) / pi)); else tmp = Float64(atan(fma(Float64(abs(B) / C), -0.5, Float64(0.0 / abs(B)))) * Float64(Float64(1.0 / pi) * 180.0)); end return Float64(copysign(1.0, B) * tmp) end
code[A_, B_, C_] := Block[{t$95$0 = N[(1.0 / N[Abs[B], $MachinePrecision]), $MachinePrecision]}, 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[(t$95$0 * 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], -40.0], N[(180.0 * N[(N[ArcTan[N[(t$95$0 * N[(N[(C - A), $MachinePrecision] - N[Sqrt[N[(C - A), $MachinePrecision] ^ 2 + N[Abs[B], $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(N[ArcTan[N[(N[(N[Abs[B], $MachinePrecision] / C), $MachinePrecision] * -0.5 + N[(0.0 / N[Abs[B], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(N[(1.0 / Pi), $MachinePrecision] * 180.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
t_0 := \frac{1}{\left|B\right|}\\
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;180 \cdot \frac{\tan^{-1} \left(t\_0 \cdot \left(\left(C - A\right) - \sqrt{{\left(A - C\right)}^{2} + {\left(\left|B\right|\right)}^{2}}\right)\right)}{\pi} \leq -40:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(t\_0 \cdot \left(\left(C - A\right) - \mathsf{hypot}\left(C - A, \left|B\right|\right)\right)\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\mathsf{fma}\left(\frac{\left|B\right|}{C}, -0.5, \frac{0}{\left|B\right|}\right)\right) \cdot \left(\frac{1}{\pi} \cdot 180\right)\\
\end{array}
\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))) < -40Initial program 53.2%
lift-sqrt.f64N/A
sqrt-fabs-revN/A
lift-sqrt.f64N/A
rem-sqrt-square-revN/A
lift-sqrt.f64N/A
lift-sqrt.f64N/A
rem-square-sqrtN/A
lift-+.f64N/A
lift-pow.f64N/A
unpow2N/A
lift-pow.f64N/A
unpow2N/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
Applied rewrites78.2%
if -40 < (*.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.2%
Taylor expanded in C around inf
lower-fma.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6425.6%
Applied rewrites25.6%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
mult-flipN/A
associate-*l*N/A
lower-*.f64N/A
Applied rewrites25.6%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= C 2.5e+64)
(/
(*
(atan (/ (* (- (/ (- C A) (fabs B)) 1.0) (fabs B)) (fabs B)))
180.0)
PI)
(* (atan (/ (* -0.5 (fabs B)) C)) (* (/ 1.0 PI) 180.0)))))double code(double A, double B, double C) {
double tmp;
if (C <= 2.5e+64) {
tmp = (atan((((((C - A) / fabs(B)) - 1.0) * fabs(B)) / fabs(B))) * 180.0) / ((double) M_PI);
} else {
tmp = atan(((-0.5 * fabs(B)) / C)) * ((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 (C <= 2.5e+64) {
tmp = (Math.atan((((((C - A) / Math.abs(B)) - 1.0) * Math.abs(B)) / Math.abs(B))) * 180.0) / Math.PI;
} else {
tmp = Math.atan(((-0.5 * Math.abs(B)) / C)) * ((1.0 / Math.PI) * 180.0);
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if C <= 2.5e+64: tmp = (math.atan((((((C - A) / math.fabs(B)) - 1.0) * math.fabs(B)) / math.fabs(B))) * 180.0) / math.pi else: tmp = math.atan(((-0.5 * math.fabs(B)) / C)) * ((1.0 / math.pi) * 180.0) return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (C <= 2.5e+64) tmp = Float64(Float64(atan(Float64(Float64(Float64(Float64(Float64(C - A) / abs(B)) - 1.0) * abs(B)) / abs(B))) * 180.0) / pi); else tmp = Float64(atan(Float64(Float64(-0.5 * abs(B)) / C)) * Float64(Float64(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 (C <= 2.5e+64) tmp = (atan((((((C - A) / abs(B)) - 1.0) * abs(B)) / abs(B))) * 180.0) / pi; else tmp = atan(((-0.5 * abs(B)) / C)) * ((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[C, 2.5e+64], N[(N[(N[ArcTan[N[(N[(N[(N[(N[(C - A), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision] * N[Abs[B], $MachinePrecision]), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision], N[(N[ArcTan[N[(N[(-0.5 * N[Abs[B], $MachinePrecision]), $MachinePrecision] / C), $MachinePrecision]], $MachinePrecision] * N[(N[(1.0 / Pi), $MachinePrecision] * 180.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;C \leq 2.5 \cdot 10^{+64}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{\left(\frac{C - A}{\left|B\right|} - 1\right) \cdot \left|B\right|}{\left|B\right|}\right) \cdot 180}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{-0.5 \cdot \left|B\right|}{C}\right) \cdot \left(\frac{1}{\pi} \cdot 180\right)\\
\end{array}
if C < 2.5e64Initial program 53.2%
Taylor expanded in B around inf
lower-*.f64N/A
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6448.9%
Applied rewrites48.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites50.2%
if 2.5e64 < C Initial program 53.2%
Taylor expanded in C around inf
lower-fma.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6425.6%
Applied rewrites25.6%
lift-fma.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
mult-flipN/A
lower-fma.f64N/A
lower-*.f64N/A
lower-/.f64N/A
mul-1-negN/A
lift-/.f64N/A
distribute-neg-fracN/A
lower-/.f64N/A
lift-+.f64N/A
lift-*.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-eval25.6%
Applied rewrites25.6%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
mult-flipN/A
associate-*l*N/A
lower-*.f64N/A
Applied rewrites25.6%
(FPCore (A B C) :precision binary64 (* (copysign 1.0 B) (if (<= C 2.5e+64) (/ (* (atan (- (/ (- C A) (fabs B)) 1.0)) 180.0) PI) (* (atan (/ (* -0.5 (fabs B)) C)) (* (/ 1.0 PI) 180.0)))))
double code(double A, double B, double C) {
double tmp;
if (C <= 2.5e+64) {
tmp = (atan((((C - A) / fabs(B)) - 1.0)) * 180.0) / ((double) M_PI);
} else {
tmp = atan(((-0.5 * fabs(B)) / C)) * ((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 (C <= 2.5e+64) {
tmp = (Math.atan((((C - A) / Math.abs(B)) - 1.0)) * 180.0) / Math.PI;
} else {
tmp = Math.atan(((-0.5 * Math.abs(B)) / C)) * ((1.0 / Math.PI) * 180.0);
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if C <= 2.5e+64: tmp = (math.atan((((C - A) / math.fabs(B)) - 1.0)) * 180.0) / math.pi else: tmp = math.atan(((-0.5 * math.fabs(B)) / C)) * ((1.0 / math.pi) * 180.0) return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (C <= 2.5e+64) tmp = Float64(Float64(atan(Float64(Float64(Float64(C - A) / abs(B)) - 1.0)) * 180.0) / pi); else tmp = Float64(atan(Float64(Float64(-0.5 * abs(B)) / C)) * Float64(Float64(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 (C <= 2.5e+64) tmp = (atan((((C - A) / abs(B)) - 1.0)) * 180.0) / pi; else tmp = atan(((-0.5 * abs(B)) / C)) * ((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[C, 2.5e+64], N[(N[(N[ArcTan[N[(N[(N[(C - A), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision], N[(N[ArcTan[N[(N[(-0.5 * N[Abs[B], $MachinePrecision]), $MachinePrecision] / C), $MachinePrecision]], $MachinePrecision] * N[(N[(1.0 / Pi), $MachinePrecision] * 180.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;C \leq 2.5 \cdot 10^{+64}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{C - A}{\left|B\right|} - 1\right) \cdot 180}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{-0.5 \cdot \left|B\right|}{C}\right) \cdot \left(\frac{1}{\pi} \cdot 180\right)\\
\end{array}
if C < 2.5e64Initial program 53.2%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6448.9%
Applied rewrites48.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites50.2%
if 2.5e64 < C Initial program 53.2%
Taylor expanded in C around inf
lower-fma.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6425.6%
Applied rewrites25.6%
lift-fma.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
mult-flipN/A
lower-fma.f64N/A
lower-*.f64N/A
lower-/.f64N/A
mul-1-negN/A
lift-/.f64N/A
distribute-neg-fracN/A
lower-/.f64N/A
lift-+.f64N/A
lift-*.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-eval25.6%
Applied rewrites25.6%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
mult-flipN/A
associate-*l*N/A
lower-*.f64N/A
Applied rewrites25.6%
(FPCore (A B C) :precision binary64 (* (copysign 1.0 B) (if (<= C 2.5e+64) (/ (* (atan (- (/ (- C A) (fabs B)) 1.0)) 180.0) PI) (/ (* (atan (/ (* -0.5 (fabs B)) C)) 180.0) PI))))
double code(double A, double B, double C) {
double tmp;
if (C <= 2.5e+64) {
tmp = (atan((((C - A) / fabs(B)) - 1.0)) * 180.0) / ((double) M_PI);
} else {
tmp = (atan(((-0.5 * fabs(B)) / C)) * 180.0) / ((double) M_PI);
}
return copysign(1.0, B) * tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (C <= 2.5e+64) {
tmp = (Math.atan((((C - A) / Math.abs(B)) - 1.0)) * 180.0) / Math.PI;
} else {
tmp = (Math.atan(((-0.5 * Math.abs(B)) / C)) * 180.0) / Math.PI;
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if C <= 2.5e+64: tmp = (math.atan((((C - A) / math.fabs(B)) - 1.0)) * 180.0) / math.pi else: tmp = (math.atan(((-0.5 * math.fabs(B)) / C)) * 180.0) / math.pi return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (C <= 2.5e+64) tmp = Float64(Float64(atan(Float64(Float64(Float64(C - A) / abs(B)) - 1.0)) * 180.0) / pi); else tmp = Float64(Float64(atan(Float64(Float64(-0.5 * abs(B)) / C)) * 180.0) / pi); end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (C <= 2.5e+64) tmp = (atan((((C - A) / abs(B)) - 1.0)) * 180.0) / pi; else tmp = (atan(((-0.5 * abs(B)) / C)) * 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[C, 2.5e+64], N[(N[(N[ArcTan[N[(N[(N[(C - A), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision], N[(N[(N[ArcTan[N[(N[(-0.5 * N[Abs[B], $MachinePrecision]), $MachinePrecision] / C), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;C \leq 2.5 \cdot 10^{+64}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{C - A}{\left|B\right|} - 1\right) \cdot 180}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{-0.5 \cdot \left|B\right|}{C}\right) \cdot 180}{\pi}\\
\end{array}
if C < 2.5e64Initial program 53.2%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6448.9%
Applied rewrites48.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites50.2%
if 2.5e64 < C Initial program 53.2%
Taylor expanded in C around inf
lower-fma.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6425.6%
Applied rewrites25.6%
lift-fma.f64N/A
+-commutativeN/A
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
mult-flipN/A
lower-fma.f64N/A
lower-*.f64N/A
lower-/.f64N/A
mul-1-negN/A
lift-/.f64N/A
distribute-neg-fracN/A
lower-/.f64N/A
lift-+.f64N/A
lift-*.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-eval25.6%
Applied rewrites25.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites25.6%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= A -6.2e+148)
(* 180.0 (/ (atan (* 0.5 (/ (fabs B) A))) PI))
(if (<= A 2e-105)
-45.0
(* 180.0 (/ (atan (/ (- C A) (fabs B))) PI))))))double code(double A, double B, double C) {
double tmp;
if (A <= -6.2e+148) {
tmp = 180.0 * (atan((0.5 * (fabs(B) / A))) / ((double) M_PI));
} else if (A <= 2e-105) {
tmp = -45.0;
} else {
tmp = 180.0 * (atan(((C - A) / fabs(B))) / ((double) M_PI));
}
return copysign(1.0, B) * tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (A <= -6.2e+148) {
tmp = 180.0 * (Math.atan((0.5 * (Math.abs(B) / A))) / Math.PI);
} else if (A <= 2e-105) {
tmp = -45.0;
} else {
tmp = 180.0 * (Math.atan(((C - A) / Math.abs(B))) / Math.PI);
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if A <= -6.2e+148: tmp = 180.0 * (math.atan((0.5 * (math.fabs(B) / A))) / math.pi) elif A <= 2e-105: tmp = -45.0 else: tmp = 180.0 * (math.atan(((C - A) / math.fabs(B))) / math.pi) return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (A <= -6.2e+148) tmp = Float64(180.0 * Float64(atan(Float64(0.5 * Float64(abs(B) / A))) / pi)); elseif (A <= 2e-105) tmp = -45.0; else tmp = Float64(180.0 * Float64(atan(Float64(Float64(C - A) / abs(B))) / pi)); end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (A <= -6.2e+148) tmp = 180.0 * (atan((0.5 * (abs(B) / A))) / pi); elseif (A <= 2e-105) tmp = -45.0; else tmp = 180.0 * (atan(((C - A) / abs(B))) / 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, -6.2e+148], N[(180.0 * N[(N[ArcTan[N[(0.5 * N[(N[Abs[B], $MachinePrecision] / A), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], If[LessEqual[A, 2e-105], -45.0, N[(180.0 * N[(N[ArcTan[N[(N[(C - A), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;A \leq -6.2 \cdot 10^{+148}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(0.5 \cdot \frac{\left|B\right|}{A}\right)}{\pi}\\
\mathbf{elif}\;A \leq 2 \cdot 10^{-105}:\\
\;\;\;\;-45\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{C - A}{\left|B\right|}\right)}{\pi}\\
\end{array}
if A < -6.1999999999999995e148Initial program 53.2%
Taylor expanded in A around -inf
lower-*.f64N/A
lower-/.f6426.1%
Applied rewrites26.1%
if -6.1999999999999995e148 < A < 1.9999999999999999e-105Initial program 53.2%
Taylor expanded in B around inf
Applied rewrites21.4%
Evaluated real constant21.4%
if 1.9999999999999999e-105 < A Initial program 53.2%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6448.9%
Applied rewrites48.9%
Taylor expanded in B around 0
lower-/.f64N/A
lower--.f6434.3%
Applied rewrites34.3%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= C -1.6e+194)
(/ (* (atan (/ (+ C C) (fabs B))) 180.0) PI)
(if (<= C 6e+64)
-45.0
(* 180.0 (/ (atan (* -0.5 (/ (fabs B) C))) PI))))))double code(double A, double B, double C) {
double tmp;
if (C <= -1.6e+194) {
tmp = (atan(((C + C) / fabs(B))) * 180.0) / ((double) M_PI);
} else if (C <= 6e+64) {
tmp = -45.0;
} else {
tmp = 180.0 * (atan((-0.5 * (fabs(B) / C))) / ((double) M_PI));
}
return copysign(1.0, B) * tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (C <= -1.6e+194) {
tmp = (Math.atan(((C + C) / Math.abs(B))) * 180.0) / Math.PI;
} else if (C <= 6e+64) {
tmp = -45.0;
} else {
tmp = 180.0 * (Math.atan((-0.5 * (Math.abs(B) / C))) / Math.PI);
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if C <= -1.6e+194: tmp = (math.atan(((C + C) / math.fabs(B))) * 180.0) / math.pi elif C <= 6e+64: tmp = -45.0 else: tmp = 180.0 * (math.atan((-0.5 * (math.fabs(B) / C))) / math.pi) return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (C <= -1.6e+194) tmp = Float64(Float64(atan(Float64(Float64(C + C) / abs(B))) * 180.0) / pi); elseif (C <= 6e+64) tmp = -45.0; else tmp = Float64(180.0 * Float64(atan(Float64(-0.5 * Float64(abs(B) / C))) / pi)); end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (C <= -1.6e+194) tmp = (atan(((C + C) / abs(B))) * 180.0) / pi; elseif (C <= 6e+64) tmp = -45.0; else tmp = 180.0 * (atan((-0.5 * (abs(B) / C))) / 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, -1.6e+194], N[(N[(N[ArcTan[N[(N[(C + C), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision], If[LessEqual[C, 6e+64], -45.0, N[(180.0 * N[(N[ArcTan[N[(-0.5 * N[(N[Abs[B], $MachinePrecision] / C), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;C \leq -1.6 \cdot 10^{+194}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{C + C}{\left|B\right|}\right) \cdot 180}{\pi}\\
\mathbf{elif}\;C \leq 6 \cdot 10^{+64}:\\
\;\;\;\;-45\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(-0.5 \cdot \frac{\left|B\right|}{C}\right)}{\pi}\\
\end{array}
if C < -1.6000000000000001e194Initial program 53.2%
Taylor expanded in C around -inf
lower-*.f64N/A
lower-/.f6422.6%
Applied rewrites22.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites22.6%
if -1.6000000000000001e194 < C < 6.0000000000000004e64Initial program 53.2%
Taylor expanded in B around inf
Applied rewrites21.4%
Evaluated real constant21.4%
if 6.0000000000000004e64 < C Initial program 53.2%
Taylor expanded in C around inf
lower-fma.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6425.6%
Applied rewrites25.6%
Taylor expanded in A around 0
lower-*.f64N/A
lower-/.f6425.6%
Applied rewrites25.6%
(FPCore (A B C) :precision binary64 (* (copysign 1.0 B) (if (<= (fabs B) 6e+76) (* 180.0 (/ (atan (/ (- C A) (fabs B))) PI)) -45.0)))
double code(double A, double B, double C) {
double tmp;
if (fabs(B) <= 6e+76) {
tmp = 180.0 * (atan(((C - A) / fabs(B))) / ((double) M_PI));
} else {
tmp = -45.0;
}
return copysign(1.0, B) * tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (Math.abs(B) <= 6e+76) {
tmp = 180.0 * (Math.atan(((C - A) / Math.abs(B))) / Math.PI);
} else {
tmp = -45.0;
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if math.fabs(B) <= 6e+76: tmp = 180.0 * (math.atan(((C - A) / math.fabs(B))) / math.pi) else: tmp = -45.0 return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (abs(B) <= 6e+76) tmp = Float64(180.0 * Float64(atan(Float64(Float64(C - A) / abs(B))) / pi)); else tmp = -45.0; end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (abs(B) <= 6e+76) tmp = 180.0 * (atan(((C - A) / abs(B))) / pi); else tmp = -45.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[N[Abs[B], $MachinePrecision], 6e+76], N[(180.0 * N[(N[ArcTan[N[(N[(C - A), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], -45.0]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;\left|B\right| \leq 6 \cdot 10^{+76}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{C - A}{\left|B\right|}\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;-45\\
\end{array}
if B < 5.9999999999999996e76Initial program 53.2%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6448.9%
Applied rewrites48.9%
Taylor expanded in B around 0
lower-/.f64N/A
lower--.f6434.3%
Applied rewrites34.3%
if 5.9999999999999996e76 < B Initial program 53.2%
Taylor expanded in B around inf
Applied rewrites21.4%
Evaluated real constant21.4%
(FPCore (A B C) :precision binary64 (* (copysign 1.0 B) -45.0))
double code(double A, double B, double C) {
return copysign(1.0, B) * -45.0;
}
public static double code(double A, double B, double C) {
return Math.copySign(1.0, B) * -45.0;
}
def code(A, B, C): return math.copysign(1.0, B) * -45.0
function code(A, B, C) return Float64(copysign(1.0, B) * -45.0) end
function tmp = code(A, B, C) tmp = (sign(B) * abs(1.0)) * -45.0; end
code[A_, B_, C_] := N[(N[With[{TMP1 = Abs[1.0], TMP2 = Sign[B]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision] * -45.0), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot -45
Initial program 53.2%
Taylor expanded in B around inf
Applied rewrites21.4%
Evaluated real constant21.4%
herbie shell --seed 2025212
(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)))