
(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 9 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 (if (<= C 6.2e+163) (* 180.0 (/ (atan (* (/ 1.0 B) (- (- C A) (hypot (- C A) B)))) PI)) (* (atan (* (/ B C) -0.5)) (/ 180.0 PI))))
double code(double A, double B, double C) {
double tmp;
if (C <= 6.2e+163) {
tmp = 180.0 * (atan(((1.0 / B) * ((C - A) - hypot((C - A), B)))) / ((double) M_PI));
} else {
tmp = atan(((B / C) * -0.5)) * (180.0 / ((double) M_PI));
}
return tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (C <= 6.2e+163) {
tmp = 180.0 * (Math.atan(((1.0 / B) * ((C - A) - Math.hypot((C - A), B)))) / Math.PI);
} else {
tmp = Math.atan(((B / C) * -0.5)) * (180.0 / Math.PI);
}
return tmp;
}
def code(A, B, C): tmp = 0 if C <= 6.2e+163: tmp = 180.0 * (math.atan(((1.0 / B) * ((C - A) - math.hypot((C - A), B)))) / math.pi) else: tmp = math.atan(((B / C) * -0.5)) * (180.0 / math.pi) return tmp
function code(A, B, C) tmp = 0.0 if (C <= 6.2e+163) tmp = Float64(180.0 * Float64(atan(Float64(Float64(1.0 / B) * Float64(Float64(C - A) - hypot(Float64(C - A), B)))) / pi)); else tmp = Float64(atan(Float64(Float64(B / C) * -0.5)) * Float64(180.0 / pi)); end return tmp end
function tmp_2 = code(A, B, C) tmp = 0.0; if (C <= 6.2e+163) tmp = 180.0 * (atan(((1.0 / B) * ((C - A) - hypot((C - A), B)))) / pi); else tmp = atan(((B / C) * -0.5)) * (180.0 / pi); end tmp_2 = tmp; end
code[A_, B_, C_] := If[LessEqual[C, 6.2e+163], N[(180.0 * N[(N[ArcTan[N[(N[(1.0 / B), $MachinePrecision] * N[(N[(C - A), $MachinePrecision] - N[Sqrt[N[(C - A), $MachinePrecision] ^ 2 + B ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(N[ArcTan[N[(N[(B / C), $MachinePrecision] * -0.5), $MachinePrecision]], $MachinePrecision] * N[(180.0 / Pi), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\mathbf{if}\;C \leq 6.2 \cdot 10^{+163}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{1}{B} \cdot \left(\left(C - A\right) - \mathsf{hypot}\left(C - A, B\right)\right)\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{B}{C} \cdot -0.5\right) \cdot \frac{180}{\pi}\\
\end{array}
if C < 6.20000000000000057e163Initial program 53.4%
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 rewrites77.9%
if 6.20000000000000057e163 < C Initial program 53.4%
Taylor expanded in C around inf
lower-fma.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6426.2%
Applied rewrites26.2%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
mult-flipN/A
associate-*l*N/A
lower-*.f64N/A
Applied rewrites26.3%
lift-fma.f64N/A
lift-/.f64N/A
div0N/A
+-rgt-identityN/A
lower-*.f6426.3%
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites26.3%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= C 2.6e+118)
(/ (* (atan (/ (- C (+ A (fabs B))) (fabs B))) 180.0) PI)
(* (atan (* (/ (fabs B) C) -0.5)) (/ 180.0 PI)))))double code(double A, double B, double C) {
double tmp;
if (C <= 2.6e+118) {
tmp = (atan(((C - (A + fabs(B))) / fabs(B))) * 180.0) / ((double) M_PI);
} else {
tmp = atan(((fabs(B) / C) * -0.5)) * (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.6e+118) {
tmp = (Math.atan(((C - (A + Math.abs(B))) / Math.abs(B))) * 180.0) / Math.PI;
} else {
tmp = Math.atan(((Math.abs(B) / C) * -0.5)) * (180.0 / Math.PI);
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if C <= 2.6e+118: tmp = (math.atan(((C - (A + math.fabs(B))) / math.fabs(B))) * 180.0) / math.pi else: tmp = math.atan(((math.fabs(B) / C) * -0.5)) * (180.0 / math.pi) return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (C <= 2.6e+118) tmp = Float64(Float64(atan(Float64(Float64(C - Float64(A + abs(B))) / abs(B))) * 180.0) / pi); else tmp = Float64(atan(Float64(Float64(abs(B) / C) * -0.5)) * Float64(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.6e+118) tmp = (atan(((C - (A + abs(B))) / abs(B))) * 180.0) / pi; else tmp = atan(((abs(B) / C) * -0.5)) * (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.6e+118], N[(N[(N[ArcTan[N[(N[(C - N[(A + N[Abs[B], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] / Pi), $MachinePrecision], N[(N[ArcTan[N[(N[(N[Abs[B], $MachinePrecision] / C), $MachinePrecision] * -0.5), $MachinePrecision]], $MachinePrecision] * N[(180.0 / Pi), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;C \leq 2.6 \cdot 10^{+118}:\\
\;\;\;\;\frac{\tan^{-1} \left(\frac{C - \left(A + \left|B\right|\right)}{\left|B\right|}\right) \cdot 180}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{\left|B\right|}{C} \cdot -0.5\right) \cdot \frac{180}{\pi}\\
\end{array}
if C < 2.60000000000000016e118Initial program 53.4%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6450.2%
Applied rewrites50.2%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites51.1%
if 2.60000000000000016e118 < C Initial program 53.4%
Taylor expanded in C around inf
lower-fma.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6426.2%
Applied rewrites26.2%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
mult-flipN/A
associate-*l*N/A
lower-*.f64N/A
Applied rewrites26.3%
lift-fma.f64N/A
lift-/.f64N/A
div0N/A
+-rgt-identityN/A
lower-*.f6426.3%
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites26.3%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= C 2.6e+118)
(* (* (atan (/ (- C (+ A (fabs B))) (fabs B))) 180.0) 0.3183098861837907)
(* (atan (* (/ (fabs B) C) -0.5)) (/ 180.0 PI)))))double code(double A, double B, double C) {
double tmp;
if (C <= 2.6e+118) {
tmp = (atan(((C - (A + fabs(B))) / fabs(B))) * 180.0) * 0.3183098861837907;
} else {
tmp = atan(((fabs(B) / C) * -0.5)) * (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.6e+118) {
tmp = (Math.atan(((C - (A + Math.abs(B))) / Math.abs(B))) * 180.0) * 0.3183098861837907;
} else {
tmp = Math.atan(((Math.abs(B) / C) * -0.5)) * (180.0 / Math.PI);
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if C <= 2.6e+118: tmp = (math.atan(((C - (A + math.fabs(B))) / math.fabs(B))) * 180.0) * 0.3183098861837907 else: tmp = math.atan(((math.fabs(B) / C) * -0.5)) * (180.0 / math.pi) return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (C <= 2.6e+118) tmp = Float64(Float64(atan(Float64(Float64(C - Float64(A + abs(B))) / abs(B))) * 180.0) * 0.3183098861837907); else tmp = Float64(atan(Float64(Float64(abs(B) / C) * -0.5)) * Float64(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.6e+118) tmp = (atan(((C - (A + abs(B))) / abs(B))) * 180.0) * 0.3183098861837907; else tmp = atan(((abs(B) / C) * -0.5)) * (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.6e+118], N[(N[(N[ArcTan[N[(N[(C - N[(A + N[Abs[B], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Abs[B], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * 180.0), $MachinePrecision] * 0.3183098861837907), $MachinePrecision], N[(N[ArcTan[N[(N[(N[Abs[B], $MachinePrecision] / C), $MachinePrecision] * -0.5), $MachinePrecision]], $MachinePrecision] * N[(180.0 / Pi), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;C \leq 2.6 \cdot 10^{+118}:\\
\;\;\;\;\left(\tan^{-1} \left(\frac{C - \left(A + \left|B\right|\right)}{\left|B\right|}\right) \cdot 180\right) \cdot 0.3183098861837907\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{\left|B\right|}{C} \cdot -0.5\right) \cdot \frac{180}{\pi}\\
\end{array}
if C < 2.60000000000000016e118Initial program 53.4%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6450.2%
Applied rewrites50.2%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
mult-flipN/A
lift-/.f64N/A
lower-*.f64N/A
Applied rewrites51.1%
Evaluated real constant51.1%
if 2.60000000000000016e118 < C Initial program 53.4%
Taylor expanded in C around inf
lower-fma.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6426.2%
Applied rewrites26.2%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
mult-flipN/A
associate-*l*N/A
lower-*.f64N/A
Applied rewrites26.3%
lift-fma.f64N/A
lift-/.f64N/A
div0N/A
+-rgt-identityN/A
lower-*.f6426.3%
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites26.3%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= C 820000.0)
(* 180.0 (/ (atan (- (/ C (fabs B)) 1.0)) PI))
(* (atan (* (/ (fabs B) C) -0.5)) (/ 180.0 PI)))))double code(double A, double B, double C) {
double tmp;
if (C <= 820000.0) {
tmp = 180.0 * (atan(((C / fabs(B)) - 1.0)) / ((double) M_PI));
} else {
tmp = atan(((fabs(B) / C) * -0.5)) * (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 <= 820000.0) {
tmp = 180.0 * (Math.atan(((C / Math.abs(B)) - 1.0)) / Math.PI);
} else {
tmp = Math.atan(((Math.abs(B) / C) * -0.5)) * (180.0 / Math.PI);
}
return Math.copySign(1.0, B) * tmp;
}
def code(A, B, C): tmp = 0 if C <= 820000.0: tmp = 180.0 * (math.atan(((C / math.fabs(B)) - 1.0)) / math.pi) else: tmp = math.atan(((math.fabs(B) / C) * -0.5)) * (180.0 / math.pi) return math.copysign(1.0, B) * tmp
function code(A, B, C) tmp = 0.0 if (C <= 820000.0) tmp = Float64(180.0 * Float64(atan(Float64(Float64(C / abs(B)) - 1.0)) / pi)); else tmp = Float64(atan(Float64(Float64(abs(B) / C) * -0.5)) * Float64(180.0 / pi)); end return Float64(copysign(1.0, B) * tmp) end
function tmp_2 = code(A, B, C) tmp = 0.0; if (C <= 820000.0) tmp = 180.0 * (atan(((C / abs(B)) - 1.0)) / pi); else tmp = atan(((abs(B) / C) * -0.5)) * (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, 820000.0], N[(180.0 * N[(N[ArcTan[N[(N[(C / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(N[ArcTan[N[(N[(N[Abs[B], $MachinePrecision] / C), $MachinePrecision] * -0.5), $MachinePrecision]], $MachinePrecision] * N[(180.0 / Pi), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\mathsf{copysign}\left(1, B\right) \cdot \begin{array}{l}
\mathbf{if}\;C \leq 820000:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{C}{\left|B\right|} - 1\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{\left|B\right|}{C} \cdot -0.5\right) \cdot \frac{180}{\pi}\\
\end{array}
if C < 8.2e5Initial program 53.4%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6450.2%
Applied rewrites50.2%
Taylor expanded in A around 0
Applied rewrites39.4%
if 8.2e5 < C Initial program 53.4%
Taylor expanded in C around inf
lower-fma.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6426.2%
Applied rewrites26.2%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
mult-flipN/A
associate-*l*N/A
lower-*.f64N/A
Applied rewrites26.3%
lift-fma.f64N/A
lift-/.f64N/A
div0N/A
+-rgt-identityN/A
lower-*.f6426.3%
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f64N/A
Applied rewrites26.3%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= A -4.2e+161)
(* (atan (* (/ (fabs B) A) 0.5)) 57.29577951308232)
(if (<= A 3e+91)
(* 180.0 (/ (atan (- (/ C (fabs B)) 1.0)) PI))
(* 180.0 (/ (atan (/ (- C A) (fabs B))) PI))))))double code(double A, double B, double C) {
double tmp;
if (A <= -4.2e+161) {
tmp = atan(((fabs(B) / A) * 0.5)) * 57.29577951308232;
} else if (A <= 3e+91) {
tmp = 180.0 * (atan(((C / fabs(B)) - 1.0)) / ((double) M_PI));
} 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 <= -4.2e+161) {
tmp = Math.atan(((Math.abs(B) / A) * 0.5)) * 57.29577951308232;
} else if (A <= 3e+91) {
tmp = 180.0 * (Math.atan(((C / Math.abs(B)) - 1.0)) / Math.PI);
} 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 <= -4.2e+161: tmp = math.atan(((math.fabs(B) / A) * 0.5)) * 57.29577951308232 elif A <= 3e+91: tmp = 180.0 * (math.atan(((C / math.fabs(B)) - 1.0)) / math.pi) 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 <= -4.2e+161) tmp = Float64(atan(Float64(Float64(abs(B) / A) * 0.5)) * 57.29577951308232); elseif (A <= 3e+91) tmp = Float64(180.0 * Float64(atan(Float64(Float64(C / abs(B)) - 1.0)) / pi)); 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 <= -4.2e+161) tmp = atan(((abs(B) / A) * 0.5)) * 57.29577951308232; elseif (A <= 3e+91) tmp = 180.0 * (atan(((C / abs(B)) - 1.0)) / pi); 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, -4.2e+161], N[(N[ArcTan[N[(N[(N[Abs[B], $MachinePrecision] / A), $MachinePrecision] * 0.5), $MachinePrecision]], $MachinePrecision] * 57.29577951308232), $MachinePrecision], If[LessEqual[A, 3e+91], N[(180.0 * N[(N[ArcTan[N[(N[(C / N[Abs[B], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], 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 -4.2 \cdot 10^{+161}:\\
\;\;\;\;\tan^{-1} \left(\frac{\left|B\right|}{A} \cdot 0.5\right) \cdot 57.29577951308232\\
\mathbf{elif}\;A \leq 3 \cdot 10^{+91}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{C}{\left|B\right|} - 1\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{C - A}{\left|B\right|}\right)}{\pi}\\
\end{array}
if A < -4.2e161Initial program 53.4%
Taylor expanded in A around -inf
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.7%
Evaluated real constant25.7%
if -4.2e161 < A < 3.00000000000000006e91Initial program 53.4%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6450.2%
Applied rewrites50.2%
Taylor expanded in A around 0
Applied rewrites39.4%
if 3.00000000000000006e91 < A Initial program 53.4%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6450.2%
Applied rewrites50.2%
Taylor expanded in B around 0
lower-/.f64N/A
lower--.f6435.2%
Applied rewrites35.2%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= A -5.8e+81)
(* (atan (* (/ (fabs B) A) 0.5)) 57.29577951308232)
(if (<= A 4.5e-81)
(* 180.0 (/ (atan -1.0) PI))
(* 180.0 (/ (atan (/ (- C A) (fabs B))) PI))))))double code(double A, double B, double C) {
double tmp;
if (A <= -5.8e+81) {
tmp = atan(((fabs(B) / A) * 0.5)) * 57.29577951308232;
} else if (A <= 4.5e-81) {
tmp = 180.0 * (atan(-1.0) / ((double) M_PI));
} 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 <= -5.8e+81) {
tmp = Math.atan(((Math.abs(B) / A) * 0.5)) * 57.29577951308232;
} else if (A <= 4.5e-81) {
tmp = 180.0 * (Math.atan(-1.0) / Math.PI);
} 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 <= -5.8e+81: tmp = math.atan(((math.fabs(B) / A) * 0.5)) * 57.29577951308232 elif A <= 4.5e-81: tmp = 180.0 * (math.atan(-1.0) / math.pi) 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 <= -5.8e+81) tmp = Float64(atan(Float64(Float64(abs(B) / A) * 0.5)) * 57.29577951308232); elseif (A <= 4.5e-81) tmp = Float64(180.0 * Float64(atan(-1.0) / pi)); 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 <= -5.8e+81) tmp = atan(((abs(B) / A) * 0.5)) * 57.29577951308232; elseif (A <= 4.5e-81) tmp = 180.0 * (atan(-1.0) / pi); 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, -5.8e+81], N[(N[ArcTan[N[(N[(N[Abs[B], $MachinePrecision] / A), $MachinePrecision] * 0.5), $MachinePrecision]], $MachinePrecision] * 57.29577951308232), $MachinePrecision], If[LessEqual[A, 4.5e-81], N[(180.0 * N[(N[ArcTan[-1.0], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], 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 -5.8 \cdot 10^{+81}:\\
\;\;\;\;\tan^{-1} \left(\frac{\left|B\right|}{A} \cdot 0.5\right) \cdot 57.29577951308232\\
\mathbf{elif}\;A \leq 4.5 \cdot 10^{-81}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} -1}{\pi}\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{C - A}{\left|B\right|}\right)}{\pi}\\
\end{array}
if A < -5.7999999999999999e81Initial program 53.4%
Taylor expanded in A around -inf
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.7%
Evaluated real constant25.7%
if -5.7999999999999999e81 < A < 4.5e-81Initial program 53.4%
Taylor expanded in B around inf
Applied rewrites21.6%
if 4.5e-81 < A Initial program 53.4%
Taylor expanded in B around inf
lower--.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-/.f6450.2%
Applied rewrites50.2%
Taylor expanded in B around 0
lower-/.f64N/A
lower--.f6435.2%
Applied rewrites35.2%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= A -5.8e+81)
(* (atan (* (/ (fabs B) A) 0.5)) 57.29577951308232)
(* 180.0 (/ (atan -1.0) PI)))))double code(double A, double B, double C) {
double tmp;
if (A <= -5.8e+81) {
tmp = atan(((fabs(B) / A) * 0.5)) * 57.29577951308232;
} 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 (A <= -5.8e+81) {
tmp = Math.atan(((Math.abs(B) / A) * 0.5)) * 57.29577951308232;
} 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 A <= -5.8e+81: tmp = math.atan(((math.fabs(B) / A) * 0.5)) * 57.29577951308232 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 (A <= -5.8e+81) tmp = Float64(atan(Float64(Float64(abs(B) / A) * 0.5)) * 57.29577951308232); 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 (A <= -5.8e+81) tmp = atan(((abs(B) / A) * 0.5)) * 57.29577951308232; 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[A, -5.8e+81], N[(N[ArcTan[N[(N[(N[Abs[B], $MachinePrecision] / A), $MachinePrecision] * 0.5), $MachinePrecision]], $MachinePrecision] * 57.29577951308232), $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}\;A \leq -5.8 \cdot 10^{+81}:\\
\;\;\;\;\tan^{-1} \left(\frac{\left|B\right|}{A} \cdot 0.5\right) \cdot 57.29577951308232\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} -1}{\pi}\\
\end{array}
if A < -5.7999999999999999e81Initial program 53.4%
Taylor expanded in A around -inf
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.7%
Evaluated real constant25.7%
if -5.7999999999999999e81 < A Initial program 53.4%
Taylor expanded in B around inf
Applied rewrites21.6%
(FPCore (A B C)
:precision binary64
(*
(copysign 1.0 B)
(if (<= (fabs B) 1.15e-138)
(* 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) <= 1.15e-138) {
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) <= 1.15e-138) {
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) <= 1.15e-138: 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) <= 1.15e-138) 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) <= 1.15e-138) 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], 1.15e-138], 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 1.15 \cdot 10^{-138}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} 0}{\pi}\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} -1}{\pi}\\
\end{array}
if B < 1.14999999999999995e-138Initial program 53.4%
Taylor expanded in C around inf
lower-*.f64N/A
lower-/.f64N/A
lower-+.f64N/A
lower-*.f6413.0%
Applied rewrites13.0%
Taylor expanded in A around 0
Applied rewrites13.0%
if 1.14999999999999995e-138 < B Initial program 53.4%
Taylor expanded in B around inf
Applied rewrites21.6%
(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.4%
Taylor expanded in B around inf
Applied rewrites21.6%
herbie shell --seed 2025188
(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)))