#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 r23283 = -2;
        float r23284 = J;
        float r23285 = r23283 * r23284;
        float r23286 = K;
        float r23287 = 2;
        float r23288 = r23286 / r23287;
        float r23289 = cos(r23288);
        float r23290 = r23285 * r23289;
        float r23291 = 1;
        float r23292 = U;
        float r23293 = r23287 * r23284;
        float r23294 = r23293 * r23289;
        float r23295 = r23292 / r23294;
        float r23296 = pow(r23295, r23287);
        float r23297 = r23291 + r23296;
        float r23298 = sqrt(r23297);
        float r23299 = r23290 * r23298;
        return r23299;
}

double f_id(double J, double K, double U) {
        double r23300 = -2;
        double r23301 = J;
        double r23302 = r23300 * r23301;
        double r23303 = K;
        double r23304 = 2;
        double r23305 = r23303 / r23304;
        double r23306 = cos(r23305);
        double r23307 = r23302 * r23306;
        double r23308 = 1;
        double r23309 = U;
        double r23310 = r23304 * r23301;
        double r23311 = r23310 * r23306;
        double r23312 = r23309 / r23311;
        double r23313 = pow(r23312, r23304);
        double r23314 = r23308 + r23313;
        double r23315 = sqrt(r23314);
        double r23316 = r23307 * r23315;
        return r23316;
}


double f_of(float J, float K, float U) {
        float r23317 = 1;
        float r23318 = U;
        float r23319 = J;
        float r23320 = r23318 / r23319;
        float r23321 = 2;
        float r23322 = r23320 / r23321;
        float r23323 = K;
        float r23324 = r23323 / r23321;
        float r23325 = cos(r23324);
        float r23326 = r23322 / r23325;
        float r23327 = hypot(r23317, r23326);
        float r23328 = -2;
        float r23329 = r23328 * r23319;
        float r23330 = r23329 * r23325;
        float r23331 = r23327 * r23330;
        float r23332 = -inf.0;
        bool r23333 = r23331 <= r23332;
        float r23334 = r23328 * r23318;
        float r23335 = 1/2;
        float r23336 = r23334 * r23335;
        float r23337 = r23319 * r23328;
        float r23338 = r23318 / r23321;
        float r23339 = r23338 / r23319;
        float r23340 = r23339 / r23325;
        float r23341 = hypot(r23317, r23340);
        float r23342 = r23325 * r23341;
        float r23343 = r23337 * r23342;
        float r23344 = r23333 ? r23336 : r23343;
        return r23344;
}

double f_od(double J, double K, double U) {
        double r23345 = 1;
        double r23346 = U;
        double r23347 = J;
        double r23348 = r23346 / r23347;
        double r23349 = 2;
        double r23350 = r23348 / r23349;
        double r23351 = K;
        double r23352 = r23351 / r23349;
        double r23353 = cos(r23352);
        double r23354 = r23350 / r23353;
        double r23355 = hypot(r23345, r23354);
        double r23356 = -2;
        double r23357 = r23356 * r23347;
        double r23358 = r23357 * r23353;
        double r23359 = r23355 * r23358;
        double r23360 = -inf.0;
        bool r23361 = r23359 <= r23360;
        double r23362 = r23356 * r23346;
        double r23363 = 1/2;
        double r23364 = r23362 * r23363;
        double r23365 = r23347 * r23356;
        double r23366 = r23346 / r23349;
        double r23367 = r23366 / r23347;
        double r23368 = r23367 / r23353;
        double r23369 = hypot(r23345, r23368);
        double r23370 = r23353 * r23369;
        double r23371 = r23365 * r23370;
        double r23372 = r23361 ? r23364 : r23371;
        return r23372;
}

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 r23373, r23374, r23375, r23376, r23377, r23378, r23379, r23380, r23381, r23382, r23383, r23384, r23385, r23386, r23387, r23388, r23389;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r23373, "-2", 10, MPFR_RNDN);
        mpfr_init(r23374);
        mpfr_init(r23375);
        mpfr_init(r23376);
        mpfr_init_set_str(r23377, "2", 10, MPFR_RNDN);
        mpfr_init(r23378);
        mpfr_init(r23379);
        mpfr_init(r23380);
        mpfr_init_set_str(r23381, "1", 10, MPFR_RNDN);
        mpfr_init(r23382);
        mpfr_init(r23383);
        mpfr_init(r23384);
        mpfr_init(r23385);
        mpfr_init(r23386);
        mpfr_init(r23387);
        mpfr_init(r23388);
        mpfr_init(r23389);
}

double f_im(double J, double K, double U) {
        ;
        mpfr_set_d(r23374, J, MPFR_RNDN);
        mpfr_mul(r23375, r23373, r23374, MPFR_RNDN);
        mpfr_set_d(r23376, K, MPFR_RNDN);
        ;
        mpfr_div(r23378, r23376, r23377, MPFR_RNDN);
        mpfr_cos(r23379, r23378, MPFR_RNDN);
        mpfr_mul(r23380, r23375, r23379, MPFR_RNDN);
        ;
        mpfr_set_d(r23382, U, MPFR_RNDN);
        mpfr_mul(r23383, r23377, r23374, MPFR_RNDN);
        mpfr_mul(r23384, r23383, r23379, MPFR_RNDN);
        mpfr_div(r23385, r23382, r23384, MPFR_RNDN);
        mpfr_pow(r23386, r23385, r23377, MPFR_RNDN);
        mpfr_add(r23387, r23381, r23386, MPFR_RNDN);
        mpfr_sqrt(r23388, r23387, MPFR_RNDN);
        mpfr_mul(r23389, r23380, r23388, MPFR_RNDN);
        return mpfr_get_d(r23389, MPFR_RNDN);
}

static mpfr_t r23390, r23391, r23392, r23393, r23394, r23395, r23396, r23397, r23398, r23399, r23400, r23401, r23402, r23403, r23404, r23405, r23406, r23407, r23408, r23409, r23410, r23411, r23412, r23413, r23414, r23415, r23416, r23417;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r23390, "1", 10, MPFR_RNDN);
        mpfr_init(r23391);
        mpfr_init(r23392);
        mpfr_init(r23393);
        mpfr_init_set_str(r23394, "2", 10, MPFR_RNDN);
        mpfr_init(r23395);
        mpfr_init(r23396);
        mpfr_init(r23397);
        mpfr_init(r23398);
        mpfr_init(r23399);
        mpfr_init(r23400);
        mpfr_init_set_str(r23401, "-2", 10, MPFR_RNDN);
        mpfr_init(r23402);
        mpfr_init(r23403);
        mpfr_init(r23404);
        mpfr_init_set_str(r23405, "-inf.0", 10, MPFR_RNDN);
        mpfr_init(r23406);
        mpfr_init(r23407);
        mpfr_init_set_str(r23408, "1/2", 10, MPFR_RNDN);
        mpfr_init(r23409);
        mpfr_init(r23410);
        mpfr_init(r23411);
        mpfr_init(r23412);
        mpfr_init(r23413);
        mpfr_init(r23414);
        mpfr_init(r23415);
        mpfr_init(r23416);
        mpfr_init(r23417);
}

double f_fm(double J, double K, double U) {
        ;
        mpfr_set_d(r23391, U, MPFR_RNDN);
        mpfr_set_d(r23392, J, MPFR_RNDN);
        mpfr_div(r23393, r23391, r23392, MPFR_RNDN);
        ;
        mpfr_div(r23395, r23393, r23394, MPFR_RNDN);
        mpfr_set_d(r23396, K, MPFR_RNDN);
        mpfr_div(r23397, r23396, r23394, MPFR_RNDN);
        mpfr_cos(r23398, r23397, MPFR_RNDN);
        mpfr_div(r23399, r23395, r23398, MPFR_RNDN);
        mpfr_hypot(r23400, r23390, r23399, MPFR_RNDN);
        ;
        mpfr_mul(r23402, r23401, r23392, MPFR_RNDN);
        mpfr_mul(r23403, r23402, r23398, MPFR_RNDN);
        mpfr_mul(r23404, r23400, r23403, MPFR_RNDN);
        ;
        mpfr_set_si(r23406, mpfr_cmp(r23404, r23405) <= 0, MPFR_RNDN);
        mpfr_mul(r23407, r23401, r23391, MPFR_RNDN);
        ;
        mpfr_mul(r23409, r23407, r23408, MPFR_RNDN);
        mpfr_mul(r23410, r23392, r23401, MPFR_RNDN);
        mpfr_div(r23411, r23391, r23394, MPFR_RNDN);
        mpfr_div(r23412, r23411, r23392, MPFR_RNDN);
        mpfr_div(r23413, r23412, r23398, MPFR_RNDN);
        mpfr_hypot(r23414, r23390, r23413, MPFR_RNDN);
        mpfr_mul(r23415, r23398, r23414, MPFR_RNDN);
        mpfr_mul(r23416, r23410, r23415, MPFR_RNDN);
        if (mpfr_get_si(r23406, MPFR_RNDN)) { mpfr_set(r23417, r23409, MPFR_RNDN); } else { mpfr_set(r23417, r23416, MPFR_RNDN); };
        return mpfr_get_d(r23417, MPFR_RNDN);
}

static mpfr_t r23418, r23419, r23420, r23421, r23422, r23423, r23424, r23425, r23426, r23427, r23428, r23429, r23430, r23431, r23432, r23433, r23434, r23435, r23436, r23437, r23438, r23439, r23440, r23441, r23442, r23443, r23444, r23445;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r23418, "1", 10, MPFR_RNDN);
        mpfr_init(r23419);
        mpfr_init(r23420);
        mpfr_init(r23421);
        mpfr_init_set_str(r23422, "2", 10, MPFR_RNDN);
        mpfr_init(r23423);
        mpfr_init(r23424);
        mpfr_init(r23425);
        mpfr_init(r23426);
        mpfr_init(r23427);
        mpfr_init(r23428);
        mpfr_init_set_str(r23429, "-2", 10, MPFR_RNDN);
        mpfr_init(r23430);
        mpfr_init(r23431);
        mpfr_init(r23432);
        mpfr_init_set_str(r23433, "-inf.0", 10, MPFR_RNDN);
        mpfr_init(r23434);
        mpfr_init(r23435);
        mpfr_init_set_str(r23436, "1/2", 10, MPFR_RNDN);
        mpfr_init(r23437);
        mpfr_init(r23438);
        mpfr_init(r23439);
        mpfr_init(r23440);
        mpfr_init(r23441);
        mpfr_init(r23442);
        mpfr_init(r23443);
        mpfr_init(r23444);
        mpfr_init(r23445);
}

double f_dm(double J, double K, double U) {
        ;
        mpfr_set_d(r23419, U, MPFR_RNDN);
        mpfr_set_d(r23420, J, MPFR_RNDN);
        mpfr_div(r23421, r23419, r23420, MPFR_RNDN);
        ;
        mpfr_div(r23423, r23421, r23422, MPFR_RNDN);
        mpfr_set_d(r23424, K, MPFR_RNDN);
        mpfr_div(r23425, r23424, r23422, MPFR_RNDN);
        mpfr_cos(r23426, r23425, MPFR_RNDN);
        mpfr_div(r23427, r23423, r23426, MPFR_RNDN);
        mpfr_hypot(r23428, r23418, r23427, MPFR_RNDN);
        ;
        mpfr_mul(r23430, r23429, r23420, MPFR_RNDN);
        mpfr_mul(r23431, r23430, r23426, MPFR_RNDN);
        mpfr_mul(r23432, r23428, r23431, MPFR_RNDN);
        ;
        mpfr_set_si(r23434, mpfr_cmp(r23432, r23433) <= 0, MPFR_RNDN);
        mpfr_mul(r23435, r23429, r23419, MPFR_RNDN);
        ;
        mpfr_mul(r23437, r23435, r23436, MPFR_RNDN);
        mpfr_mul(r23438, r23420, r23429, MPFR_RNDN);
        mpfr_div(r23439, r23419, r23422, MPFR_RNDN);
        mpfr_div(r23440, r23439, r23420, MPFR_RNDN);
        mpfr_div(r23441, r23440, r23426, MPFR_RNDN);
        mpfr_hypot(r23442, r23418, r23441, MPFR_RNDN);
        mpfr_mul(r23443, r23426, r23442, MPFR_RNDN);
        mpfr_mul(r23444, r23438, r23443, MPFR_RNDN);
        if (mpfr_get_si(r23434, MPFR_RNDN)) { mpfr_set(r23445, r23437, MPFR_RNDN); } else { mpfr_set(r23445, r23444, MPFR_RNDN); };
        return mpfr_get_d(r23445, MPFR_RNDN);
}