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

char *name = "x / (x^2 + 1)";

double f_if(float x) {
        float r22598 = x;
        float r22599 = r22598 * r22598;
        float r22600 = 1;
        float r22601 = r22599 + r22600;
        float r22602 = r22598 / r22601;
        return r22602;
}

double f_id(double x) {
        double r22603 = x;
        double r22604 = r22603 * r22603;
        double r22605 = 1;
        double r22606 = r22604 + r22605;
        double r22607 = r22603 / r22606;
        return r22607;
}


double f_of(float x) {
        float r22608 = 1;
        float r22609 = x;
        float r22610 = 5;
        float r22611 = pow(r22609, r22610);
        float r22612 = r22608 / r22611;
        float r22613 = r22608 / r22609;
        float r22614 = r22612 + r22613;
        float r22615 = 3;
        float r22616 = pow(r22609, r22615);
        float r22617 = r22608 / r22616;
        float r22618 = r22614 - r22617;
        float r22619 = -8.98475646888388e-17;
        bool r22620 = r22618 <= r22619;
        float r22621 = hypot(r22608, r22609);
        float r22622 = r22609 / r22621;
        float r22623 = r22609 * r22609;
        float r22624 = r22623 + r22608;
        float r22625 = sqrt(r22624);
        float r22626 = r22622 / r22625;
        float r22627 = 3.9064990173073113e-25;
        bool r22628 = r22618 <= r22627;
        float r22629 = fma(r22609, r22609, r22608);
        float r22630 = r22608 / r22629;
        float r22631 = r22609 * r22630;
        float r22632 = r22628 ? r22618 : r22631;
        float r22633 = r22620 ? r22626 : r22632;
        return r22633;
}

double f_od(double x) {
        double r22634 = 1;
        double r22635 = x;
        double r22636 = 5;
        double r22637 = pow(r22635, r22636);
        double r22638 = r22634 / r22637;
        double r22639 = r22634 / r22635;
        double r22640 = r22638 + r22639;
        double r22641 = 3;
        double r22642 = pow(r22635, r22641);
        double r22643 = r22634 / r22642;
        double r22644 = r22640 - r22643;
        double r22645 = -8.98475646888388e-17;
        bool r22646 = r22644 <= r22645;
        double r22647 = hypot(r22634, r22635);
        double r22648 = r22635 / r22647;
        double r22649 = r22635 * r22635;
        double r22650 = r22649 + r22634;
        double r22651 = sqrt(r22650);
        double r22652 = r22648 / r22651;
        double r22653 = 3.9064990173073113e-25;
        bool r22654 = r22644 <= r22653;
        double r22655 = fma(r22635, r22635, r22634);
        double r22656 = r22634 / r22655;
        double r22657 = r22635 * r22656;
        double r22658 = r22654 ? r22644 : r22657;
        double r22659 = r22646 ? r22652 : r22658;
        return r22659;
}

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 r22660, r22661, r22662, r22663, r22664;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r22660);
        mpfr_init(r22661);
        mpfr_init_set_str(r22662, "1", 10, MPFR_RNDN);
        mpfr_init(r22663);
        mpfr_init(r22664);
}

double f_im(double x) {
        mpfr_set_d(r22660, x, MPFR_RNDN);
        mpfr_mul(r22661, r22660, r22660, MPFR_RNDN);
        ;
        mpfr_add(r22663, r22661, r22662, MPFR_RNDN);
        mpfr_div(r22664, r22660, r22663, MPFR_RNDN);
        return mpfr_get_d(r22664, MPFR_RNDN);
}

static mpfr_t 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;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r22665, "1", 10, MPFR_RNDN);
        mpfr_init(r22666);
        mpfr_init_set_str(r22667, "5", 10, MPFR_RNDN);
        mpfr_init(r22668);
        mpfr_init(r22669);
        mpfr_init(r22670);
        mpfr_init(r22671);
        mpfr_init_set_str(r22672, "3", 10, MPFR_RNDN);
        mpfr_init(r22673);
        mpfr_init(r22674);
        mpfr_init(r22675);
        mpfr_init_set_str(r22676, "-8.98475646888388e-17", 10, MPFR_RNDN);
        mpfr_init(r22677);
        mpfr_init(r22678);
        mpfr_init(r22679);
        mpfr_init(r22680);
        mpfr_init(r22681);
        mpfr_init(r22682);
        mpfr_init(r22683);
        mpfr_init_set_str(r22684, "3.9064990173073113e-25", 10, MPFR_RNDN);
        mpfr_init(r22685);
        mpfr_init(r22686);
        mpfr_init(r22687);
        mpfr_init(r22688);
        mpfr_init(r22689);
        mpfr_init(r22690);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r22666, x, MPFR_RNDN);
        ;
        mpfr_pow(r22668, r22666, r22667, MPFR_RNDN);
        mpfr_div(r22669, r22665, r22668, MPFR_RNDN);
        mpfr_div(r22670, r22665, r22666, MPFR_RNDN);
        mpfr_add(r22671, r22669, r22670, MPFR_RNDN);
        ;
        mpfr_pow(r22673, r22666, r22672, MPFR_RNDN);
        mpfr_div(r22674, r22665, r22673, MPFR_RNDN);
        mpfr_sub(r22675, r22671, r22674, MPFR_RNDN);
        ;
        mpfr_set_si(r22677, mpfr_cmp(r22675, r22676) <= 0, MPFR_RNDN);
        mpfr_hypot(r22678, r22665, r22666, MPFR_RNDN);
        mpfr_div(r22679, r22666, r22678, MPFR_RNDN);
        mpfr_mul(r22680, r22666, r22666, MPFR_RNDN);
        mpfr_add(r22681, r22680, r22665, MPFR_RNDN);
        mpfr_sqrt(r22682, r22681, MPFR_RNDN);
        mpfr_div(r22683, r22679, r22682, MPFR_RNDN);
        ;
        mpfr_set_si(r22685, mpfr_cmp(r22675, r22684) <= 0, MPFR_RNDN);
        mpfr_fma(r22686, r22666, r22666, r22665, MPFR_RNDN);
        mpfr_div(r22687, r22665, r22686, MPFR_RNDN);
        mpfr_mul(r22688, r22666, r22687, MPFR_RNDN);
        if (mpfr_get_si(r22685, MPFR_RNDN)) { mpfr_set(r22689, r22675, MPFR_RNDN); } else { mpfr_set(r22689, r22688, MPFR_RNDN); };
        if (mpfr_get_si(r22677, MPFR_RNDN)) { mpfr_set(r22690, r22683, MPFR_RNDN); } else { mpfr_set(r22690, r22689, MPFR_RNDN); };
        return mpfr_get_d(r22690, MPFR_RNDN);
}

static mpfr_t r22691, r22692, r22693, r22694, r22695, r22696, r22697, r22698, r22699, r22700, r22701, r22702, r22703, r22704, r22705, r22706, r22707, r22708, r22709, r22710, r22711, r22712, r22713, r22714, r22715, r22716;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r22691, "1", 10, MPFR_RNDN);
        mpfr_init(r22692);
        mpfr_init_set_str(r22693, "5", 10, MPFR_RNDN);
        mpfr_init(r22694);
        mpfr_init(r22695);
        mpfr_init(r22696);
        mpfr_init(r22697);
        mpfr_init_set_str(r22698, "3", 10, MPFR_RNDN);
        mpfr_init(r22699);
        mpfr_init(r22700);
        mpfr_init(r22701);
        mpfr_init_set_str(r22702, "-8.98475646888388e-17", 10, MPFR_RNDN);
        mpfr_init(r22703);
        mpfr_init(r22704);
        mpfr_init(r22705);
        mpfr_init(r22706);
        mpfr_init(r22707);
        mpfr_init(r22708);
        mpfr_init(r22709);
        mpfr_init_set_str(r22710, "3.9064990173073113e-25", 10, MPFR_RNDN);
        mpfr_init(r22711);
        mpfr_init(r22712);
        mpfr_init(r22713);
        mpfr_init(r22714);
        mpfr_init(r22715);
        mpfr_init(r22716);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r22692, x, MPFR_RNDN);
        ;
        mpfr_pow(r22694, r22692, r22693, MPFR_RNDN);
        mpfr_div(r22695, r22691, r22694, MPFR_RNDN);
        mpfr_div(r22696, r22691, r22692, MPFR_RNDN);
        mpfr_add(r22697, r22695, r22696, MPFR_RNDN);
        ;
        mpfr_pow(r22699, r22692, r22698, MPFR_RNDN);
        mpfr_div(r22700, r22691, r22699, MPFR_RNDN);
        mpfr_sub(r22701, r22697, r22700, MPFR_RNDN);
        ;
        mpfr_set_si(r22703, mpfr_cmp(r22701, r22702) <= 0, MPFR_RNDN);
        mpfr_hypot(r22704, r22691, r22692, MPFR_RNDN);
        mpfr_div(r22705, r22692, r22704, MPFR_RNDN);
        mpfr_mul(r22706, r22692, r22692, MPFR_RNDN);
        mpfr_add(r22707, r22706, r22691, MPFR_RNDN);
        mpfr_sqrt(r22708, r22707, MPFR_RNDN);
        mpfr_div(r22709, r22705, r22708, MPFR_RNDN);
        ;
        mpfr_set_si(r22711, mpfr_cmp(r22701, r22710) <= 0, MPFR_RNDN);
        mpfr_fma(r22712, r22692, r22692, r22691, MPFR_RNDN);
        mpfr_div(r22713, r22691, r22712, MPFR_RNDN);
        mpfr_mul(r22714, r22692, r22713, MPFR_RNDN);
        if (mpfr_get_si(r22711, MPFR_RNDN)) { mpfr_set(r22715, r22701, MPFR_RNDN); } else { mpfr_set(r22715, r22714, MPFR_RNDN); };
        if (mpfr_get_si(r22703, MPFR_RNDN)) { mpfr_set(r22716, r22709, MPFR_RNDN); } else { mpfr_set(r22716, r22715, MPFR_RNDN); };
        return mpfr_get_d(r22716, MPFR_RNDN);
}

