\frac{\cos th}{\sqrt{2}} \cdot \left(a1 \cdot a1\right) + \frac{\cos th}{\sqrt{2}} \cdot \left(a2 \cdot a2\right)\frac{\frac{\left(a1 \cdot a1 + a2 \cdot a2\right) \cdot \cos th}{\left|\sqrt[3]{\sqrt{2}}\right|}}{\sqrt{\sqrt[3]{\sqrt{2}}} \cdot \sqrt{\sqrt{2}}}double f(double a1, double a2, double th) {
double r115088 = th;
double r115089 = cos(r115088);
double r115090 = 2.0;
double r115091 = sqrt(r115090);
double r115092 = r115089 / r115091;
double r115093 = a1;
double r115094 = r115093 * r115093;
double r115095 = r115092 * r115094;
double r115096 = a2;
double r115097 = r115096 * r115096;
double r115098 = r115092 * r115097;
double r115099 = r115095 + r115098;
return r115099;
}
double f(double a1, double a2, double th) {
double r115100 = a1;
double r115101 = r115100 * r115100;
double r115102 = a2;
double r115103 = r115102 * r115102;
double r115104 = r115101 + r115103;
double r115105 = th;
double r115106 = cos(r115105);
double r115107 = r115104 * r115106;
double r115108 = 2.0;
double r115109 = sqrt(r115108);
double r115110 = cbrt(r115109);
double r115111 = fabs(r115110);
double r115112 = r115107 / r115111;
double r115113 = sqrt(r115110);
double r115114 = sqrt(r115109);
double r115115 = r115113 * r115114;
double r115116 = r115112 / r115115;
return r115116;
}



Bits error versus a1



Bits error versus a2



Bits error versus th
Results
Initial program 0.5
Simplified0.5
rmApplied add-sqr-sqrt0.5
Applied sqrt-prod0.5
Applied associate-/r*0.5
rmApplied add-cube-cbrt0.5
Applied sqrt-prod0.6
Applied div-inv0.5
Applied times-frac0.5
Applied associate-*l*0.4
rmApplied associate-*r*0.5
Simplified0.5
Final simplification0.5
herbie shell --seed 2019304
(FPCore (a1 a2 th)
:name "Migdal et al, Equation (64)"
:precision binary64
(+ (* (/ (cos th) (sqrt 2)) (* a1 a1)) (* (/ (cos th) (sqrt 2)) (* a2 a2))))