#include #include #include #include #include char *name = "2-ancestry mixing, negative discriminant"; double f_if(float g, float h) { float r27054 = 2; float r27055 = atan2(1.0, 0.0); float r27056 = r27054 * r27055; float r27057 = 3; float r27058 = r27056 / r27057; float r27059 = g; float r27060 = -r27059; float r27061 = h; float r27062 = r27060 / r27061; float r27063 = acos(r27062); float r27064 = r27063 / r27057; float r27065 = r27058 + r27064; float r27066 = cos(r27065); float r27067 = r27054 * r27066; return r27067; } double f_id(double g, double h) { double r27068 = 2; double r27069 = atan2(1.0, 0.0); double r27070 = r27068 * r27069; double r27071 = 3; double r27072 = r27070 / r27071; double r27073 = g; double r27074 = -r27073; double r27075 = h; double r27076 = r27074 / r27075; double r27077 = acos(r27076); double r27078 = r27077 / r27071; double r27079 = r27072 + r27078; double r27080 = cos(r27079); double r27081 = r27068 * r27080; return r27081; } double f_of(float g, float h) { float r27082 = 2; float r27083 = atan2(1.0, 0.0); float r27084 = r27082 * r27083; float r27085 = 3; float r27086 = r27084 / r27085; float r27087 = cos(r27086); float r27088 = cbrt(r27087); float r27089 = r27088 * r27088; float r27090 = r27089 * r27088; float r27091 = g; float r27092 = -r27091; float r27093 = h; float r27094 = r27092 / r27093; float r27095 = acos(r27094); float r27096 = r27095 / r27085; float r27097 = cos(r27096); float r27098 = r27090 * r27097; float r27099 = sin(r27086); float r27100 = sin(r27096); float r27101 = r27099 * r27100; float r27102 = r27098 - r27101; float r27103 = r27082 * r27102; return r27103; } double f_od(double g, double h) { double r27104 = 2; double r27105 = atan2(1.0, 0.0); double r27106 = r27104 * r27105; double r27107 = 3; double r27108 = r27106 / r27107; double r27109 = cos(r27108); double r27110 = cbrt(r27109); double r27111 = r27110 * r27110; double r27112 = r27111 * r27110; double r27113 = g; double r27114 = -r27113; double r27115 = h; double r27116 = r27114 / r27115; double r27117 = acos(r27116); double r27118 = r27117 / r27107; double r27119 = cos(r27118); double r27120 = r27112 * r27119; double r27121 = sin(r27108); double r27122 = sin(r27118); double r27123 = r27121 * r27122; double r27124 = r27120 - r27123; double r27125 = r27104 * r27124; return r27125; } 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 r27126, r27127, r27128, r27129, r27130, r27131, r27132, r27133, r27134, r27135, r27136, r27137, r27138, r27139; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r27126, "2", 10, MPFR_RNDN); mpfr_init(r27127); mpfr_init(r27128); mpfr_init_set_str(r27129, "3", 10, MPFR_RNDN); mpfr_init(r27130); mpfr_init(r27131); mpfr_init(r27132); mpfr_init(r27133); mpfr_init(r27134); mpfr_init(r27135); mpfr_init(r27136); mpfr_init(r27137); mpfr_init(r27138); mpfr_init(r27139); } double f_im(double g, double h) { ; mpfr_const_pi(r27127, MPFR_RNDN); mpfr_mul(r27128, r27126, r27127, MPFR_RNDN); ; mpfr_div(r27130, r27128, r27129, MPFR_RNDN); mpfr_set_d(r27131, g, MPFR_RNDN); mpfr_neg(r27132, r27131, MPFR_RNDN); mpfr_set_d(r27133, h, MPFR_RNDN); mpfr_div(r27134, r27132, r27133, MPFR_RNDN); mpfr_acos(r27135, r27134, MPFR_RNDN); mpfr_div(r27136, r27135, r27129, MPFR_RNDN); mpfr_add(r27137, r27130, r27136, MPFR_RNDN); mpfr_cos(r27138, r27137, MPFR_RNDN); mpfr_mul(r27139, r27126, r27138, MPFR_RNDN); return mpfr_get_d(r27139, MPFR_RNDN); } static mpfr_t r27140, r27141, r27142, r27143, r27144, r27145, r27146, r27147, r27148, r27149, r27150, r27151, r27152, r27153, r27154, r27155, r27156, r27157, r27158, r27159, r27160, r27161; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27140, "2", 10, MPFR_RNDN); mpfr_init(r27141); mpfr_init(r27142); mpfr_init_set_str(r27143, "3", 10, MPFR_RNDN); mpfr_init(r27144); mpfr_init(r27145); mpfr_init(r27146); mpfr_init(r27147); mpfr_init(r27148); mpfr_init(r27149); mpfr_init(r27150); mpfr_init(r27151); mpfr_init(r27152); mpfr_init(r27153); mpfr_init(r27154); mpfr_init(r27155); mpfr_init(r27156); mpfr_init(r27157); mpfr_init(r27158); mpfr_init(r27159); mpfr_init(r27160); mpfr_init(r27161); } double f_fm(double g, double h) { ; mpfr_const_pi(r27141, MPFR_RNDN); mpfr_mul(r27142, r27140, r27141, MPFR_RNDN); ; mpfr_div(r27144, r27142, r27143, MPFR_RNDN); mpfr_cos(r27145, r27144, MPFR_RNDN); mpfr_cbrt(r27146, r27145, MPFR_RNDN); mpfr_mul(r27147, r27146, r27146, MPFR_RNDN); mpfr_mul(r27148, r27147, r27146, MPFR_RNDN); mpfr_set_d(r27149, g, MPFR_RNDN); mpfr_neg(r27150, r27149, MPFR_RNDN); mpfr_set_d(r27151, h, MPFR_RNDN); mpfr_div(r27152, r27150, r27151, MPFR_RNDN); mpfr_acos(r27153, r27152, MPFR_RNDN); mpfr_div(r27154, r27153, r27143, MPFR_RNDN); mpfr_cos(r27155, r27154, MPFR_RNDN); mpfr_mul(r27156, r27148, r27155, MPFR_RNDN); mpfr_sin(r27157, r27144, MPFR_RNDN); mpfr_sin(r27158, r27154, MPFR_RNDN); mpfr_mul(r27159, r27157, r27158, MPFR_RNDN); mpfr_sub(r27160, r27156, r27159, MPFR_RNDN); mpfr_mul(r27161, r27140, r27160, MPFR_RNDN); return mpfr_get_d(r27161, MPFR_RNDN); } static mpfr_t r27162, r27163, r27164, r27165, r27166, r27167, r27168, r27169, r27170, r27171, r27172, r27173, r27174, r27175, r27176, r27177, r27178, r27179, r27180, r27181, r27182, r27183; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27162, "2", 10, MPFR_RNDN); mpfr_init(r27163); mpfr_init(r27164); mpfr_init_set_str(r27165, "3", 10, MPFR_RNDN); 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); mpfr_init(r27175); mpfr_init(r27176); mpfr_init(r27177); mpfr_init(r27178); mpfr_init(r27179); mpfr_init(r27180); mpfr_init(r27181); mpfr_init(r27182); mpfr_init(r27183); } double f_dm(double g, double h) { ; mpfr_const_pi(r27163, MPFR_RNDN); mpfr_mul(r27164, r27162, r27163, MPFR_RNDN); ; mpfr_div(r27166, r27164, r27165, MPFR_RNDN); mpfr_cos(r27167, r27166, MPFR_RNDN); mpfr_cbrt(r27168, r27167, MPFR_RNDN); mpfr_mul(r27169, r27168, r27168, MPFR_RNDN); mpfr_mul(r27170, r27169, r27168, MPFR_RNDN); mpfr_set_d(r27171, g, MPFR_RNDN); mpfr_neg(r27172, r27171, MPFR_RNDN); mpfr_set_d(r27173, h, MPFR_RNDN); mpfr_div(r27174, r27172, r27173, MPFR_RNDN); mpfr_acos(r27175, r27174, MPFR_RNDN); mpfr_div(r27176, r27175, r27165, MPFR_RNDN); mpfr_cos(r27177, r27176, MPFR_RNDN); mpfr_mul(r27178, r27170, r27177, MPFR_RNDN); mpfr_sin(r27179, r27166, MPFR_RNDN); mpfr_sin(r27180, r27176, MPFR_RNDN); mpfr_mul(r27181, r27179, r27180, MPFR_RNDN); mpfr_sub(r27182, r27178, r27181, MPFR_RNDN); mpfr_mul(r27183, r27162, r27182, MPFR_RNDN); return mpfr_get_d(r27183, MPFR_RNDN); }