\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin th\frac{\sin th}{\frac{\mathsf{hypot}\left({\left(\sin kx\right)}^{\left(\frac{2}{2}\right)}, {\left(\sin ky\right)}^{\left(\frac{2}{2}\right)}\right)}{\sin ky}}double code(double kx, double ky, double th) {
return ((sin(ky) / sqrt((pow(sin(kx), 2.0) + pow(sin(ky), 2.0)))) * sin(th));
}
double code(double kx, double ky, double th) {
return (sin(th) / (hypot(pow(sin(kx), (2.0 / 2.0)), pow(sin(ky), (2.0 / 2.0))) / sin(ky)));
}



Bits error versus kx



Bits error versus ky



Bits error versus th
Results
Initial program 3.8
rmApplied sqr-pow3.8
Applied sqr-pow3.8
Applied hypot-def0.2
rmApplied clear-num0.3
rmApplied associate-*l/0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2020075 +o rules:numerics
(FPCore (kx ky th)
:name "Toniolo and Linder, Equation (3b), real"
:precision binary64
(* (/ (sin ky) (sqrt (+ (pow (sin kx) 2) (pow (sin ky) 2)))) (sin th)))