#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 r26017 = K;
        float r26018 = m;
        float r26019 = n;
        float r26020 = r26018 + r26019;
        float r26021 = r26017 * r26020;
        float r26022 = 2;
        float r26023 = r26021 / r26022;
        float r26024 = M;
        float r26025 = r26023 - r26024;
        float r26026 = cos(r26025);
        float r26027 = r26020 / r26022;
        float r26028 = r26027 - r26024;
        float r26029 = pow(r26028, r26022);
        float r26030 = -r26029;
        float r26031 = l;
        float r26032 = r26018 - r26019;
        float r26033 = fabs(r26032);
        float r26034 = r26031 - r26033;
        float r26035 = r26030 - r26034;
        float r26036 = exp(r26035);
        float r26037 = r26026 * r26036;
        return r26037;
}

double f_id(double K, double m, double n, double M, double l) {
        double r26038 = K;
        double r26039 = m;
        double r26040 = n;
        double r26041 = r26039 + r26040;
        double r26042 = r26038 * r26041;
        double r26043 = 2;
        double r26044 = r26042 / r26043;
        double r26045 = M;
        double r26046 = r26044 - r26045;
        double r26047 = cos(r26046);
        double r26048 = r26041 / r26043;
        double r26049 = r26048 - r26045;
        double r26050 = pow(r26049, r26043);
        double r26051 = -r26050;
        double r26052 = l;
        double r26053 = r26039 - r26040;
        double r26054 = fabs(r26053);
        double r26055 = r26052 - r26054;
        double r26056 = r26051 - r26055;
        double r26057 = exp(r26056);
        double r26058 = r26047 * r26057;
        return r26058;
}


double f_of(float __attribute__((unused)) K, float m, float n, float M, float l) {
        float r26059 = m;
        float r26060 = n;
        float r26061 = r26059 + r26060;
        float r26062 = 2;
        float r26063 = r26061 / r26062;
        float r26064 = M;
        float r26065 = r26063 - r26064;
        float r26066 = pow(r26065, r26062);
        float r26067 = -r26066;
        float r26068 = l;
        float r26069 = r26059 - r26060;
        float r26070 = fabs(r26069);
        float r26071 = r26068 - r26070;
        float r26072 = r26067 - r26071;
        float r26073 = exp(r26072);
        return r26073;
}

double f_od(double __attribute__((unused)) K, double m, double n, double M, double l) {
        double r26074 = m;
        double r26075 = n;
        double r26076 = r26074 + r26075;
        double r26077 = 2;
        double r26078 = r26076 / r26077;
        double r26079 = M;
        double r26080 = r26078 - r26079;
        double r26081 = pow(r26080, r26077);
        double r26082 = -r26081;
        double r26083 = l;
        double r26084 = r26074 - r26075;
        double r26085 = fabs(r26084);
        double r26086 = r26083 - r26085;
        double r26087 = r26082 - r26086;
        double r26088 = exp(r26087);
        return r26088;
}

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 r26089, r26090, r26091, r26092, r26093, r26094, r26095, r26096, r26097, r26098, r26099, r26100, r26101, r26102, r26103, r26104, r26105, r26106, r26107, r26108, r26109;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26089);
        mpfr_init(r26090);
        mpfr_init(r26091);
        mpfr_init(r26092);
        mpfr_init(r26093);
        mpfr_init_set_str(r26094, "2", 10, MPFR_RNDN);
        mpfr_init(r26095);
        mpfr_init(r26096);
        mpfr_init(r26097);
        mpfr_init(r26098);
        mpfr_init(r26099);
        mpfr_init(r26100);
        mpfr_init(r26101);
        mpfr_init(r26102);
        mpfr_init(r26103);
        mpfr_init(r26104);
        mpfr_init(r26105);
        mpfr_init(r26106);
        mpfr_init(r26107);
        mpfr_init(r26108);
        mpfr_init(r26109);
}

double f_im(double K, double m, double n, double M, double l) {
        mpfr_set_d(r26089, K, MPFR_RNDN);
        mpfr_set_d(r26090, m, MPFR_RNDN);
        mpfr_set_d(r26091, n, MPFR_RNDN);
        mpfr_add(r26092, r26090, r26091, MPFR_RNDN);
        mpfr_mul(r26093, r26089, r26092, MPFR_RNDN);
        ;
        mpfr_div(r26095, r26093, r26094, MPFR_RNDN);
        mpfr_set_d(r26096, M, MPFR_RNDN);
        mpfr_sub(r26097, r26095, r26096, MPFR_RNDN);
        mpfr_cos(r26098, r26097, MPFR_RNDN);
        mpfr_div(r26099, r26092, r26094, MPFR_RNDN);
        mpfr_sub(r26100, r26099, r26096, MPFR_RNDN);
        mpfr_pow(r26101, r26100, r26094, MPFR_RNDN);
        mpfr_neg(r26102, r26101, MPFR_RNDN);
        mpfr_set_d(r26103, l, MPFR_RNDN);
        mpfr_sub(r26104, r26090, r26091, MPFR_RNDN);
        mpfr_abs(r26105, r26104, MPFR_RNDN);
        mpfr_sub(r26106, r26103, r26105, MPFR_RNDN);
        mpfr_sub(r26107, r26102, r26106, MPFR_RNDN);
        mpfr_exp(r26108, r26107, MPFR_RNDN);
        mpfr_mul(r26109, r26098, r26108, MPFR_RNDN);
        return mpfr_get_d(r26109, MPFR_RNDN);
}

static mpfr_t r26110, r26111, r26112, r26113, r26114, r26115, r26116, r26117, r26118, r26119, r26120, r26121, r26122, r26123, r26124;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26110);
        mpfr_init(r26111);
        mpfr_init(r26112);
        mpfr_init_set_str(r26113, "2", 10, MPFR_RNDN);
        mpfr_init(r26114);
        mpfr_init(r26115);
        mpfr_init(r26116);
        mpfr_init(r26117);
        mpfr_init(r26118);
        mpfr_init(r26119);
        mpfr_init(r26120);
        mpfr_init(r26121);
        mpfr_init(r26122);
        mpfr_init(r26123);
        mpfr_init(r26124);
}

double f_fm(double __attribute__((unused)) K, double m, double n, double M, double l) {
        mpfr_set_d(r26110, m, MPFR_RNDN);
        mpfr_set_d(r26111, n, MPFR_RNDN);
        mpfr_add(r26112, r26110, r26111, MPFR_RNDN);
        ;
        mpfr_div(r26114, r26112, r26113, MPFR_RNDN);
        mpfr_set_d(r26115, M, MPFR_RNDN);
        mpfr_sub(r26116, r26114, r26115, MPFR_RNDN);
        mpfr_pow(r26117, r26116, r26113, MPFR_RNDN);
        mpfr_neg(r26118, r26117, MPFR_RNDN);
        mpfr_set_d(r26119, l, MPFR_RNDN);
        mpfr_sub(r26120, r26110, r26111, MPFR_RNDN);
        mpfr_abs(r26121, r26120, MPFR_RNDN);
        mpfr_sub(r26122, r26119, r26121, MPFR_RNDN);
        mpfr_sub(r26123, r26118, r26122, MPFR_RNDN);
        mpfr_exp(r26124, r26123, MPFR_RNDN);
        return mpfr_get_d(r26124, MPFR_RNDN);
}

static mpfr_t r26125, r26126, r26127, r26128, r26129, r26130, r26131, r26132, r26133, r26134, r26135, r26136, r26137, r26138, r26139;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26125);
        mpfr_init(r26126);
        mpfr_init(r26127);
        mpfr_init_set_str(r26128, "2", 10, MPFR_RNDN);
        mpfr_init(r26129);
        mpfr_init(r26130);
        mpfr_init(r26131);
        mpfr_init(r26132);
        mpfr_init(r26133);
        mpfr_init(r26134);
        mpfr_init(r26135);
        mpfr_init(r26136);
        mpfr_init(r26137);
        mpfr_init(r26138);
        mpfr_init(r26139);
}

double f_dm(double __attribute__((unused)) K, double m, double n, double M, double l) {
        mpfr_set_d(r26125, m, MPFR_RNDN);
        mpfr_set_d(r26126, n, MPFR_RNDN);
        mpfr_add(r26127, r26125, r26126, MPFR_RNDN);
        ;
        mpfr_div(r26129, r26127, r26128, MPFR_RNDN);
        mpfr_set_d(r26130, M, MPFR_RNDN);
        mpfr_sub(r26131, r26129, r26130, MPFR_RNDN);
        mpfr_pow(r26132, r26131, r26128, MPFR_RNDN);
        mpfr_neg(r26133, r26132, MPFR_RNDN);
        mpfr_set_d(r26134, l, MPFR_RNDN);
        mpfr_sub(r26135, r26125, r26126, MPFR_RNDN);
        mpfr_abs(r26136, r26135, MPFR_RNDN);
        mpfr_sub(r26137, r26134, r26136, MPFR_RNDN);
        mpfr_sub(r26138, r26133, r26137, MPFR_RNDN);
        mpfr_exp(r26139, r26138, MPFR_RNDN);
        return mpfr_get_d(r26139, MPFR_RNDN);
}