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

char *name = "Harley's example";

double f_if(float c_p, float c_n, float t, float s) {
        float r27376 = 1;
        float r27377 = s;
        float r27378 = -r27377;
        float r27379 = exp(r27378);
        float r27380 = r27376 + r27379;
        float r27381 = r27376 / r27380;
        float r27382 = c_p;
        float r27383 = pow(r27381, r27382);
        float r27384 = r27376 - r27381;
        float r27385 = c_n;
        float r27386 = pow(r27384, r27385);
        float r27387 = r27383 * r27386;
        float r27388 = t;
        float r27389 = -r27388;
        float r27390 = exp(r27389);
        float r27391 = r27376 + r27390;
        float r27392 = r27376 / r27391;
        float r27393 = pow(r27392, r27382);
        float r27394 = r27376 - r27392;
        float r27395 = pow(r27394, r27385);
        float r27396 = r27393 * r27395;
        float r27397 = r27387 / r27396;
        return r27397;
}

double f_id(double c_p, double c_n, double t, double s) {
        double r27398 = 1;
        double r27399 = s;
        double r27400 = -r27399;
        double r27401 = exp(r27400);
        double r27402 = r27398 + r27401;
        double r27403 = r27398 / r27402;
        double r27404 = c_p;
        double r27405 = pow(r27403, r27404);
        double r27406 = r27398 - r27403;
        double r27407 = c_n;
        double r27408 = pow(r27406, r27407);
        double r27409 = r27405 * r27408;
        double r27410 = t;
        double r27411 = -r27410;
        double r27412 = exp(r27411);
        double r27413 = r27398 + r27412;
        double r27414 = r27398 / r27413;
        double r27415 = pow(r27414, r27404);
        double r27416 = r27398 - r27414;
        double r27417 = pow(r27416, r27407);
        double r27418 = r27415 * r27417;
        double r27419 = r27409 / r27418;
        return r27419;
}


double f_of(float c_p, float c_n, float t, float s) {
        float r27420 = 1;
        float r27421 = s;
        float r27422 = -r27421;
        float r27423 = exp(r27422);
        float r27424 = r27420 + r27423;
        float r27425 = r27420 / r27424;
        float r27426 = 0.49183472926145494;
        bool r27427 = r27425 <= r27426;
        float r27428 = 0.5000000000000013;
        bool r27429 = r27425 <= r27428;
        float r27430 = !r27429;
        bool r27431 = r27427 || r27430;
        float r27432 = c_p;
        float r27433 = pow(r27425, r27432);
        float r27434 = t;
        float r27435 = 1/2;
        float r27436 = r27434 * r27435;
        float r27437 = log(r27435);
        float r27438 = r27436 + r27437;
        float r27439 = r27432 * r27438;
        float r27440 = r27439 + r27420;
        float r27441 = r27433 / r27440;
        float r27442 = r27420 - r27425;
        float r27443 = c_n;
        float r27444 = pow(r27442, r27443);
        float r27445 = -r27434;
        float r27446 = exp(r27445);
        float r27447 = r27420 + r27446;
        float r27448 = r27420 / r27447;
        float r27449 = r27420 - r27448;
        float r27450 = pow(r27449, r27443);
        float r27451 = r27444 / r27450;
        float r27452 = r27441 * r27451;
        float r27453 = r27432 - r27443;
        float r27454 = r27435 * r27421;
        float r27455 = r27453 * r27454;
        float r27456 = r27420 + r27455;
        float r27457 = r27431 ? r27452 : r27456;
        return r27457;
}

double f_od(double c_p, double c_n, double t, double s) {
        double r27458 = 1;
        double r27459 = s;
        double r27460 = -r27459;
        double r27461 = exp(r27460);
        double r27462 = r27458 + r27461;
        double r27463 = r27458 / r27462;
        double r27464 = 0.49183472926145494;
        bool r27465 = r27463 <= r27464;
        double r27466 = 0.5000000000000013;
        bool r27467 = r27463 <= r27466;
        double r27468 = !r27467;
        bool r27469 = r27465 || r27468;
        double r27470 = c_p;
        double r27471 = pow(r27463, r27470);
        double r27472 = t;
        double r27473 = 1/2;
        double r27474 = r27472 * r27473;
        double r27475 = log(r27473);
        double r27476 = r27474 + r27475;
        double r27477 = r27470 * r27476;
        double r27478 = r27477 + r27458;
        double r27479 = r27471 / r27478;
        double r27480 = r27458 - r27463;
        double r27481 = c_n;
        double r27482 = pow(r27480, r27481);
        double r27483 = -r27472;
        double r27484 = exp(r27483);
        double r27485 = r27458 + r27484;
        double r27486 = r27458 / r27485;
        double r27487 = r27458 - r27486;
        double r27488 = pow(r27487, r27481);
        double r27489 = r27482 / r27488;
        double r27490 = r27479 * r27489;
        double r27491 = r27470 - r27481;
        double r27492 = r27473 * r27459;
        double r27493 = r27491 * r27492;
        double r27494 = r27458 + r27493;
        double r27495 = r27469 ? r27490 : r27494;
        return r27495;
}

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 r27496, r27497, r27498, r27499, r27500, r27501, r27502, r27503, r27504, r27505, r27506, r27507, r27508, r27509, r27510, r27511, r27512, r27513, r27514, r27515, r27516, r27517;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r27496, "1", 10, MPFR_RNDN);
        mpfr_init(r27497);
        mpfr_init(r27498);
        mpfr_init(r27499);
        mpfr_init(r27500);
        mpfr_init(r27501);
        mpfr_init(r27502);
        mpfr_init(r27503);
        mpfr_init(r27504);
        mpfr_init(r27505);
        mpfr_init(r27506);
        mpfr_init(r27507);
        mpfr_init(r27508);
        mpfr_init(r27509);
        mpfr_init(r27510);
        mpfr_init(r27511);
        mpfr_init(r27512);
        mpfr_init(r27513);
        mpfr_init(r27514);
        mpfr_init(r27515);
        mpfr_init(r27516);
        mpfr_init(r27517);
}

double f_im(double c_p, double c_n, double t, double s) {
        ;
        mpfr_set_d(r27497, s, MPFR_RNDN);
        mpfr_neg(r27498, r27497, MPFR_RNDN);
        mpfr_exp(r27499, r27498, MPFR_RNDN);
        mpfr_add(r27500, r27496, r27499, MPFR_RNDN);
        mpfr_div(r27501, r27496, r27500, MPFR_RNDN);
        mpfr_set_d(r27502, c_p, MPFR_RNDN);
        mpfr_pow(r27503, r27501, r27502, MPFR_RNDN);
        mpfr_sub(r27504, r27496, r27501, MPFR_RNDN);
        mpfr_set_d(r27505, c_n, MPFR_RNDN);
        mpfr_pow(r27506, r27504, r27505, MPFR_RNDN);
        mpfr_mul(r27507, r27503, r27506, MPFR_RNDN);
        mpfr_set_d(r27508, t, MPFR_RNDN);
        mpfr_neg(r27509, r27508, MPFR_RNDN);
        mpfr_exp(r27510, r27509, MPFR_RNDN);
        mpfr_add(r27511, r27496, r27510, MPFR_RNDN);
        mpfr_div(r27512, r27496, r27511, MPFR_RNDN);
        mpfr_pow(r27513, r27512, r27502, MPFR_RNDN);
        mpfr_sub(r27514, r27496, r27512, MPFR_RNDN);
        mpfr_pow(r27515, r27514, r27505, MPFR_RNDN);
        mpfr_mul(r27516, r27513, r27515, MPFR_RNDN);
        mpfr_div(r27517, r27507, r27516, MPFR_RNDN);
        return mpfr_get_d(r27517, MPFR_RNDN);
}

static mpfr_t r27518, r27519, r27520, r27521, r27522, r27523, r27524, r27525, r27526, r27527, r27528, r27529, r27530, r27531, r27532, r27533, r27534, r27535, r27536, r27537, r27538, r27539, r27540, r27541, r27542, r27543, r27544, r27545, r27546, r27547, r27548, r27549, r27550, r27551, r27552, r27553, r27554, r27555;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r27518, "1", 10, MPFR_RNDN);
        mpfr_init(r27519);
        mpfr_init(r27520);
        mpfr_init(r27521);
        mpfr_init(r27522);
        mpfr_init(r27523);
        mpfr_init_set_str(r27524, "0.49183472926145494", 10, MPFR_RNDN);
        mpfr_init(r27525);
        mpfr_init_set_str(r27526, "0.5000000000000013", 10, MPFR_RNDN);
        mpfr_init(r27527);
        mpfr_init(r27528);
        mpfr_init(r27529);
        mpfr_init(r27530);
        mpfr_init(r27531);
        mpfr_init(r27532);
        mpfr_init_set_str(r27533, "1/2", 10, MPFR_RNDN);
        mpfr_init(r27534);
        mpfr_init(r27535);
        mpfr_init(r27536);
        mpfr_init(r27537);
        mpfr_init(r27538);
        mpfr_init(r27539);
        mpfr_init(r27540);
        mpfr_init(r27541);
        mpfr_init(r27542);
        mpfr_init(r27543);
        mpfr_init(r27544);
        mpfr_init(r27545);
        mpfr_init(r27546);
        mpfr_init(r27547);
        mpfr_init(r27548);
        mpfr_init(r27549);
        mpfr_init(r27550);
        mpfr_init(r27551);
        mpfr_init(r27552);
        mpfr_init(r27553);
        mpfr_init(r27554);
        mpfr_init(r27555);
}

double f_fm(double c_p, double c_n, double t, double s) {
        ;
        mpfr_set_d(r27519, s, MPFR_RNDN);
        mpfr_neg(r27520, r27519, MPFR_RNDN);
        mpfr_exp(r27521, r27520, MPFR_RNDN);
        mpfr_add(r27522, r27518, r27521, MPFR_RNDN);
        mpfr_div(r27523, r27518, r27522, MPFR_RNDN);
        ;
        mpfr_set_si(r27525, mpfr_cmp(r27523, r27524) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r27527, mpfr_cmp(r27523, r27526) <= 0, MPFR_RNDN);
        mpfr_set_si(r27528, !mpfr_get_si(r27527, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r27529, mpfr_get_si(r27525, MPFR_RNDN) || mpfr_get_si(r27528, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_d(r27530, c_p, MPFR_RNDN);
        mpfr_pow(r27531, r27523, r27530, MPFR_RNDN);
        mpfr_set_d(r27532, t, MPFR_RNDN);
        ;
        mpfr_mul(r27534, r27532, r27533, MPFR_RNDN);
        mpfr_log(r27535, r27533, MPFR_RNDN);
        mpfr_add(r27536, r27534, r27535, MPFR_RNDN);
        mpfr_mul(r27537, r27530, r27536, MPFR_RNDN);
        mpfr_add(r27538, r27537, r27518, MPFR_RNDN);
        mpfr_div(r27539, r27531, r27538, MPFR_RNDN);
        mpfr_sub(r27540, r27518, r27523, MPFR_RNDN);
        mpfr_set_d(r27541, c_n, MPFR_RNDN);
        mpfr_pow(r27542, r27540, r27541, MPFR_RNDN);
        mpfr_neg(r27543, r27532, MPFR_RNDN);
        mpfr_exp(r27544, r27543, MPFR_RNDN);
        mpfr_add(r27545, r27518, r27544, MPFR_RNDN);
        mpfr_div(r27546, r27518, r27545, MPFR_RNDN);
        mpfr_sub(r27547, r27518, r27546, MPFR_RNDN);
        mpfr_pow(r27548, r27547, r27541, MPFR_RNDN);
        mpfr_div(r27549, r27542, r27548, MPFR_RNDN);
        mpfr_mul(r27550, r27539, r27549, MPFR_RNDN);
        mpfr_sub(r27551, r27530, r27541, MPFR_RNDN);
        mpfr_mul(r27552, r27533, r27519, MPFR_RNDN);
        mpfr_mul(r27553, r27551, r27552, MPFR_RNDN);
        mpfr_add(r27554, r27518, r27553, MPFR_RNDN);
        if (mpfr_get_si(r27529, MPFR_RNDN)) { mpfr_set(r27555, r27550, MPFR_RNDN); } else { mpfr_set(r27555, r27554, MPFR_RNDN); };
        return mpfr_get_d(r27555, MPFR_RNDN);
}

static mpfr_t r27556, r27557, r27558, r27559, r27560, r27561, r27562, r27563, r27564, r27565, r27566, r27567, r27568, r27569, r27570, r27571, r27572, r27573, r27574, r27575, r27576, r27577, r27578, r27579, r27580, r27581, r27582, r27583, r27584, r27585, r27586, r27587, r27588, r27589, r27590, r27591, r27592, r27593;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r27556, "1", 10, MPFR_RNDN);
        mpfr_init(r27557);
        mpfr_init(r27558);
        mpfr_init(r27559);
        mpfr_init(r27560);
        mpfr_init(r27561);
        mpfr_init_set_str(r27562, "0.49183472926145494", 10, MPFR_RNDN);
        mpfr_init(r27563);
        mpfr_init_set_str(r27564, "0.5000000000000013", 10, MPFR_RNDN);
        mpfr_init(r27565);
        mpfr_init(r27566);
        mpfr_init(r27567);
        mpfr_init(r27568);
        mpfr_init(r27569);
        mpfr_init(r27570);
        mpfr_init_set_str(r27571, "1/2", 10, MPFR_RNDN);
        mpfr_init(r27572);
        mpfr_init(r27573);
        mpfr_init(r27574);
        mpfr_init(r27575);
        mpfr_init(r27576);
        mpfr_init(r27577);
        mpfr_init(r27578);
        mpfr_init(r27579);
        mpfr_init(r27580);
        mpfr_init(r27581);
        mpfr_init(r27582);
        mpfr_init(r27583);
        mpfr_init(r27584);
        mpfr_init(r27585);
        mpfr_init(r27586);
        mpfr_init(r27587);
        mpfr_init(r27588);
        mpfr_init(r27589);
        mpfr_init(r27590);
        mpfr_init(r27591);
        mpfr_init(r27592);
        mpfr_init(r27593);
}

double f_dm(double c_p, double c_n, double t, double s) {
        ;
        mpfr_set_d(r27557, s, MPFR_RNDN);
        mpfr_neg(r27558, r27557, MPFR_RNDN);
        mpfr_exp(r27559, r27558, MPFR_RNDN);
        mpfr_add(r27560, r27556, r27559, MPFR_RNDN);
        mpfr_div(r27561, r27556, r27560, MPFR_RNDN);
        ;
        mpfr_set_si(r27563, mpfr_cmp(r27561, r27562) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r27565, mpfr_cmp(r27561, r27564) <= 0, MPFR_RNDN);
        mpfr_set_si(r27566, !mpfr_get_si(r27565, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r27567, mpfr_get_si(r27563, MPFR_RNDN) || mpfr_get_si(r27566, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_d(r27568, c_p, MPFR_RNDN);
        mpfr_pow(r27569, r27561, r27568, MPFR_RNDN);
        mpfr_set_d(r27570, t, MPFR_RNDN);
        ;
        mpfr_mul(r27572, r27570, r27571, MPFR_RNDN);
        mpfr_log(r27573, r27571, MPFR_RNDN);
        mpfr_add(r27574, r27572, r27573, MPFR_RNDN);
        mpfr_mul(r27575, r27568, r27574, MPFR_RNDN);
        mpfr_add(r27576, r27575, r27556, MPFR_RNDN);
        mpfr_div(r27577, r27569, r27576, MPFR_RNDN);
        mpfr_sub(r27578, r27556, r27561, MPFR_RNDN);
        mpfr_set_d(r27579, c_n, MPFR_RNDN);
        mpfr_pow(r27580, r27578, r27579, MPFR_RNDN);
        mpfr_neg(r27581, r27570, MPFR_RNDN);
        mpfr_exp(r27582, r27581, MPFR_RNDN);
        mpfr_add(r27583, r27556, r27582, MPFR_RNDN);
        mpfr_div(r27584, r27556, r27583, MPFR_RNDN);
        mpfr_sub(r27585, r27556, r27584, MPFR_RNDN);
        mpfr_pow(r27586, r27585, r27579, MPFR_RNDN);
        mpfr_div(r27587, r27580, r27586, MPFR_RNDN);
        mpfr_mul(r27588, r27577, r27587, MPFR_RNDN);
        mpfr_sub(r27589, r27568, r27579, MPFR_RNDN);
        mpfr_mul(r27590, r27571, r27557, MPFR_RNDN);
        mpfr_mul(r27591, r27589, r27590, MPFR_RNDN);
        mpfr_add(r27592, r27556, r27591, MPFR_RNDN);
        if (mpfr_get_si(r27567, MPFR_RNDN)) { mpfr_set(r27593, r27588, MPFR_RNDN); } else { mpfr_set(r27593, r27592, MPFR_RNDN); };
        return mpfr_get_d(r27593, MPFR_RNDN);
}