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

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

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

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


double f_of(float J, float K, float U) {
        float r26032 = -2;
        float r26033 = J;
        float r26034 = r26032 * r26033;
        float r26035 = K;
        float r26036 = 2;
        float r26037 = r26035 / r26036;
        float r26038 = cos(r26037);
        float r26039 = r26034 * r26038;
        float r26040 = 1;
        float r26041 = U;
        float r26042 = r26036 * r26033;
        float r26043 = r26042 * r26038;
        float r26044 = r26041 / r26043;
        float r26045 = pow(r26044, r26036);
        float r26046 = r26040 + r26045;
        float r26047 = sqrt(r26046);
        float r26048 = r26039 * r26047;
        return r26048;
}

double f_od(double J, double K, double U) {
        double r26049 = -2;
        double r26050 = J;
        double r26051 = r26049 * r26050;
        double r26052 = K;
        double r26053 = 2;
        double r26054 = r26052 / r26053;
        double r26055 = cos(r26054);
        double r26056 = r26051 * r26055;
        double r26057 = 1;
        double r26058 = U;
        double r26059 = r26053 * r26050;
        double r26060 = r26059 * r26055;
        double r26061 = r26058 / r26060;
        double r26062 = pow(r26061, r26053);
        double r26063 = r26057 + r26062;
        double r26064 = sqrt(r26063);
        double r26065 = r26056 * r26064;
        return r26065;
}

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 r26066, r26067, r26068, r26069, r26070, r26071, r26072, r26073, r26074, r26075, r26076, r26077, r26078, r26079, r26080, r26081, r26082;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r26066, "-2", 10, MPFR_RNDN);
        mpfr_init(r26067);
        mpfr_init(r26068);
        mpfr_init(r26069);
        mpfr_init_set_str(r26070, "2", 10, MPFR_RNDN);
        mpfr_init(r26071);
        mpfr_init(r26072);
        mpfr_init(r26073);
        mpfr_init_set_str(r26074, "1", 10, MPFR_RNDN);
        mpfr_init(r26075);
        mpfr_init(r26076);
        mpfr_init(r26077);
        mpfr_init(r26078);
        mpfr_init(r26079);
        mpfr_init(r26080);
        mpfr_init(r26081);
        mpfr_init(r26082);
}

double f_im(double J, double K, double U) {
        ;
        mpfr_set_d(r26067, J, MPFR_RNDN);
        mpfr_mul(r26068, r26066, r26067, MPFR_RNDN);
        mpfr_set_d(r26069, K, MPFR_RNDN);
        ;
        mpfr_div(r26071, r26069, r26070, MPFR_RNDN);
        mpfr_cos(r26072, r26071, MPFR_RNDN);
        mpfr_mul(r26073, r26068, r26072, MPFR_RNDN);
        ;
        mpfr_set_d(r26075, U, MPFR_RNDN);
        mpfr_mul(r26076, r26070, r26067, MPFR_RNDN);
        mpfr_mul(r26077, r26076, r26072, MPFR_RNDN);
        mpfr_div(r26078, r26075, r26077, MPFR_RNDN);
        mpfr_pow(r26079, r26078, r26070, MPFR_RNDN);
        mpfr_add(r26080, r26074, r26079, MPFR_RNDN);
        mpfr_sqrt(r26081, r26080, MPFR_RNDN);
        mpfr_mul(r26082, r26073, r26081, MPFR_RNDN);
        return mpfr_get_d(r26082, MPFR_RNDN);
}

static mpfr_t r26083, r26084, r26085, r26086, r26087, r26088, r26089, r26090, r26091, r26092, r26093, r26094, r26095, r26096, r26097, r26098, r26099;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r26083, "-2", 10, MPFR_RNDN);
        mpfr_init(r26084);
        mpfr_init(r26085);
        mpfr_init(r26086);
        mpfr_init_set_str(r26087, "2", 10, MPFR_RNDN);
        mpfr_init(r26088);
        mpfr_init(r26089);
        mpfr_init(r26090);
        mpfr_init_set_str(r26091, "1", 10, MPFR_RNDN);
        mpfr_init(r26092);
        mpfr_init(r26093);
        mpfr_init(r26094);
        mpfr_init(r26095);
        mpfr_init(r26096);
        mpfr_init(r26097);
        mpfr_init(r26098);
        mpfr_init(r26099);
}

double f_fm(double J, double K, double U) {
        ;
        mpfr_set_d(r26084, J, MPFR_RNDN);
        mpfr_mul(r26085, r26083, r26084, MPFR_RNDN);
        mpfr_set_d(r26086, K, MPFR_RNDN);
        ;
        mpfr_div(r26088, r26086, r26087, MPFR_RNDN);
        mpfr_cos(r26089, r26088, MPFR_RNDN);
        mpfr_mul(r26090, r26085, r26089, MPFR_RNDN);
        ;
        mpfr_set_d(r26092, U, MPFR_RNDN);
        mpfr_mul(r26093, r26087, r26084, MPFR_RNDN);
        mpfr_mul(r26094, r26093, r26089, MPFR_RNDN);
        mpfr_div(r26095, r26092, r26094, MPFR_RNDN);
        mpfr_pow(r26096, r26095, r26087, MPFR_RNDN);
        mpfr_add(r26097, r26091, r26096, MPFR_RNDN);
        mpfr_sqrt(r26098, r26097, MPFR_RNDN);
        mpfr_mul(r26099, r26090, r26098, MPFR_RNDN);
        return mpfr_get_d(r26099, MPFR_RNDN);
}

static mpfr_t r26100, r26101, r26102, r26103, r26104, r26105, r26106, r26107, r26108, r26109, r26110, r26111, r26112, r26113, r26114, r26115, r26116;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r26100, "-2", 10, MPFR_RNDN);
        mpfr_init(r26101);
        mpfr_init(r26102);
        mpfr_init(r26103);
        mpfr_init_set_str(r26104, "2", 10, MPFR_RNDN);
        mpfr_init(r26105);
        mpfr_init(r26106);
        mpfr_init(r26107);
        mpfr_init_set_str(r26108, "1", 10, MPFR_RNDN);
        mpfr_init(r26109);
        mpfr_init(r26110);
        mpfr_init(r26111);
        mpfr_init(r26112);
        mpfr_init(r26113);
        mpfr_init(r26114);
        mpfr_init(r26115);
        mpfr_init(r26116);
}

double f_dm(double J, double K, double U) {
        ;
        mpfr_set_d(r26101, J, MPFR_RNDN);
        mpfr_mul(r26102, r26100, r26101, MPFR_RNDN);
        mpfr_set_d(r26103, K, MPFR_RNDN);
        ;
        mpfr_div(r26105, r26103, r26104, MPFR_RNDN);
        mpfr_cos(r26106, r26105, MPFR_RNDN);
        mpfr_mul(r26107, r26102, r26106, MPFR_RNDN);
        ;
        mpfr_set_d(r26109, U, MPFR_RNDN);
        mpfr_mul(r26110, r26104, r26101, MPFR_RNDN);
        mpfr_mul(r26111, r26110, r26106, MPFR_RNDN);
        mpfr_div(r26112, r26109, r26111, MPFR_RNDN);
        mpfr_pow(r26113, r26112, r26104, MPFR_RNDN);
        mpfr_add(r26114, r26108, r26113, MPFR_RNDN);
        mpfr_sqrt(r26115, r26114, MPFR_RNDN);
        mpfr_mul(r26116, r26107, r26115, MPFR_RNDN);
        return mpfr_get_d(r26116, MPFR_RNDN);
}