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

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

double f_if() {
        float r22538 = 1;
        float r22539 = -r22538;
        float r22540 = 2;
        float r22541 = r22539 / r22540;
        float r22542 = 3;
        float r22543 = sqrt(r22542);
        float r22544 = r22543 / r22540;
        float r22545 = /* ERROR: no complex support in C */;
        float r22546 = r22545 * r22545;
        float r22547 = r22546 * r22545;
        float r22548 = r22547 * r22545;
        float r22549 = -r22540;
        float r22550 = 0;
        float r22551 = /* ERROR: no complex support in C */;
        float r22552 = r22551 * r22545;
        float r22553 = r22552 * r22545;
        float r22554 = r22553 * r22545;
        float r22555 = r22548 + r22554;
        float r22556 = 5;
        float r22557 = /* ERROR: no complex support in C */;
        float r22558 = r22557 * r22545;
        float r22559 = r22558 * r22545;
        float r22560 = r22555 + r22559;
        float r22561 = 4;
        float r22562 = /* ERROR: no complex support in C */;
        float r22563 = r22562 * r22545;
        float r22564 = r22560 + r22563;
        float r22565 = 7;
        float r22566 = /* ERROR: no complex support in C */;
        float r22567 = r22564 + r22566;
        float r22568 = /* ERROR: no complex support in C */;
        return r22568;
}

double f_id() {
        double r22569 = 1;
        double r22570 = -r22569;
        double r22571 = 2;
        double r22572 = r22570 / r22571;
        double r22573 = 3;
        double r22574 = sqrt(r22573);
        double r22575 = r22574 / r22571;
        double r22576 = /* ERROR: no complex support in C */;
        double r22577 = r22576 * r22576;
        double r22578 = r22577 * r22576;
        double r22579 = r22578 * r22576;
        double r22580 = -r22571;
        double r22581 = 0;
        double r22582 = /* ERROR: no complex support in C */;
        double r22583 = r22582 * r22576;
        double r22584 = r22583 * r22576;
        double r22585 = r22584 * r22576;
        double r22586 = r22579 + r22585;
        double r22587 = 5;
        double r22588 = /* ERROR: no complex support in C */;
        double r22589 = r22588 * r22576;
        double r22590 = r22589 * r22576;
        double r22591 = r22586 + r22590;
        double r22592 = 4;
        double r22593 = /* ERROR: no complex support in C */;
        double r22594 = r22593 * r22576;
        double r22595 = r22591 + r22594;
        double r22596 = 7;
        double r22597 = /* ERROR: no complex support in C */;
        double r22598 = r22595 + r22597;
        double r22599 = /* ERROR: no complex support in C */;
        return r22599;
}


double f_of() {
        float r22600 = 1;
        float r22601 = -r22600;
        float r22602 = 2;
        float r22603 = r22601 / r22602;
        float r22604 = 3;
        float r22605 = sqrt(r22604);
        float r22606 = r22605 / r22602;
        float r22607 = /* ERROR: no complex support in C */;
        float r22608 = r22607 * r22607;
        float r22609 = r22608 * r22607;
        float r22610 = r22609 * r22607;
        float r22611 = -r22602;
        float r22612 = 0;
        float r22613 = /* ERROR: no complex support in C */;
        float r22614 = r22613 * r22607;
        float r22615 = r22614 * r22607;
        float r22616 = r22615 * r22607;
        float r22617 = r22610 + r22616;
        float r22618 = 5;
        float r22619 = /* ERROR: no complex support in C */;
        float r22620 = r22619 * r22607;
        float r22621 = r22620 * r22607;
        float r22622 = r22617 + r22621;
        float r22623 = 4;
        float r22624 = /* ERROR: no complex support in C */;
        float r22625 = r22624 * r22607;
        float r22626 = r22622 + r22625;
        float r22627 = 7;
        float r22628 = /* ERROR: no complex support in C */;
        float r22629 = r22626 + r22628;
        float r22630 = /* ERROR: no complex support in C */;
        return r22630;
}

double f_od() {
        double r22631 = 1;
        double r22632 = -r22631;
        double r22633 = 2;
        double r22634 = r22632 / r22633;
        double r22635 = 3;
        double r22636 = sqrt(r22635);
        double r22637 = r22636 / r22633;
        double r22638 = /* ERROR: no complex support in C */;
        double r22639 = r22638 * r22638;
        double r22640 = r22639 * r22638;
        double r22641 = r22640 * r22638;
        double r22642 = -r22633;
        double r22643 = 0;
        double r22644 = /* ERROR: no complex support in C */;
        double r22645 = r22644 * r22638;
        double r22646 = r22645 * r22638;
        double r22647 = r22646 * r22638;
        double r22648 = r22641 + r22647;
        double r22649 = 5;
        double r22650 = /* ERROR: no complex support in C */;
        double r22651 = r22650 * r22638;
        double r22652 = r22651 * r22638;
        double r22653 = r22648 + r22652;
        double r22654 = 4;
        double r22655 = /* ERROR: no complex support in C */;
        double r22656 = r22655 * r22638;
        double r22657 = r22653 + r22656;
        double r22658 = 7;
        double r22659 = /* ERROR: no complex support in C */;
        double r22660 = r22657 + r22659;
        double r22661 = /* ERROR: no complex support in C */;
        return r22661;
}

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 r22662, r22663, r22664, r22665, r22666, r22667, r22668, r22669, r22670, r22671, r22672, r22673, r22674, r22675, r22676, r22677, r22678, r22679, r22680, r22681, r22682, r22683, r22684, r22685, r22686, r22687, r22688, r22689, r22690, r22691, r22692;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r22662, "1", 10, MPFR_RNDN);
        mpfr_init(r22663);
        mpfr_init_set_str(r22664, "2", 10, MPFR_RNDN);
        mpfr_init(r22665);
        mpfr_init_set_str(r22666, "3", 10, MPFR_RNDN);
        mpfr_init(r22667);
        mpfr_init(r22668);
        mpfr_init(r22669);
        mpfr_init(r22670);
        mpfr_init(r22671);
        mpfr_init(r22672);
        mpfr_init(r22673);
        mpfr_init_set_str(r22674, "0", 10, MPFR_RNDN);
        mpfr_init(r22675);
        mpfr_init(r22676);
        mpfr_init(r22677);
        mpfr_init(r22678);
        mpfr_init(r22679);
        mpfr_init_set_str(r22680, "5", 10, MPFR_RNDN);
        mpfr_init(r22681);
        mpfr_init(r22682);
        mpfr_init(r22683);
        mpfr_init(r22684);
        mpfr_init_set_str(r22685, "4", 10, MPFR_RNDN);
        mpfr_init(r22686);
        mpfr_init(r22687);
        mpfr_init(r22688);
        mpfr_init_set_str(r22689, "7", 10, MPFR_RNDN);
        mpfr_init(r22690);
        mpfr_init(r22691);
        mpfr_init(r22692);
}

double f_im() {
        ;
        mpfr_neg(r22663, r22662, MPFR_RNDN);
        ;
        mpfr_div(r22665, r22663, r22664, MPFR_RNDN);
        ;
        mpfr_sqrt(r22667, r22666, MPFR_RNDN);
        mpfr_div(r22668, r22667, r22664, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_mul(r22670, r22669, r22669, MPFR_RNDN);
        mpfr_mul(r22671, r22670, r22669, MPFR_RNDN);
        mpfr_mul(r22672, r22671, r22669, MPFR_RNDN);
        mpfr_neg(r22673, r22664, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r22676, r22675, r22669, MPFR_RNDN);
        mpfr_mul(r22677, r22676, r22669, MPFR_RNDN);
        mpfr_mul(r22678, r22677, r22669, MPFR_RNDN);
        mpfr_add(r22679, r22672, r22678, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r22682, r22681, r22669, MPFR_RNDN);
        mpfr_mul(r22683, r22682, r22669, MPFR_RNDN);
        mpfr_add(r22684, r22679, r22683, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r22687, r22686, r22669, MPFR_RNDN);
        mpfr_add(r22688, r22684, r22687, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_add(r22691, r22688, r22690, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        return mpfr_get_d(r22692, MPFR_RNDN);
}

static mpfr_t r22693, r22694, r22695, r22696, r22697, r22698, r22699, r22700, r22701, r22702, r22703, r22704, r22705, r22706, r22707, r22708, r22709, r22710, r22711, r22712, r22713, r22714, r22715, r22716, r22717, r22718, r22719, r22720, r22721, r22722, r22723;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r22693, "1", 10, MPFR_RNDN);
        mpfr_init(r22694);
        mpfr_init_set_str(r22695, "2", 10, MPFR_RNDN);
        mpfr_init(r22696);
        mpfr_init_set_str(r22697, "3", 10, MPFR_RNDN);
        mpfr_init(r22698);
        mpfr_init(r22699);
        mpfr_init(r22700);
        mpfr_init(r22701);
        mpfr_init(r22702);
        mpfr_init(r22703);
        mpfr_init(r22704);
        mpfr_init_set_str(r22705, "0", 10, MPFR_RNDN);
        mpfr_init(r22706);
        mpfr_init(r22707);
        mpfr_init(r22708);
        mpfr_init(r22709);
        mpfr_init(r22710);
        mpfr_init_set_str(r22711, "5", 10, MPFR_RNDN);
        mpfr_init(r22712);
        mpfr_init(r22713);
        mpfr_init(r22714);
        mpfr_init(r22715);
        mpfr_init_set_str(r22716, "4", 10, MPFR_RNDN);
        mpfr_init(r22717);
        mpfr_init(r22718);
        mpfr_init(r22719);
        mpfr_init_set_str(r22720, "7", 10, MPFR_RNDN);
        mpfr_init(r22721);
        mpfr_init(r22722);
        mpfr_init(r22723);
}

double f_fm() {
        ;
        mpfr_neg(r22694, r22693, MPFR_RNDN);
        ;
        mpfr_div(r22696, r22694, r22695, MPFR_RNDN);
        ;
        mpfr_sqrt(r22698, r22697, MPFR_RNDN);
        mpfr_div(r22699, r22698, r22695, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_mul(r22701, r22700, r22700, MPFR_RNDN);
        mpfr_mul(r22702, r22701, r22700, MPFR_RNDN);
        mpfr_mul(r22703, r22702, r22700, MPFR_RNDN);
        mpfr_neg(r22704, r22695, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r22707, r22706, r22700, MPFR_RNDN);
        mpfr_mul(r22708, r22707, r22700, MPFR_RNDN);
        mpfr_mul(r22709, r22708, r22700, MPFR_RNDN);
        mpfr_add(r22710, r22703, r22709, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r22713, r22712, r22700, MPFR_RNDN);
        mpfr_mul(r22714, r22713, r22700, MPFR_RNDN);
        mpfr_add(r22715, r22710, r22714, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r22718, r22717, r22700, MPFR_RNDN);
        mpfr_add(r22719, r22715, r22718, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_add(r22722, r22719, r22721, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        return mpfr_get_d(r22723, MPFR_RNDN);
}

static mpfr_t r22724, r22725, r22726, r22727, r22728, r22729, r22730, r22731, r22732, r22733, r22734, r22735, r22736, r22737, r22738, r22739, r22740, r22741, r22742, r22743, r22744, r22745, r22746, r22747, r22748, r22749, r22750, r22751, r22752, r22753, r22754;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r22724, "1", 10, MPFR_RNDN);
        mpfr_init(r22725);
        mpfr_init_set_str(r22726, "2", 10, MPFR_RNDN);
        mpfr_init(r22727);
        mpfr_init_set_str(r22728, "3", 10, MPFR_RNDN);
        mpfr_init(r22729);
        mpfr_init(r22730);
        mpfr_init(r22731);
        mpfr_init(r22732);
        mpfr_init(r22733);
        mpfr_init(r22734);
        mpfr_init(r22735);
        mpfr_init_set_str(r22736, "0", 10, MPFR_RNDN);
        mpfr_init(r22737);
        mpfr_init(r22738);
        mpfr_init(r22739);
        mpfr_init(r22740);
        mpfr_init(r22741);
        mpfr_init_set_str(r22742, "5", 10, MPFR_RNDN);
        mpfr_init(r22743);
        mpfr_init(r22744);
        mpfr_init(r22745);
        mpfr_init(r22746);
        mpfr_init_set_str(r22747, "4", 10, MPFR_RNDN);
        mpfr_init(r22748);
        mpfr_init(r22749);
        mpfr_init(r22750);
        mpfr_init_set_str(r22751, "7", 10, MPFR_RNDN);
        mpfr_init(r22752);
        mpfr_init(r22753);
        mpfr_init(r22754);
}

double f_dm() {
        ;
        mpfr_neg(r22725, r22724, MPFR_RNDN);
        ;
        mpfr_div(r22727, r22725, r22726, MPFR_RNDN);
        ;
        mpfr_sqrt(r22729, r22728, MPFR_RNDN);
        mpfr_div(r22730, r22729, r22726, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_mul(r22732, r22731, r22731, MPFR_RNDN);
        mpfr_mul(r22733, r22732, r22731, MPFR_RNDN);
        mpfr_mul(r22734, r22733, r22731, MPFR_RNDN);
        mpfr_neg(r22735, r22726, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r22738, r22737, r22731, MPFR_RNDN);
        mpfr_mul(r22739, r22738, r22731, MPFR_RNDN);
        mpfr_mul(r22740, r22739, r22731, MPFR_RNDN);
        mpfr_add(r22741, r22734, r22740, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r22744, r22743, r22731, MPFR_RNDN);
        mpfr_mul(r22745, r22744, r22731, MPFR_RNDN);
        mpfr_add(r22746, r22741, r22745, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r22749, r22748, r22731, MPFR_RNDN);
        mpfr_add(r22750, r22746, r22749, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_add(r22753, r22750, r22752, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        return mpfr_get_d(r22754, MPFR_RNDN);
}