\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(\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}}}} \cdot \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 r1179543 = ky;
double r1179544 = sin(r1179543);
double r1179545 = kx;
double r1179546 = sin(r1179545);
double r1179547 = 2.0;
double r1179548 = pow(r1179546, r1179547);
double r1179549 = pow(r1179544, r1179547);
double r1179550 = r1179548 + r1179549;
double r1179551 = sqrt(r1179550);
double r1179552 = r1179544 / r1179551;
double r1179553 = th;
double r1179554 = sin(r1179553);
double r1179555 = r1179552 * r1179554;
return r1179555;
}
double f(double kx, double ky, double th) {
double r1179556 = ky;
double r1179557 = sin(r1179556);
double r1179558 = kx;
double r1179559 = sin(r1179558);
double r1179560 = 2.0;
double r1179561 = pow(r1179559, r1179560);
double r1179562 = pow(r1179557, r1179560);
double r1179563 = r1179561 + r1179562;
double r1179564 = sqrt(r1179563);
double r1179565 = r1179557 / r1179564;
double r1179566 = cbrt(r1179565);
double r1179567 = r1179566 * r1179566;
double r1179568 = cbrt(r1179566);
double r1179569 = cbrt(r1179567);
double r1179570 = r1179568 * r1179569;
double r1179571 = th;
double r1179572 = sin(r1179571);
double r1179573 = r1179570 * r1179572;
double r1179574 = r1179567 * r1179573;
return r1179574;
}



Bits error versus kx



Bits error versus ky



Bits error versus th
Results
Initial program 12.5
rmApplied add-cube-cbrt12.8
Applied associate-*l*12.8
rmApplied add-cube-cbrt12.9
Applied cbrt-prod12.9
Final simplification12.9
herbie shell --seed 2019171
(FPCore (kx ky th)
:name "Toniolo and Linder, Equation (3b), real"
(* (/ (sin ky) (sqrt (+ (pow (sin kx) 2.0) (pow (sin ky) 2.0)))) (sin th)))