
(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]
\begin{array}{l}
\\
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}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 7 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]
\begin{array}{l}
\\
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}
\end{array}
(FPCore (A B C) :precision binary64 (if (<= A -9.5e+99) (/ (* 180.0 (atan (/ (* (fma (/ C A) B B) -0.5) (- A)))) PI) (* 180.0 (/ (atan (/ (- (- C A) (hypot (- A C) B)) B)) PI))))
double code(double A, double B, double C) {
double tmp;
if (A <= -9.5e+99) {
tmp = (180.0 * atan(((fma((C / A), B, B) * -0.5) / -A))) / ((double) M_PI);
} else {
tmp = 180.0 * (atan((((C - A) - hypot((A - C), B)) / B)) / ((double) M_PI));
}
return tmp;
}
function code(A, B, C) tmp = 0.0 if (A <= -9.5e+99) tmp = Float64(Float64(180.0 * atan(Float64(Float64(fma(Float64(C / A), B, B) * -0.5) / Float64(-A)))) / pi); else tmp = Float64(180.0 * Float64(atan(Float64(Float64(Float64(C - A) - hypot(Float64(A - C), B)) / B)) / pi)); end return tmp end
code[A_, B_, C_] := If[LessEqual[A, -9.5e+99], N[(N[(180.0 * N[ArcTan[N[(N[(N[(N[(C / A), $MachinePrecision] * B + B), $MachinePrecision] * -0.5), $MachinePrecision] / (-A)), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / Pi), $MachinePrecision], N[(180.0 * N[(N[ArcTan[N[(N[(N[(C - A), $MachinePrecision] - N[Sqrt[N[(A - C), $MachinePrecision] ^ 2 + B ^ 2], $MachinePrecision]), $MachinePrecision] / B), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;A \leq -9.5 \cdot 10^{+99}:\\
\;\;\;\;\frac{180 \cdot \tan^{-1} \left(\frac{\mathsf{fma}\left(\frac{C}{A}, B, B\right) \cdot -0.5}{-A}\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{\left(C - A\right) - \mathsf{hypot}\left(A - C, B\right)}{B}\right)}{\pi}\\
\end{array}
\end{array}
if A < -9.49999999999999908e99Initial program 11.7%
Taylor expanded in A around -inf
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6487.0
Applied rewrites87.0%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
Applied rewrites87.2%
if -9.49999999999999908e99 < A Initial program 60.2%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
lift--.f64N/A
lift-sqrt.f64N/A
lift-+.f64N/A
lift--.f64N/A
lift-pow.f64N/A
lift-pow.f64N/A
associate-*l/N/A
lower-/.f64N/A
Applied rewrites84.8%
Final simplification85.2%
(FPCore (A B C) :precision binary64 (if (<= A -9.6e+99) (* 180.0 (/ (atan (/ (* -0.5 (+ B (* B (/ C A)))) (- A))) PI)) (* 180.0 (/ (atan (/ (- (- C A) (hypot (- A C) B)) B)) PI))))
double code(double A, double B, double C) {
double tmp;
if (A <= -9.6e+99) {
tmp = 180.0 * (atan(((-0.5 * (B + (B * (C / A)))) / -A)) / ((double) M_PI));
} else {
tmp = 180.0 * (atan((((C - A) - hypot((A - C), B)) / B)) / ((double) M_PI));
}
return tmp;
}
public static double code(double A, double B, double C) {
double tmp;
if (A <= -9.6e+99) {
tmp = 180.0 * (Math.atan(((-0.5 * (B + (B * (C / A)))) / -A)) / Math.PI);
} else {
tmp = 180.0 * (Math.atan((((C - A) - Math.hypot((A - C), B)) / B)) / Math.PI);
}
return tmp;
}
def code(A, B, C): tmp = 0 if A <= -9.6e+99: tmp = 180.0 * (math.atan(((-0.5 * (B + (B * (C / A)))) / -A)) / math.pi) else: tmp = 180.0 * (math.atan((((C - A) - math.hypot((A - C), B)) / B)) / math.pi) return tmp
function code(A, B, C) tmp = 0.0 if (A <= -9.6e+99) tmp = Float64(180.0 * Float64(atan(Float64(Float64(-0.5 * Float64(B + Float64(B * Float64(C / A)))) / Float64(-A))) / pi)); else tmp = Float64(180.0 * Float64(atan(Float64(Float64(Float64(C - A) - hypot(Float64(A - C), B)) / B)) / pi)); end return tmp end
function tmp_2 = code(A, B, C) tmp = 0.0; if (A <= -9.6e+99) tmp = 180.0 * (atan(((-0.5 * (B + (B * (C / A)))) / -A)) / pi); else tmp = 180.0 * (atan((((C - A) - hypot((A - C), B)) / B)) / pi); end tmp_2 = tmp; end
code[A_, B_, C_] := If[LessEqual[A, -9.6e+99], N[(180.0 * N[(N[ArcTan[N[(N[(-0.5 * N[(B + N[(B * N[(C / A), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / (-A)), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(180.0 * N[(N[ArcTan[N[(N[(N[(C - A), $MachinePrecision] - N[Sqrt[N[(A - C), $MachinePrecision] ^ 2 + B ^ 2], $MachinePrecision]), $MachinePrecision] / B), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;A \leq -9.6 \cdot 10^{+99}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{-0.5 \cdot \left(B + B \cdot \frac{C}{A}\right)}{-A}\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{\left(C - A\right) - \mathsf{hypot}\left(A - C, B\right)}{B}\right)}{\pi}\\
\end{array}
\end{array}
if A < -9.6000000000000005e99Initial program 11.7%
Taylor expanded in A around -inf
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6487.0
Applied rewrites87.0%
if -9.6000000000000005e99 < A Initial program 60.2%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
lift--.f64N/A
lift-sqrt.f64N/A
lift-+.f64N/A
lift--.f64N/A
lift-pow.f64N/A
lift-pow.f64N/A
associate-*l/N/A
lower-/.f64N/A
Applied rewrites84.8%
Final simplification85.2%
(FPCore (A B C)
:precision binary64
(let* ((t_0 (pow (hypot C B) -1.0)))
(if (<= C -1.25e+145)
(*
180.0
(/
(atan
(* (/ 1.0 B) (- (- C A) (sqrt (+ (pow (- A C) 2.0) (pow B 2.0))))))
PI))
(if (<= C 5.9e+81)
(*
180.0
(/
(atan
(*
(/ 1.0 B)
(-
(fma
(-
(* (* -0.5 (* (- 1.0 (/ (* C C) (fma C C (* B B)))) A)) t_0)
(fma (- C) t_0 1.0))
A
C)
(hypot C B))))
PI))
(*
180.0
(/
(atan (fma -0.5 (+ (/ B C) (* A (/ B (* C C)))) (/ (* 0.0 A) B)))
PI))))))
double code(double A, double B, double C) {
double t_0 = pow(hypot(C, B), -1.0);
double tmp;
if (C <= -1.25e+145) {
tmp = 180.0 * (atan(((1.0 / B) * ((C - A) - sqrt((pow((A - C), 2.0) + pow(B, 2.0)))))) / ((double) M_PI));
} else if (C <= 5.9e+81) {
tmp = 180.0 * (atan(((1.0 / B) * (fma((((-0.5 * ((1.0 - ((C * C) / fma(C, C, (B * B)))) * A)) * t_0) - fma(-C, t_0, 1.0)), A, C) - hypot(C, B)))) / ((double) M_PI));
} else {
tmp = 180.0 * (atan(fma(-0.5, ((B / C) + (A * (B / (C * C)))), ((0.0 * A) / B))) / ((double) M_PI));
}
return tmp;
}
function code(A, B, C) t_0 = hypot(C, B) ^ -1.0 tmp = 0.0 if (C <= -1.25e+145) tmp = 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)); elseif (C <= 5.9e+81) tmp = Float64(180.0 * Float64(atan(Float64(Float64(1.0 / B) * Float64(fma(Float64(Float64(Float64(-0.5 * Float64(Float64(1.0 - Float64(Float64(C * C) / fma(C, C, Float64(B * B)))) * A)) * t_0) - fma(Float64(-C), t_0, 1.0)), A, C) - hypot(C, B)))) / pi)); else tmp = Float64(180.0 * Float64(atan(fma(-0.5, Float64(Float64(B / C) + Float64(A * Float64(B / Float64(C * C)))), Float64(Float64(0.0 * A) / B))) / pi)); end return tmp end
code[A_, B_, C_] := Block[{t$95$0 = N[Power[N[Sqrt[C ^ 2 + B ^ 2], $MachinePrecision], -1.0], $MachinePrecision]}, If[LessEqual[C, -1.25e+145], 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], If[LessEqual[C, 5.9e+81], N[(180.0 * N[(N[ArcTan[N[(N[(1.0 / B), $MachinePrecision] * N[(N[(N[(N[(N[(-0.5 * N[(N[(1.0 - N[(N[(C * C), $MachinePrecision] / N[(C * C + N[(B * B), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * A), $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision] - N[((-C) * t$95$0 + 1.0), $MachinePrecision]), $MachinePrecision] * A + C), $MachinePrecision] - N[Sqrt[C ^ 2 + B ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(180.0 * N[(N[ArcTan[N[(-0.5 * N[(N[(B / C), $MachinePrecision] + N[(A * N[(B / N[(C * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(0.0 * A), $MachinePrecision] / B), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(\mathsf{hypot}\left(C, B\right)\right)}^{-1}\\
\mathbf{if}\;C \leq -1.25 \cdot 10^{+145}:\\
\;\;\;\;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}\\
\mathbf{elif}\;C \leq 5.9 \cdot 10^{+81}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{1}{B} \cdot \left(\mathsf{fma}\left(\left(-0.5 \cdot \left(\left(1 - \frac{C \cdot C}{\mathsf{fma}\left(C, C, B \cdot B\right)}\right) \cdot A\right)\right) \cdot t\_0 - \mathsf{fma}\left(-C, t\_0, 1\right), A, C\right) - \mathsf{hypot}\left(C, B\right)\right)\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\mathsf{fma}\left(-0.5, \frac{B}{C} + A \cdot \frac{B}{C \cdot C}, \frac{0 \cdot A}{B}\right)\right)}{\pi}\\
\end{array}
\end{array}
if C < -1.24999999999999992e145Initial program 77.4%
if -1.24999999999999992e145 < C < 5.9000000000000004e81Initial program 53.7%
Taylor expanded in A around 0
Applied rewrites70.0%
if 5.9000000000000004e81 < C Initial program 26.8%
Taylor expanded in C around inf
+-commutativeN/A
distribute-lft-outN/A
lower-fma.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
lower-*.f6481.2
Applied rewrites81.2%
Final simplification73.2%
(FPCore (A B C)
:precision binary64
(let* ((t_0 (pow (hypot C B) -1.0)))
(if (<= C -1.26e+132)
(* 180.0 (/ (atan (/ (* -0.5 (+ B (* B (/ C A)))) (- A))) PI))
(if (<= C 5.9e+81)
(*
180.0
(/
(atan
(*
(/ 1.0 B)
(-
(fma
(-
(* (* -0.5 (* (- 1.0 (/ (* C C) (fma C C (* B B)))) A)) t_0)
(fma (- C) t_0 1.0))
A
C)
(hypot C B))))
PI))
(*
180.0
(/
(atan (fma -0.5 (+ (/ B C) (* A (/ B (* C C)))) (/ (* 0.0 A) B)))
PI))))))
double code(double A, double B, double C) {
double t_0 = pow(hypot(C, B), -1.0);
double tmp;
if (C <= -1.26e+132) {
tmp = 180.0 * (atan(((-0.5 * (B + (B * (C / A)))) / -A)) / ((double) M_PI));
} else if (C <= 5.9e+81) {
tmp = 180.0 * (atan(((1.0 / B) * (fma((((-0.5 * ((1.0 - ((C * C) / fma(C, C, (B * B)))) * A)) * t_0) - fma(-C, t_0, 1.0)), A, C) - hypot(C, B)))) / ((double) M_PI));
} else {
tmp = 180.0 * (atan(fma(-0.5, ((B / C) + (A * (B / (C * C)))), ((0.0 * A) / B))) / ((double) M_PI));
}
return tmp;
}
function code(A, B, C) t_0 = hypot(C, B) ^ -1.0 tmp = 0.0 if (C <= -1.26e+132) tmp = Float64(180.0 * Float64(atan(Float64(Float64(-0.5 * Float64(B + Float64(B * Float64(C / A)))) / Float64(-A))) / pi)); elseif (C <= 5.9e+81) tmp = Float64(180.0 * Float64(atan(Float64(Float64(1.0 / B) * Float64(fma(Float64(Float64(Float64(-0.5 * Float64(Float64(1.0 - Float64(Float64(C * C) / fma(C, C, Float64(B * B)))) * A)) * t_0) - fma(Float64(-C), t_0, 1.0)), A, C) - hypot(C, B)))) / pi)); else tmp = Float64(180.0 * Float64(atan(fma(-0.5, Float64(Float64(B / C) + Float64(A * Float64(B / Float64(C * C)))), Float64(Float64(0.0 * A) / B))) / pi)); end return tmp end
code[A_, B_, C_] := Block[{t$95$0 = N[Power[N[Sqrt[C ^ 2 + B ^ 2], $MachinePrecision], -1.0], $MachinePrecision]}, If[LessEqual[C, -1.26e+132], N[(180.0 * N[(N[ArcTan[N[(N[(-0.5 * N[(B + N[(B * N[(C / A), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / (-A)), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], If[LessEqual[C, 5.9e+81], N[(180.0 * N[(N[ArcTan[N[(N[(1.0 / B), $MachinePrecision] * N[(N[(N[(N[(N[(-0.5 * N[(N[(1.0 - N[(N[(C * C), $MachinePrecision] / N[(C * C + N[(B * B), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * A), $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision] - N[((-C) * t$95$0 + 1.0), $MachinePrecision]), $MachinePrecision] * A + C), $MachinePrecision] - N[Sqrt[C ^ 2 + B ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(180.0 * N[(N[ArcTan[N[(-0.5 * N[(N[(B / C), $MachinePrecision] + N[(A * N[(B / N[(C * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(0.0 * A), $MachinePrecision] / B), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(\mathsf{hypot}\left(C, B\right)\right)}^{-1}\\
\mathbf{if}\;C \leq -1.26 \cdot 10^{+132}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{-0.5 \cdot \left(B + B \cdot \frac{C}{A}\right)}{-A}\right)}{\pi}\\
\mathbf{elif}\;C \leq 5.9 \cdot 10^{+81}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{1}{B} \cdot \left(\mathsf{fma}\left(\left(-0.5 \cdot \left(\left(1 - \frac{C \cdot C}{\mathsf{fma}\left(C, C, B \cdot B\right)}\right) \cdot A\right)\right) \cdot t\_0 - \mathsf{fma}\left(-C, t\_0, 1\right), A, C\right) - \mathsf{hypot}\left(C, B\right)\right)\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\mathsf{fma}\left(-0.5, \frac{B}{C} + A \cdot \frac{B}{C \cdot C}, \frac{0 \cdot A}{B}\right)\right)}{\pi}\\
\end{array}
\end{array}
if C < -1.25999999999999999e132Initial program 74.5%
Taylor expanded in A around -inf
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6473.9
Applied rewrites73.9%
if -1.25999999999999999e132 < C < 5.9000000000000004e81Initial program 54.0%
Taylor expanded in A around 0
Applied rewrites70.1%
if 5.9000000000000004e81 < C Initial program 26.8%
Taylor expanded in C around inf
+-commutativeN/A
distribute-lft-outN/A
lower-fma.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
lower-*.f6481.2
Applied rewrites81.2%
Final simplification72.8%
(FPCore (A B C)
:precision binary64
(if (<= C 1.05e+51)
(* 180.0 (/ (atan (/ (* -0.5 (+ B (* B (/ C A)))) (- A))) PI))
(*
180.0
(/ (atan (fma -0.5 (+ (/ B C) (* A (/ B (* C C)))) (/ (* 0.0 A) B))) PI))))
double code(double A, double B, double C) {
double tmp;
if (C <= 1.05e+51) {
tmp = 180.0 * (atan(((-0.5 * (B + (B * (C / A)))) / -A)) / ((double) M_PI));
} else {
tmp = 180.0 * (atan(fma(-0.5, ((B / C) + (A * (B / (C * C)))), ((0.0 * A) / B))) / ((double) M_PI));
}
return tmp;
}
function code(A, B, C) tmp = 0.0 if (C <= 1.05e+51) tmp = Float64(180.0 * Float64(atan(Float64(Float64(-0.5 * Float64(B + Float64(B * Float64(C / A)))) / Float64(-A))) / pi)); else tmp = Float64(180.0 * Float64(atan(fma(-0.5, Float64(Float64(B / C) + Float64(A * Float64(B / Float64(C * C)))), Float64(Float64(0.0 * A) / B))) / pi)); end return tmp end
code[A_, B_, C_] := If[LessEqual[C, 1.05e+51], N[(180.0 * N[(N[ArcTan[N[(N[(-0.5 * N[(B + N[(B * N[(C / A), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / (-A)), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(180.0 * N[(N[ArcTan[N[(-0.5 * N[(N[(B / C), $MachinePrecision] + N[(A * N[(B / N[(C * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(0.0 * A), $MachinePrecision] / B), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;C \leq 1.05 \cdot 10^{+51}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{-0.5 \cdot \left(B + B \cdot \frac{C}{A}\right)}{-A}\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\mathsf{fma}\left(-0.5, \frac{B}{C} + A \cdot \frac{B}{C \cdot C}, \frac{0 \cdot A}{B}\right)\right)}{\pi}\\
\end{array}
\end{array}
if C < 1.0500000000000001e51Initial program 57.9%
Taylor expanded in A around -inf
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6443.7
Applied rewrites43.7%
if 1.0500000000000001e51 < C Initial program 29.0%
Taylor expanded in C around inf
+-commutativeN/A
distribute-lft-outN/A
lower-fma.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
lower-*.f6477.0
Applied rewrites77.0%
Final simplification50.9%
(FPCore (A B C)
:precision binary64
(if (<= C 1.05e+51)
(* 180.0 (/ (atan (/ (* -0.5 (+ B (* B (/ C A)))) (- A))) PI))
(*
180.0
(/
(atan (fma -0.5 (* (+ (/ (/ B C) C) (/ (/ B A) C)) A) (/ (* 0.0 A) B)))
PI))))
double code(double A, double B, double C) {
double tmp;
if (C <= 1.05e+51) {
tmp = 180.0 * (atan(((-0.5 * (B + (B * (C / A)))) / -A)) / ((double) M_PI));
} else {
tmp = 180.0 * (atan(fma(-0.5, ((((B / C) / C) + ((B / A) / C)) * A), ((0.0 * A) / B))) / ((double) M_PI));
}
return tmp;
}
function code(A, B, C) tmp = 0.0 if (C <= 1.05e+51) tmp = Float64(180.0 * Float64(atan(Float64(Float64(-0.5 * Float64(B + Float64(B * Float64(C / A)))) / Float64(-A))) / pi)); else tmp = Float64(180.0 * Float64(atan(fma(-0.5, Float64(Float64(Float64(Float64(B / C) / C) + Float64(Float64(B / A) / C)) * A), Float64(Float64(0.0 * A) / B))) / pi)); end return tmp end
code[A_, B_, C_] := If[LessEqual[C, 1.05e+51], N[(180.0 * N[(N[ArcTan[N[(N[(-0.5 * N[(B + N[(B * N[(C / A), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / (-A)), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(180.0 * N[(N[ArcTan[N[(-0.5 * N[(N[(N[(N[(B / C), $MachinePrecision] / C), $MachinePrecision] + N[(N[(B / A), $MachinePrecision] / C), $MachinePrecision]), $MachinePrecision] * A), $MachinePrecision] + N[(N[(0.0 * A), $MachinePrecision] / B), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;C \leq 1.05 \cdot 10^{+51}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{-0.5 \cdot \left(B + B \cdot \frac{C}{A}\right)}{-A}\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\mathsf{fma}\left(-0.5, \left(\frac{\frac{B}{C}}{C} + \frac{\frac{B}{A}}{C}\right) \cdot A, \frac{0 \cdot A}{B}\right)\right)}{\pi}\\
\end{array}
\end{array}
if C < 1.0500000000000001e51Initial program 57.9%
Taylor expanded in A around -inf
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6443.7
Applied rewrites43.7%
if 1.0500000000000001e51 < C Initial program 29.0%
Taylor expanded in C around inf
+-commutativeN/A
distribute-lft-outN/A
lower-fma.f64N/A
lower-+.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
distribute-rgt1-inN/A
metadata-evalN/A
lower-*.f6477.0
Applied rewrites77.0%
Taylor expanded in A around inf
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f64N/A
pow2N/A
associate-/r*N/A
lower-/.f64N/A
lift-/.f64N/A
associate-/r*N/A
lower-/.f64N/A
lift-/.f6473.8
Applied rewrites73.8%
Final simplification50.2%
(FPCore (A B C) :precision binary64 (* 180.0 (/ (atan (/ (* -0.5 (+ B (* B (/ C A)))) (- A))) PI)))
double code(double A, double B, double C) {
return 180.0 * (atan(((-0.5 * (B + (B * (C / A)))) / -A)) / ((double) M_PI));
}
public static double code(double A, double B, double C) {
return 180.0 * (Math.atan(((-0.5 * (B + (B * (C / A)))) / -A)) / Math.PI);
}
def code(A, B, C): return 180.0 * (math.atan(((-0.5 * (B + (B * (C / A)))) / -A)) / math.pi)
function code(A, B, C) return Float64(180.0 * Float64(atan(Float64(Float64(-0.5 * Float64(B + Float64(B * Float64(C / A)))) / Float64(-A))) / pi)) end
function tmp = code(A, B, C) tmp = 180.0 * (atan(((-0.5 * (B + (B * (C / A)))) / -A)) / pi); end
code[A_, B_, C_] := N[(180.0 * N[(N[ArcTan[N[(N[(-0.5 * N[(B + N[(B * N[(C / A), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / (-A)), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
180 \cdot \frac{\tan^{-1} \left(\frac{-0.5 \cdot \left(B + B \cdot \frac{C}{A}\right)}{-A}\right)}{\pi}
\end{array}
Initial program 51.7%
Taylor expanded in A around -inf
mul-1-negN/A
lower-neg.f64N/A
lower-/.f64N/A
distribute-lft-outN/A
lower-*.f64N/A
lower-+.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6436.0
Applied rewrites36.0%
Final simplification36.0%
herbie shell --seed 2025057
(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)))