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

char *name = "3.9.2 real part (p56)";

double f_if() {
        float r8384 = 1;
        float r8385 = -r8384;
        float r8386 = 2;
        float r8387 = r8385 / r8386;
        float r8388 = 3;
        float r8389 = sqrt(r8388);
        float r8390 = r8389 / r8386;
        float r8391 = /* ERROR: no complex support in C */;
        float r8392 = r8391 * r8391;
        float r8393 = r8392 * r8391;
        float r8394 = r8393 * r8391;
        float r8395 = -r8386;
        float r8396 = 0;
        float r8397 = /* ERROR: no complex support in C */;
        float r8398 = r8397 * r8391;
        float r8399 = r8398 * r8391;
        float r8400 = r8399 * r8391;
        float r8401 = r8394 + r8400;
        float r8402 = 5;
        float r8403 = /* ERROR: no complex support in C */;
        float r8404 = r8403 * r8391;
        float r8405 = r8404 * r8391;
        float r8406 = r8401 + r8405;
        float r8407 = 4;
        float r8408 = /* ERROR: no complex support in C */;
        float r8409 = r8408 * r8391;
        float r8410 = r8406 + r8409;
        float r8411 = 7;
        float r8412 = /* ERROR: no complex support in C */;
        float r8413 = r8410 + r8412;
        float r8414 = /* ERROR: no complex support in C */;
        return r8414;
}

double f_id() {
        double r8415 = 1;
        double r8416 = -r8415;
        double r8417 = 2;
        double r8418 = r8416 / r8417;
        double r8419 = 3;
        double r8420 = sqrt(r8419);
        double r8421 = r8420 / r8417;
        double r8422 = /* ERROR: no complex support in C */;
        double r8423 = r8422 * r8422;
        double r8424 = r8423 * r8422;
        double r8425 = r8424 * r8422;
        double r8426 = -r8417;
        double r8427 = 0;
        double r8428 = /* ERROR: no complex support in C */;
        double r8429 = r8428 * r8422;
        double r8430 = r8429 * r8422;
        double r8431 = r8430 * r8422;
        double r8432 = r8425 + r8431;
        double r8433 = 5;
        double r8434 = /* ERROR: no complex support in C */;
        double r8435 = r8434 * r8422;
        double r8436 = r8435 * r8422;
        double r8437 = r8432 + r8436;
        double r8438 = 4;
        double r8439 = /* ERROR: no complex support in C */;
        double r8440 = r8439 * r8422;
        double r8441 = r8437 + r8440;
        double r8442 = 7;
        double r8443 = /* ERROR: no complex support in C */;
        double r8444 = r8441 + r8443;
        double r8445 = /* ERROR: no complex support in C */;
        return r8445;
}


double f_of() {
        float r8446 = 5;
        float r8447 = -r8446;
        float r8448 = 1;
        float r8449 = 2;
        float r8450 = r8448 / r8449;
        float r8451 = r8447 * r8450;
        float r8452 = 4;
        float r8453 = r8451 + r8452;
        float r8454 = 3;
        float r8455 = sqrt(r8454);
        float r8456 = r8455 / r8449;
        float r8457 = r8446 * r8456;
        float r8458 = /* ERROR: no complex support in C */;
        float r8459 = -r8448;
        float r8460 = r8459 / r8449;
        float r8461 = /* ERROR: no complex support in C */;
        float r8462 = r8458 * r8461;
        float r8463 = 7;
        float r8464 = 0;
        float r8465 = /* ERROR: no complex support in C */;
        float r8466 = r8462 + r8465;
        float r8467 = r8449 * r8449;
        float r8468 = r8448 / r8467;
        float r8469 = r8454 / r8467;
        float r8470 = r8468 - r8469;
        float r8471 = r8459 / r8467;
        float r8472 = r8455 + r8455;
        float r8473 = r8471 * r8472;
        float r8474 = /* ERROR: no complex support in C */;
        float r8475 = -r8449;
        float r8476 = /* ERROR: no complex support in C */;
        float r8477 = r8476 + r8461;
        float r8478 = r8474 * r8477;
        float r8479 = r8478 * r8461;
        float r8480 = r8466 + r8479;
        float r8481 = /* ERROR: no complex support in C */;
        return r8481;
}

double f_od() {
        double r8482 = 5;
        double r8483 = -r8482;
        double r8484 = 1;
        double r8485 = 2;
        double r8486 = r8484 / r8485;
        double r8487 = r8483 * r8486;
        double r8488 = 4;
        double r8489 = r8487 + r8488;
        double r8490 = 3;
        double r8491 = sqrt(r8490);
        double r8492 = r8491 / r8485;
        double r8493 = r8482 * r8492;
        double r8494 = /* ERROR: no complex support in C */;
        double r8495 = -r8484;
        double r8496 = r8495 / r8485;
        double r8497 = /* ERROR: no complex support in C */;
        double r8498 = r8494 * r8497;
        double r8499 = 7;
        double r8500 = 0;
        double r8501 = /* ERROR: no complex support in C */;
        double r8502 = r8498 + r8501;
        double r8503 = r8485 * r8485;
        double r8504 = r8484 / r8503;
        double r8505 = r8490 / r8503;
        double r8506 = r8504 - r8505;
        double r8507 = r8495 / r8503;
        double r8508 = r8491 + r8491;
        double r8509 = r8507 * r8508;
        double r8510 = /* ERROR: no complex support in C */;
        double r8511 = -r8485;
        double r8512 = /* ERROR: no complex support in C */;
        double r8513 = r8512 + r8497;
        double r8514 = r8510 * r8513;
        double r8515 = r8514 * r8497;
        double r8516 = r8502 + r8515;
        double r8517 = /* ERROR: no complex support in C */;
        return r8517;
}

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 r8518, r8519, r8520, r8521, r8522, r8523, r8524, r8525, r8526, r8527, r8528, r8529, r8530, r8531, r8532, r8533, r8534, r8535, r8536, r8537, r8538, r8539, r8540, r8541, r8542, r8543, r8544, r8545, r8546, r8547, r8548;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r8518, "1", 10, MPFR_RNDN);
        mpfr_init(r8519);
        mpfr_init_set_str(r8520, "2", 10, MPFR_RNDN);
        mpfr_init(r8521);
        mpfr_init_set_str(r8522, "3", 10, MPFR_RNDN);
        mpfr_init(r8523);
        mpfr_init(r8524);
        mpfr_init(r8525);
        mpfr_init(r8526);
        mpfr_init(r8527);
        mpfr_init(r8528);
        mpfr_init(r8529);
        mpfr_init_set_str(r8530, "0", 10, MPFR_RNDN);
        mpfr_init(r8531);
        mpfr_init(r8532);
        mpfr_init(r8533);
        mpfr_init(r8534);
        mpfr_init(r8535);
        mpfr_init_set_str(r8536, "5", 10, MPFR_RNDN);
        mpfr_init(r8537);
        mpfr_init(r8538);
        mpfr_init(r8539);
        mpfr_init(r8540);
        mpfr_init_set_str(r8541, "4", 10, MPFR_RNDN);
        mpfr_init(r8542);
        mpfr_init(r8543);
        mpfr_init(r8544);
        mpfr_init_set_str(r8545, "7", 10, MPFR_RNDN);
        mpfr_init(r8546);
        mpfr_init(r8547);
        mpfr_init(r8548);
}

double f_im() {
        ;
        mpfr_neg(r8519, r8518, MPFR_RNDN);
        ;
        mpfr_div(r8521, r8519, r8520, MPFR_RNDN);
        ;
        mpfr_sqrt(r8523, r8522, MPFR_RNDN);
        mpfr_div(r8524, r8523, r8520, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_mul(r8526, r8525, r8525, MPFR_RNDN);
        mpfr_mul(r8527, r8526, r8525, MPFR_RNDN);
        mpfr_mul(r8528, r8527, r8525, MPFR_RNDN);
        mpfr_neg(r8529, r8520, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8532, r8531, r8525, MPFR_RNDN);
        mpfr_mul(r8533, r8532, r8525, MPFR_RNDN);
        mpfr_mul(r8534, r8533, r8525, MPFR_RNDN);
        mpfr_add(r8535, r8528, r8534, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8538, r8537, r8525, MPFR_RNDN);
        mpfr_mul(r8539, r8538, r8525, MPFR_RNDN);
        mpfr_add(r8540, r8535, r8539, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8543, r8542, r8525, MPFR_RNDN);
        mpfr_add(r8544, r8540, r8543, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_add(r8547, r8544, r8546, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        return mpfr_get_d(r8548, MPFR_RNDN);
}

static mpfr_t r8549, r8550, r8551, r8552, r8553, r8554, r8555, r8556, r8557, r8558, r8559, r8560, r8561, r8562, r8563, r8564, r8565, r8566, r8567, r8568, r8569, r8570, r8571, r8572, r8573, r8574, r8575, r8576, r8577, r8578, r8579, r8580, r8581, r8582, r8583, r8584;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r8549, "5", 10, MPFR_RNDN);
        mpfr_init(r8550);
        mpfr_init_set_str(r8551, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r8552, "2", 10, MPFR_RNDN);
        mpfr_init(r8553);
        mpfr_init(r8554);
        mpfr_init_set_str(r8555, "4", 10, MPFR_RNDN);
        mpfr_init(r8556);
        mpfr_init_set_str(r8557, "3", 10, MPFR_RNDN);
        mpfr_init(r8558);
        mpfr_init(r8559);
        mpfr_init(r8560);
        mpfr_init(r8561);
        mpfr_init(r8562);
        mpfr_init(r8563);
        mpfr_init(r8564);
        mpfr_init(r8565);
        mpfr_init_set_str(r8566, "7", 10, MPFR_RNDN);
        mpfr_init_set_str(r8567, "0", 10, MPFR_RNDN);
        mpfr_init(r8568);
        mpfr_init(r8569);
        mpfr_init(r8570);
        mpfr_init(r8571);
        mpfr_init(r8572);
        mpfr_init(r8573);
        mpfr_init(r8574);
        mpfr_init(r8575);
        mpfr_init(r8576);
        mpfr_init(r8577);
        mpfr_init(r8578);
        mpfr_init(r8579);
        mpfr_init(r8580);
        mpfr_init(r8581);
        mpfr_init(r8582);
        mpfr_init(r8583);
        mpfr_init(r8584);
}

double f_fm() {
        ;
        mpfr_neg(r8550, r8549, MPFR_RNDN);
        ;
        ;
        mpfr_div(r8553, r8551, r8552, MPFR_RNDN);
        mpfr_mul(r8554, r8550, r8553, MPFR_RNDN);
        ;
        mpfr_add(r8556, r8554, r8555, MPFR_RNDN);
        ;
        mpfr_sqrt(r8558, r8557, MPFR_RNDN);
        mpfr_div(r8559, r8558, r8552, MPFR_RNDN);
        mpfr_mul(r8560, r8549, r8559, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_neg(r8562, r8551, MPFR_RNDN);
        mpfr_div(r8563, r8562, r8552, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_mul(r8565, r8561, r8564, MPFR_RNDN);
        ;
        ;
        /* ERROR: no complex support in C */;
        mpfr_add(r8569, r8565, r8568, MPFR_RNDN);
        mpfr_mul(r8570, r8552, r8552, MPFR_RNDN);
        mpfr_div(r8571, r8551, r8570, MPFR_RNDN);
        mpfr_div(r8572, r8557, r8570, MPFR_RNDN);
        mpfr_sub(r8573, r8571, r8572, MPFR_RNDN);
        mpfr_div(r8574, r8562, r8570, MPFR_RNDN);
        mpfr_add(r8575, r8558, r8558, MPFR_RNDN);
        mpfr_mul(r8576, r8574, r8575, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_neg(r8578, r8552, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_add(r8580, r8579, r8564, MPFR_RNDN);
        mpfr_mul(r8581, r8577, r8580, MPFR_RNDN);
        mpfr_mul(r8582, r8581, r8564, MPFR_RNDN);
        mpfr_add(r8583, r8569, r8582, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        return mpfr_get_d(r8584, MPFR_RNDN);
}

static mpfr_t r8585, r8586, r8587, r8588, r8589, r8590, r8591, r8592, r8593, r8594, r8595, r8596, r8597, r8598, r8599, r8600, r8601, r8602, r8603, r8604, r8605, r8606, r8607, r8608, r8609, r8610, r8611, r8612, r8613, r8614, r8615, r8616, r8617, r8618, r8619, r8620;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r8585, "5", 10, MPFR_RNDN);
        mpfr_init(r8586);
        mpfr_init_set_str(r8587, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r8588, "2", 10, MPFR_RNDN);
        mpfr_init(r8589);
        mpfr_init(r8590);
        mpfr_init_set_str(r8591, "4", 10, MPFR_RNDN);
        mpfr_init(r8592);
        mpfr_init_set_str(r8593, "3", 10, MPFR_RNDN);
        mpfr_init(r8594);
        mpfr_init(r8595);
        mpfr_init(r8596);
        mpfr_init(r8597);
        mpfr_init(r8598);
        mpfr_init(r8599);
        mpfr_init(r8600);
        mpfr_init(r8601);
        mpfr_init_set_str(r8602, "7", 10, MPFR_RNDN);
        mpfr_init_set_str(r8603, "0", 10, MPFR_RNDN);
        mpfr_init(r8604);
        mpfr_init(r8605);
        mpfr_init(r8606);
        mpfr_init(r8607);
        mpfr_init(r8608);
        mpfr_init(r8609);
        mpfr_init(r8610);
        mpfr_init(r8611);
        mpfr_init(r8612);
        mpfr_init(r8613);
        mpfr_init(r8614);
        mpfr_init(r8615);
        mpfr_init(r8616);
        mpfr_init(r8617);
        mpfr_init(r8618);
        mpfr_init(r8619);
        mpfr_init(r8620);
}

double f_dm() {
        ;
        mpfr_neg(r8586, r8585, MPFR_RNDN);
        ;
        ;
        mpfr_div(r8589, r8587, r8588, MPFR_RNDN);
        mpfr_mul(r8590, r8586, r8589, MPFR_RNDN);
        ;
        mpfr_add(r8592, r8590, r8591, MPFR_RNDN);
        ;
        mpfr_sqrt(r8594, r8593, MPFR_RNDN);
        mpfr_div(r8595, r8594, r8588, MPFR_RNDN);
        mpfr_mul(r8596, r8585, r8595, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_neg(r8598, r8587, MPFR_RNDN);
        mpfr_div(r8599, r8598, r8588, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_mul(r8601, r8597, r8600, MPFR_RNDN);
        ;
        ;
        /* ERROR: no complex support in C */;
        mpfr_add(r8605, r8601, r8604, MPFR_RNDN);
        mpfr_mul(r8606, r8588, r8588, MPFR_RNDN);
        mpfr_div(r8607, r8587, r8606, MPFR_RNDN);
        mpfr_div(r8608, r8593, r8606, MPFR_RNDN);
        mpfr_sub(r8609, r8607, r8608, MPFR_RNDN);
        mpfr_div(r8610, r8598, r8606, MPFR_RNDN);
        mpfr_add(r8611, r8594, r8594, MPFR_RNDN);
        mpfr_mul(r8612, r8610, r8611, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_neg(r8614, r8588, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_add(r8616, r8615, r8600, MPFR_RNDN);
        mpfr_mul(r8617, r8613, r8616, MPFR_RNDN);
        mpfr_mul(r8618, r8617, r8600, MPFR_RNDN);
        mpfr_add(r8619, r8605, r8618, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        return mpfr_get_d(r8620, MPFR_RNDN);
}