double f(double g, double h) {
double r17184198 = 2.0;
double r17184199 = atan2(1.0, 0.0);
double r17184200 = r17184198 * r17184199;
double r17184201 = 3.0;
double r17184202 = r17184200 / r17184201;
double r17184203 = g;
double r17184204 = -r17184203;
double r17184205 = h;
double r17184206 = r17184204 / r17184205;
double r17184207 = acos(r17184206);
double r17184208 = r17184207 / r17184201;
double r17184209 = r17184202 + r17184208;
double r17184210 = cos(r17184209);
double r17184211 = r17184198 * r17184210;
return r17184211;
}
double f(double g, double h) {
double r17184212 = 2.0;
double r17184213 = 0.5;
double r17184214 = atan2(1.0, 0.0);
double r17184215 = 6.0;
double r17184216 = r17184214 / r17184215;
double r17184217 = 1.5;
double r17184218 = r17184214 / r17184217;
double r17184219 = g;
double r17184220 = h;
double r17184221 = r17184219 / r17184220;
double r17184222 = asin(r17184221);
double r17184223 = 3.0;
double r17184224 = r17184222 / r17184223;
double r17184225 = r17184218 + r17184224;
double r17184226 = r17184216 - r17184225;
double r17184227 = cos(r17184226);
double r17184228 = r17184213 * r17184227;
double r17184229 = r17184214 / r17184223;
double r17184230 = sin(r17184229);
double r17184231 = acos(r17184221);
double r17184232 = r17184231 / r17184223;
double r17184233 = r17184232 - r17184218;
double r17184234 = sin(r17184233);
double r17184235 = r17184230 * r17184234;
double r17184236 = r17184228 + r17184235;
double r17184237 = r17184212 * r17184236;
return r17184237;
}
2 \cdot \cos \left(\frac{2 \cdot \pi}{3} + \frac{\cos^{-1} \left(\frac{-g}{h}\right)}{3}\right)2 \cdot \left(\frac{1}{2} \cdot \cos \left(\frac{\pi}{6} - \left(\frac{\pi}{\frac{3}{2}} + \frac{\sin^{-1} \left(\frac{g}{h}\right)}{3}\right)\right) + \sin \left(\frac{\pi}{3}\right) \cdot \sin \left(\frac{\cos^{-1} \left(\frac{g}{h}\right)}{3} - \frac{\pi}{\frac{3}{2}}\right)\right)


Bits error versus g



Bits error versus h
Initial program 1.0
Simplified1.0
rmApplied distribute-frac-neg1.0
Applied acos-neg1.0
Applied div-sub1.0
Applied associate-+l-1.0
Applied cos-diff0.1
Simplified0.1
rmApplied acos-asin0.1
Applied div-sub0.1
Applied associate--l-0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019102
(FPCore (g h)
:name "2-ancestry mixing, negative discriminant"
(* 2 (cos (+ (/ (* 2 PI) 3) (/ (acos (/ (- g) h)) 3)))))