#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 r9129 = th;
        float r9130 = cos(r9129);
        float r9131 = 2.0f;
        float r9132 = sqrt(r9131);
        float r9133 = r9130 / r9132;
        float r9134 = a1;
        float r9135 = r9134 * r9134;
        float r9136 = r9133 * r9135;
        float r9137 = a2;
        float r9138 = r9137 * r9137;
        float r9139 = r9133 * r9138;
        float r9140 = r9136 + r9139;
        return r9140;
}

double f_id(double a1, double a2, double th) {
        double r9141 = th;
        double r9142 = cos(r9141);
        double r9143 = 2.0;
        double r9144 = sqrt(r9143);
        double r9145 = r9142 / r9144;
        double r9146 = a1;
        double r9147 = r9146 * r9146;
        double r9148 = r9145 * r9147;
        double r9149 = a2;
        double r9150 = r9149 * r9149;
        double r9151 = r9145 * r9150;
        double r9152 = r9148 + r9151;
        return r9152;
}


double f_of(float a1, float a2, float th) {
        float r9153 = th;
        float r9154 = cos(r9153);
        float r9155 = a1;
        float r9156 = a2;
        float r9157 = hypot(r9155, r9156);
        float r9158 = r9154 * r9157;
        float r9159 = 2.0f;
        float r9160 = sqrt(r9159);
        float r9161 = r9160 / r9157;
        float r9162 = r9158 / r9161;
        return r9162;
}

double f_od(double a1, double a2, double th) {
        double r9163 = th;
        double r9164 = cos(r9163);
        double r9165 = a1;
        double r9166 = a2;
        double r9167 = hypot(r9165, r9166);
        double r9168 = r9164 * r9167;
        double r9169 = 2.0;
        double r9170 = sqrt(r9169);
        double r9171 = r9170 / r9167;
        double r9172 = r9168 / r9171;
        return r9172;
}

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 r9173, r9174, r9175, r9176, r9177, r9178, r9179, r9180, r9181, r9182, r9183, r9184;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9173);
        mpfr_init(r9174);
        mpfr_init_set_str(r9175, "2", 10, MPFR_RNDN);
        mpfr_init(r9176);
        mpfr_init(r9177);
        mpfr_init(r9178);
        mpfr_init(r9179);
        mpfr_init(r9180);
        mpfr_init(r9181);
        mpfr_init(r9182);
        mpfr_init(r9183);
        mpfr_init(r9184);
}

double f_im(double a1, double a2, double th) {
        mpfr_set_d(r9173, th, MPFR_RNDN);
        mpfr_cos(r9174, r9173, MPFR_RNDN);
        ;
        mpfr_sqrt(r9176, r9175, MPFR_RNDN);
        mpfr_div(r9177, r9174, r9176, MPFR_RNDN);
        mpfr_set_d(r9178, a1, MPFR_RNDN);
        mpfr_mul(r9179, r9178, r9178, MPFR_RNDN);
        mpfr_mul(r9180, r9177, r9179, MPFR_RNDN);
        mpfr_set_d(r9181, a2, MPFR_RNDN);
        mpfr_mul(r9182, r9181, r9181, MPFR_RNDN);
        mpfr_mul(r9183, r9177, r9182, MPFR_RNDN);
        mpfr_add(r9184, r9180, r9183, MPFR_RNDN);
        return mpfr_get_d(r9184, MPFR_RNDN);
}

static mpfr_t r9185, r9186, r9187, r9188, r9189, r9190, r9191, r9192, r9193, r9194;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9185);
        mpfr_init(r9186);
        mpfr_init(r9187);
        mpfr_init(r9188);
        mpfr_init(r9189);
        mpfr_init(r9190);
        mpfr_init_set_str(r9191, "2", 10, MPFR_RNDN);
        mpfr_init(r9192);
        mpfr_init(r9193);
        mpfr_init(r9194);
}

double f_fm(double a1, double a2, double th) {
        mpfr_set_d(r9185, th, MPFR_RNDN);
        mpfr_cos(r9186, r9185, MPFR_RNDN);
        mpfr_set_d(r9187, a1, MPFR_RNDN);
        mpfr_set_d(r9188, a2, MPFR_RNDN);
        mpfr_hypot(r9189, r9187, r9188, MPFR_RNDN);
        mpfr_mul(r9190, r9186, r9189, MPFR_RNDN);
        ;
        mpfr_sqrt(r9192, r9191, MPFR_RNDN);
        mpfr_div(r9193, r9192, r9189, MPFR_RNDN);
        mpfr_div(r9194, r9190, r9193, MPFR_RNDN);
        return mpfr_get_d(r9194, MPFR_RNDN);
}

static mpfr_t r9195, r9196, r9197, r9198, r9199, r9200, r9201, r9202, r9203, r9204;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9195);
        mpfr_init(r9196);
        mpfr_init(r9197);
        mpfr_init(r9198);
        mpfr_init(r9199);
        mpfr_init(r9200);
        mpfr_init_set_str(r9201, "2", 10, MPFR_RNDN);
        mpfr_init(r9202);
        mpfr_init(r9203);
        mpfr_init(r9204);
}

double f_dm(double a1, double a2, double th) {
        mpfr_set_d(r9195, th, MPFR_RNDN);
        mpfr_cos(r9196, r9195, MPFR_RNDN);
        mpfr_set_d(r9197, a1, MPFR_RNDN);
        mpfr_set_d(r9198, a2, MPFR_RNDN);
        mpfr_hypot(r9199, r9197, r9198, MPFR_RNDN);
        mpfr_mul(r9200, r9196, r9199, MPFR_RNDN);
        ;
        mpfr_sqrt(r9202, r9201, MPFR_RNDN);
        mpfr_div(r9203, r9202, r9199, MPFR_RNDN);
        mpfr_div(r9204, r9200, r9203, MPFR_RNDN);
        return mpfr_get_d(r9204, MPFR_RNDN);
}