#include #include #include #include #include char *name = "2-ancestry mixing, negative discriminant"; double f_if(float g, float h) { float r21099 = 2; float r21100 = atan2(1.0, 0.0); float r21101 = r21099 * r21100; float r21102 = 3; float r21103 = r21101 / r21102; float r21104 = g; float r21105 = -r21104; float r21106 = h; float r21107 = r21105 / r21106; float r21108 = acos(r21107); float r21109 = r21108 / r21102; float r21110 = r21103 + r21109; float r21111 = cos(r21110); float r21112 = r21099 * r21111; return r21112; } double f_id(double g, double h) { double r21113 = 2; double r21114 = atan2(1.0, 0.0); double r21115 = r21113 * r21114; double r21116 = 3; double r21117 = r21115 / r21116; double r21118 = g; double r21119 = -r21118; double r21120 = h; double r21121 = r21119 / r21120; double r21122 = acos(r21121); double r21123 = r21122 / r21116; double r21124 = r21117 + r21123; double r21125 = cos(r21124); double r21126 = r21113 * r21125; return r21126; } double f_of(float g, float h) { float r21127 = 2; float r21128 = atan2(1.0, 0.0); float r21129 = r21128 + r21128; float r21130 = 3; float r21131 = r21129 / r21130; float r21132 = g; float r21133 = -r21132; float r21134 = h; float r21135 = r21133 / r21134; float r21136 = acos(r21135); float r21137 = r21136 / r21130; float r21138 = r21131 + r21137; float r21139 = cos(r21138); float r21140 = r21139 * r21139; float r21141 = cbrt(r21140); float r21142 = r21128 / r21130; float r21143 = r21142 + r21142; float r21144 = r21132 / r21134; float r21145 = -r21144; float r21146 = acos(r21145); float r21147 = r21146 / r21130; float r21148 = r21143 + r21147; float r21149 = cos(r21148); float r21150 = cbrt(r21149); float r21151 = r21141 * r21150; float r21152 = r21127 * r21151; return r21152; } double f_od(double g, double h) { double r21153 = 2; double r21154 = atan2(1.0, 0.0); double r21155 = r21154 + r21154; double r21156 = 3; double r21157 = r21155 / r21156; double r21158 = g; double r21159 = -r21158; double r21160 = h; double r21161 = r21159 / r21160; double r21162 = acos(r21161); double r21163 = r21162 / r21156; double r21164 = r21157 + r21163; double r21165 = cos(r21164); double r21166 = r21165 * r21165; double r21167 = cbrt(r21166); double r21168 = r21154 / r21156; double r21169 = r21168 + r21168; double r21170 = r21158 / r21160; double r21171 = -r21170; double r21172 = acos(r21171); double r21173 = r21172 / r21156; double r21174 = r21169 + r21173; double r21175 = cos(r21174); double r21176 = cbrt(r21175); double r21177 = r21167 * r21176; double r21178 = r21153 * r21177; return r21178; } void mpfr_fmod2(mpfr_t r, mpfr_t n, mpfr_t d, mpfr_rnd_t rmd) { mpfr_fmod(r, n, d, rmd); if (mpfr_cmp_ui(r, 0) < 0) mpfr_add(r, r, d, rmd); } static mpfr_t r21179, r21180, r21181, r21182, r21183, r21184, r21185, r21186, r21187, r21188, r21189, r21190, r21191, r21192; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21179, "2", 10, MPFR_RNDN); mpfr_init(r21180); mpfr_init(r21181); mpfr_init_set_str(r21182, "3", 10, MPFR_RNDN); mpfr_init(r21183); mpfr_init(r21184); mpfr_init(r21185); mpfr_init(r21186); mpfr_init(r21187); mpfr_init(r21188); mpfr_init(r21189); mpfr_init(r21190); mpfr_init(r21191); mpfr_init(r21192); } double f_im(double g, double h) { ; mpfr_const_pi(r21180, MPFR_RNDN); mpfr_mul(r21181, r21179, r21180, MPFR_RNDN); ; mpfr_div(r21183, r21181, r21182, MPFR_RNDN); mpfr_set_d(r21184, g, MPFR_RNDN); mpfr_neg(r21185, r21184, MPFR_RNDN); mpfr_set_d(r21186, h, MPFR_RNDN); mpfr_div(r21187, r21185, r21186, MPFR_RNDN); mpfr_acos(r21188, r21187, MPFR_RNDN); mpfr_div(r21189, r21188, r21182, MPFR_RNDN); mpfr_add(r21190, r21183, r21189, MPFR_RNDN); mpfr_cos(r21191, r21190, MPFR_RNDN); mpfr_mul(r21192, r21179, r21191, MPFR_RNDN); return mpfr_get_d(r21192, MPFR_RNDN); } static mpfr_t r21193, r21194, r21195, r21196, r21197, r21198, r21199, r21200, r21201, r21202, r21203, r21204, r21205, r21206, r21207, r21208, r21209, r21210, r21211, r21212, r21213, r21214, r21215, r21216, r21217, r21218; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21193, "2", 10, MPFR_RNDN); mpfr_init(r21194); mpfr_init(r21195); mpfr_init_set_str(r21196, "3", 10, MPFR_RNDN); mpfr_init(r21197); mpfr_init(r21198); mpfr_init(r21199); mpfr_init(r21200); mpfr_init(r21201); mpfr_init(r21202); mpfr_init(r21203); mpfr_init(r21204); mpfr_init(r21205); mpfr_init(r21206); mpfr_init(r21207); mpfr_init(r21208); mpfr_init(r21209); mpfr_init(r21210); mpfr_init(r21211); mpfr_init(r21212); mpfr_init(r21213); mpfr_init(r21214); mpfr_init(r21215); mpfr_init(r21216); mpfr_init(r21217); mpfr_init(r21218); } double f_fm(double g, double h) { ; mpfr_const_pi(r21194, MPFR_RNDN); mpfr_add(r21195, r21194, r21194, MPFR_RNDN); ; mpfr_div(r21197, r21195, r21196, MPFR_RNDN); mpfr_set_d(r21198, g, MPFR_RNDN); mpfr_neg(r21199, r21198, MPFR_RNDN); mpfr_set_d(r21200, h, MPFR_RNDN); mpfr_div(r21201, r21199, r21200, MPFR_RNDN); mpfr_acos(r21202, r21201, MPFR_RNDN); mpfr_div(r21203, r21202, r21196, MPFR_RNDN); mpfr_add(r21204, r21197, r21203, MPFR_RNDN); mpfr_cos(r21205, r21204, MPFR_RNDN); mpfr_mul(r21206, r21205, r21205, MPFR_RNDN); mpfr_cbrt(r21207, r21206, MPFR_RNDN); mpfr_div(r21208, r21194, r21196, MPFR_RNDN); mpfr_add(r21209, r21208, r21208, MPFR_RNDN); mpfr_div(r21210, r21198, r21200, MPFR_RNDN); mpfr_neg(r21211, r21210, MPFR_RNDN); mpfr_acos(r21212, r21211, MPFR_RNDN); mpfr_div(r21213, r21212, r21196, MPFR_RNDN); mpfr_add(r21214, r21209, r21213, MPFR_RNDN); mpfr_cos(r21215, r21214, MPFR_RNDN); mpfr_cbrt(r21216, r21215, MPFR_RNDN); mpfr_mul(r21217, r21207, r21216, MPFR_RNDN); mpfr_mul(r21218, r21193, r21217, MPFR_RNDN); return mpfr_get_d(r21218, MPFR_RNDN); } static mpfr_t r21219, r21220, r21221, r21222, r21223, r21224, r21225, r21226, r21227, r21228, r21229, r21230, r21231, r21232, r21233, r21234, r21235, r21236, r21237, r21238, r21239, r21240, r21241, r21242, r21243, r21244; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21219, "2", 10, MPFR_RNDN); mpfr_init(r21220); mpfr_init(r21221); mpfr_init_set_str(r21222, "3", 10, MPFR_RNDN); mpfr_init(r21223); mpfr_init(r21224); mpfr_init(r21225); mpfr_init(r21226); mpfr_init(r21227); mpfr_init(r21228); mpfr_init(r21229); mpfr_init(r21230); mpfr_init(r21231); mpfr_init(r21232); mpfr_init(r21233); mpfr_init(r21234); mpfr_init(r21235); mpfr_init(r21236); mpfr_init(r21237); mpfr_init(r21238); mpfr_init(r21239); mpfr_init(r21240); mpfr_init(r21241); mpfr_init(r21242); mpfr_init(r21243); mpfr_init(r21244); } double f_dm(double g, double h) { ; mpfr_const_pi(r21220, MPFR_RNDN); mpfr_add(r21221, r21220, r21220, MPFR_RNDN); ; mpfr_div(r21223, r21221, r21222, MPFR_RNDN); mpfr_set_d(r21224, g, MPFR_RNDN); mpfr_neg(r21225, r21224, MPFR_RNDN); mpfr_set_d(r21226, h, MPFR_RNDN); mpfr_div(r21227, r21225, r21226, MPFR_RNDN); mpfr_acos(r21228, r21227, MPFR_RNDN); mpfr_div(r21229, r21228, r21222, MPFR_RNDN); mpfr_add(r21230, r21223, r21229, MPFR_RNDN); mpfr_cos(r21231, r21230, MPFR_RNDN); mpfr_mul(r21232, r21231, r21231, MPFR_RNDN); mpfr_cbrt(r21233, r21232, MPFR_RNDN); mpfr_div(r21234, r21220, r21222, MPFR_RNDN); mpfr_add(r21235, r21234, r21234, MPFR_RNDN); mpfr_div(r21236, r21224, r21226, MPFR_RNDN); mpfr_neg(r21237, r21236, MPFR_RNDN); mpfr_acos(r21238, r21237, MPFR_RNDN); mpfr_div(r21239, r21238, r21222, MPFR_RNDN); mpfr_add(r21240, r21235, r21239, MPFR_RNDN); mpfr_cos(r21241, r21240, MPFR_RNDN); mpfr_cbrt(r21242, r21241, MPFR_RNDN); mpfr_mul(r21243, r21233, r21242, MPFR_RNDN); mpfr_mul(r21244, r21219, r21243, MPFR_RNDN); return mpfr_get_d(r21244, MPFR_RNDN); }