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

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

double f_if(float k, float n) {
        float r9102 = 1.0f;
        float r9103 = k;
        float r9104 = sqrt(r9103);
        float r9105 = r9102 / r9104;
        float r9106 = 2.0f;
        float r9107 = atan2(1.0, 0.0);
        float r9108 = r9106 * r9107;
        float r9109 = n;
        float r9110 = r9108 * r9109;
        float r9111 = r9102 - r9103;
        float r9112 = r9111 / r9106;
        float r9113 = pow(r9110, r9112);
        float r9114 = r9105 * r9113;
        return r9114;
}

double f_id(double k, double n) {
        double r9115 = 1.0;
        double r9116 = k;
        double r9117 = sqrt(r9116);
        double r9118 = r9115 / r9117;
        double r9119 = 2.0;
        double r9120 = atan2(1.0, 0.0);
        double r9121 = r9119 * r9120;
        double r9122 = n;
        double r9123 = r9121 * r9122;
        double r9124 = r9115 - r9116;
        double r9125 = r9124 / r9119;
        double r9126 = pow(r9123, r9125);
        double r9127 = r9118 * r9126;
        return r9127;
}


double f_of(float k, float n) {
        float r9128 = n;
        float r9129 = 2.0f;
        float r9130 = r9128 * r9129;
        float r9131 = atan2(1.0, 0.0);
        float r9132 = r9130 * r9131;
        float r9133 = 1.0f;
        float r9134 = r9133 / r9129;
        float r9135 = pow(r9132, r9134);
        float r9136 = k;
        float r9137 = sqrt(r9136);
        float r9138 = r9136 / r9129;
        float r9139 = pow(r9130, r9138);
        float r9140 = pow(r9131, r9138);
        float r9141 = r9139 * r9140;
        float r9142 = r9137 * r9141;
        float r9143 = r9135 / r9142;
        return r9143;
}

double f_od(double k, double n) {
        double r9144 = n;
        double r9145 = 2.0;
        double r9146 = r9144 * r9145;
        double r9147 = atan2(1.0, 0.0);
        double r9148 = r9146 * r9147;
        double r9149 = 1.0;
        double r9150 = r9149 / r9145;
        double r9151 = pow(r9148, r9150);
        double r9152 = k;
        double r9153 = sqrt(r9152);
        double r9154 = r9152 / r9145;
        double r9155 = pow(r9146, r9154);
        double r9156 = pow(r9147, r9154);
        double r9157 = r9155 * r9156;
        double r9158 = r9153 * r9157;
        double r9159 = r9151 / r9158;
        return r9159;
}

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 r9160, r9161, r9162, r9163, r9164, r9165, r9166, r9167, r9168, r9169, r9170, r9171, r9172;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r9160, "1", 10, MPFR_RNDN);
        mpfr_init(r9161);
        mpfr_init(r9162);
        mpfr_init(r9163);
        mpfr_init_set_str(r9164, "2", 10, MPFR_RNDN);
        mpfr_init(r9165);
        mpfr_init(r9166);
        mpfr_init(r9167);
        mpfr_init(r9168);
        mpfr_init(r9169);
        mpfr_init(r9170);
        mpfr_init(r9171);
        mpfr_init(r9172);
}

double f_im(double k, double n) {
        ;
        mpfr_set_d(r9161, k, MPFR_RNDN);
        mpfr_sqrt(r9162, r9161, MPFR_RNDN);
        mpfr_div(r9163, r9160, r9162, MPFR_RNDN);
        ;
        mpfr_const_pi(r9165, MPFR_RNDN);
        mpfr_mul(r9166, r9164, r9165, MPFR_RNDN);
        mpfr_set_d(r9167, n, MPFR_RNDN);
        mpfr_mul(r9168, r9166, r9167, MPFR_RNDN);
        mpfr_sub(r9169, r9160, r9161, MPFR_RNDN);
        mpfr_div(r9170, r9169, r9164, MPFR_RNDN);
        mpfr_pow(r9171, r9168, r9170, MPFR_RNDN);
        mpfr_mul(r9172, r9163, r9171, MPFR_RNDN);
        return mpfr_get_d(r9172, MPFR_RNDN);
}

static mpfr_t r9173, r9174, r9175, r9176, r9177, r9178, r9179, r9180, r9181, r9182, r9183, r9184, r9185, r9186, r9187, r9188;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9173);
        mpfr_init_set_str(r9174, "2", 10, MPFR_RNDN);
        mpfr_init(r9175);
        mpfr_init(r9176);
        mpfr_init(r9177);
        mpfr_init_set_str(r9178, "1", 10, MPFR_RNDN);
        mpfr_init(r9179);
        mpfr_init(r9180);
        mpfr_init(r9181);
        mpfr_init(r9182);
        mpfr_init(r9183);
        mpfr_init(r9184);
        mpfr_init(r9185);
        mpfr_init(r9186);
        mpfr_init(r9187);
        mpfr_init(r9188);
}

double f_fm(double k, double n) {
        mpfr_set_d(r9173, n, MPFR_RNDN);
        ;
        mpfr_mul(r9175, r9173, r9174, MPFR_RNDN);
        mpfr_const_pi(r9176, MPFR_RNDN);
        mpfr_mul(r9177, r9175, r9176, MPFR_RNDN);
        ;
        mpfr_div(r9179, r9178, r9174, MPFR_RNDN);
        mpfr_pow(r9180, r9177, r9179, MPFR_RNDN);
        mpfr_set_d(r9181, k, MPFR_RNDN);
        mpfr_sqrt(r9182, r9181, MPFR_RNDN);
        mpfr_div(r9183, r9181, r9174, MPFR_RNDN);
        mpfr_pow(r9184, r9175, r9183, MPFR_RNDN);
        mpfr_pow(r9185, r9176, r9183, MPFR_RNDN);
        mpfr_mul(r9186, r9184, r9185, MPFR_RNDN);
        mpfr_mul(r9187, r9182, r9186, MPFR_RNDN);
        mpfr_div(r9188, r9180, r9187, MPFR_RNDN);
        return mpfr_get_d(r9188, MPFR_RNDN);
}

static mpfr_t r9189, r9190, r9191, r9192, r9193, r9194, r9195, r9196, r9197, r9198, r9199, r9200, r9201, r9202, r9203, r9204;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9189);
        mpfr_init_set_str(r9190, "2", 10, MPFR_RNDN);
        mpfr_init(r9191);
        mpfr_init(r9192);
        mpfr_init(r9193);
        mpfr_init_set_str(r9194, "1", 10, MPFR_RNDN);
        mpfr_init(r9195);
        mpfr_init(r9196);
        mpfr_init(r9197);
        mpfr_init(r9198);
        mpfr_init(r9199);
        mpfr_init(r9200);
        mpfr_init(r9201);
        mpfr_init(r9202);
        mpfr_init(r9203);
        mpfr_init(r9204);
}

double f_dm(double k, double n) {
        mpfr_set_d(r9189, n, MPFR_RNDN);
        ;
        mpfr_mul(r9191, r9189, r9190, MPFR_RNDN);
        mpfr_const_pi(r9192, MPFR_RNDN);
        mpfr_mul(r9193, r9191, r9192, MPFR_RNDN);
        ;
        mpfr_div(r9195, r9194, r9190, MPFR_RNDN);
        mpfr_pow(r9196, r9193, r9195, MPFR_RNDN);
        mpfr_set_d(r9197, k, MPFR_RNDN);
        mpfr_sqrt(r9198, r9197, MPFR_RNDN);
        mpfr_div(r9199, r9197, r9190, MPFR_RNDN);
        mpfr_pow(r9200, r9191, r9199, MPFR_RNDN);
        mpfr_pow(r9201, r9192, r9199, MPFR_RNDN);
        mpfr_mul(r9202, r9200, r9201, MPFR_RNDN);
        mpfr_mul(r9203, r9198, r9202, MPFR_RNDN);
        mpfr_div(r9204, r9196, r9203, MPFR_RNDN);
        return mpfr_get_d(r9204, MPFR_RNDN);
}