#include #include #include #include #include char *name = "2-ancestry mixing, positive discriminant"; double f_if(float g, float h, float a) { float r21150 = 1; float r21151 = 2; float r21152 = a; float r21153 = r21151 * r21152; float r21154 = r21150 / r21153; float r21155 = g; float r21156 = -r21155; float r21157 = r21155 * r21155; float r21158 = h; float r21159 = r21158 * r21158; float r21160 = r21157 - r21159; float r21161 = sqrt(r21160); float r21162 = r21156 + r21161; float r21163 = r21154 * r21162; float r21164 = cbrt(r21163); float r21165 = r21156 - r21161; float r21166 = r21154 * r21165; float r21167 = cbrt(r21166); float r21168 = r21164 + r21167; return r21168; } double f_id(double g, double h, double a) { double r21169 = 1; double r21170 = 2; double r21171 = a; double r21172 = r21170 * r21171; double r21173 = r21169 / r21172; double r21174 = g; double r21175 = -r21174; double r21176 = r21174 * r21174; double r21177 = h; double r21178 = r21177 * r21177; double r21179 = r21176 - r21178; double r21180 = sqrt(r21179); double r21181 = r21175 + r21180; double r21182 = r21173 * r21181; double r21183 = cbrt(r21182); double r21184 = r21175 - r21180; double r21185 = r21173 * r21184; double r21186 = cbrt(r21185); double r21187 = r21183 + r21186; return r21187; } double f_of(float g, float h, float a) { float r21188 = g; float r21189 = 2.6166369468200776e-178; bool r21190 = r21188 <= r21189; float r21191 = r21188 + r21188; float r21192 = -r21191; float r21193 = a; float r21194 = r21193 + r21193; float r21195 = cbrt(r21194); float r21196 = r21192 / r21195; float r21197 = cbrt(r21196); float r21198 = 1; float r21199 = r21195 * r21195; float r21200 = r21198 / r21199; float r21201 = cbrt(r21200); float r21202 = r21197 * r21201; float r21203 = -r21188; float r21204 = h; float r21205 = r21204 + r21188; float r21206 = r21188 - r21204; float r21207 = r21205 * r21206; float r21208 = sqrt(r21207); float r21209 = r21203 - r21208; float r21210 = r21209 / r21194; float r21211 = cbrt(r21210); float r21212 = r21202 + r21211; float r21213 = r21188 + r21204; float r21214 = r21213 * r21206; float r21215 = sqrt(r21214); float r21216 = r21203 + r21215; float r21217 = r21216 / r21194; float r21218 = cbrt(r21217); float r21219 = r21203 - r21215; float r21220 = cbrt(r21219); float r21221 = r21198 / r21194; float r21222 = cbrt(r21221); float r21223 = r21220 * r21222; float r21224 = r21218 + r21223; float r21225 = r21190 ? r21212 : r21224; return r21225; } double f_od(double g, double h, double a) { double r21226 = g; double r21227 = 2.6166369468200776e-178; bool r21228 = r21226 <= r21227; double r21229 = r21226 + r21226; double r21230 = -r21229; double r21231 = a; double r21232 = r21231 + r21231; double r21233 = cbrt(r21232); double r21234 = r21230 / r21233; double r21235 = cbrt(r21234); double r21236 = 1; double r21237 = r21233 * r21233; double r21238 = r21236 / r21237; double r21239 = cbrt(r21238); double r21240 = r21235 * r21239; double r21241 = -r21226; double r21242 = h; double r21243 = r21242 + r21226; double r21244 = r21226 - r21242; double r21245 = r21243 * r21244; double r21246 = sqrt(r21245); double r21247 = r21241 - r21246; double r21248 = r21247 / r21232; double r21249 = cbrt(r21248); double r21250 = r21240 + r21249; double r21251 = r21226 + r21242; double r21252 = r21251 * r21244; double r21253 = sqrt(r21252); double r21254 = r21241 + r21253; double r21255 = r21254 / r21232; double r21256 = cbrt(r21255); double r21257 = r21241 - r21253; double r21258 = cbrt(r21257); double r21259 = r21236 / r21232; double r21260 = cbrt(r21259); double r21261 = r21258 * r21260; double r21262 = r21256 + r21261; double r21263 = r21228 ? r21250 : r21262; return r21263; } 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 r21264, r21265, r21266, r21267, r21268, r21269, r21270, r21271, r21272, r21273, r21274, r21275, r21276, r21277, r21278, r21279, r21280, r21281, r21282; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init_set_str(r21264, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21265, "2", 10, MPFR_RNDN); mpfr_init(r21266); mpfr_init(r21267); mpfr_init(r21268); mpfr_init(r21269); mpfr_init(r21270); mpfr_init(r21271); mpfr_init(r21272); mpfr_init(r21273); mpfr_init(r21274); mpfr_init(r21275); mpfr_init(r21276); mpfr_init(r21277); mpfr_init(r21278); mpfr_init(r21279); mpfr_init(r21280); mpfr_init(r21281); mpfr_init(r21282); } double f_im(double g, double h, double a) { ; ; mpfr_set_d(r21266, a, MPFR_RNDN); mpfr_mul(r21267, r21265, r21266, MPFR_RNDN); mpfr_div(r21268, r21264, r21267, MPFR_RNDN); mpfr_set_d(r21269, g, MPFR_RNDN); mpfr_neg(r21270, r21269, MPFR_RNDN); mpfr_mul(r21271, r21269, r21269, MPFR_RNDN); mpfr_set_d(r21272, h, MPFR_RNDN); mpfr_mul(r21273, r21272, r21272, MPFR_RNDN); mpfr_sub(r21274, r21271, r21273, MPFR_RNDN); mpfr_sqrt(r21275, r21274, MPFR_RNDN); mpfr_add(r21276, r21270, r21275, MPFR_RNDN); mpfr_mul(r21277, r21268, r21276, MPFR_RNDN); mpfr_cbrt(r21278, r21277, MPFR_RNDN); mpfr_sub(r21279, r21270, r21275, MPFR_RNDN); mpfr_mul(r21280, r21268, r21279, MPFR_RNDN); mpfr_cbrt(r21281, r21280, MPFR_RNDN); mpfr_add(r21282, r21278, r21281, MPFR_RNDN); return mpfr_get_d(r21282, MPFR_RNDN); } static mpfr_t r21283, r21284, r21285, r21286, r21287, r21288, r21289, r21290, r21291, r21292, r21293, r21294, r21295, r21296, r21297, r21298, r21299, r21300, r21301, r21302, r21303, r21304, r21305, r21306, r21307, r21308, r21309, r21310, r21311, r21312, r21313, r21314, r21315, r21316, r21317, r21318, r21319, r21320; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init(r21283); mpfr_init_set_str(r21284, "2.6166369468200776e-178", 10, MPFR_RNDN); mpfr_init(r21285); mpfr_init(r21286); mpfr_init(r21287); mpfr_init(r21288); mpfr_init(r21289); mpfr_init(r21290); mpfr_init(r21291); mpfr_init(r21292); mpfr_init_set_str(r21293, "1", 10, MPFR_RNDN); mpfr_init(r21294); mpfr_init(r21295); mpfr_init(r21296); mpfr_init(r21297); mpfr_init(r21298); mpfr_init(r21299); mpfr_init(r21300); mpfr_init(r21301); mpfr_init(r21302); mpfr_init(r21303); mpfr_init(r21304); mpfr_init(r21305); mpfr_init(r21306); mpfr_init(r21307); mpfr_init(r21308); mpfr_init(r21309); mpfr_init(r21310); mpfr_init(r21311); mpfr_init(r21312); mpfr_init(r21313); mpfr_init(r21314); mpfr_init(r21315); mpfr_init(r21316); mpfr_init(r21317); mpfr_init(r21318); mpfr_init(r21319); mpfr_init(r21320); } double f_fm(double g, double h, double a) { mpfr_set_d(r21283, g, MPFR_RNDN); ; mpfr_set_si(r21285, mpfr_cmp(r21283, r21284) <= 0, MPFR_RNDN); mpfr_add(r21286, r21283, r21283, MPFR_RNDN); mpfr_neg(r21287, r21286, MPFR_RNDN); mpfr_set_d(r21288, a, MPFR_RNDN); mpfr_add(r21289, r21288, r21288, MPFR_RNDN); mpfr_cbrt(r21290, r21289, MPFR_RNDN); mpfr_div(r21291, r21287, r21290, MPFR_RNDN); mpfr_cbrt(r21292, r21291, MPFR_RNDN); ; mpfr_mul(r21294, r21290, r21290, MPFR_RNDN); mpfr_div(r21295, r21293, r21294, MPFR_RNDN); mpfr_cbrt(r21296, r21295, MPFR_RNDN); mpfr_mul(r21297, r21292, r21296, MPFR_RNDN); mpfr_neg(r21298, r21283, MPFR_RNDN); mpfr_set_d(r21299, h, MPFR_RNDN); mpfr_add(r21300, r21299, r21283, MPFR_RNDN); mpfr_sub(r21301, r21283, r21299, MPFR_RNDN); mpfr_mul(r21302, r21300, r21301, MPFR_RNDN); mpfr_sqrt(r21303, r21302, MPFR_RNDN); mpfr_sub(r21304, r21298, r21303, MPFR_RNDN); mpfr_div(r21305, r21304, r21289, MPFR_RNDN); mpfr_cbrt(r21306, r21305, MPFR_RNDN); mpfr_add(r21307, r21297, r21306, MPFR_RNDN); mpfr_add(r21308, r21283, r21299, MPFR_RNDN); mpfr_mul(r21309, r21308, r21301, MPFR_RNDN); mpfr_sqrt(r21310, r21309, MPFR_RNDN); mpfr_add(r21311, r21298, r21310, MPFR_RNDN); mpfr_div(r21312, r21311, r21289, MPFR_RNDN); mpfr_cbrt(r21313, r21312, MPFR_RNDN); mpfr_sub(r21314, r21298, r21310, MPFR_RNDN); mpfr_cbrt(r21315, r21314, MPFR_RNDN); mpfr_div(r21316, r21293, r21289, MPFR_RNDN); mpfr_cbrt(r21317, r21316, MPFR_RNDN); mpfr_mul(r21318, r21315, r21317, MPFR_RNDN); mpfr_add(r21319, r21313, r21318, MPFR_RNDN); if (mpfr_get_si(r21285, MPFR_RNDN)) { mpfr_set(r21320, r21307, MPFR_RNDN); } else { mpfr_set(r21320, r21319, MPFR_RNDN); }; return mpfr_get_d(r21320, MPFR_RNDN); } static mpfr_t r21321, r21322, r21323, r21324, r21325, r21326, r21327, r21328, r21329, r21330, r21331, r21332, r21333, r21334, r21335, r21336, r21337, r21338, r21339, r21340, r21341, r21342, r21343, r21344, r21345, r21346, r21347, r21348, r21349, r21350, r21351, r21352, r21353, r21354, r21355, r21356, r21357, r21358; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init(r21321); mpfr_init_set_str(r21322, "2.6166369468200776e-178", 10, MPFR_RNDN); mpfr_init(r21323); mpfr_init(r21324); mpfr_init(r21325); mpfr_init(r21326); mpfr_init(r21327); mpfr_init(r21328); mpfr_init(r21329); mpfr_init(r21330); mpfr_init_set_str(r21331, "1", 10, MPFR_RNDN); mpfr_init(r21332); mpfr_init(r21333); mpfr_init(r21334); mpfr_init(r21335); mpfr_init(r21336); mpfr_init(r21337); mpfr_init(r21338); mpfr_init(r21339); mpfr_init(r21340); mpfr_init(r21341); mpfr_init(r21342); mpfr_init(r21343); mpfr_init(r21344); mpfr_init(r21345); mpfr_init(r21346); mpfr_init(r21347); mpfr_init(r21348); mpfr_init(r21349); mpfr_init(r21350); mpfr_init(r21351); mpfr_init(r21352); mpfr_init(r21353); mpfr_init(r21354); mpfr_init(r21355); mpfr_init(r21356); mpfr_init(r21357); mpfr_init(r21358); } double f_dm(double g, double h, double a) { mpfr_set_d(r21321, g, MPFR_RNDN); ; mpfr_set_si(r21323, mpfr_cmp(r21321, r21322) <= 0, MPFR_RNDN); mpfr_add(r21324, r21321, r21321, MPFR_RNDN); mpfr_neg(r21325, r21324, MPFR_RNDN); mpfr_set_d(r21326, a, MPFR_RNDN); mpfr_add(r21327, r21326, r21326, MPFR_RNDN); mpfr_cbrt(r21328, r21327, MPFR_RNDN); mpfr_div(r21329, r21325, r21328, MPFR_RNDN); mpfr_cbrt(r21330, r21329, MPFR_RNDN); ; mpfr_mul(r21332, r21328, r21328, MPFR_RNDN); mpfr_div(r21333, r21331, r21332, MPFR_RNDN); mpfr_cbrt(r21334, r21333, MPFR_RNDN); mpfr_mul(r21335, r21330, r21334, MPFR_RNDN); mpfr_neg(r21336, r21321, MPFR_RNDN); mpfr_set_d(r21337, h, MPFR_RNDN); mpfr_add(r21338, r21337, r21321, MPFR_RNDN); mpfr_sub(r21339, r21321, r21337, MPFR_RNDN); mpfr_mul(r21340, r21338, r21339, MPFR_RNDN); mpfr_sqrt(r21341, r21340, MPFR_RNDN); mpfr_sub(r21342, r21336, r21341, MPFR_RNDN); mpfr_div(r21343, r21342, r21327, MPFR_RNDN); mpfr_cbrt(r21344, r21343, MPFR_RNDN); mpfr_add(r21345, r21335, r21344, MPFR_RNDN); mpfr_add(r21346, r21321, r21337, MPFR_RNDN); mpfr_mul(r21347, r21346, r21339, MPFR_RNDN); mpfr_sqrt(r21348, r21347, MPFR_RNDN); mpfr_add(r21349, r21336, r21348, MPFR_RNDN); mpfr_div(r21350, r21349, r21327, MPFR_RNDN); mpfr_cbrt(r21351, r21350, MPFR_RNDN); mpfr_sub(r21352, r21336, r21348, MPFR_RNDN); mpfr_cbrt(r21353, r21352, MPFR_RNDN); mpfr_div(r21354, r21331, r21327, MPFR_RNDN); mpfr_cbrt(r21355, r21354, MPFR_RNDN); mpfr_mul(r21356, r21353, r21355, MPFR_RNDN); mpfr_add(r21357, r21351, r21356, MPFR_RNDN); if (mpfr_get_si(r21323, MPFR_RNDN)) { mpfr_set(r21358, r21345, MPFR_RNDN); } else { mpfr_set(r21358, r21357, MPFR_RNDN); }; return mpfr_get_d(r21358, MPFR_RNDN); }