\frac{\sin ky}{\sqrt{{\sin kx}^{2} + {\sin ky}^{2}}} \cdot \sin th\sin th \cdot \frac{\sin ky}{\sqrt{{\sin kx}^{2} + {\sin ky}^{2}}}(FPCore (kx ky th) :precision binary64 (* (/ (sin ky) (sqrt (+ (pow (sin kx) 2.0) (pow (sin ky) 2.0)))) (sin th)))
(FPCore (kx ky th) :precision binary64 (* (sin th) (/ (sin ky) (sqrt (+ (pow (sin kx) 2.0) (pow (sin ky) 2.0))))))
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) * (sin(ky) / sqrt(pow(sin(kx), 2.0) + pow(sin(ky), 2.0)));
}



Bits error versus kx



Bits error versus ky



Bits error versus th
Results
Initial program 3.8
rmApplied *-un-lft-identity_binary643.8
Applied associate-*l*_binary643.8
Simplified3.8
rmApplied pow1_binary643.8
Final simplification3.8
herbie shell --seed 2021110
(FPCore (kx ky th)
:name "Toniolo and Linder, Equation (3b), real"
:precision binary64
(* (/ (sin ky) (sqrt (+ (pow (sin kx) 2.0) (pow (sin ky) 2.0)))) (sin th)))