\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 f(double kx, double ky, double th) {
double r35978 = ky;
double r35979 = sin(r35978);
double r35980 = kx;
double r35981 = sin(r35980);
double r35982 = 2.0;
double r35983 = pow(r35981, r35982);
double r35984 = pow(r35979, r35982);
double r35985 = r35983 + r35984;
double r35986 = sqrt(r35985);
double r35987 = r35979 / r35986;
double r35988 = th;
double r35989 = sin(r35988);
double r35990 = r35987 * r35989;
return r35990;
}
double f(double kx, double ky, double th) {
double r35991 = ky;
double r35992 = sin(r35991);
double r35993 = cbrt(r35992);
double r35994 = r35993 * r35993;
double r35995 = kx;
double r35996 = sin(r35995);
double r35997 = 2.0;
double r35998 = pow(r35996, r35997);
double r35999 = pow(r35992, r35997);
double r36000 = r35998 + r35999;
double r36001 = sqrt(r36000);
double r36002 = cbrt(r36001);
double r36003 = r36002 * r36002;
double r36004 = r35994 / r36003;
double r36005 = r35993 / r36002;
double r36006 = th;
double r36007 = sin(r36006);
double r36008 = r36005 * r36007;
double r36009 = r36004 * r36008;
return r36009;
}



Bits error versus kx



Bits error versus ky



Bits error versus th
Results
Initial program 12.3
rmApplied add-cube-cbrt13.1
Applied add-cube-cbrt12.7
Applied times-frac12.7
Applied associate-*l*12.7
Final simplification12.7
herbie shell --seed 2020001
(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)))