#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 r26424 = k;
        float r26425 = sqrt(r26424);
        float r26426 = 1.0/r26425;
        float r26427 = 2.0f;
        float r26428 = atan2(1.0, 0.0);
        float r26429 = r26427 * r26428;
        float r26430 = n;
        float r26431 = r26429 * r26430;
        float r26432 = 1.0f;
        float r26433 = r26432 - r26424;
        float r26434 = r26433 / r26427;
        float r26435 = pow(r26431, r26434);
        float r26436 = r26426 * r26435;
        return r26436;
}

double f_id(double k, double n) {
        double r26437 = k;
        double r26438 = sqrt(r26437);
        double r26439 = 1.0/r26438;
        double r26440 = 2.0;
        double r26441 = atan2(1.0, 0.0);
        double r26442 = r26440 * r26441;
        double r26443 = n;
        double r26444 = r26442 * r26443;
        double r26445 = 1.0;
        double r26446 = r26445 - r26437;
        double r26447 = r26446 / r26440;
        double r26448 = pow(r26444, r26447);
        double r26449 = r26439 * r26448;
        return r26449;
}


double f_of(float k, float n) {
        float r26450 = k;
        float r26451 = sqrt(r26450);
        float r26452 = 1.0/r26451;
        float r26453 = 2.0f;
        float r26454 = atan2(1.0, 0.0);
        float r26455 = r26453 * r26454;
        float r26456 = 1.0f;
        float r26457 = r26456 - r26450;
        float r26458 = r26457 / r26453;
        float r26459 = pow(r26455, r26458);
        float r26460 = sqrt(r26459);
        float r26461 = r26460 * r26460;
        float r26462 = r26452 * r26461;
        float r26463 = n;
        float r26464 = r26456 / r26453;
        float r26465 = pow(r26463, r26464);
        float r26466 = r26450 / r26453;
        float r26467 = pow(r26463, r26466);
        float r26468 = r26465 / r26467;
        float r26469 = r26462 * r26468;
        return r26469;
}

double f_od(double k, double n) {
        double r26470 = k;
        double r26471 = sqrt(r26470);
        double r26472 = 1.0/r26471;
        double r26473 = 2.0;
        double r26474 = atan2(1.0, 0.0);
        double r26475 = r26473 * r26474;
        double r26476 = 1.0;
        double r26477 = r26476 - r26470;
        double r26478 = r26477 / r26473;
        double r26479 = pow(r26475, r26478);
        double r26480 = sqrt(r26479);
        double r26481 = r26480 * r26480;
        double r26482 = r26472 * r26481;
        double r26483 = n;
        double r26484 = r26476 / r26473;
        double r26485 = pow(r26483, r26484);
        double r26486 = r26470 / r26473;
        double r26487 = pow(r26483, r26486);
        double r26488 = r26485 / r26487;
        double r26489 = r26482 * r26488;
        return r26489;
}

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 r26490, r26491, r26492, r26493, r26494, r26495, r26496, r26497, r26498, r26499, r26500, r26501, r26502;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r26490);
        mpfr_init(r26491);
        mpfr_init(r26492);
        mpfr_init_set_str(r26493, "2", 10, MPFR_RNDN);
        mpfr_init(r26494);
        mpfr_init(r26495);
        mpfr_init(r26496);
        mpfr_init(r26497);
        mpfr_init_set_str(r26498, "1", 10, MPFR_RNDN);
        mpfr_init(r26499);
        mpfr_init(r26500);
        mpfr_init(r26501);
        mpfr_init(r26502);
}

double f_im(double k, double n) {
        mpfr_set_d(r26490, k, MPFR_RNDN);
        mpfr_sqrt(r26491, r26490, MPFR_RNDN);
        mpfr_ui_div(r26492, 1, r26491, MPFR_RNDN);
        ;
        mpfr_const_pi(r26494, MPFR_RNDN);
        mpfr_mul(r26495, r26493, r26494, MPFR_RNDN);
        mpfr_set_d(r26496, n, MPFR_RNDN);
        mpfr_mul(r26497, r26495, r26496, MPFR_RNDN);
        ;
        mpfr_sub(r26499, r26498, r26490, MPFR_RNDN);
        mpfr_div(r26500, r26499, r26493, MPFR_RNDN);
        mpfr_pow(r26501, r26497, r26500, MPFR_RNDN);
        mpfr_mul(r26502, r26492, r26501, MPFR_RNDN);
        return mpfr_get_d(r26502, MPFR_RNDN);
}

static mpfr_t r26503, r26504, r26505, r26506, r26507, r26508, r26509, r26510, r26511, r26512, r26513, r26514, r26515, r26516, r26517, r26518, r26519, r26520, r26521, r26522;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r26503);
        mpfr_init(r26504);
        mpfr_init(r26505);
        mpfr_init_set_str(r26506, "2", 10, MPFR_RNDN);
        mpfr_init(r26507);
        mpfr_init(r26508);
        mpfr_init_set_str(r26509, "1", 10, MPFR_RNDN);
        mpfr_init(r26510);
        mpfr_init(r26511);
        mpfr_init(r26512);
        mpfr_init(r26513);
        mpfr_init(r26514);
        mpfr_init(r26515);
        mpfr_init(r26516);
        mpfr_init(r26517);
        mpfr_init(r26518);
        mpfr_init(r26519);
        mpfr_init(r26520);
        mpfr_init(r26521);
        mpfr_init(r26522);
}

double f_fm(double k, double n) {
        mpfr_set_d(r26503, k, MPFR_RNDN);
        mpfr_sqrt(r26504, r26503, MPFR_RNDN);
        mpfr_ui_div(r26505, 1, r26504, MPFR_RNDN);
        ;
        mpfr_const_pi(r26507, MPFR_RNDN);
        mpfr_mul(r26508, r26506, r26507, MPFR_RNDN);
        ;
        mpfr_sub(r26510, r26509, r26503, MPFR_RNDN);
        mpfr_div(r26511, r26510, r26506, MPFR_RNDN);
        mpfr_pow(r26512, r26508, r26511, MPFR_RNDN);
        mpfr_sqrt(r26513, r26512, MPFR_RNDN);
        mpfr_sqr(r26514, r26513, MPFR_RNDN);
        mpfr_mul(r26515, r26505, r26514, MPFR_RNDN);
        mpfr_set_d(r26516, n, MPFR_RNDN);
        mpfr_div(r26517, r26509, r26506, MPFR_RNDN);
        mpfr_pow(r26518, r26516, r26517, MPFR_RNDN);
        mpfr_div(r26519, r26503, r26506, MPFR_RNDN);
        mpfr_pow(r26520, r26516, r26519, MPFR_RNDN);
        mpfr_div(r26521, r26518, r26520, MPFR_RNDN);
        mpfr_mul(r26522, r26515, r26521, MPFR_RNDN);
        return mpfr_get_d(r26522, MPFR_RNDN);
}

static mpfr_t r26523, r26524, r26525, r26526, r26527, r26528, r26529, r26530, r26531, r26532, r26533, r26534, r26535, r26536, r26537, r26538, r26539, r26540, r26541, r26542;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r26523);
        mpfr_init(r26524);
        mpfr_init(r26525);
        mpfr_init_set_str(r26526, "2", 10, MPFR_RNDN);
        mpfr_init(r26527);
        mpfr_init(r26528);
        mpfr_init_set_str(r26529, "1", 10, MPFR_RNDN);
        mpfr_init(r26530);
        mpfr_init(r26531);
        mpfr_init(r26532);
        mpfr_init(r26533);
        mpfr_init(r26534);
        mpfr_init(r26535);
        mpfr_init(r26536);
        mpfr_init(r26537);
        mpfr_init(r26538);
        mpfr_init(r26539);
        mpfr_init(r26540);
        mpfr_init(r26541);
        mpfr_init(r26542);
}

double f_dm(double k, double n) {
        mpfr_set_d(r26523, k, MPFR_RNDN);
        mpfr_sqrt(r26524, r26523, MPFR_RNDN);
        mpfr_ui_div(r26525, 1, r26524, MPFR_RNDN);
        ;
        mpfr_const_pi(r26527, MPFR_RNDN);
        mpfr_mul(r26528, r26526, r26527, MPFR_RNDN);
        ;
        mpfr_sub(r26530, r26529, r26523, MPFR_RNDN);
        mpfr_div(r26531, r26530, r26526, MPFR_RNDN);
        mpfr_pow(r26532, r26528, r26531, MPFR_RNDN);
        mpfr_sqrt(r26533, r26532, MPFR_RNDN);
        mpfr_sqr(r26534, r26533, MPFR_RNDN);
        mpfr_mul(r26535, r26525, r26534, MPFR_RNDN);
        mpfr_set_d(r26536, n, MPFR_RNDN);
        mpfr_div(r26537, r26529, r26526, MPFR_RNDN);
        mpfr_pow(r26538, r26536, r26537, MPFR_RNDN);
        mpfr_div(r26539, r26523, r26526, MPFR_RNDN);
        mpfr_pow(r26540, r26536, r26539, MPFR_RNDN);
        mpfr_div(r26541, r26538, r26540, MPFR_RNDN);
        mpfr_mul(r26542, r26535, r26541, MPFR_RNDN);
        return mpfr_get_d(r26542, MPFR_RNDN);
}