#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 r25890 = -2;
        float r25891 = J;
        float r25892 = r25890 * r25891;
        float r25893 = K;
        float r25894 = 2;
        float r25895 = r25893 / r25894;
        float r25896 = cos(r25895);
        float r25897 = r25892 * r25896;
        float r25898 = 1;
        float r25899 = U;
        float r25900 = r25894 * r25891;
        float r25901 = r25900 * r25896;
        float r25902 = r25899 / r25901;
        float r25903 = pow(r25902, r25894);
        float r25904 = r25898 + r25903;
        float r25905 = sqrt(r25904);
        float r25906 = r25897 * r25905;
        return r25906;
}

double f_id(double J, double K, double U) {
        double r25907 = -2;
        double r25908 = J;
        double r25909 = r25907 * r25908;
        double r25910 = K;
        double r25911 = 2;
        double r25912 = r25910 / r25911;
        double r25913 = cos(r25912);
        double r25914 = r25909 * r25913;
        double r25915 = 1;
        double r25916 = U;
        double r25917 = r25911 * r25908;
        double r25918 = r25917 * r25913;
        double r25919 = r25916 / r25918;
        double r25920 = pow(r25919, r25911);
        double r25921 = r25915 + r25920;
        double r25922 = sqrt(r25921);
        double r25923 = r25914 * r25922;
        return r25923;
}


double f_of(float J, float K, float U) {
        float r25924 = J;
        float r25925 = -2;
        float r25926 = r25924 * r25925;
        float r25927 = K;
        float r25928 = 2;
        float r25929 = r25927 / r25928;
        float r25930 = cos(r25929);
        float r25931 = r25926 * r25930;
        float r25932 = cbrt(r25931);
        float r25933 = r25932 * r25932;
        float r25934 = r25933 * r25932;
        float r25935 = 1;
        float r25936 = U;
        float r25937 = r25936 / r25928;
        float r25938 = r25937 / r25924;
        float r25939 = r25938 / r25930;
        float r25940 = hypot(r25935, r25939);
        float r25941 = r25934 * r25940;
        float r25942 = -inf.0;
        bool r25943 = r25941 <= r25942;
        float r25944 = 1/2;
        float r25945 = r25944 * r25936;
        float r25946 = r25925 * r25945;
        float r25947 = 1.7976931348623157e+308;
        bool r25948 = r25941 <= r25947;
        float r25949 = r25930 * r25940;
        float r25950 = r25926 * r25949;
        float r25951 = r25948 ? r25950 : r25946;
        float r25952 = r25943 ? r25946 : r25951;
        return r25952;
}

double f_od(double J, double K, double U) {
        double r25953 = J;
        double r25954 = -2;
        double r25955 = r25953 * r25954;
        double r25956 = K;
        double r25957 = 2;
        double r25958 = r25956 / r25957;
        double r25959 = cos(r25958);
        double r25960 = r25955 * r25959;
        double r25961 = cbrt(r25960);
        double r25962 = r25961 * r25961;
        double r25963 = r25962 * r25961;
        double r25964 = 1;
        double r25965 = U;
        double r25966 = r25965 / r25957;
        double r25967 = r25966 / r25953;
        double r25968 = r25967 / r25959;
        double r25969 = hypot(r25964, r25968);
        double r25970 = r25963 * r25969;
        double r25971 = -inf.0;
        bool r25972 = r25970 <= r25971;
        double r25973 = 1/2;
        double r25974 = r25973 * r25965;
        double r25975 = r25954 * r25974;
        double r25976 = 1.7976931348623157e+308;
        bool r25977 = r25970 <= r25976;
        double r25978 = r25959 * r25969;
        double r25979 = r25955 * r25978;
        double r25980 = r25977 ? r25979 : r25975;
        double r25981 = r25972 ? r25975 : r25980;
        return r25981;
}

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 r25982, r25983, r25984, r25985, r25986, r25987, r25988, r25989, r25990, r25991, r25992, r25993, r25994, r25995, r25996, r25997, r25998;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r25982, "-2", 10, MPFR_RNDN);
        mpfr_init(r25983);
        mpfr_init(r25984);
        mpfr_init(r25985);
        mpfr_init_set_str(r25986, "2", 10, MPFR_RNDN);
        mpfr_init(r25987);
        mpfr_init(r25988);
        mpfr_init(r25989);
        mpfr_init_set_str(r25990, "1", 10, MPFR_RNDN);
        mpfr_init(r25991);
        mpfr_init(r25992);
        mpfr_init(r25993);
        mpfr_init(r25994);
        mpfr_init(r25995);
        mpfr_init(r25996);
        mpfr_init(r25997);
        mpfr_init(r25998);
}

double f_im(double J, double K, double U) {
        ;
        mpfr_set_d(r25983, J, MPFR_RNDN);
        mpfr_mul(r25984, r25982, r25983, MPFR_RNDN);
        mpfr_set_d(r25985, K, MPFR_RNDN);
        ;
        mpfr_div(r25987, r25985, r25986, MPFR_RNDN);
        mpfr_cos(r25988, r25987, MPFR_RNDN);
        mpfr_mul(r25989, r25984, r25988, MPFR_RNDN);
        ;
        mpfr_set_d(r25991, U, MPFR_RNDN);
        mpfr_mul(r25992, r25986, r25983, MPFR_RNDN);
        mpfr_mul(r25993, r25992, r25988, MPFR_RNDN);
        mpfr_div(r25994, r25991, r25993, MPFR_RNDN);
        mpfr_pow(r25995, r25994, r25986, MPFR_RNDN);
        mpfr_add(r25996, r25990, r25995, MPFR_RNDN);
        mpfr_sqrt(r25997, r25996, MPFR_RNDN);
        mpfr_mul(r25998, r25989, r25997, MPFR_RNDN);
        return mpfr_get_d(r25998, MPFR_RNDN);
}

static mpfr_t r25999, r26000, r26001, r26002, r26003, r26004, r26005, r26006, r26007, r26008, r26009, r26010, r26011, r26012, r26013, r26014, r26015, r26016, r26017, r26018, r26019, r26020, r26021, r26022, r26023, r26024, r26025, r26026, r26027;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25999);
        mpfr_init_set_str(r26000, "-2", 10, MPFR_RNDN);
        mpfr_init(r26001);
        mpfr_init(r26002);
        mpfr_init_set_str(r26003, "2", 10, MPFR_RNDN);
        mpfr_init(r26004);
        mpfr_init(r26005);
        mpfr_init(r26006);
        mpfr_init(r26007);
        mpfr_init(r26008);
        mpfr_init(r26009);
        mpfr_init_set_str(r26010, "1", 10, MPFR_RNDN);
        mpfr_init(r26011);
        mpfr_init(r26012);
        mpfr_init(r26013);
        mpfr_init(r26014);
        mpfr_init(r26015);
        mpfr_init(r26016);
        mpfr_init_set_str(r26017, "-inf.0", 10, MPFR_RNDN);
        mpfr_init(r26018);
        mpfr_init_set_str(r26019, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26020);
        mpfr_init(r26021);
        mpfr_init_set_str(r26022, "1.7976931348623157e+308", 10, MPFR_RNDN);
        mpfr_init(r26023);
        mpfr_init(r26024);
        mpfr_init(r26025);
        mpfr_init(r26026);
        mpfr_init(r26027);
}

double f_fm(double J, double K, double U) {
        mpfr_set_d(r25999, J, MPFR_RNDN);
        ;
        mpfr_mul(r26001, r25999, r26000, MPFR_RNDN);
        mpfr_set_d(r26002, K, MPFR_RNDN);
        ;
        mpfr_div(r26004, r26002, r26003, MPFR_RNDN);
        mpfr_cos(r26005, r26004, MPFR_RNDN);
        mpfr_mul(r26006, r26001, r26005, MPFR_RNDN);
        mpfr_cbrt(r26007, r26006, MPFR_RNDN);
        mpfr_mul(r26008, r26007, r26007, MPFR_RNDN);
        mpfr_mul(r26009, r26008, r26007, MPFR_RNDN);
        ;
        mpfr_set_d(r26011, U, MPFR_RNDN);
        mpfr_div(r26012, r26011, r26003, MPFR_RNDN);
        mpfr_div(r26013, r26012, r25999, MPFR_RNDN);
        mpfr_div(r26014, r26013, r26005, MPFR_RNDN);
        mpfr_hypot(r26015, r26010, r26014, MPFR_RNDN);
        mpfr_mul(r26016, r26009, r26015, MPFR_RNDN);
        ;
        mpfr_set_si(r26018, mpfr_cmp(r26016, r26017) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r26020, r26019, r26011, MPFR_RNDN);
        mpfr_mul(r26021, r26000, r26020, MPFR_RNDN);
        ;
        mpfr_set_si(r26023, mpfr_cmp(r26016, r26022) <= 0, MPFR_RNDN);
        mpfr_mul(r26024, r26005, r26015, MPFR_RNDN);
        mpfr_mul(r26025, r26001, r26024, MPFR_RNDN);
        if (mpfr_get_si(r26023, MPFR_RNDN)) { mpfr_set(r26026, r26025, MPFR_RNDN); } else { mpfr_set(r26026, r26021, MPFR_RNDN); };
        if (mpfr_get_si(r26018, MPFR_RNDN)) { mpfr_set(r26027, r26021, MPFR_RNDN); } else { mpfr_set(r26027, r26026, MPFR_RNDN); };
        return mpfr_get_d(r26027, MPFR_RNDN);
}

static mpfr_t r26028, r26029, r26030, r26031, r26032, r26033, r26034, r26035, r26036, r26037, r26038, r26039, r26040, r26041, r26042, r26043, r26044, r26045, r26046, r26047, r26048, r26049, r26050, r26051, r26052, r26053, r26054, r26055, r26056;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26028);
        mpfr_init_set_str(r26029, "-2", 10, MPFR_RNDN);
        mpfr_init(r26030);
        mpfr_init(r26031);
        mpfr_init_set_str(r26032, "2", 10, MPFR_RNDN);
        mpfr_init(r26033);
        mpfr_init(r26034);
        mpfr_init(r26035);
        mpfr_init(r26036);
        mpfr_init(r26037);
        mpfr_init(r26038);
        mpfr_init_set_str(r26039, "1", 10, MPFR_RNDN);
        mpfr_init(r26040);
        mpfr_init(r26041);
        mpfr_init(r26042);
        mpfr_init(r26043);
        mpfr_init(r26044);
        mpfr_init(r26045);
        mpfr_init_set_str(r26046, "-inf.0", 10, MPFR_RNDN);
        mpfr_init(r26047);
        mpfr_init_set_str(r26048, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26049);
        mpfr_init(r26050);
        mpfr_init_set_str(r26051, "1.7976931348623157e+308", 10, MPFR_RNDN);
        mpfr_init(r26052);
        mpfr_init(r26053);
        mpfr_init(r26054);
        mpfr_init(r26055);
        mpfr_init(r26056);
}

double f_dm(double J, double K, double U) {
        mpfr_set_d(r26028, J, MPFR_RNDN);
        ;
        mpfr_mul(r26030, r26028, r26029, MPFR_RNDN);
        mpfr_set_d(r26031, K, MPFR_RNDN);
        ;
        mpfr_div(r26033, r26031, r26032, MPFR_RNDN);
        mpfr_cos(r26034, r26033, MPFR_RNDN);
        mpfr_mul(r26035, r26030, r26034, MPFR_RNDN);
        mpfr_cbrt(r26036, r26035, MPFR_RNDN);
        mpfr_mul(r26037, r26036, r26036, MPFR_RNDN);
        mpfr_mul(r26038, r26037, r26036, MPFR_RNDN);
        ;
        mpfr_set_d(r26040, U, MPFR_RNDN);
        mpfr_div(r26041, r26040, r26032, MPFR_RNDN);
        mpfr_div(r26042, r26041, r26028, MPFR_RNDN);
        mpfr_div(r26043, r26042, r26034, MPFR_RNDN);
        mpfr_hypot(r26044, r26039, r26043, MPFR_RNDN);
        mpfr_mul(r26045, r26038, r26044, MPFR_RNDN);
        ;
        mpfr_set_si(r26047, mpfr_cmp(r26045, r26046) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r26049, r26048, r26040, MPFR_RNDN);
        mpfr_mul(r26050, r26029, r26049, MPFR_RNDN);
        ;
        mpfr_set_si(r26052, mpfr_cmp(r26045, r26051) <= 0, MPFR_RNDN);
        mpfr_mul(r26053, r26034, r26044, MPFR_RNDN);
        mpfr_mul(r26054, r26030, r26053, MPFR_RNDN);
        if (mpfr_get_si(r26052, MPFR_RNDN)) { mpfr_set(r26055, r26054, MPFR_RNDN); } else { mpfr_set(r26055, r26050, MPFR_RNDN); };
        if (mpfr_get_si(r26047, MPFR_RNDN)) { mpfr_set(r26056, r26050, MPFR_RNDN); } else { mpfr_set(r26056, r26055, MPFR_RNDN); };
        return mpfr_get_d(r26056, MPFR_RNDN);
}