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

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

double f_if(float J, float l, float K, float U) {
        float r26001 = J;
        float r26002 = l;
        float r26003 = exp(r26002);
        float r26004 = -r26002;
        float r26005 = exp(r26004);
        float r26006 = r26003 - r26005;
        float r26007 = r26001 * r26006;
        float r26008 = K;
        float r26009 = 2;
        float r26010 = r26008 / r26009;
        float r26011 = cos(r26010);
        float r26012 = r26007 * r26011;
        float r26013 = U;
        float r26014 = r26012 + r26013;
        return r26014;
}

double f_id(double J, double l, double K, double U) {
        double r26015 = J;
        double r26016 = l;
        double r26017 = exp(r26016);
        double r26018 = -r26016;
        double r26019 = exp(r26018);
        double r26020 = r26017 - r26019;
        double r26021 = r26015 * r26020;
        double r26022 = K;
        double r26023 = 2;
        double r26024 = r26022 / r26023;
        double r26025 = cos(r26024);
        double r26026 = r26021 * r26025;
        double r26027 = U;
        double r26028 = r26026 + r26027;
        return r26028;
}


double f_of(float J, float l, float K, float U) {
        float r26029 = K;
        float r26030 = 2;
        float r26031 = r26029 / r26030;
        float r26032 = cos(r26031);
        float r26033 = cbrt(r26032);
        float r26034 = r26033 * r26033;
        float r26035 = exp(r26033);
        float r26036 = log(r26035);
        float r26037 = J;
        float r26038 = r26036 * r26037;
        float r26039 = r26034 * r26038;
        float r26040 = l;
        float r26041 = 1/3;
        float r26042 = r26040 * r26041;
        float r26043 = fma(r26042, r26040, r26030);
        float r26044 = 5;
        float r26045 = pow(r26040, r26044);
        float r26046 = 1/60;
        float r26047 = r26045 * r26046;
        float r26048 = fma(r26040, r26043, r26047);
        float r26049 = U;
        float r26050 = fma(r26039, r26048, r26049);
        return r26050;
}

double f_od(double J, double l, double K, double U) {
        double r26051 = K;
        double r26052 = 2;
        double r26053 = r26051 / r26052;
        double r26054 = cos(r26053);
        double r26055 = cbrt(r26054);
        double r26056 = r26055 * r26055;
        double r26057 = exp(r26055);
        double r26058 = log(r26057);
        double r26059 = J;
        double r26060 = r26058 * r26059;
        double r26061 = r26056 * r26060;
        double r26062 = l;
        double r26063 = 1/3;
        double r26064 = r26062 * r26063;
        double r26065 = fma(r26064, r26062, r26052);
        double r26066 = 5;
        double r26067 = pow(r26062, r26066);
        double r26068 = 1/60;
        double r26069 = r26067 * r26068;
        double r26070 = fma(r26062, r26065, r26069);
        double r26071 = U;
        double r26072 = fma(r26061, r26070, r26071);
        return r26072;
}

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 r26073, r26074, r26075, r26076, r26077, r26078, r26079, r26080, r26081, r26082, r26083, r26084, r26085, r26086;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26073);
        mpfr_init(r26074);
        mpfr_init(r26075);
        mpfr_init(r26076);
        mpfr_init(r26077);
        mpfr_init(r26078);
        mpfr_init(r26079);
        mpfr_init(r26080);
        mpfr_init_set_str(r26081, "2", 10, MPFR_RNDN);
        mpfr_init(r26082);
        mpfr_init(r26083);
        mpfr_init(r26084);
        mpfr_init(r26085);
        mpfr_init(r26086);
}

double f_im(double J, double l, double K, double U) {
        mpfr_set_d(r26073, J, MPFR_RNDN);
        mpfr_set_d(r26074, l, MPFR_RNDN);
        mpfr_exp(r26075, r26074, MPFR_RNDN);
        mpfr_neg(r26076, r26074, MPFR_RNDN);
        mpfr_exp(r26077, r26076, MPFR_RNDN);
        mpfr_sub(r26078, r26075, r26077, MPFR_RNDN);
        mpfr_mul(r26079, r26073, r26078, MPFR_RNDN);
        mpfr_set_d(r26080, K, MPFR_RNDN);
        ;
        mpfr_div(r26082, r26080, r26081, MPFR_RNDN);
        mpfr_cos(r26083, r26082, MPFR_RNDN);
        mpfr_mul(r26084, r26079, r26083, MPFR_RNDN);
        mpfr_set_d(r26085, U, MPFR_RNDN);
        mpfr_add(r26086, r26084, r26085, MPFR_RNDN);
        return mpfr_get_d(r26086, MPFR_RNDN);
}

static mpfr_t r26087, r26088, r26089, r26090, r26091, r26092, r26093, r26094, r26095, r26096, r26097, r26098, r26099, r26100, r26101, r26102, r26103, r26104, r26105, r26106, r26107, r26108;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26087);
        mpfr_init_set_str(r26088, "2", 10, MPFR_RNDN);
        mpfr_init(r26089);
        mpfr_init(r26090);
        mpfr_init(r26091);
        mpfr_init(r26092);
        mpfr_init(r26093);
        mpfr_init(r26094);
        mpfr_init(r26095);
        mpfr_init(r26096);
        mpfr_init(r26097);
        mpfr_init(r26098);
        mpfr_init_set_str(r26099, "1/3", 10, MPFR_RNDN);
        mpfr_init(r26100);
        mpfr_init(r26101);
        mpfr_init_set_str(r26102, "5", 10, MPFR_RNDN);
        mpfr_init(r26103);
        mpfr_init_set_str(r26104, "1/60", 10, MPFR_RNDN);
        mpfr_init(r26105);
        mpfr_init(r26106);
        mpfr_init(r26107);
        mpfr_init(r26108);
}

double f_fm(double J, double l, double K, double U) {
        mpfr_set_d(r26087, K, MPFR_RNDN);
        ;
        mpfr_div(r26089, r26087, r26088, MPFR_RNDN);
        mpfr_cos(r26090, r26089, MPFR_RNDN);
        mpfr_cbrt(r26091, r26090, MPFR_RNDN);
        mpfr_mul(r26092, r26091, r26091, MPFR_RNDN);
        mpfr_exp(r26093, r26091, MPFR_RNDN);
        mpfr_log(r26094, r26093, MPFR_RNDN);
        mpfr_set_d(r26095, J, MPFR_RNDN);
        mpfr_mul(r26096, r26094, r26095, MPFR_RNDN);
        mpfr_mul(r26097, r26092, r26096, MPFR_RNDN);
        mpfr_set_d(r26098, l, MPFR_RNDN);
        ;
        mpfr_mul(r26100, r26098, r26099, MPFR_RNDN);
        mpfr_fma(r26101, r26100, r26098, r26088, MPFR_RNDN);
        ;
        mpfr_pow(r26103, r26098, r26102, MPFR_RNDN);
        ;
        mpfr_mul(r26105, r26103, r26104, MPFR_RNDN);
        mpfr_fma(r26106, r26098, r26101, r26105, MPFR_RNDN);
        mpfr_set_d(r26107, U, MPFR_RNDN);
        mpfr_fma(r26108, r26097, r26106, r26107, MPFR_RNDN);
        return mpfr_get_d(r26108, MPFR_RNDN);
}

static mpfr_t r26109, r26110, r26111, r26112, r26113, r26114, r26115, r26116, r26117, r26118, r26119, r26120, r26121, r26122, r26123, r26124, r26125, r26126, r26127, r26128, r26129, r26130;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26109);
        mpfr_init_set_str(r26110, "2", 10, MPFR_RNDN);
        mpfr_init(r26111);
        mpfr_init(r26112);
        mpfr_init(r26113);
        mpfr_init(r26114);
        mpfr_init(r26115);
        mpfr_init(r26116);
        mpfr_init(r26117);
        mpfr_init(r26118);
        mpfr_init(r26119);
        mpfr_init(r26120);
        mpfr_init_set_str(r26121, "1/3", 10, MPFR_RNDN);
        mpfr_init(r26122);
        mpfr_init(r26123);
        mpfr_init_set_str(r26124, "5", 10, MPFR_RNDN);
        mpfr_init(r26125);
        mpfr_init_set_str(r26126, "1/60", 10, MPFR_RNDN);
        mpfr_init(r26127);
        mpfr_init(r26128);
        mpfr_init(r26129);
        mpfr_init(r26130);
}

double f_dm(double J, double l, double K, double U) {
        mpfr_set_d(r26109, K, MPFR_RNDN);
        ;
        mpfr_div(r26111, r26109, r26110, MPFR_RNDN);
        mpfr_cos(r26112, r26111, MPFR_RNDN);
        mpfr_cbrt(r26113, r26112, MPFR_RNDN);
        mpfr_mul(r26114, r26113, r26113, MPFR_RNDN);
        mpfr_exp(r26115, r26113, MPFR_RNDN);
        mpfr_log(r26116, r26115, MPFR_RNDN);
        mpfr_set_d(r26117, J, MPFR_RNDN);
        mpfr_mul(r26118, r26116, r26117, MPFR_RNDN);
        mpfr_mul(r26119, r26114, r26118, MPFR_RNDN);
        mpfr_set_d(r26120, l, MPFR_RNDN);
        ;
        mpfr_mul(r26122, r26120, r26121, MPFR_RNDN);
        mpfr_fma(r26123, r26122, r26120, r26110, MPFR_RNDN);
        ;
        mpfr_pow(r26125, r26120, r26124, MPFR_RNDN);
        ;
        mpfr_mul(r26127, r26125, r26126, MPFR_RNDN);
        mpfr_fma(r26128, r26120, r26123, r26127, MPFR_RNDN);
        mpfr_set_d(r26129, U, MPFR_RNDN);
        mpfr_fma(r26130, r26119, r26128, r26129, MPFR_RNDN);
        return mpfr_get_d(r26130, MPFR_RNDN);
}