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

char *name = "From Warwick Tucker's Validated Numerics";

double f_if() {
        float r8361 = 333.75f;
        float r8362 = 33096.0f;
        float r8363 = 6.0f;
        float r8364 = pow(r8362, r8363);
        float r8365 = r8361 * r8364;
        float r8366 = 77617.0f;
        float r8367 = r8366 * r8366;
        float r8368 = 11.0f;
        float r8369 = r8368 * r8367;
        float r8370 = r8362 * r8362;
        float r8371 = r8369 * r8370;
        float r8372 = -r8364;
        float r8373 = r8371 + r8372;
        float r8374 = -121.0f;
        float r8375 = 4.0f;
        float r8376 = pow(r8362, r8375);
        float r8377 = r8374 * r8376;
        float r8378 = r8373 + r8377;
        float r8379 = -2.0f;
        float r8380 = r8378 + r8379;
        float r8381 = r8367 * r8380;
        float r8382 = r8365 + r8381;
        float r8383 = 5.5f;
        float r8384 = 8.0f;
        float r8385 = pow(r8362, r8384);
        float r8386 = r8383 * r8385;
        float r8387 = r8382 + r8386;
        float r8388 = 2.0f;
        float r8389 = r8388 * r8362;
        float r8390 = r8366 / r8389;
        float r8391 = r8387 + r8390;
        return r8391;
}

double f_id() {
        double r8392 = 333.75;
        double r8393 = 33096.0;
        double r8394 = 6.0;
        double r8395 = pow(r8393, r8394);
        double r8396 = r8392 * r8395;
        double r8397 = 77617.0;
        double r8398 = r8397 * r8397;
        double r8399 = 11.0;
        double r8400 = r8399 * r8398;
        double r8401 = r8393 * r8393;
        double r8402 = r8400 * r8401;
        double r8403 = -r8395;
        double r8404 = r8402 + r8403;
        double r8405 = -121.0;
        double r8406 = 4.0;
        double r8407 = pow(r8393, r8406);
        double r8408 = r8405 * r8407;
        double r8409 = r8404 + r8408;
        double r8410 = -2.0;
        double r8411 = r8409 + r8410;
        double r8412 = r8398 * r8411;
        double r8413 = r8396 + r8412;
        double r8414 = 5.5;
        double r8415 = 8.0;
        double r8416 = pow(r8393, r8415);
        double r8417 = r8414 * r8416;
        double r8418 = r8413 + r8417;
        double r8419 = 2.0;
        double r8420 = r8419 * r8393;
        double r8421 = r8397 / r8420;
        double r8422 = r8418 + r8421;
        return r8422;
}


double f_of() {
        float r8423 = 77617.0f;
        float r8424 = r8423 * r8423;
        float r8425 = 33096.0f;
        float r8426 = r8425 * r8425;
        float r8427 = r8424 * r8426;
        float r8428 = 11.0f;
        float r8429 = r8427 * r8428;
        float r8430 = 6.0f;
        float r8431 = pow(r8425, r8430);
        float r8432 = -121.0f;
        float r8433 = 4.0f;
        float r8434 = pow(r8425, r8433);
        float r8435 = -2.0f;
        float r8436 = fma(r8432, r8434, r8435);
        float r8437 = r8431 - r8436;
        float r8438 = r8429 - r8437;
        float r8439 = 333.75f;
        float r8440 = 8.0f;
        float r8441 = pow(r8425, r8440);
        float r8442 = 5.5f;
        float r8443 = 2.0f;
        float r8444 = r8423 / r8443;
        float r8445 = r8444 / r8425;
        float r8446 = fma(r8441, r8442, r8445);
        float r8447 = fma(r8431, r8439, r8446);
        float r8448 = fma(r8438, r8424, r8447);
        return r8448;
}

double f_od() {
        double r8449 = 77617.0;
        double r8450 = r8449 * r8449;
        double r8451 = 33096.0;
        double r8452 = r8451 * r8451;
        double r8453 = r8450 * r8452;
        double r8454 = 11.0;
        double r8455 = r8453 * r8454;
        double r8456 = 6.0;
        double r8457 = pow(r8451, r8456);
        double r8458 = -121.0;
        double r8459 = 4.0;
        double r8460 = pow(r8451, r8459);
        double r8461 = -2.0;
        double r8462 = fma(r8458, r8460, r8461);
        double r8463 = r8457 - r8462;
        double r8464 = r8455 - r8463;
        double r8465 = 333.75;
        double r8466 = 8.0;
        double r8467 = pow(r8451, r8466);
        double r8468 = 5.5;
        double r8469 = 2.0;
        double r8470 = r8449 / r8469;
        double r8471 = r8470 / r8451;
        double r8472 = fma(r8467, r8468, r8471);
        double r8473 = fma(r8457, r8465, r8472);
        double r8474 = fma(r8464, r8450, r8473);
        return r8474;
}

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 r8475, r8476, r8477, r8478, r8479, r8480, r8481, r8482, r8483, r8484, r8485, r8486, r8487, r8488, r8489, r8490, r8491, r8492, r8493, r8494, r8495, r8496, r8497, r8498, r8499, r8500, r8501, r8502, r8503, r8504, r8505;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8475, "333.75", 10, MPFR_RNDN);
        mpfr_init_set_str(r8476, "33096", 10, MPFR_RNDN);
        mpfr_init_set_str(r8477, "6", 10, MPFR_RNDN);
        mpfr_init(r8478);
        mpfr_init(r8479);
        mpfr_init_set_str(r8480, "77617", 10, MPFR_RNDN);
        mpfr_init(r8481);
        mpfr_init_set_str(r8482, "11", 10, MPFR_RNDN);
        mpfr_init(r8483);
        mpfr_init(r8484);
        mpfr_init(r8485);
        mpfr_init(r8486);
        mpfr_init(r8487);
        mpfr_init_set_str(r8488, "-121", 10, MPFR_RNDN);
        mpfr_init_set_str(r8489, "4", 10, MPFR_RNDN);
        mpfr_init(r8490);
        mpfr_init(r8491);
        mpfr_init(r8492);
        mpfr_init_set_str(r8493, "-2", 10, MPFR_RNDN);
        mpfr_init(r8494);
        mpfr_init(r8495);
        mpfr_init(r8496);
        mpfr_init_set_str(r8497, "5.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r8498, "8", 10, MPFR_RNDN);
        mpfr_init(r8499);
        mpfr_init(r8500);
        mpfr_init(r8501);
        mpfr_init_set_str(r8502, "2", 10, MPFR_RNDN);
        mpfr_init(r8503);
        mpfr_init(r8504);
        mpfr_init(r8505);
}

double f_im() {
        ;
        ;
        ;
        mpfr_pow(r8478, r8476, r8477, MPFR_RNDN);
        mpfr_mul(r8479, r8475, r8478, MPFR_RNDN);
        ;
        mpfr_mul(r8481, r8480, r8480, MPFR_RNDN);
        ;
        mpfr_mul(r8483, r8482, r8481, MPFR_RNDN);
        mpfr_mul(r8484, r8476, r8476, MPFR_RNDN);
        mpfr_mul(r8485, r8483, r8484, MPFR_RNDN);
        mpfr_neg(r8486, r8478, MPFR_RNDN);
        mpfr_add(r8487, r8485, r8486, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r8490, r8476, r8489, MPFR_RNDN);
        mpfr_mul(r8491, r8488, r8490, MPFR_RNDN);
        mpfr_add(r8492, r8487, r8491, MPFR_RNDN);
        ;
        mpfr_add(r8494, r8492, r8493, MPFR_RNDN);
        mpfr_mul(r8495, r8481, r8494, MPFR_RNDN);
        mpfr_add(r8496, r8479, r8495, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r8499, r8476, r8498, MPFR_RNDN);
        mpfr_mul(r8500, r8497, r8499, MPFR_RNDN);
        mpfr_add(r8501, r8496, r8500, MPFR_RNDN);
        ;
        mpfr_mul(r8503, r8502, r8476, MPFR_RNDN);
        mpfr_div(r8504, r8480, r8503, MPFR_RNDN);
        mpfr_add(r8505, r8501, r8504, MPFR_RNDN);
        return mpfr_get_d(r8505, MPFR_RNDN);
}

static mpfr_t r8506, r8507, r8508, r8509, r8510, r8511, r8512, r8513, r8514, r8515, r8516, r8517, r8518, r8519, r8520, r8521, r8522, r8523, r8524, r8525, r8526, r8527, r8528, r8529, r8530, r8531;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8506, "77617", 10, MPFR_RNDN);
        mpfr_init(r8507);
        mpfr_init_set_str(r8508, "33096", 10, MPFR_RNDN);
        mpfr_init(r8509);
        mpfr_init(r8510);
        mpfr_init_set_str(r8511, "11", 10, MPFR_RNDN);
        mpfr_init(r8512);
        mpfr_init_set_str(r8513, "6", 10, MPFR_RNDN);
        mpfr_init(r8514);
        mpfr_init_set_str(r8515, "-121", 10, MPFR_RNDN);
        mpfr_init_set_str(r8516, "4", 10, MPFR_RNDN);
        mpfr_init(r8517);
        mpfr_init_set_str(r8518, "-2", 10, MPFR_RNDN);
        mpfr_init(r8519);
        mpfr_init(r8520);
        mpfr_init(r8521);
        mpfr_init_set_str(r8522, "333.75", 10, MPFR_RNDN);
        mpfr_init_set_str(r8523, "8", 10, MPFR_RNDN);
        mpfr_init(r8524);
        mpfr_init_set_str(r8525, "5.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r8526, "2", 10, MPFR_RNDN);
        mpfr_init(r8527);
        mpfr_init(r8528);
        mpfr_init(r8529);
        mpfr_init(r8530);
        mpfr_init(r8531);
}

double f_fm() {
        ;
        mpfr_mul(r8507, r8506, r8506, MPFR_RNDN);
        ;
        mpfr_mul(r8509, r8508, r8508, MPFR_RNDN);
        mpfr_mul(r8510, r8507, r8509, MPFR_RNDN);
        ;
        mpfr_mul(r8512, r8510, r8511, MPFR_RNDN);
        ;
        mpfr_pow(r8514, r8508, r8513, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r8517, r8508, r8516, MPFR_RNDN);
        ;
        mpfr_fma(r8519, r8515, r8517, r8518, MPFR_RNDN);
        mpfr_sub(r8520, r8514, r8519, MPFR_RNDN);
        mpfr_sub(r8521, r8512, r8520, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r8524, r8508, r8523, MPFR_RNDN);
        ;
        ;
        mpfr_div(r8527, r8506, r8526, MPFR_RNDN);
        mpfr_div(r8528, r8527, r8508, MPFR_RNDN);
        mpfr_fma(r8529, r8524, r8525, r8528, MPFR_RNDN);
        mpfr_fma(r8530, r8514, r8522, r8529, MPFR_RNDN);
        mpfr_fma(r8531, r8521, r8507, r8530, MPFR_RNDN);
        return mpfr_get_d(r8531, MPFR_RNDN);
}

static mpfr_t r8532, r8533, r8534, r8535, r8536, r8537, r8538, r8539, r8540, r8541, r8542, r8543, r8544, r8545, r8546, r8547, r8548, r8549, r8550, r8551, r8552, r8553, r8554, r8555, r8556, r8557;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8532, "77617", 10, MPFR_RNDN);
        mpfr_init(r8533);
        mpfr_init_set_str(r8534, "33096", 10, MPFR_RNDN);
        mpfr_init(r8535);
        mpfr_init(r8536);
        mpfr_init_set_str(r8537, "11", 10, MPFR_RNDN);
        mpfr_init(r8538);
        mpfr_init_set_str(r8539, "6", 10, MPFR_RNDN);
        mpfr_init(r8540);
        mpfr_init_set_str(r8541, "-121", 10, MPFR_RNDN);
        mpfr_init_set_str(r8542, "4", 10, MPFR_RNDN);
        mpfr_init(r8543);
        mpfr_init_set_str(r8544, "-2", 10, MPFR_RNDN);
        mpfr_init(r8545);
        mpfr_init(r8546);
        mpfr_init(r8547);
        mpfr_init_set_str(r8548, "333.75", 10, MPFR_RNDN);
        mpfr_init_set_str(r8549, "8", 10, MPFR_RNDN);
        mpfr_init(r8550);
        mpfr_init_set_str(r8551, "5.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r8552, "2", 10, MPFR_RNDN);
        mpfr_init(r8553);
        mpfr_init(r8554);
        mpfr_init(r8555);
        mpfr_init(r8556);
        mpfr_init(r8557);
}

double f_dm() {
        ;
        mpfr_mul(r8533, r8532, r8532, MPFR_RNDN);
        ;
        mpfr_mul(r8535, r8534, r8534, MPFR_RNDN);
        mpfr_mul(r8536, r8533, r8535, MPFR_RNDN);
        ;
        mpfr_mul(r8538, r8536, r8537, MPFR_RNDN);
        ;
        mpfr_pow(r8540, r8534, r8539, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r8543, r8534, r8542, MPFR_RNDN);
        ;
        mpfr_fma(r8545, r8541, r8543, r8544, MPFR_RNDN);
        mpfr_sub(r8546, r8540, r8545, MPFR_RNDN);
        mpfr_sub(r8547, r8538, r8546, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r8550, r8534, r8549, MPFR_RNDN);
        ;
        ;
        mpfr_div(r8553, r8532, r8552, MPFR_RNDN);
        mpfr_div(r8554, r8553, r8534, MPFR_RNDN);
        mpfr_fma(r8555, r8550, r8551, r8554, MPFR_RNDN);
        mpfr_fma(r8556, r8540, r8548, r8555, MPFR_RNDN);
        mpfr_fma(r8557, r8547, r8533, r8556, MPFR_RNDN);
        return mpfr_get_d(r8557, MPFR_RNDN);
}