#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (26)"; double f_if(float a, float b) { float r25060 = a; float r25061 = r25060 * r25060; float r25062 = b; float r25063 = r25062 * r25062; float r25064 = r25061 + r25063; float r25065 = 2; float r25066 = pow(r25064, r25065); float r25067 = 4; float r25068 = r25067 * r25063; float r25069 = r25066 + r25068; float r25070 = 1; float r25071 = r25069 - r25070; return r25071; } double f_id(double a, double b) { double r25072 = a; double r25073 = r25072 * r25072; double r25074 = b; double r25075 = r25074 * r25074; double r25076 = r25073 + r25075; double r25077 = 2; double r25078 = pow(r25076, r25077); double r25079 = 4; double r25080 = r25079 * r25075; double r25081 = r25078 + r25080; double r25082 = 1; double r25083 = r25081 - r25082; return r25083; } double f_of(float a, float b) { float r25084 = b; float r25085 = r25084 * r25084; float r25086 = a; float r25087 = r25086 * r25086; float r25088 = r25085 + r25087; float r25089 = cbrt(r25088); float r25090 = r25088 * r25089; float r25091 = r25087 + r25085; float r25092 = cbrt(r25091); float r25093 = 2; float r25094 = pow(r25092, r25093); float r25095 = r25090 * r25094; float r25096 = 4; float r25097 = r25096 * r25085; float r25098 = r25095 + r25097; float r25099 = 1; float r25100 = r25098 - r25099; return r25100; } double f_od(double a, double b) { double r25101 = b; double r25102 = r25101 * r25101; double r25103 = a; double r25104 = r25103 * r25103; double r25105 = r25102 + r25104; double r25106 = cbrt(r25105); double r25107 = r25105 * r25106; double r25108 = r25104 + r25102; double r25109 = cbrt(r25108); double r25110 = 2; double r25111 = pow(r25109, r25110); double r25112 = r25107 * r25111; double r25113 = 4; double r25114 = r25113 * r25102; double r25115 = r25112 + r25114; double r25116 = 1; double r25117 = r25115 - r25116; return r25117; } 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 r25118, r25119, r25120, r25121, r25122, r25123, r25124, r25125, r25126, r25127, r25128, r25129; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r25118); mpfr_init(r25119); mpfr_init(r25120); mpfr_init(r25121); mpfr_init(r25122); mpfr_init_set_str(r25123, "2", 10, MPFR_RNDN); mpfr_init(r25124); mpfr_init_set_str(r25125, "4", 10, MPFR_RNDN); mpfr_init(r25126); mpfr_init(r25127); mpfr_init_set_str(r25128, "1", 10, MPFR_RNDN); mpfr_init(r25129); } double f_im(double a, double b) { mpfr_set_d(r25118, a, MPFR_RNDN); mpfr_mul(r25119, r25118, r25118, MPFR_RNDN); mpfr_set_d(r25120, b, MPFR_RNDN); mpfr_mul(r25121, r25120, r25120, MPFR_RNDN); mpfr_add(r25122, r25119, r25121, MPFR_RNDN); ; mpfr_pow(r25124, r25122, r25123, MPFR_RNDN); ; mpfr_mul(r25126, r25125, r25121, MPFR_RNDN); mpfr_add(r25127, r25124, r25126, MPFR_RNDN); ; mpfr_sub(r25129, r25127, r25128, MPFR_RNDN); return mpfr_get_d(r25129, MPFR_RNDN); } static mpfr_t 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(400); mpfr_init(r25130); mpfr_init(r25131); mpfr_init(r25132); mpfr_init(r25133); mpfr_init(r25134); mpfr_init(r25135); mpfr_init(r25136); mpfr_init(r25137); mpfr_init(r25138); mpfr_init_set_str(r25139, "2", 10, MPFR_RNDN); mpfr_init(r25140); mpfr_init(r25141); mpfr_init_set_str(r25142, "4", 10, MPFR_RNDN); mpfr_init(r25143); mpfr_init(r25144); mpfr_init_set_str(r25145, "1", 10, MPFR_RNDN); mpfr_init(r25146); } double f_fm(double a, double b) { mpfr_set_d(r25130, b, MPFR_RNDN); mpfr_mul(r25131, r25130, r25130, MPFR_RNDN); mpfr_set_d(r25132, a, MPFR_RNDN); mpfr_mul(r25133, r25132, r25132, MPFR_RNDN); mpfr_add(r25134, r25131, r25133, MPFR_RNDN); mpfr_cbrt(r25135, r25134, MPFR_RNDN); mpfr_mul(r25136, r25134, r25135, MPFR_RNDN); mpfr_add(r25137, r25133, r25131, MPFR_RNDN); mpfr_cbrt(r25138, r25137, MPFR_RNDN); ; mpfr_pow(r25140, r25138, r25139, MPFR_RNDN); mpfr_mul(r25141, r25136, r25140, MPFR_RNDN); ; mpfr_mul(r25143, r25142, r25131, MPFR_RNDN); mpfr_add(r25144, r25141, r25143, MPFR_RNDN); ; mpfr_sub(r25146, r25144, 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; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r25147); 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_set_str(r25156, "2", 10, MPFR_RNDN); mpfr_init(r25157); mpfr_init(r25158); mpfr_init_set_str(r25159, "4", 10, MPFR_RNDN); mpfr_init(r25160); mpfr_init(r25161); mpfr_init_set_str(r25162, "1", 10, MPFR_RNDN); mpfr_init(r25163); } double f_dm(double a, double b) { mpfr_set_d(r25147, b, MPFR_RNDN); mpfr_mul(r25148, r25147, r25147, MPFR_RNDN); mpfr_set_d(r25149, a, MPFR_RNDN); mpfr_mul(r25150, r25149, r25149, MPFR_RNDN); mpfr_add(r25151, r25148, r25150, MPFR_RNDN); mpfr_cbrt(r25152, r25151, MPFR_RNDN); mpfr_mul(r25153, r25151, r25152, MPFR_RNDN); mpfr_add(r25154, r25150, r25148, MPFR_RNDN); mpfr_cbrt(r25155, r25154, MPFR_RNDN); ; mpfr_pow(r25157, r25155, r25156, MPFR_RNDN); mpfr_mul(r25158, r25153, r25157, MPFR_RNDN); ; mpfr_mul(r25160, r25159, r25148, MPFR_RNDN); mpfr_add(r25161, r25158, r25160, MPFR_RNDN); ; mpfr_sub(r25163, r25161, r25162, MPFR_RNDN); return mpfr_get_d(r25163, MPFR_RNDN); }