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