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

char *name = "Maksimov and Kolovsky, Equation (32)";

double f_if(float K, float m, float n, float M, float l) {
        float r9319 = K;
        float r9320 = m;
        float r9321 = n;
        float r9322 = r9320 + r9321;
        float r9323 = r9319 * r9322;
        float r9324 = 2.0f;
        float r9325 = r9323 / r9324;
        float r9326 = M;
        float r9327 = r9325 - r9326;
        float r9328 = cos(r9327);
        float r9329 = r9322 / r9324;
        float r9330 = r9329 - r9326;
        float r9331 = pow(r9330, r9324);
        float r9332 = -r9331;
        float r9333 = l;
        float r9334 = r9320 - r9321;
        float r9335 = fabs(r9334);
        float r9336 = r9333 - r9335;
        float r9337 = r9332 - r9336;
        float r9338 = exp(r9337);
        float r9339 = r9328 * r9338;
        return r9339;
}

double f_id(double K, double m, double n, double M, double l) {
        double r9340 = K;
        double r9341 = m;
        double r9342 = n;
        double r9343 = r9341 + r9342;
        double r9344 = r9340 * r9343;
        double r9345 = 2.0;
        double r9346 = r9344 / r9345;
        double r9347 = M;
        double r9348 = r9346 - r9347;
        double r9349 = cos(r9348);
        double r9350 = r9343 / r9345;
        double r9351 = r9350 - r9347;
        double r9352 = pow(r9351, r9345);
        double r9353 = -r9352;
        double r9354 = l;
        double r9355 = r9341 - r9342;
        double r9356 = fabs(r9355);
        double r9357 = r9354 - r9356;
        double r9358 = r9353 - r9357;
        double r9359 = exp(r9358);
        double r9360 = r9349 * r9359;
        return r9360;
}


double f_of(float __attribute__((unused)) K, float m, float n, float M, float l) {
        float r9361 = 1.0f;
        float r9362 = m;
        float r9363 = n;
        float r9364 = r9362 + r9363;
        float r9365 = 2.0f;
        float r9366 = r9364 / r9365;
        float r9367 = M;
        float r9368 = r9366 - r9367;
        float r9369 = l;
        float r9370 = r9362 - r9363;
        float r9371 = fabs(r9370);
        float r9372 = r9369 - r9371;
        float r9373 = fma(r9368, r9368, r9372);
        float r9374 = exp(r9373);
        float r9375 = r9361 / r9374;
        return r9375;
}

double f_od(double __attribute__((unused)) K, double m, double n, double M, double l) {
        double r9376 = 1.0;
        double r9377 = m;
        double r9378 = n;
        double r9379 = r9377 + r9378;
        double r9380 = 2.0;
        double r9381 = r9379 / r9380;
        double r9382 = M;
        double r9383 = r9381 - r9382;
        double r9384 = l;
        double r9385 = r9377 - r9378;
        double r9386 = fabs(r9385);
        double r9387 = r9384 - r9386;
        double r9388 = fma(r9383, r9383, r9387);
        double r9389 = exp(r9388);
        double r9390 = r9376 / r9389;
        return r9390;
}

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 r9391, r9392, r9393, r9394, r9395, r9396, r9397, r9398, r9399, r9400, r9401, r9402, r9403, r9404, r9405, r9406, r9407, r9408, r9409, r9410, r9411;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1616);
        mpfr_init(r9391);
        mpfr_init(r9392);
        mpfr_init(r9393);
        mpfr_init(r9394);
        mpfr_init(r9395);
        mpfr_init_set_str(r9396, "2", 10, MPFR_RNDN);
        mpfr_init(r9397);
        mpfr_init(r9398);
        mpfr_init(r9399);
        mpfr_init(r9400);
        mpfr_init(r9401);
        mpfr_init(r9402);
        mpfr_init(r9403);
        mpfr_init(r9404);
        mpfr_init(r9405);
        mpfr_init(r9406);
        mpfr_init(r9407);
        mpfr_init(r9408);
        mpfr_init(r9409);
        mpfr_init(r9410);
        mpfr_init(r9411);
}

double f_im(double K, double m, double n, double M, double l) {
        mpfr_set_d(r9391, K, MPFR_RNDN);
        mpfr_set_d(r9392, m, MPFR_RNDN);
        mpfr_set_d(r9393, n, MPFR_RNDN);
        mpfr_add(r9394, r9392, r9393, MPFR_RNDN);
        mpfr_mul(r9395, r9391, r9394, MPFR_RNDN);
        ;
        mpfr_div(r9397, r9395, r9396, MPFR_RNDN);
        mpfr_set_d(r9398, M, MPFR_RNDN);
        mpfr_sub(r9399, r9397, r9398, MPFR_RNDN);
        mpfr_cos(r9400, r9399, MPFR_RNDN);
        mpfr_div(r9401, r9394, r9396, MPFR_RNDN);
        mpfr_sub(r9402, r9401, r9398, MPFR_RNDN);
        mpfr_pow(r9403, r9402, r9396, MPFR_RNDN);
        mpfr_neg(r9404, r9403, MPFR_RNDN);
        mpfr_set_d(r9405, l, MPFR_RNDN);
        mpfr_sub(r9406, r9392, r9393, MPFR_RNDN);
        mpfr_abs(r9407, r9406, MPFR_RNDN);
        mpfr_sub(r9408, r9405, r9407, MPFR_RNDN);
        mpfr_sub(r9409, r9404, r9408, MPFR_RNDN);
        mpfr_exp(r9410, r9409, MPFR_RNDN);
        mpfr_mul(r9411, r9400, r9410, MPFR_RNDN);
        return mpfr_get_d(r9411, MPFR_RNDN);
}

static mpfr_t r9412, r9413, r9414, r9415, r9416, r9417, r9418, r9419, r9420, r9421, r9422, r9423, r9424, r9425, r9426;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1616);
        mpfr_init_set_str(r9412, "1", 10, MPFR_RNDN);
        mpfr_init(r9413);
        mpfr_init(r9414);
        mpfr_init(r9415);
        mpfr_init_set_str(r9416, "2", 10, MPFR_RNDN);
        mpfr_init(r9417);
        mpfr_init(r9418);
        mpfr_init(r9419);
        mpfr_init(r9420);
        mpfr_init(r9421);
        mpfr_init(r9422);
        mpfr_init(r9423);
        mpfr_init(r9424);
        mpfr_init(r9425);
        mpfr_init(r9426);
}

double f_fm(double __attribute__((unused)) K, double m, double n, double M, double l) {
        ;
        mpfr_set_d(r9413, m, MPFR_RNDN);
        mpfr_set_d(r9414, n, MPFR_RNDN);
        mpfr_add(r9415, r9413, r9414, MPFR_RNDN);
        ;
        mpfr_div(r9417, r9415, r9416, MPFR_RNDN);
        mpfr_set_d(r9418, M, MPFR_RNDN);
        mpfr_sub(r9419, r9417, r9418, MPFR_RNDN);
        mpfr_set_d(r9420, l, MPFR_RNDN);
        mpfr_sub(r9421, r9413, r9414, MPFR_RNDN);
        mpfr_abs(r9422, r9421, MPFR_RNDN);
        mpfr_sub(r9423, r9420, r9422, MPFR_RNDN);
        mpfr_fma(r9424, r9419, r9419, r9423, MPFR_RNDN);
        mpfr_exp(r9425, r9424, MPFR_RNDN);
        mpfr_div(r9426, r9412, r9425, MPFR_RNDN);
        return mpfr_get_d(r9426, MPFR_RNDN);
}

static mpfr_t r9427, r9428, r9429, r9430, r9431, r9432, r9433, r9434, r9435, r9436, r9437, r9438, r9439, r9440, r9441;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1616);
        mpfr_init_set_str(r9427, "1", 10, MPFR_RNDN);
        mpfr_init(r9428);
        mpfr_init(r9429);
        mpfr_init(r9430);
        mpfr_init_set_str(r9431, "2", 10, MPFR_RNDN);
        mpfr_init(r9432);
        mpfr_init(r9433);
        mpfr_init(r9434);
        mpfr_init(r9435);
        mpfr_init(r9436);
        mpfr_init(r9437);
        mpfr_init(r9438);
        mpfr_init(r9439);
        mpfr_init(r9440);
        mpfr_init(r9441);
}

double f_dm(double __attribute__((unused)) K, double m, double n, double M, double l) {
        ;
        mpfr_set_d(r9428, m, MPFR_RNDN);
        mpfr_set_d(r9429, n, MPFR_RNDN);
        mpfr_add(r9430, r9428, r9429, MPFR_RNDN);
        ;
        mpfr_div(r9432, r9430, r9431, MPFR_RNDN);
        mpfr_set_d(r9433, M, MPFR_RNDN);
        mpfr_sub(r9434, r9432, r9433, MPFR_RNDN);
        mpfr_set_d(r9435, l, MPFR_RNDN);
        mpfr_sub(r9436, r9428, r9429, MPFR_RNDN);
        mpfr_abs(r9437, r9436, MPFR_RNDN);
        mpfr_sub(r9438, r9435, r9437, MPFR_RNDN);
        mpfr_fma(r9439, r9434, r9434, r9438, MPFR_RNDN);
        mpfr_exp(r9440, r9439, MPFR_RNDN);
        mpfr_div(r9441, r9427, r9440, MPFR_RNDN);
        return mpfr_get_d(r9441, MPFR_RNDN);
}