\frac{\cos th}{\sqrt{2}} \cdot \left(a1 \cdot a1\right) + \frac{\cos th}{\sqrt{2}} \cdot \left(a2 \cdot a2\right)\cos th \cdot \frac{a2 \cdot a2}{\sqrt{2}} + \frac{\frac{\cos th}{\sqrt{\sqrt{2}}}}{\sqrt{\sqrt{2}}} \cdot \left(a1 \cdot a1\right)double f(double a1, double a2, double th) {
double r2398636 = th;
double r2398637 = cos(r2398636);
double r2398638 = 2.0;
double r2398639 = sqrt(r2398638);
double r2398640 = r2398637 / r2398639;
double r2398641 = a1;
double r2398642 = r2398641 * r2398641;
double r2398643 = r2398640 * r2398642;
double r2398644 = a2;
double r2398645 = r2398644 * r2398644;
double r2398646 = r2398640 * r2398645;
double r2398647 = r2398643 + r2398646;
return r2398647;
}
double f(double a1, double a2, double th) {
double r2398648 = th;
double r2398649 = cos(r2398648);
double r2398650 = a2;
double r2398651 = r2398650 * r2398650;
double r2398652 = 2.0;
double r2398653 = sqrt(r2398652);
double r2398654 = r2398651 / r2398653;
double r2398655 = r2398649 * r2398654;
double r2398656 = sqrt(r2398653);
double r2398657 = r2398649 / r2398656;
double r2398658 = r2398657 / r2398656;
double r2398659 = a1;
double r2398660 = r2398659 * r2398659;
double r2398661 = r2398658 * r2398660;
double r2398662 = r2398655 + r2398661;
return r2398662;
}



Bits error versus a1



Bits error versus a2



Bits error versus th
Results
Initial program 0.5
rmApplied div-inv0.5
Applied associate-*l*0.5
Simplified0.5
rmApplied add-sqr-sqrt0.5
Applied sqrt-prod0.5
Applied associate-/r*0.5
Final simplification0.5
herbie shell --seed 2019172
(FPCore (a1 a2 th)
:name "Migdal et al, Equation (64)"
(+ (* (/ (cos th) (sqrt 2.0)) (* a1 a1)) (* (/ (cos th) (sqrt 2.0)) (* a2 a2))))