#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 r26753 = g;
        float r26754 = 2;
        float r26755 = a;
        float r26756 = r26754 * r26755;
        float r26757 = r26753 / r26756;
        float r26758 = cbrt(r26757);
        return r26758;
}

double f_id(double g, double a) {
        double r26759 = g;
        double r26760 = 2;
        double r26761 = a;
        double r26762 = r26760 * r26761;
        double r26763 = r26759 / r26762;
        double r26764 = cbrt(r26763);
        return r26764;
}


double f_of(float g, float a) {
        float r26765 = g;
        float r26766 = cbrt(r26765);
        float r26767 = 1;
        float r26768 = 2;
        float r26769 = a;
        float r26770 = r26768 * r26769;
        float r26771 = r26767 / r26770;
        float r26772 = cbrt(r26771);
        float r26773 = r26766 * r26772;
        return r26773;
}

double f_od(double g, double a) {
        double r26774 = g;
        double r26775 = cbrt(r26774);
        double r26776 = 1;
        double r26777 = 2;
        double r26778 = a;
        double r26779 = r26777 * r26778;
        double r26780 = r26776 / r26779;
        double r26781 = cbrt(r26780);
        double r26782 = r26775 * r26781;
        return r26782;
}

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 r26783, r26784, r26785, r26786, r26787, r26788;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26783);
        mpfr_init_set_str(r26784, "2", 10, MPFR_RNDN);
        mpfr_init(r26785);
        mpfr_init(r26786);
        mpfr_init(r26787);
        mpfr_init(r26788);
}

double f_im(double g, double a) {
        mpfr_set_d(r26783, g, MPFR_RNDN);
        ;
        mpfr_set_d(r26785, a, MPFR_RNDN);
        mpfr_mul(r26786, r26784, r26785, MPFR_RNDN);
        mpfr_div(r26787, r26783, r26786, MPFR_RNDN);
        mpfr_cbrt(r26788, r26787, MPFR_RNDN);
        return mpfr_get_d(r26788, MPFR_RNDN);
}

static mpfr_t r26789, r26790, r26791, r26792, r26793, r26794, r26795, r26796, r26797;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26789);
        mpfr_init(r26790);
        mpfr_init_set_str(r26791, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r26792, "2", 10, MPFR_RNDN);
        mpfr_init(r26793);
        mpfr_init(r26794);
        mpfr_init(r26795);
        mpfr_init(r26796);
        mpfr_init(r26797);
}

double f_fm(double g, double a) {
        mpfr_set_d(r26789, g, MPFR_RNDN);
        mpfr_cbrt(r26790, r26789, MPFR_RNDN);
        ;
        ;
        mpfr_set_d(r26793, a, MPFR_RNDN);
        mpfr_mul(r26794, r26792, r26793, MPFR_RNDN);
        mpfr_div(r26795, r26791, r26794, MPFR_RNDN);
        mpfr_cbrt(r26796, r26795, MPFR_RNDN);
        mpfr_mul(r26797, r26790, r26796, MPFR_RNDN);
        return mpfr_get_d(r26797, MPFR_RNDN);
}

static mpfr_t r26798, r26799, r26800, r26801, r26802, r26803, r26804, r26805, r26806;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26798);
        mpfr_init(r26799);
        mpfr_init_set_str(r26800, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r26801, "2", 10, MPFR_RNDN);
        mpfr_init(r26802);
        mpfr_init(r26803);
        mpfr_init(r26804);
        mpfr_init(r26805);
        mpfr_init(r26806);
}

double f_dm(double g, double a) {
        mpfr_set_d(r26798, g, MPFR_RNDN);
        mpfr_cbrt(r26799, r26798, MPFR_RNDN);
        ;
        ;
        mpfr_set_d(r26802, a, MPFR_RNDN);
        mpfr_mul(r26803, r26801, r26802, MPFR_RNDN);
        mpfr_div(r26804, r26800, r26803, MPFR_RNDN);
        mpfr_cbrt(r26805, r26804, MPFR_RNDN);
        mpfr_mul(r26806, r26799, r26805, MPFR_RNDN);
        return mpfr_get_d(r26806, MPFR_RNDN);
}