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

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

double f_if(float g, float a) {
        float r26972 = g;
        float r26973 = 2;
        float r26974 = a;
        float r26975 = r26973 * r26974;
        float r26976 = r26972 / r26975;
        float r26977 = cbrt(r26976);
        return r26977;
}

double f_id(double g, double a) {
        double r26978 = g;
        double r26979 = 2;
        double r26980 = a;
        double r26981 = r26979 * r26980;
        double r26982 = r26978 / r26981;
        double r26983 = cbrt(r26982);
        return r26983;
}


double f_of(float g, float a) {
        float r26984 = g;
        float r26985 = cbrt(r26984);
        float r26986 = 1;
        float r26987 = 2;
        float r26988 = a;
        float r26989 = r26987 * r26988;
        float r26990 = r26986 / r26989;
        float r26991 = cbrt(r26990);
        float r26992 = r26985 * r26991;
        return r26992;
}

double f_od(double g, double a) {
        double r26993 = g;
        double r26994 = cbrt(r26993);
        double r26995 = 1;
        double r26996 = 2;
        double r26997 = a;
        double r26998 = r26996 * r26997;
        double r26999 = r26995 / r26998;
        double r27000 = cbrt(r26999);
        double r27001 = r26994 * r27000;
        return r27001;
}

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 r27002, r27003, r27004, r27005, r27006, r27007;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r27002);
        mpfr_init_set_str(r27003, "2", 10, MPFR_RNDN);
        mpfr_init(r27004);
        mpfr_init(r27005);
        mpfr_init(r27006);
        mpfr_init(r27007);
}

double f_im(double g, double a) {
        mpfr_set_d(r27002, g, MPFR_RNDN);
        ;
        mpfr_set_d(r27004, a, MPFR_RNDN);
        mpfr_mul(r27005, r27003, r27004, MPFR_RNDN);
        mpfr_div(r27006, r27002, r27005, MPFR_RNDN);
        mpfr_cbrt(r27007, r27006, MPFR_RNDN);
        return mpfr_get_d(r27007, MPFR_RNDN);
}

static mpfr_t r27008, r27009, r27010, r27011, r27012, r27013, r27014, r27015, r27016;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27008);
        mpfr_init(r27009);
        mpfr_init_set_str(r27010, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27011, "2", 10, MPFR_RNDN);
        mpfr_init(r27012);
        mpfr_init(r27013);
        mpfr_init(r27014);
        mpfr_init(r27015);
        mpfr_init(r27016);
}

double f_fm(double g, double a) {
        mpfr_set_d(r27008, g, MPFR_RNDN);
        mpfr_cbrt(r27009, r27008, MPFR_RNDN);
        ;
        ;
        mpfr_set_d(r27012, a, MPFR_RNDN);
        mpfr_mul(r27013, r27011, r27012, MPFR_RNDN);
        mpfr_div(r27014, r27010, r27013, MPFR_RNDN);
        mpfr_cbrt(r27015, r27014, MPFR_RNDN);
        mpfr_mul(r27016, r27009, r27015, MPFR_RNDN);
        return mpfr_get_d(r27016, MPFR_RNDN);
}

static mpfr_t r27017, r27018, r27019, r27020, r27021, r27022, r27023, r27024, r27025;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27017);
        mpfr_init(r27018);
        mpfr_init_set_str(r27019, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27020, "2", 10, MPFR_RNDN);
        mpfr_init(r27021);
        mpfr_init(r27022);
        mpfr_init(r27023);
        mpfr_init(r27024);
        mpfr_init(r27025);
}

double f_dm(double g, double a) {
        mpfr_set_d(r27017, g, MPFR_RNDN);
        mpfr_cbrt(r27018, r27017, MPFR_RNDN);
        ;
        ;
        mpfr_set_d(r27021, a, MPFR_RNDN);
        mpfr_mul(r27022, r27020, r27021, MPFR_RNDN);
        mpfr_div(r27023, r27019, r27022, MPFR_RNDN);
        mpfr_cbrt(r27024, r27023, MPFR_RNDN);
        mpfr_mul(r27025, r27018, r27024, MPFR_RNDN);
        return mpfr_get_d(r27025, MPFR_RNDN);
}