#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 r25047 = x;
        float r25048 = r25047 * r25047;
        float r25049 = 1;
        float r25050 = r25048 + r25049;
        float r25051 = r25047 / r25050;
        return r25051;
}

double f_id(double x) {
        double r25052 = x;
        double r25053 = r25052 * r25052;
        double r25054 = 1;
        double r25055 = r25053 + r25054;
        double r25056 = r25052 / r25055;
        return r25056;
}


double f_of(float x) {
        float r25057 = x;
        float r25058 = -578.1237044981676;
        bool r25059 = r25057 <= r25058;
        float r25060 = 402.7494417162069;
        bool r25061 = r25057 <= r25060;
        float r25062 = !r25061;
        bool r25063 = r25059 || r25062;
        float r25064 = 1;
        float r25065 = 5;
        float r25066 = pow(r25057, r25065);
        float r25067 = r25064 / r25066;
        float r25068 = r25064 / r25057;
        float r25069 = r25067 + r25068;
        float r25070 = 3;
        float r25071 = pow(r25057, r25070);
        float r25072 = r25064 / r25071;
        float r25073 = r25069 - r25072;
        float r25074 = hypot(r25064, r25057);
        float r25075 = r25057 / r25074;
        float r25076 = r25057 * r25057;
        float r25077 = r25064 + r25076;
        float r25078 = sqrt(r25077);
        float r25079 = r25075 / r25078;
        float r25080 = r25063 ? r25073 : r25079;
        return r25080;
}

double f_od(double x) {
        double r25081 = x;
        double r25082 = -578.1237044981676;
        bool r25083 = r25081 <= r25082;
        double r25084 = 402.7494417162069;
        bool r25085 = r25081 <= r25084;
        double r25086 = !r25085;
        bool r25087 = r25083 || r25086;
        double r25088 = 1;
        double r25089 = 5;
        double r25090 = pow(r25081, r25089);
        double r25091 = r25088 / r25090;
        double r25092 = r25088 / r25081;
        double r25093 = r25091 + r25092;
        double r25094 = 3;
        double r25095 = pow(r25081, r25094);
        double r25096 = r25088 / r25095;
        double r25097 = r25093 - r25096;
        double r25098 = hypot(r25088, r25081);
        double r25099 = r25081 / r25098;
        double r25100 = r25081 * r25081;
        double r25101 = r25088 + r25100;
        double r25102 = sqrt(r25101);
        double r25103 = r25099 / r25102;
        double r25104 = r25087 ? r25097 : r25103;
        return r25104;
}

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 r25105, r25106, r25107, r25108, r25109;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r25105);
        mpfr_init(r25106);
        mpfr_init_set_str(r25107, "1", 10, MPFR_RNDN);
        mpfr_init(r25108);
        mpfr_init(r25109);
}

double f_im(double x) {
        mpfr_set_d(r25105, x, MPFR_RNDN);
        mpfr_mul(r25106, r25105, r25105, MPFR_RNDN);
        ;
        mpfr_add(r25108, r25106, r25107, MPFR_RNDN);
        mpfr_div(r25109, r25105, r25108, MPFR_RNDN);
        return mpfr_get_d(r25109, MPFR_RNDN);
}

static mpfr_t r25110, r25111, r25112, r25113, r25114, r25115, r25116, r25117, r25118, r25119, r25120, r25121, r25122, r25123, r25124, r25125, r25126, r25127, r25128, r25129, r25130, r25131, r25132, r25133;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25110);
        mpfr_init_set_str(r25111, "-578.1237044981676", 10, MPFR_RNDN);
        mpfr_init(r25112);
        mpfr_init_set_str(r25113, "402.7494417162069", 10, MPFR_RNDN);
        mpfr_init(r25114);
        mpfr_init(r25115);
        mpfr_init(r25116);
        mpfr_init_set_str(r25117, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r25118, "5", 10, MPFR_RNDN);
        mpfr_init(r25119);
        mpfr_init(r25120);
        mpfr_init(r25121);
        mpfr_init(r25122);
        mpfr_init_set_str(r25123, "3", 10, MPFR_RNDN);
        mpfr_init(r25124);
        mpfr_init(r25125);
        mpfr_init(r25126);
        mpfr_init(r25127);
        mpfr_init(r25128);
        mpfr_init(r25129);
        mpfr_init(r25130);
        mpfr_init(r25131);
        mpfr_init(r25132);
        mpfr_init(r25133);
}

double f_fm(double x) {
        mpfr_set_d(r25110, x, MPFR_RNDN);
        ;
        mpfr_set_si(r25112, mpfr_cmp(r25110, r25111) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r25114, mpfr_cmp(r25110, r25113) <= 0, MPFR_RNDN);
        mpfr_set_si(r25115, !mpfr_get_si(r25114, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r25116, mpfr_get_si(r25112, MPFR_RNDN) || mpfr_get_si(r25115, MPFR_RNDN), MPFR_RNDN);
        ;
        ;
        mpfr_pow(r25119, r25110, r25118, MPFR_RNDN);
        mpfr_div(r25120, r25117, r25119, MPFR_RNDN);
        mpfr_div(r25121, r25117, r25110, MPFR_RNDN);
        mpfr_add(r25122, r25120, r25121, MPFR_RNDN);
        ;
        mpfr_pow(r25124, r25110, r25123, MPFR_RNDN);
        mpfr_div(r25125, r25117, r25124, MPFR_RNDN);
        mpfr_sub(r25126, r25122, r25125, MPFR_RNDN);
        mpfr_hypot(r25127, r25117, r25110, MPFR_RNDN);
        mpfr_div(r25128, r25110, r25127, MPFR_RNDN);
        mpfr_mul(r25129, r25110, r25110, MPFR_RNDN);
        mpfr_add(r25130, r25117, r25129, MPFR_RNDN);
        mpfr_sqrt(r25131, r25130, MPFR_RNDN);
        mpfr_div(r25132, r25128, r25131, MPFR_RNDN);
        if (mpfr_get_si(r25116, MPFR_RNDN)) { mpfr_set(r25133, r25126, MPFR_RNDN); } else { mpfr_set(r25133, r25132, MPFR_RNDN); };
        return mpfr_get_d(r25133, MPFR_RNDN);
}

static mpfr_t r25134, r25135, r25136, r25137, r25138, r25139, r25140, r25141, r25142, r25143, r25144, r25145, r25146, r25147, r25148, r25149, r25150, r25151, r25152, r25153, r25154, r25155, r25156, r25157;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25134);
        mpfr_init_set_str(r25135, "-578.1237044981676", 10, MPFR_RNDN);
        mpfr_init(r25136);
        mpfr_init_set_str(r25137, "402.7494417162069", 10, MPFR_RNDN);
        mpfr_init(r25138);
        mpfr_init(r25139);
        mpfr_init(r25140);
        mpfr_init_set_str(r25141, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r25142, "5", 10, MPFR_RNDN);
        mpfr_init(r25143);
        mpfr_init(r25144);
        mpfr_init(r25145);
        mpfr_init(r25146);
        mpfr_init_set_str(r25147, "3", 10, MPFR_RNDN);
        mpfr_init(r25148);
        mpfr_init(r25149);
        mpfr_init(r25150);
        mpfr_init(r25151);
        mpfr_init(r25152);
        mpfr_init(r25153);
        mpfr_init(r25154);
        mpfr_init(r25155);
        mpfr_init(r25156);
        mpfr_init(r25157);
}

double f_dm(double x) {
        mpfr_set_d(r25134, x, MPFR_RNDN);
        ;
        mpfr_set_si(r25136, mpfr_cmp(r25134, r25135) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r25138, mpfr_cmp(r25134, r25137) <= 0, MPFR_RNDN);
        mpfr_set_si(r25139, !mpfr_get_si(r25138, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r25140, mpfr_get_si(r25136, MPFR_RNDN) || mpfr_get_si(r25139, MPFR_RNDN), MPFR_RNDN);
        ;
        ;
        mpfr_pow(r25143, r25134, r25142, MPFR_RNDN);
        mpfr_div(r25144, r25141, r25143, MPFR_RNDN);
        mpfr_div(r25145, r25141, r25134, MPFR_RNDN);
        mpfr_add(r25146, r25144, r25145, MPFR_RNDN);
        ;
        mpfr_pow(r25148, r25134, r25147, MPFR_RNDN);
        mpfr_div(r25149, r25141, r25148, MPFR_RNDN);
        mpfr_sub(r25150, r25146, r25149, MPFR_RNDN);
        mpfr_hypot(r25151, r25141, r25134, MPFR_RNDN);
        mpfr_div(r25152, r25134, r25151, MPFR_RNDN);
        mpfr_mul(r25153, r25134, r25134, MPFR_RNDN);
        mpfr_add(r25154, r25141, r25153, MPFR_RNDN);
        mpfr_sqrt(r25155, r25154, MPFR_RNDN);
        mpfr_div(r25156, r25152, r25155, MPFR_RNDN);
        if (mpfr_get_si(r25140, MPFR_RNDN)) { mpfr_set(r25157, r25150, MPFR_RNDN); } else { mpfr_set(r25157, r25156, MPFR_RNDN); };
        return mpfr_get_d(r25157, MPFR_RNDN);
}