\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin th\left(\sqrt[3]{\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}}} \cdot \sqrt[3]{\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}}}\right) \cdot \left(\left(\left(\sqrt[3]{\sqrt[3]{\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}}}} \cdot \sqrt[3]{\sqrt[3]{\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}}}}\right) \cdot \sqrt[3]{\sqrt[3]{\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}}}}\right) \cdot \sin th\right)double f(double kx, double ky, double th) {
double r36571 = ky;
double r36572 = sin(r36571);
double r36573 = kx;
double r36574 = sin(r36573);
double r36575 = 2.0;
double r36576 = pow(r36574, r36575);
double r36577 = pow(r36572, r36575);
double r36578 = r36576 + r36577;
double r36579 = sqrt(r36578);
double r36580 = r36572 / r36579;
double r36581 = th;
double r36582 = sin(r36581);
double r36583 = r36580 * r36582;
return r36583;
}
double f(double kx, double ky, double th) {
double r36584 = ky;
double r36585 = sin(r36584);
double r36586 = kx;
double r36587 = sin(r36586);
double r36588 = 2.0;
double r36589 = pow(r36587, r36588);
double r36590 = pow(r36585, r36588);
double r36591 = r36589 + r36590;
double r36592 = sqrt(r36591);
double r36593 = r36585 / r36592;
double r36594 = cbrt(r36593);
double r36595 = r36594 * r36594;
double r36596 = cbrt(r36594);
double r36597 = r36596 * r36596;
double r36598 = r36597 * r36596;
double r36599 = th;
double r36600 = sin(r36599);
double r36601 = r36598 * r36600;
double r36602 = r36595 * r36601;
return r36602;
}



Bits error versus kx



Bits error versus ky



Bits error versus th
Results
Initial program 12.2
rmApplied add-cube-cbrt12.6
Applied associate-*l*12.6
rmApplied add-cube-cbrt12.7
Final simplification12.7
herbie shell --seed 2020056
(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)))