#include #include #include #include #include char *name = "2-ancestry mixing, negative discriminant"; double f_if(float g, float h) { float r23800 = 2; float r23801 = atan2(1.0, 0.0); float r23802 = r23800 * r23801; float r23803 = 3; float r23804 = r23802 / r23803; float r23805 = g; float r23806 = -r23805; float r23807 = h; float r23808 = r23806 / r23807; float r23809 = acos(r23808); float r23810 = r23809 / r23803; float r23811 = r23804 + r23810; float r23812 = cos(r23811); float r23813 = r23800 * r23812; return r23813; } double f_id(double g, double h) { double r23814 = 2; double r23815 = atan2(1.0, 0.0); double r23816 = r23814 * r23815; double r23817 = 3; double r23818 = r23816 / r23817; double r23819 = g; double r23820 = -r23819; double r23821 = h; double r23822 = r23820 / r23821; double r23823 = acos(r23822); double r23824 = r23823 / r23817; double r23825 = r23818 + r23824; double r23826 = cos(r23825); double r23827 = r23814 * r23826; return r23827; } double f_of(float g, float h) { float r23828 = 2; float r23829 = 3; float r23830 = r23828 / r23829; float r23831 = atan2(1.0, 0.0); float r23832 = g; float r23833 = h; float r23834 = r23832 / r23833; float r23835 = -r23834; float r23836 = acos(r23835); float r23837 = r23836 / r23829; float r23838 = fma(r23830, r23831, r23837); float r23839 = cos(r23838); float r23840 = r23839 * r23839; float r23841 = cbrt(r23840); float r23842 = cbrt(r23839); float r23843 = r23841 * r23842; float r23844 = r23828 * r23843; return r23844; } double f_od(double g, double h) { double r23845 = 2; double r23846 = 3; double r23847 = r23845 / r23846; double r23848 = atan2(1.0, 0.0); double r23849 = g; double r23850 = h; double r23851 = r23849 / r23850; double r23852 = -r23851; double r23853 = acos(r23852); double r23854 = r23853 / r23846; double r23855 = fma(r23847, r23848, r23854); double r23856 = cos(r23855); double r23857 = r23856 * r23856; double r23858 = cbrt(r23857); double r23859 = cbrt(r23856); double r23860 = r23858 * r23859; double r23861 = r23845 * r23860; return r23861; } 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 r23862, r23863, r23864, r23865, r23866, r23867, r23868, r23869, r23870, r23871, r23872, r23873, r23874, r23875; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r23862, "2", 10, MPFR_RNDN); mpfr_init(r23863); mpfr_init(r23864); mpfr_init_set_str(r23865, "3", 10, MPFR_RNDN); mpfr_init(r23866); mpfr_init(r23867); mpfr_init(r23868); mpfr_init(r23869); mpfr_init(r23870); mpfr_init(r23871); mpfr_init(r23872); mpfr_init(r23873); mpfr_init(r23874); mpfr_init(r23875); } double f_im(double g, double h) { ; mpfr_const_pi(r23863, MPFR_RNDN); mpfr_mul(r23864, r23862, r23863, MPFR_RNDN); ; mpfr_div(r23866, r23864, r23865, MPFR_RNDN); mpfr_set_d(r23867, g, MPFR_RNDN); mpfr_neg(r23868, r23867, MPFR_RNDN); mpfr_set_d(r23869, h, MPFR_RNDN); mpfr_div(r23870, r23868, r23869, MPFR_RNDN); mpfr_acos(r23871, r23870, MPFR_RNDN); mpfr_div(r23872, r23871, r23865, MPFR_RNDN); mpfr_add(r23873, r23866, r23872, MPFR_RNDN); mpfr_cos(r23874, r23873, MPFR_RNDN); mpfr_mul(r23875, r23862, r23874, MPFR_RNDN); return mpfr_get_d(r23875, MPFR_RNDN); } static mpfr_t r23876, r23877, r23878, r23879, r23880, r23881, r23882, r23883, r23884, r23885, r23886, r23887, r23888, r23889, r23890, r23891, r23892; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23876, "2", 10, MPFR_RNDN); mpfr_init_set_str(r23877, "3", 10, MPFR_RNDN); mpfr_init(r23878); mpfr_init(r23879); mpfr_init(r23880); mpfr_init(r23881); mpfr_init(r23882); mpfr_init(r23883); mpfr_init(r23884); mpfr_init(r23885); mpfr_init(r23886); mpfr_init(r23887); mpfr_init(r23888); mpfr_init(r23889); mpfr_init(r23890); mpfr_init(r23891); mpfr_init(r23892); } double f_fm(double g, double h) { ; ; mpfr_div(r23878, r23876, r23877, MPFR_RNDN); mpfr_const_pi(r23879, MPFR_RNDN); mpfr_set_d(r23880, g, MPFR_RNDN); mpfr_set_d(r23881, h, MPFR_RNDN); mpfr_div(r23882, r23880, r23881, MPFR_RNDN); mpfr_neg(r23883, r23882, MPFR_RNDN); mpfr_acos(r23884, r23883, MPFR_RNDN); mpfr_div(r23885, r23884, r23877, MPFR_RNDN); mpfr_fma(r23886, r23878, r23879, r23885, MPFR_RNDN); mpfr_cos(r23887, r23886, MPFR_RNDN); mpfr_mul(r23888, r23887, r23887, MPFR_RNDN); mpfr_cbrt(r23889, r23888, MPFR_RNDN); mpfr_cbrt(r23890, r23887, MPFR_RNDN); mpfr_mul(r23891, r23889, r23890, MPFR_RNDN); mpfr_mul(r23892, r23876, r23891, MPFR_RNDN); return mpfr_get_d(r23892, MPFR_RNDN); } static mpfr_t r23893, r23894, r23895, r23896, r23897, r23898, r23899, r23900, r23901, r23902, r23903, r23904, r23905, r23906, r23907, r23908, r23909; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23893, "2", 10, MPFR_RNDN); mpfr_init_set_str(r23894, "3", 10, MPFR_RNDN); mpfr_init(r23895); mpfr_init(r23896); mpfr_init(r23897); mpfr_init(r23898); mpfr_init(r23899); mpfr_init(r23900); mpfr_init(r23901); mpfr_init(r23902); mpfr_init(r23903); mpfr_init(r23904); mpfr_init(r23905); mpfr_init(r23906); mpfr_init(r23907); mpfr_init(r23908); mpfr_init(r23909); } double f_dm(double g, double h) { ; ; mpfr_div(r23895, r23893, r23894, MPFR_RNDN); mpfr_const_pi(r23896, MPFR_RNDN); mpfr_set_d(r23897, g, MPFR_RNDN); mpfr_set_d(r23898, h, MPFR_RNDN); mpfr_div(r23899, r23897, r23898, MPFR_RNDN); mpfr_neg(r23900, r23899, MPFR_RNDN); mpfr_acos(r23901, r23900, MPFR_RNDN); mpfr_div(r23902, r23901, r23894, MPFR_RNDN); mpfr_fma(r23903, r23895, r23896, r23902, MPFR_RNDN); mpfr_cos(r23904, r23903, MPFR_RNDN); mpfr_mul(r23905, r23904, r23904, MPFR_RNDN); mpfr_cbrt(r23906, r23905, MPFR_RNDN); mpfr_cbrt(r23907, r23904, MPFR_RNDN); mpfr_mul(r23908, r23906, r23907, MPFR_RNDN); mpfr_mul(r23909, r23893, r23908, MPFR_RNDN); return mpfr_get_d(r23909, MPFR_RNDN); }