#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "2-ancestry mixing, negative discriminant";

double f_if(float g, float h) {
        float r26862 = 2;
        float r26863 = atan2(1.0, 0.0);
        float r26864 = r26862 * r26863;
        float r26865 = 3;
        float r26866 = r26864 / r26865;
        float r26867 = g;
        float r26868 = -r26867;
        float r26869 = h;
        float r26870 = r26868 / r26869;
        float r26871 = acos(r26870);
        float r26872 = r26871 / r26865;
        float r26873 = r26866 + r26872;
        float r26874 = cos(r26873);
        float r26875 = r26862 * r26874;
        return r26875;
}

double f_id(double g, double h) {
        double r26876 = 2;
        double r26877 = atan2(1.0, 0.0);
        double r26878 = r26876 * r26877;
        double r26879 = 3;
        double r26880 = r26878 / r26879;
        double r26881 = g;
        double r26882 = -r26881;
        double r26883 = h;
        double r26884 = r26882 / r26883;
        double r26885 = acos(r26884);
        double r26886 = r26885 / r26879;
        double r26887 = r26880 + r26886;
        double r26888 = cos(r26887);
        double r26889 = r26876 * r26888;
        return r26889;
}


double f_of(float g, float h) {
        float r26890 = 2;
        float r26891 = g;
        float r26892 = h;
        float r26893 = r26891 / r26892;
        float r26894 = acos(r26893);
        float r26895 = 3;
        float r26896 = r26894 / r26895;
        float r26897 = cos(r26896);
        float r26898 = 1;
        float r26899 = r26890 + r26898;
        float r26900 = atan2(1.0, 0.0);
        float r26901 = r26900 / r26895;
        float r26902 = r26899 * r26901;
        float r26903 = cos(r26902);
        float r26904 = r26897 * r26903;
        float r26905 = sin(r26896);
        float r26906 = sin(r26902);
        float r26907 = r26905 * r26906;
        float r26908 = r26904 + r26907;
        float r26909 = r26890 * r26908;
        return r26909;
}

double f_od(double g, double h) {
        double r26910 = 2;
        double r26911 = g;
        double r26912 = h;
        double r26913 = r26911 / r26912;
        double r26914 = acos(r26913);
        double r26915 = 3;
        double r26916 = r26914 / r26915;
        double r26917 = cos(r26916);
        double r26918 = 1;
        double r26919 = r26910 + r26918;
        double r26920 = atan2(1.0, 0.0);
        double r26921 = r26920 / r26915;
        double r26922 = r26919 * r26921;
        double r26923 = cos(r26922);
        double r26924 = r26917 * r26923;
        double r26925 = sin(r26916);
        double r26926 = sin(r26922);
        double r26927 = r26925 * r26926;
        double r26928 = r26924 + r26927;
        double r26929 = r26910 * r26928;
        return r26929;
}

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 r26930, r26931, r26932, r26933, r26934, r26935, r26936, r26937, r26938, r26939, r26940, r26941, r26942, r26943;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r26930, "2", 10, MPFR_RNDN);
        mpfr_init(r26931);
        mpfr_init(r26932);
        mpfr_init_set_str(r26933, "3", 10, MPFR_RNDN);
        mpfr_init(r26934);
        mpfr_init(r26935);
        mpfr_init(r26936);
        mpfr_init(r26937);
        mpfr_init(r26938);
        mpfr_init(r26939);
        mpfr_init(r26940);
        mpfr_init(r26941);
        mpfr_init(r26942);
        mpfr_init(r26943);
}

double f_im(double g, double h) {
        ;
        mpfr_const_pi(r26931, MPFR_RNDN);
        mpfr_mul(r26932, r26930, r26931, MPFR_RNDN);
        ;
        mpfr_div(r26934, r26932, r26933, MPFR_RNDN);
        mpfr_set_d(r26935, g, MPFR_RNDN);
        mpfr_neg(r26936, r26935, MPFR_RNDN);
        mpfr_set_d(r26937, h, MPFR_RNDN);
        mpfr_div(r26938, r26936, r26937, MPFR_RNDN);
        mpfr_acos(r26939, r26938, MPFR_RNDN);
        mpfr_div(r26940, r26939, r26933, MPFR_RNDN);
        mpfr_add(r26941, r26934, r26940, MPFR_RNDN);
        mpfr_cos(r26942, r26941, MPFR_RNDN);
        mpfr_mul(r26943, r26930, r26942, MPFR_RNDN);
        return mpfr_get_d(r26943, MPFR_RNDN);
}

static mpfr_t r26944, r26945, r26946, r26947, r26948, r26949, r26950, r26951, r26952, r26953, r26954, r26955, r26956, r26957, r26958, r26959, r26960, r26961, r26962, r26963;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r26944, "2", 10, MPFR_RNDN);
        mpfr_init(r26945);
        mpfr_init(r26946);
        mpfr_init(r26947);
        mpfr_init(r26948);
        mpfr_init_set_str(r26949, "3", 10, MPFR_RNDN);
        mpfr_init(r26950);
        mpfr_init(r26951);
        mpfr_init_set_str(r26952, "1", 10, MPFR_RNDN);
        mpfr_init(r26953);
        mpfr_init(r26954);
        mpfr_init(r26955);
        mpfr_init(r26956);
        mpfr_init(r26957);
        mpfr_init(r26958);
        mpfr_init(r26959);
        mpfr_init(r26960);
        mpfr_init(r26961);
        mpfr_init(r26962);
        mpfr_init(r26963);
}

double f_fm(double g, double h) {
        ;
        mpfr_set_d(r26945, g, MPFR_RNDN);
        mpfr_set_d(r26946, h, MPFR_RNDN);
        mpfr_div(r26947, r26945, r26946, MPFR_RNDN);
        mpfr_acos(r26948, r26947, MPFR_RNDN);
        ;
        mpfr_div(r26950, r26948, r26949, MPFR_RNDN);
        mpfr_cos(r26951, r26950, MPFR_RNDN);
        ;
        mpfr_add(r26953, r26944, r26952, MPFR_RNDN);
        mpfr_const_pi(r26954, MPFR_RNDN);
        mpfr_div(r26955, r26954, r26949, MPFR_RNDN);
        mpfr_mul(r26956, r26953, r26955, MPFR_RNDN);
        mpfr_cos(r26957, r26956, MPFR_RNDN);
        mpfr_mul(r26958, r26951, r26957, MPFR_RNDN);
        mpfr_sin(r26959, r26950, MPFR_RNDN);
        mpfr_sin(r26960, r26956, MPFR_RNDN);
        mpfr_mul(r26961, r26959, r26960, MPFR_RNDN);
        mpfr_add(r26962, r26958, r26961, MPFR_RNDN);
        mpfr_mul(r26963, r26944, r26962, MPFR_RNDN);
        return mpfr_get_d(r26963, MPFR_RNDN);
}

static mpfr_t r26964, r26965, r26966, r26967, r26968, r26969, r26970, r26971, r26972, r26973, r26974, r26975, r26976, r26977, r26978, r26979, r26980, r26981, r26982, r26983;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r26964, "2", 10, MPFR_RNDN);
        mpfr_init(r26965);
        mpfr_init(r26966);
        mpfr_init(r26967);
        mpfr_init(r26968);
        mpfr_init_set_str(r26969, "3", 10, MPFR_RNDN);
        mpfr_init(r26970);
        mpfr_init(r26971);
        mpfr_init_set_str(r26972, "1", 10, MPFR_RNDN);
        mpfr_init(r26973);
        mpfr_init(r26974);
        mpfr_init(r26975);
        mpfr_init(r26976);
        mpfr_init(r26977);
        mpfr_init(r26978);
        mpfr_init(r26979);
        mpfr_init(r26980);
        mpfr_init(r26981);
        mpfr_init(r26982);
        mpfr_init(r26983);
}

double f_dm(double g, double h) {
        ;
        mpfr_set_d(r26965, g, MPFR_RNDN);
        mpfr_set_d(r26966, h, MPFR_RNDN);
        mpfr_div(r26967, r26965, r26966, MPFR_RNDN);
        mpfr_acos(r26968, r26967, MPFR_RNDN);
        ;
        mpfr_div(r26970, r26968, r26969, MPFR_RNDN);
        mpfr_cos(r26971, r26970, MPFR_RNDN);
        ;
        mpfr_add(r26973, r26964, r26972, MPFR_RNDN);
        mpfr_const_pi(r26974, MPFR_RNDN);
        mpfr_div(r26975, r26974, r26969, MPFR_RNDN);
        mpfr_mul(r26976, r26973, r26975, MPFR_RNDN);
        mpfr_cos(r26977, r26976, MPFR_RNDN);
        mpfr_mul(r26978, r26971, r26977, MPFR_RNDN);
        mpfr_sin(r26979, r26970, MPFR_RNDN);
        mpfr_sin(r26980, r26976, MPFR_RNDN);
        mpfr_mul(r26981, r26979, r26980, MPFR_RNDN);
        mpfr_add(r26982, r26978, r26981, MPFR_RNDN);
        mpfr_mul(r26983, r26964, r26982, MPFR_RNDN);
        return mpfr_get_d(r26983, MPFR_RNDN);
}