\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 r35639 = ky;
double r35640 = sin(r35639);
double r35641 = kx;
double r35642 = sin(r35641);
double r35643 = 2.0;
double r35644 = pow(r35642, r35643);
double r35645 = pow(r35640, r35643);
double r35646 = r35644 + r35645;
double r35647 = sqrt(r35646);
double r35648 = r35640 / r35647;
double r35649 = th;
double r35650 = sin(r35649);
double r35651 = r35648 * r35650;
return r35651;
}
double f(double kx, double ky, double th) {
double r35652 = ky;
double r35653 = sin(r35652);
double r35654 = kx;
double r35655 = sin(r35654);
double r35656 = 2.0;
double r35657 = pow(r35655, r35656);
double r35658 = pow(r35653, r35656);
double r35659 = r35657 + r35658;
double r35660 = sqrt(r35659);
double r35661 = r35653 / r35660;
double r35662 = cbrt(r35661);
double r35663 = r35662 * r35662;
double r35664 = cbrt(r35662);
double r35665 = r35664 * r35664;
double r35666 = r35665 * r35664;
double r35667 = th;
double r35668 = sin(r35667);
double r35669 = r35666 * r35668;
double r35670 = r35663 * r35669;
return r35670;
}



Bits error versus kx



Bits error versus ky



Bits error versus th
Results
Initial program 12.7
rmApplied add-cube-cbrt13.0
Applied associate-*l*13.0
rmApplied add-cube-cbrt13.1
Final simplification13.1
herbie shell --seed 2020081
(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)))