\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin th\frac{\sin ky}{\sqrt{{\left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right)}^{2} \cdot {\left(\sqrt[3]{\sin kx}\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin thdouble 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(ky) / sqrt(((pow((cbrt(sin(kx)) * cbrt(sin(kx))), 2.0) * pow(cbrt(sin(kx)), 2.0)) + pow(sin(ky), 2.0)))) * sin(th));
}



Bits error versus kx



Bits error versus ky



Bits error versus th
Results
Initial program 4.0
rmApplied add-cube-cbrt4.3
Applied unpow-prod-down4.3
Final simplification4.3
herbie shell --seed 2020103
(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)))