\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin th\frac{\sqrt[3]{\sin ky} \cdot \sqrt[3]{\sin ky}}{\sqrt[3]{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sqrt[3]{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}}} \cdot \left(\frac{\sqrt[3]{\sin ky}}{\sqrt[3]{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}}} \cdot \sin th\right)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 (((cbrt(sin(ky)) * cbrt(sin(ky))) / (cbrt(sqrt((pow(sin(kx), 2.0) + pow(sin(ky), 2.0)))) * cbrt(sqrt((pow(sin(kx), 2.0) + pow(sin(ky), 2.0)))))) * ((cbrt(sin(ky)) / cbrt(sqrt((pow(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 12.5
rmApplied add-cube-cbrt13.3
Applied add-cube-cbrt12.9
Applied times-frac12.9
Applied associate-*l*12.9
Final simplification12.9
herbie shell --seed 2020049
(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)))