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

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

double f_if(float g, float h, float a) {
        float r10516 = 1.0f;
        float r10517 = 2.0f;
        float r10518 = a;
        float r10519 = r10517 * r10518;
        float r10520 = r10516 / r10519;
        float r10521 = g;
        float r10522 = -r10521;
        float r10523 = r10521 * r10521;
        float r10524 = h;
        float r10525 = r10524 * r10524;
        float r10526 = r10523 - r10525;
        float r10527 = sqrt(r10526);
        float r10528 = r10522 + r10527;
        float r10529 = r10520 * r10528;
        float r10530 = cbrt(r10529);
        float r10531 = r10522 - r10527;
        float r10532 = r10520 * r10531;
        float r10533 = cbrt(r10532);
        float r10534 = r10530 + r10533;
        return r10534;
}

double f_id(double g, double h, double a) {
        double r10535 = 1.0;
        double r10536 = 2.0;
        double r10537 = a;
        double r10538 = r10536 * r10537;
        double r10539 = r10535 / r10538;
        double r10540 = g;
        double r10541 = -r10540;
        double r10542 = r10540 * r10540;
        double r10543 = h;
        double r10544 = r10543 * r10543;
        double r10545 = r10542 - r10544;
        double r10546 = sqrt(r10545);
        double r10547 = r10541 + r10546;
        double r10548 = r10539 * r10547;
        double r10549 = cbrt(r10548);
        double r10550 = r10541 - r10546;
        double r10551 = r10539 * r10550;
        double r10552 = cbrt(r10551);
        double r10553 = r10549 + r10552;
        return r10553;
}


double f_of(float g, float h, float a) {
        float r10554 = 1.0f;
        float r10555 = a;
        float r10556 = 2.0f;
        float r10557 = r10555 * r10556;
        float r10558 = r10554 / r10557;
        float r10559 = g;
        float r10560 = -r10559;
        float r10561 = r10559 * r10559;
        float r10562 = h;
        float r10563 = r10562 * r10562;
        float r10564 = r10561 - r10563;
        float r10565 = sqrt(r10564);
        float r10566 = r10560 - r10565;
        float r10567 = r10558 * r10566;
        float r10568 = cbrt(r10567);
        float r10569 = cbrt(r10564);
        float r10570 = fabs(r10569);
        float r10571 = sqrt(r10569);
        float r10572 = r10570 * r10571;
        float r10573 = r10560 + r10572;
        float r10574 = cbrt(r10573);
        float r10575 = r10560 + r10565;
        float r10576 = cbrt(r10575);
        float r10577 = r10576 * r10574;
        float r10578 = r10574 * r10577;
        float r10579 = r10578 * r10558;
        float r10580 = cbrt(r10579);
        float r10581 = r10568 + r10580;
        return r10581;
}

double f_od(double g, double h, double a) {
        double r10582 = 1.0;
        double r10583 = a;
        double r10584 = 2.0;
        double r10585 = r10583 * r10584;
        double r10586 = r10582 / r10585;
        double r10587 = g;
        double r10588 = -r10587;
        double r10589 = r10587 * r10587;
        double r10590 = h;
        double r10591 = r10590 * r10590;
        double r10592 = r10589 - r10591;
        double r10593 = sqrt(r10592);
        double r10594 = r10588 - r10593;
        double r10595 = r10586 * r10594;
        double r10596 = cbrt(r10595);
        double r10597 = cbrt(r10592);
        double r10598 = fabs(r10597);
        double r10599 = sqrt(r10597);
        double r10600 = r10598 * r10599;
        double r10601 = r10588 + r10600;
        double r10602 = cbrt(r10601);
        double r10603 = r10588 + r10593;
        double r10604 = cbrt(r10603);
        double r10605 = r10604 * r10602;
        double r10606 = r10602 * r10605;
        double r10607 = r10606 * r10586;
        double r10608 = cbrt(r10607);
        double r10609 = r10596 + r10608;
        return r10609;
}

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 r10610, r10611, r10612, r10613, r10614, r10615, r10616, r10617, r10618, r10619, r10620, r10621, r10622, r10623, r10624, r10625, r10626, r10627, r10628;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r10610, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r10611, "2", 10, MPFR_RNDN);
        mpfr_init(r10612);
        mpfr_init(r10613);
        mpfr_init(r10614);
        mpfr_init(r10615);
        mpfr_init(r10616);
        mpfr_init(r10617);
        mpfr_init(r10618);
        mpfr_init(r10619);
        mpfr_init(r10620);
        mpfr_init(r10621);
        mpfr_init(r10622);
        mpfr_init(r10623);
        mpfr_init(r10624);
        mpfr_init(r10625);
        mpfr_init(r10626);
        mpfr_init(r10627);
        mpfr_init(r10628);
}

double f_im(double g, double h, double a) {
        ;
        ;
        mpfr_set_d(r10612, a, MPFR_RNDN);
        mpfr_mul(r10613, r10611, r10612, MPFR_RNDN);
        mpfr_div(r10614, r10610, r10613, MPFR_RNDN);
        mpfr_set_d(r10615, g, MPFR_RNDN);
        mpfr_neg(r10616, r10615, MPFR_RNDN);
        mpfr_mul(r10617, r10615, r10615, MPFR_RNDN);
        mpfr_set_d(r10618, h, MPFR_RNDN);
        mpfr_mul(r10619, r10618, r10618, MPFR_RNDN);
        mpfr_sub(r10620, r10617, r10619, MPFR_RNDN);
        mpfr_sqrt(r10621, r10620, MPFR_RNDN);
        mpfr_add(r10622, r10616, r10621, MPFR_RNDN);
        mpfr_mul(r10623, r10614, r10622, MPFR_RNDN);
        mpfr_cbrt(r10624, r10623, MPFR_RNDN);
        mpfr_sub(r10625, r10616, r10621, MPFR_RNDN);
        mpfr_mul(r10626, r10614, r10625, MPFR_RNDN);
        mpfr_cbrt(r10627, r10626, MPFR_RNDN);
        mpfr_add(r10628, r10624, r10627, MPFR_RNDN);
        return mpfr_get_d(r10628, MPFR_RNDN);
}

static mpfr_t r10629, r10630, r10631, r10632, r10633, r10634, r10635, r10636, r10637, r10638, r10639, r10640, r10641, r10642, r10643, r10644, r10645, r10646, r10647, r10648, r10649, r10650, r10651, r10652, r10653, r10654, r10655, r10656;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r10629, "1", 10, MPFR_RNDN);
        mpfr_init(r10630);
        mpfr_init_set_str(r10631, "2", 10, MPFR_RNDN);
        mpfr_init(r10632);
        mpfr_init(r10633);
        mpfr_init(r10634);
        mpfr_init(r10635);
        mpfr_init(r10636);
        mpfr_init(r10637);
        mpfr_init(r10638);
        mpfr_init(r10639);
        mpfr_init(r10640);
        mpfr_init(r10641);
        mpfr_init(r10642);
        mpfr_init(r10643);
        mpfr_init(r10644);
        mpfr_init(r10645);
        mpfr_init(r10646);
        mpfr_init(r10647);
        mpfr_init(r10648);
        mpfr_init(r10649);
        mpfr_init(r10650);
        mpfr_init(r10651);
        mpfr_init(r10652);
        mpfr_init(r10653);
        mpfr_init(r10654);
        mpfr_init(r10655);
        mpfr_init(r10656);
}

double f_fm(double g, double h, double a) {
        ;
        mpfr_set_d(r10630, a, MPFR_RNDN);
        ;
        mpfr_mul(r10632, r10630, r10631, MPFR_RNDN);
        mpfr_div(r10633, r10629, r10632, MPFR_RNDN);
        mpfr_set_d(r10634, g, MPFR_RNDN);
        mpfr_neg(r10635, r10634, MPFR_RNDN);
        mpfr_mul(r10636, r10634, r10634, MPFR_RNDN);
        mpfr_set_d(r10637, h, MPFR_RNDN);
        mpfr_mul(r10638, r10637, r10637, MPFR_RNDN);
        mpfr_sub(r10639, r10636, r10638, MPFR_RNDN);
        mpfr_sqrt(r10640, r10639, MPFR_RNDN);
        mpfr_sub(r10641, r10635, r10640, MPFR_RNDN);
        mpfr_mul(r10642, r10633, r10641, MPFR_RNDN);
        mpfr_cbrt(r10643, r10642, MPFR_RNDN);
        mpfr_cbrt(r10644, r10639, MPFR_RNDN);
        mpfr_abs(r10645, r10644, MPFR_RNDN);
        mpfr_sqrt(r10646, r10644, MPFR_RNDN);
        mpfr_mul(r10647, r10645, r10646, MPFR_RNDN);
        mpfr_add(r10648, r10635, r10647, MPFR_RNDN);
        mpfr_cbrt(r10649, r10648, MPFR_RNDN);
        mpfr_add(r10650, r10635, r10640, MPFR_RNDN);
        mpfr_cbrt(r10651, r10650, MPFR_RNDN);
        mpfr_mul(r10652, r10651, r10649, MPFR_RNDN);
        mpfr_mul(r10653, r10649, r10652, MPFR_RNDN);
        mpfr_mul(r10654, r10653, r10633, MPFR_RNDN);
        mpfr_cbrt(r10655, r10654, MPFR_RNDN);
        mpfr_add(r10656, r10643, r10655, MPFR_RNDN);
        return mpfr_get_d(r10656, MPFR_RNDN);
}

static mpfr_t r10657, r10658, r10659, r10660, r10661, r10662, r10663, r10664, r10665, r10666, r10667, r10668, r10669, r10670, r10671, r10672, r10673, r10674, r10675, r10676, r10677, r10678, r10679, r10680, r10681, r10682, r10683, r10684;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r10657, "1", 10, MPFR_RNDN);
        mpfr_init(r10658);
        mpfr_init_set_str(r10659, "2", 10, MPFR_RNDN);
        mpfr_init(r10660);
        mpfr_init(r10661);
        mpfr_init(r10662);
        mpfr_init(r10663);
        mpfr_init(r10664);
        mpfr_init(r10665);
        mpfr_init(r10666);
        mpfr_init(r10667);
        mpfr_init(r10668);
        mpfr_init(r10669);
        mpfr_init(r10670);
        mpfr_init(r10671);
        mpfr_init(r10672);
        mpfr_init(r10673);
        mpfr_init(r10674);
        mpfr_init(r10675);
        mpfr_init(r10676);
        mpfr_init(r10677);
        mpfr_init(r10678);
        mpfr_init(r10679);
        mpfr_init(r10680);
        mpfr_init(r10681);
        mpfr_init(r10682);
        mpfr_init(r10683);
        mpfr_init(r10684);
}

double f_dm(double g, double h, double a) {
        ;
        mpfr_set_d(r10658, a, MPFR_RNDN);
        ;
        mpfr_mul(r10660, r10658, r10659, MPFR_RNDN);
        mpfr_div(r10661, r10657, r10660, MPFR_RNDN);
        mpfr_set_d(r10662, g, MPFR_RNDN);
        mpfr_neg(r10663, r10662, MPFR_RNDN);
        mpfr_mul(r10664, r10662, r10662, MPFR_RNDN);
        mpfr_set_d(r10665, h, MPFR_RNDN);
        mpfr_mul(r10666, r10665, r10665, MPFR_RNDN);
        mpfr_sub(r10667, r10664, r10666, MPFR_RNDN);
        mpfr_sqrt(r10668, r10667, MPFR_RNDN);
        mpfr_sub(r10669, r10663, r10668, MPFR_RNDN);
        mpfr_mul(r10670, r10661, r10669, MPFR_RNDN);
        mpfr_cbrt(r10671, r10670, MPFR_RNDN);
        mpfr_cbrt(r10672, r10667, MPFR_RNDN);
        mpfr_abs(r10673, r10672, MPFR_RNDN);
        mpfr_sqrt(r10674, r10672, MPFR_RNDN);
        mpfr_mul(r10675, r10673, r10674, MPFR_RNDN);
        mpfr_add(r10676, r10663, r10675, MPFR_RNDN);
        mpfr_cbrt(r10677, r10676, MPFR_RNDN);
        mpfr_add(r10678, r10663, r10668, MPFR_RNDN);
        mpfr_cbrt(r10679, r10678, MPFR_RNDN);
        mpfr_mul(r10680, r10679, r10677, MPFR_RNDN);
        mpfr_mul(r10681, r10677, r10680, MPFR_RNDN);
        mpfr_mul(r10682, r10681, r10661, MPFR_RNDN);
        mpfr_cbrt(r10683, r10682, MPFR_RNDN);
        mpfr_add(r10684, r10671, r10683, MPFR_RNDN);
        return mpfr_get_d(r10684, MPFR_RNDN);
}