\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin th\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin thdouble f(double kx, double ky, double th) {
double r36280 = ky;
double r36281 = sin(r36280);
double r36282 = kx;
double r36283 = sin(r36282);
double r36284 = 2.0;
double r36285 = pow(r36283, r36284);
double r36286 = pow(r36281, r36284);
double r36287 = r36285 + r36286;
double r36288 = sqrt(r36287);
double r36289 = r36281 / r36288;
double r36290 = th;
double r36291 = sin(r36290);
double r36292 = r36289 * r36291;
return r36292;
}
double f(double kx, double ky, double th) {
double r36293 = ky;
double r36294 = sin(r36293);
double r36295 = kx;
double r36296 = sin(r36295);
double r36297 = 2.0;
double r36298 = pow(r36296, r36297);
double r36299 = pow(r36294, r36297);
double r36300 = r36298 + r36299;
double r36301 = sqrt(r36300);
double r36302 = r36294 / r36301;
double r36303 = th;
double r36304 = sin(r36303);
double r36305 = r36302 * r36304;
return r36305;
}



Bits error versus kx



Bits error versus ky



Bits error versus th
Results
Initial program 12.4
rmApplied *-un-lft-identity12.4
Applied *-un-lft-identity12.4
Applied times-frac12.4
Applied associate-*l*12.4
Final simplification12.4
herbie shell --seed 2020035
(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)))