#include #include #include #include #include char *name = "2-ancestry mixing, negative discriminant"; double f_if(float g, float h) { float r27099 = 2; float r27100 = atan2(1.0, 0.0); float r27101 = r27099 * r27100; float r27102 = 3; float r27103 = r27101 / r27102; float r27104 = g; float r27105 = -r27104; float r27106 = h; float r27107 = r27105 / r27106; float r27108 = acos(r27107); float r27109 = r27108 / r27102; float r27110 = r27103 + r27109; float r27111 = cos(r27110); float r27112 = r27099 * r27111; return r27112; } double f_id(double g, double h) { double r27113 = 2; double r27114 = atan2(1.0, 0.0); double r27115 = r27113 * r27114; double r27116 = 3; double r27117 = r27115 / r27116; double r27118 = g; double r27119 = -r27118; double r27120 = h; double r27121 = r27119 / r27120; double r27122 = acos(r27121); double r27123 = r27122 / r27116; double r27124 = r27117 + r27123; double r27125 = cos(r27124); double r27126 = r27113 * r27125; return r27126; } double f_of(float g, float h) { float r27127 = 2; float r27128 = 3; float r27129 = r27127 / r27128; float r27130 = atan2(1.0, 0.0); float r27131 = g; float r27132 = h; float r27133 = r27131 / r27132; float r27134 = -r27133; float r27135 = acos(r27134); float r27136 = r27135 / r27128; float r27137 = fma(r27129, r27130, r27136); float r27138 = cos(r27137); float r27139 = r27138 * r27138; float r27140 = cbrt(r27139); float r27141 = cbrt(r27138); float r27142 = r27140 * r27141; float r27143 = r27127 * r27142; return r27143; } double f_od(double g, double h) { double r27144 = 2; double r27145 = 3; double r27146 = r27144 / r27145; double r27147 = atan2(1.0, 0.0); double r27148 = g; double r27149 = h; double r27150 = r27148 / r27149; double r27151 = -r27150; double r27152 = acos(r27151); double r27153 = r27152 / r27145; double r27154 = fma(r27146, r27147, r27153); double r27155 = cos(r27154); double r27156 = r27155 * r27155; double r27157 = cbrt(r27156); double r27158 = cbrt(r27155); double r27159 = r27157 * r27158; double r27160 = r27144 * r27159; return r27160; } 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 r27161, r27162, r27163, r27164, r27165, r27166, r27167, r27168, r27169, r27170, r27171, r27172, r27173, r27174; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r27161, "2", 10, MPFR_RNDN); mpfr_init(r27162); mpfr_init(r27163); mpfr_init_set_str(r27164, "3", 10, MPFR_RNDN); mpfr_init(r27165); mpfr_init(r27166); mpfr_init(r27167); mpfr_init(r27168); mpfr_init(r27169); mpfr_init(r27170); mpfr_init(r27171); mpfr_init(r27172); mpfr_init(r27173); mpfr_init(r27174); } double f_im(double g, double h) { ; mpfr_const_pi(r27162, MPFR_RNDN); mpfr_mul(r27163, r27161, r27162, MPFR_RNDN); ; mpfr_div(r27165, r27163, r27164, MPFR_RNDN); mpfr_set_d(r27166, g, MPFR_RNDN); mpfr_neg(r27167, r27166, MPFR_RNDN); mpfr_set_d(r27168, h, MPFR_RNDN); mpfr_div(r27169, r27167, r27168, MPFR_RNDN); mpfr_acos(r27170, r27169, MPFR_RNDN); mpfr_div(r27171, r27170, r27164, MPFR_RNDN); mpfr_add(r27172, r27165, r27171, MPFR_RNDN); mpfr_cos(r27173, r27172, MPFR_RNDN); mpfr_mul(r27174, r27161, r27173, MPFR_RNDN); return mpfr_get_d(r27174, MPFR_RNDN); } static mpfr_t r27175, r27176, r27177, r27178, r27179, r27180, r27181, r27182, r27183, r27184, r27185, r27186, r27187, r27188, r27189, r27190, r27191; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27175, "2", 10, MPFR_RNDN); mpfr_init_set_str(r27176, "3", 10, MPFR_RNDN); mpfr_init(r27177); mpfr_init(r27178); mpfr_init(r27179); mpfr_init(r27180); mpfr_init(r27181); mpfr_init(r27182); mpfr_init(r27183); mpfr_init(r27184); mpfr_init(r27185); mpfr_init(r27186); mpfr_init(r27187); mpfr_init(r27188); mpfr_init(r27189); mpfr_init(r27190); mpfr_init(r27191); } double f_fm(double g, double h) { ; ; mpfr_div(r27177, r27175, r27176, MPFR_RNDN); mpfr_const_pi(r27178, MPFR_RNDN); mpfr_set_d(r27179, g, MPFR_RNDN); mpfr_set_d(r27180, h, MPFR_RNDN); mpfr_div(r27181, r27179, r27180, MPFR_RNDN); mpfr_neg(r27182, r27181, MPFR_RNDN); mpfr_acos(r27183, r27182, MPFR_RNDN); mpfr_div(r27184, r27183, r27176, MPFR_RNDN); mpfr_fma(r27185, r27177, r27178, r27184, MPFR_RNDN); mpfr_cos(r27186, r27185, MPFR_RNDN); mpfr_mul(r27187, r27186, r27186, MPFR_RNDN); mpfr_cbrt(r27188, r27187, MPFR_RNDN); mpfr_cbrt(r27189, r27186, MPFR_RNDN); mpfr_mul(r27190, r27188, r27189, MPFR_RNDN); mpfr_mul(r27191, r27175, r27190, MPFR_RNDN); return mpfr_get_d(r27191, MPFR_RNDN); } static mpfr_t r27192, r27193, r27194, r27195, r27196, r27197, r27198, r27199, r27200, r27201, r27202, r27203, r27204, r27205, r27206, r27207, r27208; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27192, "2", 10, MPFR_RNDN); mpfr_init_set_str(r27193, "3", 10, MPFR_RNDN); mpfr_init(r27194); mpfr_init(r27195); mpfr_init(r27196); mpfr_init(r27197); mpfr_init(r27198); mpfr_init(r27199); mpfr_init(r27200); mpfr_init(r27201); mpfr_init(r27202); mpfr_init(r27203); mpfr_init(r27204); mpfr_init(r27205); mpfr_init(r27206); mpfr_init(r27207); mpfr_init(r27208); } double f_dm(double g, double h) { ; ; mpfr_div(r27194, r27192, r27193, MPFR_RNDN); mpfr_const_pi(r27195, MPFR_RNDN); mpfr_set_d(r27196, g, MPFR_RNDN); mpfr_set_d(r27197, h, MPFR_RNDN); mpfr_div(r27198, r27196, r27197, MPFR_RNDN); mpfr_neg(r27199, r27198, MPFR_RNDN); mpfr_acos(r27200, r27199, MPFR_RNDN); mpfr_div(r27201, r27200, r27193, MPFR_RNDN); mpfr_fma(r27202, r27194, r27195, r27201, MPFR_RNDN); mpfr_cos(r27203, r27202, MPFR_RNDN); mpfr_mul(r27204, r27203, r27203, MPFR_RNDN); mpfr_cbrt(r27205, r27204, MPFR_RNDN); mpfr_cbrt(r27206, r27203, MPFR_RNDN); mpfr_mul(r27207, r27205, r27206, MPFR_RNDN); mpfr_mul(r27208, r27192, r27207, MPFR_RNDN); return mpfr_get_d(r27208, MPFR_RNDN); }