#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 r9334 = 1;
        float r9335 = k;
        float r9336 = sqrt(r9335);
        float r9337 = r9334 / r9336;
        float r9338 = 2;
        float r9339 = atan2(1.0, 0.0);
        float r9340 = r9338 * r9339;
        float r9341 = n;
        float r9342 = r9340 * r9341;
        float r9343 = r9334 - r9335;
        float r9344 = r9343 / r9338;
        float r9345 = pow(r9342, r9344);
        float r9346 = r9337 * r9345;
        return r9346;
}

double f_id(double k, double n) {
        double r9347 = 1;
        double r9348 = k;
        double r9349 = sqrt(r9348);
        double r9350 = r9347 / r9349;
        double r9351 = 2;
        double r9352 = atan2(1.0, 0.0);
        double r9353 = r9351 * r9352;
        double r9354 = n;
        double r9355 = r9353 * r9354;
        double r9356 = r9347 - r9348;
        double r9357 = r9356 / r9351;
        double r9358 = pow(r9355, r9357);
        double r9359 = r9350 * r9358;
        return r9359;
}


double f_of(float k, float n) {
        float r9360 = n;
        float r9361 = 2;
        float r9362 = r9360 * r9361;
        float r9363 = atan2(1.0, 0.0);
        float r9364 = r9362 * r9363;
        float r9365 = 1;
        float r9366 = r9365 / r9361;
        float r9367 = pow(r9364, r9366);
        float r9368 = k;
        float r9369 = r9368 / r9361;
        float r9370 = pow(r9364, r9369);
        float r9371 = r9367 / r9370;
        float r9372 = sqrt(r9368);
        float r9373 = r9371 / r9372;
        return r9373;
}

double f_od(double k, double n) {
        double r9374 = n;
        double r9375 = 2;
        double r9376 = r9374 * r9375;
        double r9377 = atan2(1.0, 0.0);
        double r9378 = r9376 * r9377;
        double r9379 = 1;
        double r9380 = r9379 / r9375;
        double r9381 = pow(r9378, r9380);
        double r9382 = k;
        double r9383 = r9382 / r9375;
        double r9384 = pow(r9378, r9383);
        double r9385 = r9381 / r9384;
        double r9386 = sqrt(r9382);
        double r9387 = r9385 / r9386;
        return r9387;
}

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 r9388, r9389, r9390, r9391, r9392, r9393, r9394, r9395, r9396, r9397, r9398, r9399, r9400;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r9388, "1", 10, MPFR_RNDN);
        mpfr_init(r9389);
        mpfr_init(r9390);
        mpfr_init(r9391);
        mpfr_init_set_str(r9392, "2", 10, MPFR_RNDN);
        mpfr_init(r9393);
        mpfr_init(r9394);
        mpfr_init(r9395);
        mpfr_init(r9396);
        mpfr_init(r9397);
        mpfr_init(r9398);
        mpfr_init(r9399);
        mpfr_init(r9400);
}

double f_im(double k, double n) {
        ;
        mpfr_set_d(r9389, k, MPFR_RNDN);
        mpfr_sqrt(r9390, r9389, MPFR_RNDN);
        mpfr_div(r9391, r9388, r9390, MPFR_RNDN);
        ;
        mpfr_const_pi(r9393, MPFR_RNDN);
        mpfr_mul(r9394, r9392, r9393, MPFR_RNDN);
        mpfr_set_d(r9395, n, MPFR_RNDN);
        mpfr_mul(r9396, r9394, r9395, MPFR_RNDN);
        mpfr_sub(r9397, r9388, r9389, MPFR_RNDN);
        mpfr_div(r9398, r9397, r9392, MPFR_RNDN);
        mpfr_pow(r9399, r9396, r9398, MPFR_RNDN);
        mpfr_mul(r9400, r9391, r9399, MPFR_RNDN);
        return mpfr_get_d(r9400, MPFR_RNDN);
}

static mpfr_t r9401, r9402, r9403, r9404, r9405, r9406, r9407, r9408, r9409, r9410, r9411, r9412, r9413, r9414;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9401);
        mpfr_init_set_str(r9402, "2", 10, MPFR_RNDN);
        mpfr_init(r9403);
        mpfr_init(r9404);
        mpfr_init(r9405);
        mpfr_init_set_str(r9406, "1", 10, MPFR_RNDN);
        mpfr_init(r9407);
        mpfr_init(r9408);
        mpfr_init(r9409);
        mpfr_init(r9410);
        mpfr_init(r9411);
        mpfr_init(r9412);
        mpfr_init(r9413);
        mpfr_init(r9414);
}

double f_fm(double k, double n) {
        mpfr_set_d(r9401, n, MPFR_RNDN);
        ;
        mpfr_mul(r9403, r9401, r9402, MPFR_RNDN);
        mpfr_const_pi(r9404, MPFR_RNDN);
        mpfr_mul(r9405, r9403, r9404, MPFR_RNDN);
        ;
        mpfr_div(r9407, r9406, r9402, MPFR_RNDN);
        mpfr_pow(r9408, r9405, r9407, MPFR_RNDN);
        mpfr_set_d(r9409, k, MPFR_RNDN);
        mpfr_div(r9410, r9409, r9402, MPFR_RNDN);
        mpfr_pow(r9411, r9405, r9410, MPFR_RNDN);
        mpfr_div(r9412, r9408, r9411, MPFR_RNDN);
        mpfr_sqrt(r9413, r9409, MPFR_RNDN);
        mpfr_div(r9414, r9412, r9413, MPFR_RNDN);
        return mpfr_get_d(r9414, MPFR_RNDN);
}

static mpfr_t r9415, r9416, r9417, r9418, r9419, r9420, r9421, r9422, r9423, r9424, r9425, r9426, r9427, r9428;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9415);
        mpfr_init_set_str(r9416, "2", 10, MPFR_RNDN);
        mpfr_init(r9417);
        mpfr_init(r9418);
        mpfr_init(r9419);
        mpfr_init_set_str(r9420, "1", 10, MPFR_RNDN);
        mpfr_init(r9421);
        mpfr_init(r9422);
        mpfr_init(r9423);
        mpfr_init(r9424);
        mpfr_init(r9425);
        mpfr_init(r9426);
        mpfr_init(r9427);
        mpfr_init(r9428);
}

double f_dm(double k, double n) {
        mpfr_set_d(r9415, n, MPFR_RNDN);
        ;
        mpfr_mul(r9417, r9415, r9416, MPFR_RNDN);
        mpfr_const_pi(r9418, MPFR_RNDN);
        mpfr_mul(r9419, r9417, r9418, MPFR_RNDN);
        ;
        mpfr_div(r9421, r9420, r9416, MPFR_RNDN);
        mpfr_pow(r9422, r9419, r9421, MPFR_RNDN);
        mpfr_set_d(r9423, k, MPFR_RNDN);
        mpfr_div(r9424, r9423, r9416, MPFR_RNDN);
        mpfr_pow(r9425, r9419, r9424, MPFR_RNDN);
        mpfr_div(r9426, r9422, r9425, MPFR_RNDN);
        mpfr_sqrt(r9427, r9423, MPFR_RNDN);
        mpfr_div(r9428, r9426, r9427, MPFR_RNDN);
        return mpfr_get_d(r9428, MPFR_RNDN);
}