\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin th\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + \left(\sqrt[3]{{\left(\sin ky\right)}^{2}} \cdot \sqrt[3]{{\left(\sin ky\right)}^{2}}\right) \cdot \sqrt[3]{{\left(\sin ky\right)}^{2}}}} \cdot \sin thdouble code(double kx, double ky, double th) {
return ((double) (((double) (((double) sin(ky)) / ((double) sqrt(((double) (((double) pow(((double) sin(kx)), 2.0)) + ((double) pow(((double) sin(ky)), 2.0)))))))) * ((double) sin(th))));
}
double code(double kx, double ky, double th) {
return ((double) (((double) (((double) sin(ky)) / ((double) sqrt(((double) (((double) pow(((double) sin(kx)), 2.0)) + ((double) (((double) (((double) cbrt(((double) pow(((double) sin(ky)), 2.0)))) * ((double) cbrt(((double) pow(((double) sin(ky)), 2.0)))))) * ((double) cbrt(((double) pow(((double) sin(ky)), 2.0)))))))))))) * ((double) sin(th))));
}



Bits error versus kx



Bits error versus ky



Bits error versus th
Results
Initial program 4.0
rmApplied add-cube-cbrt4.4
Final simplification4.4
herbie shell --seed 2020163
(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)))