2 \cdot \cos \left(\frac{2 \cdot \pi}{3} + \frac{\cos^{-1} \left(\frac{-g}{h}\right)}{3}\right)2 \cdot \left(\cos \left(\frac{2 \cdot \pi}{3} + \frac{\pi}{3}\right) \cdot \cos \left(\frac{\cos^{-1} \left(\frac{g}{h}\right)}{3}\right) + \sin \left(\frac{2 \cdot \pi}{3} + \frac{\pi}{3}\right) \cdot \sin \left(\frac{\cos^{-1} \left(\frac{g}{h}\right)}{3}\right)\right)double f(double g, double h) {
double r169171 = 2.0;
double r169172 = atan2(1.0, 0.0);
double r169173 = r169171 * r169172;
double r169174 = 3.0;
double r169175 = r169173 / r169174;
double r169176 = g;
double r169177 = -r169176;
double r169178 = h;
double r169179 = r169177 / r169178;
double r169180 = acos(r169179);
double r169181 = r169180 / r169174;
double r169182 = r169175 + r169181;
double r169183 = cos(r169182);
double r169184 = r169171 * r169183;
return r169184;
}
double f(double g, double h) {
double r169185 = 2.0;
double r169186 = atan2(1.0, 0.0);
double r169187 = r169185 * r169186;
double r169188 = 3.0;
double r169189 = r169187 / r169188;
double r169190 = r169186 / r169188;
double r169191 = r169189 + r169190;
double r169192 = cos(r169191);
double r169193 = g;
double r169194 = h;
double r169195 = r169193 / r169194;
double r169196 = acos(r169195);
double r169197 = r169196 / r169188;
double r169198 = cos(r169197);
double r169199 = r169192 * r169198;
double r169200 = sin(r169191);
double r169201 = sin(r169197);
double r169202 = r169200 * r169201;
double r169203 = r169199 + r169202;
double r169204 = r169185 * r169203;
return r169204;
}



Bits error versus g



Bits error versus h
Results
Initial program 1.0
rmApplied distribute-frac-neg1.0
Applied acos-neg1.0
Applied div-sub1.0
Applied associate-+r-1.0
Applied cos-diff0.0
Final simplification0.0
herbie shell --seed 2020036 +o rules:numerics
(FPCore (g h)
:name "2-ancestry mixing, negative discriminant"
:precision binary64
(* 2 (cos (+ (/ (* 2 PI) 3) (/ (acos (/ (- g) h)) 3)))))