#include #include #include #include #include char *name = "2-ancestry mixing, negative discriminant"; double f_if(float g, float h) { float r21071 = 2; float r21072 = atan2(1.0, 0.0); float r21073 = r21071 * r21072; float r21074 = 3; float r21075 = r21073 / r21074; float r21076 = g; float r21077 = -r21076; float r21078 = h; float r21079 = r21077 / r21078; float r21080 = acos(r21079); float r21081 = r21080 / r21074; float r21082 = r21075 + r21081; float r21083 = cos(r21082); float r21084 = r21071 * r21083; return r21084; } double f_id(double g, double h) { double r21085 = 2; double r21086 = atan2(1.0, 0.0); double r21087 = r21085 * r21086; double r21088 = 3; double r21089 = r21087 / r21088; double r21090 = g; double r21091 = -r21090; double r21092 = h; double r21093 = r21091 / r21092; double r21094 = acos(r21093); double r21095 = r21094 / r21088; double r21096 = r21089 + r21095; double r21097 = cos(r21096); double r21098 = r21085 * r21097; return r21098; } double f_of(float g, float h) { float r21099 = 2; float r21100 = atan2(1.0, 0.0); float r21101 = r21100 + 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 = r21111 * r21111; float r21113 = cbrt(r21112); float r21114 = r21100 / r21102; float r21115 = r21114 + r21114; float r21116 = r21104 / r21106; float r21117 = -r21116; float r21118 = acos(r21117); float r21119 = r21118 / r21102; float r21120 = r21115 + r21119; float r21121 = cos(r21120); float r21122 = cbrt(r21121); float r21123 = r21113 * r21122; float r21124 = r21099 * r21123; return r21124; } double f_od(double g, double h) { double r21125 = 2; double r21126 = atan2(1.0, 0.0); double r21127 = r21126 + r21126; double r21128 = 3; double r21129 = r21127 / r21128; double r21130 = g; double r21131 = -r21130; double r21132 = h; double r21133 = r21131 / r21132; double r21134 = acos(r21133); double r21135 = r21134 / r21128; double r21136 = r21129 + r21135; double r21137 = cos(r21136); double r21138 = r21137 * r21137; double r21139 = cbrt(r21138); double r21140 = r21126 / r21128; double r21141 = r21140 + r21140; double r21142 = r21130 / r21132; double r21143 = -r21142; double r21144 = acos(r21143); double r21145 = r21144 / r21128; double r21146 = r21141 + r21145; double r21147 = cos(r21146); double r21148 = cbrt(r21147); double r21149 = r21139 * r21148; double r21150 = r21125 * r21149; return r21150; } 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 r21151, r21152, r21153, r21154, r21155, r21156, r21157, r21158, r21159, r21160, r21161, r21162, r21163, r21164; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21151, "2", 10, MPFR_RNDN); mpfr_init(r21152); mpfr_init(r21153); mpfr_init_set_str(r21154, "3", 10, MPFR_RNDN); mpfr_init(r21155); mpfr_init(r21156); mpfr_init(r21157); mpfr_init(r21158); mpfr_init(r21159); mpfr_init(r21160); mpfr_init(r21161); mpfr_init(r21162); mpfr_init(r21163); mpfr_init(r21164); } double f_im(double g, double h) { ; mpfr_const_pi(r21152, MPFR_RNDN); mpfr_mul(r21153, r21151, r21152, MPFR_RNDN); ; mpfr_div(r21155, r21153, r21154, MPFR_RNDN); mpfr_set_d(r21156, g, MPFR_RNDN); mpfr_neg(r21157, r21156, MPFR_RNDN); mpfr_set_d(r21158, h, MPFR_RNDN); mpfr_div(r21159, r21157, r21158, MPFR_RNDN); mpfr_acos(r21160, r21159, MPFR_RNDN); mpfr_div(r21161, r21160, r21154, MPFR_RNDN); mpfr_add(r21162, r21155, r21161, MPFR_RNDN); mpfr_cos(r21163, r21162, MPFR_RNDN); mpfr_mul(r21164, r21151, r21163, MPFR_RNDN); return mpfr_get_d(r21164, MPFR_RNDN); } static mpfr_t r21165, r21166, r21167, r21168, r21169, r21170, r21171, r21172, r21173, r21174, r21175, r21176, r21177, r21178, r21179, r21180, r21181, r21182, r21183, r21184, r21185, r21186, r21187, r21188, r21189, r21190; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21165, "2", 10, MPFR_RNDN); mpfr_init(r21166); mpfr_init(r21167); mpfr_init_set_str(r21168, "3", 10, MPFR_RNDN); mpfr_init(r21169); mpfr_init(r21170); mpfr_init(r21171); mpfr_init(r21172); mpfr_init(r21173); mpfr_init(r21174); mpfr_init(r21175); mpfr_init(r21176); mpfr_init(r21177); mpfr_init(r21178); mpfr_init(r21179); mpfr_init(r21180); mpfr_init(r21181); mpfr_init(r21182); mpfr_init(r21183); mpfr_init(r21184); mpfr_init(r21185); mpfr_init(r21186); mpfr_init(r21187); mpfr_init(r21188); mpfr_init(r21189); mpfr_init(r21190); } double f_fm(double g, double h) { ; mpfr_const_pi(r21166, MPFR_RNDN); mpfr_add(r21167, r21166, r21166, MPFR_RNDN); ; mpfr_div(r21169, r21167, r21168, MPFR_RNDN); mpfr_set_d(r21170, g, MPFR_RNDN); mpfr_neg(r21171, r21170, MPFR_RNDN); mpfr_set_d(r21172, h, MPFR_RNDN); mpfr_div(r21173, r21171, r21172, MPFR_RNDN); mpfr_acos(r21174, r21173, MPFR_RNDN); mpfr_div(r21175, r21174, r21168, MPFR_RNDN); mpfr_add(r21176, r21169, r21175, MPFR_RNDN); mpfr_cos(r21177, r21176, MPFR_RNDN); mpfr_mul(r21178, r21177, r21177, MPFR_RNDN); mpfr_cbrt(r21179, r21178, MPFR_RNDN); mpfr_div(r21180, r21166, r21168, MPFR_RNDN); mpfr_add(r21181, r21180, r21180, MPFR_RNDN); mpfr_div(r21182, r21170, r21172, MPFR_RNDN); mpfr_neg(r21183, r21182, MPFR_RNDN); mpfr_acos(r21184, r21183, MPFR_RNDN); mpfr_div(r21185, r21184, r21168, MPFR_RNDN); mpfr_add(r21186, r21181, r21185, MPFR_RNDN); mpfr_cos(r21187, r21186, MPFR_RNDN); mpfr_cbrt(r21188, r21187, MPFR_RNDN); mpfr_mul(r21189, r21179, r21188, MPFR_RNDN); mpfr_mul(r21190, r21165, r21189, MPFR_RNDN); return mpfr_get_d(r21190, MPFR_RNDN); } static mpfr_t r21191, r21192, r21193, r21194, r21195, r21196, r21197, r21198, r21199, r21200, r21201, r21202, r21203, r21204, r21205, r21206, r21207, r21208, r21209, r21210, r21211, r21212, r21213, r21214, r21215, r21216; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21191, "2", 10, MPFR_RNDN); mpfr_init(r21192); mpfr_init(r21193); mpfr_init_set_str(r21194, "3", 10, MPFR_RNDN); mpfr_init(r21195); mpfr_init(r21196); 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); } double f_dm(double g, double h) { ; mpfr_const_pi(r21192, MPFR_RNDN); mpfr_add(r21193, r21192, r21192, MPFR_RNDN); ; mpfr_div(r21195, r21193, r21194, MPFR_RNDN); mpfr_set_d(r21196, g, MPFR_RNDN); mpfr_neg(r21197, r21196, MPFR_RNDN); mpfr_set_d(r21198, h, MPFR_RNDN); mpfr_div(r21199, r21197, r21198, MPFR_RNDN); mpfr_acos(r21200, r21199, MPFR_RNDN); mpfr_div(r21201, r21200, r21194, MPFR_RNDN); mpfr_add(r21202, r21195, r21201, MPFR_RNDN); mpfr_cos(r21203, r21202, MPFR_RNDN); mpfr_mul(r21204, r21203, r21203, MPFR_RNDN); mpfr_cbrt(r21205, r21204, MPFR_RNDN); mpfr_div(r21206, r21192, r21194, MPFR_RNDN); mpfr_add(r21207, r21206, r21206, MPFR_RNDN); mpfr_div(r21208, r21196, r21198, MPFR_RNDN); mpfr_neg(r21209, r21208, MPFR_RNDN); mpfr_acos(r21210, r21209, MPFR_RNDN); mpfr_div(r21211, r21210, r21194, MPFR_RNDN); mpfr_add(r21212, r21207, r21211, MPFR_RNDN); mpfr_cos(r21213, r21212, MPFR_RNDN); mpfr_cbrt(r21214, r21213, MPFR_RNDN); mpfr_mul(r21215, r21205, r21214, MPFR_RNDN); mpfr_mul(r21216, r21191, r21215, MPFR_RNDN); return mpfr_get_d(r21216, MPFR_RNDN); }