#include #include #include #include #include char *name = "x / (x^2 + 1)"; double f_if(float x) { float r25069 = x; float r25070 = r25069 * r25069; float r25071 = 1; float r25072 = r25070 + r25071; float r25073 = r25069 / r25072; return r25073; } double f_id(double x) { double r25074 = x; double r25075 = r25074 * r25074; double r25076 = 1; double r25077 = r25075 + r25076; double r25078 = r25074 / r25077; return r25078; } double f_of(float x) { float r25079 = x; float r25080 = -8.07298555250867e+25; bool r25081 = r25079 <= r25080; float r25082 = 1565.0996446725549; bool r25083 = r25079 <= r25082; float r25084 = !r25083; bool r25085 = r25081 || r25084; float r25086 = 1; float r25087 = 5; float r25088 = pow(r25079, r25087); float r25089 = r25086 / r25088; float r25090 = r25086 / r25079; float r25091 = r25089 + r25090; float r25092 = 3; float r25093 = pow(r25079, r25092); float r25094 = r25086 / r25093; float r25095 = r25091 - r25094; float r25096 = r25079 * r25079; float r25097 = r25096 + r25086; float r25098 = r25079 / r25097; float r25099 = r25085 ? r25095 : r25098; return r25099; } double f_od(double x) { double r25100 = x; double r25101 = -8.07298555250867e+25; bool r25102 = r25100 <= r25101; double r25103 = 1565.0996446725549; bool r25104 = r25100 <= r25103; double r25105 = !r25104; bool r25106 = r25102 || r25105; double r25107 = 1; double r25108 = 5; double r25109 = pow(r25100, r25108); double r25110 = r25107 / r25109; double r25111 = r25107 / r25100; double r25112 = r25110 + r25111; double r25113 = 3; double r25114 = pow(r25100, r25113); double r25115 = r25107 / r25114; double r25116 = r25112 - r25115; double r25117 = r25100 * r25100; double r25118 = r25117 + r25107; double r25119 = r25100 / r25118; double r25120 = r25106 ? r25116 : r25119; return r25120; } 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 r25121, r25122, r25123, r25124, r25125; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25121); mpfr_init(r25122); mpfr_init_set_str(r25123, "1", 10, MPFR_RNDN); mpfr_init(r25124); mpfr_init(r25125); } double f_im(double x) { mpfr_set_d(r25121, x, MPFR_RNDN); mpfr_mul(r25122, r25121, r25121, MPFR_RNDN); ; mpfr_add(r25124, r25122, r25123, MPFR_RNDN); mpfr_div(r25125, r25121, r25124, MPFR_RNDN); return mpfr_get_d(r25125, MPFR_RNDN); } static mpfr_t r25126, r25127, r25128, r25129, r25130, r25131, r25132, r25133, r25134, r25135, r25136, r25137, r25138, r25139, r25140, r25141, r25142, r25143, r25144, r25145, r25146; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25126); mpfr_init_set_str(r25127, "-8.07298555250867e+25", 10, MPFR_RNDN); mpfr_init(r25128); mpfr_init_set_str(r25129, "1565.0996446725549", 10, MPFR_RNDN); mpfr_init(r25130); mpfr_init(r25131); mpfr_init(r25132); mpfr_init_set_str(r25133, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25134, "5", 10, MPFR_RNDN); mpfr_init(r25135); mpfr_init(r25136); mpfr_init(r25137); mpfr_init(r25138); mpfr_init_set_str(r25139, "3", 10, MPFR_RNDN); mpfr_init(r25140); mpfr_init(r25141); mpfr_init(r25142); mpfr_init(r25143); mpfr_init(r25144); mpfr_init(r25145); mpfr_init(r25146); } double f_fm(double x) { mpfr_set_d(r25126, x, MPFR_RNDN); ; mpfr_set_si(r25128, mpfr_cmp(r25126, r25127) <= 0, MPFR_RNDN); ; mpfr_set_si(r25130, mpfr_cmp(r25126, r25129) <= 0, MPFR_RNDN); mpfr_set_si(r25131, !mpfr_get_si(r25130, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r25132, mpfr_get_si(r25128, MPFR_RNDN) || mpfr_get_si(r25131, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r25135, r25126, r25134, MPFR_RNDN); mpfr_div(r25136, r25133, r25135, MPFR_RNDN); mpfr_div(r25137, r25133, r25126, MPFR_RNDN); mpfr_add(r25138, r25136, r25137, MPFR_RNDN); ; mpfr_pow(r25140, r25126, r25139, MPFR_RNDN); mpfr_div(r25141, r25133, r25140, MPFR_RNDN); mpfr_sub(r25142, r25138, r25141, MPFR_RNDN); mpfr_mul(r25143, r25126, r25126, MPFR_RNDN); mpfr_add(r25144, r25143, r25133, MPFR_RNDN); mpfr_div(r25145, r25126, r25144, MPFR_RNDN); if (mpfr_get_si(r25132, MPFR_RNDN)) { mpfr_set(r25146, r25142, MPFR_RNDN); } else { mpfr_set(r25146, r25145, MPFR_RNDN); }; return mpfr_get_d(r25146, MPFR_RNDN); } static mpfr_t r25147, r25148, r25149, r25150, r25151, r25152, r25153, r25154, r25155, r25156, r25157, r25158, r25159, r25160, r25161, r25162, r25163, r25164, r25165, r25166, r25167; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25147); mpfr_init_set_str(r25148, "-8.07298555250867e+25", 10, MPFR_RNDN); mpfr_init(r25149); mpfr_init_set_str(r25150, "1565.0996446725549", 10, MPFR_RNDN); mpfr_init(r25151); mpfr_init(r25152); mpfr_init(r25153); mpfr_init_set_str(r25154, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25155, "5", 10, MPFR_RNDN); mpfr_init(r25156); mpfr_init(r25157); mpfr_init(r25158); mpfr_init(r25159); mpfr_init_set_str(r25160, "3", 10, MPFR_RNDN); mpfr_init(r25161); mpfr_init(r25162); mpfr_init(r25163); mpfr_init(r25164); mpfr_init(r25165); mpfr_init(r25166); mpfr_init(r25167); } double f_dm(double x) { mpfr_set_d(r25147, x, MPFR_RNDN); ; mpfr_set_si(r25149, mpfr_cmp(r25147, r25148) <= 0, MPFR_RNDN); ; mpfr_set_si(r25151, mpfr_cmp(r25147, r25150) <= 0, MPFR_RNDN); mpfr_set_si(r25152, !mpfr_get_si(r25151, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r25153, mpfr_get_si(r25149, MPFR_RNDN) || mpfr_get_si(r25152, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r25156, r25147, r25155, MPFR_RNDN); mpfr_div(r25157, r25154, r25156, MPFR_RNDN); mpfr_div(r25158, r25154, r25147, MPFR_RNDN); mpfr_add(r25159, r25157, r25158, MPFR_RNDN); ; mpfr_pow(r25161, r25147, r25160, MPFR_RNDN); mpfr_div(r25162, r25154, r25161, MPFR_RNDN); mpfr_sub(r25163, r25159, r25162, MPFR_RNDN); mpfr_mul(r25164, r25147, r25147, MPFR_RNDN); mpfr_add(r25165, r25164, r25154, MPFR_RNDN); mpfr_div(r25166, r25147, r25165, MPFR_RNDN); if (mpfr_get_si(r25153, MPFR_RNDN)) { mpfr_set(r25167, r25163, MPFR_RNDN); } else { mpfr_set(r25167, r25166, MPFR_RNDN); }; return mpfr_get_d(r25167, MPFR_RNDN); }