\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin th\frac{\sin ky}{\sqrt{{\left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right)}^{2} \cdot {\left(\left(\sqrt[3]{\sqrt[3]{\sin kx}} \cdot \sqrt[3]{\sqrt[3]{\sin kx}}\right) \cdot \left(\sqrt[3]{\sqrt[3]{\sqrt[3]{\sin kx}} \cdot \sqrt[3]{\sqrt[3]{\sin kx}}} \cdot \sqrt[3]{\sqrt[3]{\sqrt[3]{\sqrt[3]{\sin kx}} \cdot \sqrt[3]{\sqrt[3]{\sin kx}}} \cdot \sqrt[3]{\sqrt[3]{\sqrt[3]{\sin kx}}}}\right)\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin thdouble f(double kx, double ky, double th) {
double r35437 = ky;
double r35438 = sin(r35437);
double r35439 = kx;
double r35440 = sin(r35439);
double r35441 = 2.0;
double r35442 = pow(r35440, r35441);
double r35443 = pow(r35438, r35441);
double r35444 = r35442 + r35443;
double r35445 = sqrt(r35444);
double r35446 = r35438 / r35445;
double r35447 = th;
double r35448 = sin(r35447);
double r35449 = r35446 * r35448;
return r35449;
}
double f(double kx, double ky, double th) {
double r35450 = ky;
double r35451 = sin(r35450);
double r35452 = kx;
double r35453 = sin(r35452);
double r35454 = cbrt(r35453);
double r35455 = r35454 * r35454;
double r35456 = 2.0;
double r35457 = pow(r35455, r35456);
double r35458 = cbrt(r35454);
double r35459 = r35458 * r35458;
double r35460 = cbrt(r35459);
double r35461 = cbrt(r35458);
double r35462 = r35460 * r35461;
double r35463 = cbrt(r35462);
double r35464 = r35460 * r35463;
double r35465 = r35459 * r35464;
double r35466 = pow(r35465, r35456);
double r35467 = r35457 * r35466;
double r35468 = pow(r35451, r35456);
double r35469 = r35467 + r35468;
double r35470 = sqrt(r35469);
double r35471 = r35451 / r35470;
double r35472 = th;
double r35473 = sin(r35472);
double r35474 = r35471 * r35473;
return r35474;
}



Bits error versus kx



Bits error versus ky



Bits error versus th
Results
Initial program 12.4
rmApplied add-cube-cbrt12.6
Applied unpow-prod-down12.6
rmApplied add-cube-cbrt12.7
rmApplied add-cube-cbrt12.7
Applied cbrt-prod12.7
rmApplied add-cube-cbrt12.7
Applied cbrt-prod12.7
Final simplification12.7
herbie shell --seed 2019195
(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)))