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}
\begin{array}{l}
\mathbf{if}\;A \leq -2.259626659035608 \cdot 10^{+150}:\\
\;\;\;\;\left(180 \cdot \frac{1}{\sqrt{\pi}}\right) \cdot \frac{\tan^{-1} \left(\frac{0.5 \cdot \frac{B \cdot B}{A}}{B}\right)}{\sqrt{\pi}}\\
\mathbf{else}:\\
\;\;\;\;\frac{180 \cdot \tan^{-1} \left(\frac{\left(C - A\right) - \mathsf{hypot}\left(B, C - A\right)}{B}\right)}{\pi}\\
\end{array}
(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)))
(FPCore (A B C)
:precision binary64
(if (<= A -2.259626659035608e+150)
(*
(* 180.0 (/ 1.0 (sqrt PI)))
(/ (atan (/ (* 0.5 (/ (* B B) A)) B)) (sqrt PI)))
(/ (* 180.0 (atan (/ (- (- C A) (hypot B (- C A))) B))) 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));
}
double code(double A, double B, double C) {
double tmp;
if (A <= -2.259626659035608e+150) {
tmp = (180.0 * (1.0 / sqrt((double) M_PI))) * (atan((0.5 * ((B * B) / A)) / B) / sqrt((double) M_PI));
} else {
tmp = (180.0 * atan(((C - A) - hypot(B, (C - A))) / B)) / ((double) M_PI);
}
return tmp;
}



Bits error versus A



Bits error versus B



Bits error versus C
Results
if A < -2.25962665903560801e150Initial program 56.3
Simplified28.1
Applied add-sqr-sqrt_binary6428.2
Applied *-un-lft-identity_binary6428.2
Applied times-frac_binary6428.1
Applied associate-*r*_binary6428.1
Taylor expanded in A around -inf 21.0
Simplified21.0
if -2.25962665903560801e150 < A Initial program 26.0
Simplified12.2
Applied associate-*r/_binary6412.2
Final simplification13.3
herbie shell --seed 2022067
(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)))