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

char *name = "Migdal et al, Equation (64)";

double f_if(float a1, float a2, float th) {
        float r9578 = th;
        float r9579 = cos(r9578);
        float r9580 = 2;
        float r9581 = sqrt(r9580);
        float r9582 = r9579 / r9581;
        float r9583 = a1;
        float r9584 = r9583 * r9583;
        float r9585 = r9582 * r9584;
        float r9586 = a2;
        float r9587 = r9586 * r9586;
        float r9588 = r9582 * r9587;
        float r9589 = r9585 + r9588;
        return r9589;
}

double f_id(double a1, double a2, double th) {
        double r9590 = th;
        double r9591 = cos(r9590);
        double r9592 = 2;
        double r9593 = sqrt(r9592);
        double r9594 = r9591 / r9593;
        double r9595 = a1;
        double r9596 = r9595 * r9595;
        double r9597 = r9594 * r9596;
        double r9598 = a2;
        double r9599 = r9598 * r9598;
        double r9600 = r9594 * r9599;
        double r9601 = r9597 + r9600;
        return r9601;
}


double f_of(float a1, float a2, float th) {
        float r9602 = th;
        float r9603 = cos(r9602);
        float r9604 = 2;
        float r9605 = sqrt(r9604);
        float r9606 = sqrt(r9605);
        float r9607 = r9603 / r9606;
        float r9608 = r9607 / r9606;
        float r9609 = a1;
        float r9610 = r9609 * r9609;
        float r9611 = r9608 * r9610;
        float r9612 = r9603 / r9605;
        float r9613 = a2;
        float r9614 = r9613 * r9613;
        float r9615 = r9612 * r9614;
        float r9616 = r9611 + r9615;
        return r9616;
}

double f_od(double a1, double a2, double th) {
        double r9617 = th;
        double r9618 = cos(r9617);
        double r9619 = 2;
        double r9620 = sqrt(r9619);
        double r9621 = sqrt(r9620);
        double r9622 = r9618 / r9621;
        double r9623 = r9622 / r9621;
        double r9624 = a1;
        double r9625 = r9624 * r9624;
        double r9626 = r9623 * r9625;
        double r9627 = r9618 / r9620;
        double r9628 = a2;
        double r9629 = r9628 * r9628;
        double r9630 = r9627 * r9629;
        double r9631 = r9626 + r9630;
        return r9631;
}

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 r9632, r9633, r9634, r9635, r9636, r9637, r9638, r9639, r9640, r9641, r9642, r9643;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9632);
        mpfr_init(r9633);
        mpfr_init_set_str(r9634, "2", 10, MPFR_RNDN);
        mpfr_init(r9635);
        mpfr_init(r9636);
        mpfr_init(r9637);
        mpfr_init(r9638);
        mpfr_init(r9639);
        mpfr_init(r9640);
        mpfr_init(r9641);
        mpfr_init(r9642);
        mpfr_init(r9643);
}

double f_im(double a1, double a2, double th) {
        mpfr_set_d(r9632, th, MPFR_RNDN);
        mpfr_cos(r9633, r9632, MPFR_RNDN);
        ;
        mpfr_sqrt(r9635, r9634, MPFR_RNDN);
        mpfr_div(r9636, r9633, r9635, MPFR_RNDN);
        mpfr_set_d(r9637, a1, MPFR_RNDN);
        mpfr_mul(r9638, r9637, r9637, MPFR_RNDN);
        mpfr_mul(r9639, r9636, r9638, MPFR_RNDN);
        mpfr_set_d(r9640, a2, MPFR_RNDN);
        mpfr_mul(r9641, r9640, r9640, MPFR_RNDN);
        mpfr_mul(r9642, r9636, r9641, MPFR_RNDN);
        mpfr_add(r9643, r9639, r9642, MPFR_RNDN);
        return mpfr_get_d(r9643, MPFR_RNDN);
}

static mpfr_t r9644, r9645, r9646, r9647, r9648, r9649, r9650, r9651, r9652, r9653, r9654, r9655, r9656, r9657, r9658;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9644);
        mpfr_init(r9645);
        mpfr_init_set_str(r9646, "2", 10, MPFR_RNDN);
        mpfr_init(r9647);
        mpfr_init(r9648);
        mpfr_init(r9649);
        mpfr_init(r9650);
        mpfr_init(r9651);
        mpfr_init(r9652);
        mpfr_init(r9653);
        mpfr_init(r9654);
        mpfr_init(r9655);
        mpfr_init(r9656);
        mpfr_init(r9657);
        mpfr_init(r9658);
}

double f_fm(double a1, double a2, double th) {
        mpfr_set_d(r9644, th, MPFR_RNDN);
        mpfr_cos(r9645, r9644, MPFR_RNDN);
        ;
        mpfr_sqrt(r9647, r9646, MPFR_RNDN);
        mpfr_sqrt(r9648, r9647, MPFR_RNDN);
        mpfr_div(r9649, r9645, r9648, MPFR_RNDN);
        mpfr_div(r9650, r9649, r9648, MPFR_RNDN);
        mpfr_set_d(r9651, a1, MPFR_RNDN);
        mpfr_mul(r9652, r9651, r9651, MPFR_RNDN);
        mpfr_mul(r9653, r9650, r9652, MPFR_RNDN);
        mpfr_div(r9654, r9645, r9647, MPFR_RNDN);
        mpfr_set_d(r9655, a2, MPFR_RNDN);
        mpfr_mul(r9656, r9655, r9655, MPFR_RNDN);
        mpfr_mul(r9657, r9654, r9656, MPFR_RNDN);
        mpfr_add(r9658, r9653, r9657, MPFR_RNDN);
        return mpfr_get_d(r9658, MPFR_RNDN);
}

static mpfr_t r9659, r9660, r9661, r9662, r9663, r9664, r9665, r9666, r9667, r9668, r9669, r9670, r9671, r9672, r9673;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9659);
        mpfr_init(r9660);
        mpfr_init_set_str(r9661, "2", 10, MPFR_RNDN);
        mpfr_init(r9662);
        mpfr_init(r9663);
        mpfr_init(r9664);
        mpfr_init(r9665);
        mpfr_init(r9666);
        mpfr_init(r9667);
        mpfr_init(r9668);
        mpfr_init(r9669);
        mpfr_init(r9670);
        mpfr_init(r9671);
        mpfr_init(r9672);
        mpfr_init(r9673);
}

double f_dm(double a1, double a2, double th) {
        mpfr_set_d(r9659, th, MPFR_RNDN);
        mpfr_cos(r9660, r9659, MPFR_RNDN);
        ;
        mpfr_sqrt(r9662, r9661, MPFR_RNDN);
        mpfr_sqrt(r9663, r9662, MPFR_RNDN);
        mpfr_div(r9664, r9660, r9663, MPFR_RNDN);
        mpfr_div(r9665, r9664, r9663, MPFR_RNDN);
        mpfr_set_d(r9666, a1, MPFR_RNDN);
        mpfr_mul(r9667, r9666, r9666, MPFR_RNDN);
        mpfr_mul(r9668, r9665, r9667, MPFR_RNDN);
        mpfr_div(r9669, r9660, r9662, MPFR_RNDN);
        mpfr_set_d(r9670, a2, MPFR_RNDN);
        mpfr_mul(r9671, r9670, r9670, MPFR_RNDN);
        mpfr_mul(r9672, r9669, r9671, MPFR_RNDN);
        mpfr_add(r9673, r9668, r9672, MPFR_RNDN);
        return mpfr_get_d(r9673, MPFR_RNDN);
}