\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin th\frac{\frac{\sin ky}{\sqrt{\sqrt{\sqrt{\sin kx \cdot \sin kx + \sin ky \cdot \sin ky}}}}}{\sqrt{\sqrt{\sqrt{\sqrt[3]{\sin kx \cdot \sin kx + \sin ky \cdot \sin ky}} \cdot \left|\sqrt[3]{\sin kx \cdot \sin kx + \sin ky \cdot \sin ky}\right|}}} \cdot \left(\sin th \cdot \frac{1}{\sqrt{\sqrt{\sin kx \cdot \sin kx + \sin ky \cdot \sin ky}}}\right)double f(double kx, double ky, double th) {
double r888615 = ky;
double r888616 = sin(r888615);
double r888617 = kx;
double r888618 = sin(r888617);
double r888619 = 2.0;
double r888620 = pow(r888618, r888619);
double r888621 = pow(r888616, r888619);
double r888622 = r888620 + r888621;
double r888623 = sqrt(r888622);
double r888624 = r888616 / r888623;
double r888625 = th;
double r888626 = sin(r888625);
double r888627 = r888624 * r888626;
return r888627;
}
double f(double kx, double ky, double th) {
double r888628 = ky;
double r888629 = sin(r888628);
double r888630 = kx;
double r888631 = sin(r888630);
double r888632 = r888631 * r888631;
double r888633 = r888629 * r888629;
double r888634 = r888632 + r888633;
double r888635 = sqrt(r888634);
double r888636 = sqrt(r888635);
double r888637 = sqrt(r888636);
double r888638 = r888629 / r888637;
double r888639 = cbrt(r888634);
double r888640 = sqrt(r888639);
double r888641 = fabs(r888639);
double r888642 = r888640 * r888641;
double r888643 = sqrt(r888642);
double r888644 = sqrt(r888643);
double r888645 = r888638 / r888644;
double r888646 = th;
double r888647 = sin(r888646);
double r888648 = 1.0;
double r888649 = r888648 / r888636;
double r888650 = r888647 * r888649;
double r888651 = r888645 * r888650;
return r888651;
}



Bits error versus kx



Bits error versus ky



Bits error versus th
Results
Initial program 12.5
Simplified12.5
rmApplied add-sqr-sqrt12.5
Applied sqrt-prod12.8
Applied *-un-lft-identity12.8
Applied times-frac12.8
Applied associate-*r*12.8
rmApplied add-sqr-sqrt12.9
Applied associate-/r*12.9
rmApplied add-cube-cbrt12.9
Applied sqrt-prod12.9
Simplified12.9
Final simplification12.9
herbie shell --seed 2019133
(FPCore (kx ky th)
:name "Toniolo and Linder, Equation (3b), real"
(* (/ (sin ky) (sqrt (+ (pow (sin kx) 2) (pow (sin ky) 2)))) (sin th)))