{\left(a \cdot \cos \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2} + {\left(b \cdot \sin \left(\pi \cdot \frac{angle}{180}\right)\right)}^{2}
{\left(a \cdot \cos \left(\frac{\pi \cdot angle}{180}\right)\right)}^{2} + {\left(b \cdot \mathsf{log1p}\left(\mathsf{expm1}\left(\sin \left(\sqrt{0.005555555555555556} \cdot \left(\left(\pi \cdot angle\right) \cdot \sqrt{0.005555555555555556}\right)\right)\right)\right)\right)}^{2}
(FPCore (a b angle) :precision binary64 (+ (pow (* a (cos (* PI (/ angle 180.0)))) 2.0) (pow (* b (sin (* PI (/ angle 180.0)))) 2.0)))
(FPCore (a b angle)
:precision binary64
(+
(pow (* a (cos (/ (* PI angle) 180.0))) 2.0)
(pow
(*
b
(log1p
(expm1
(sin
(*
(sqrt 0.005555555555555556)
(* (* PI angle) (sqrt 0.005555555555555556)))))))
2.0)))double code(double a, double b, double angle) {
return pow((a * cos(((double) M_PI) * (angle / 180.0))), 2.0) + pow((b * sin(((double) M_PI) * (angle / 180.0))), 2.0);
}
double code(double a, double b, double angle) {
return pow((a * cos((((double) M_PI) * angle) / 180.0)), 2.0) + pow((b * log1p(expm1(sin(sqrt(0.005555555555555556) * ((((double) M_PI) * angle) * sqrt(0.005555555555555556)))))), 2.0);
}



Bits error versus a



Bits error versus b



Bits error versus angle
Results
Initial program 12.9
Applied associate-*r/_binary6413.0
Applied associate-*r/_binary6413.0
Simplified13.0
Applied log1p-expm1-u_binary6413.0
Simplified13.0
Applied add-sqr-sqrt_binary6413.0
Applied associate-*l*_binary6413.0
Final simplification13.0
herbie shell --seed 2022067
(FPCore (a b angle)
:name "ab-angle->ABCF C"
:precision binary64
(+ (pow (* a (cos (* PI (/ angle 180.0)))) 2.0) (pow (* b (sin (* PI (/ angle 180.0)))) 2.0)))