#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r17108 = a; float r17109 = r17108 * r17108; float r17110 = b; float r17111 = r17110 * r17110; float r17112 = r17109 + r17111; float r17113 = r17112 * r17112; float r17114 = 4.0f; float r17115 = 1.0f; float r17116 = r17115 + r17108; float r17117 = r17109 * r17116; float r17118 = 3.0f; float r17119 = r17118 * r17108; float r17120 = r17115 - r17119; float r17121 = r17111 * r17120; float r17122 = r17117 + r17121; float r17123 = r17114 * r17122; float r17124 = r17113 + r17123; float r17125 = r17124 - r17115; return r17125; } double f_id(double a, double b) { double r17126 = a; double r17127 = r17126 * r17126; double r17128 = b; double r17129 = r17128 * r17128; double r17130 = r17127 + r17129; double r17131 = r17130 * r17130; double r17132 = 4.0; double r17133 = 1.0; double r17134 = r17133 + r17126; double r17135 = r17127 * r17134; double r17136 = 3.0; double r17137 = r17136 * r17126; double r17138 = r17133 - r17137; double r17139 = r17129 * r17138; double r17140 = r17135 + r17139; double r17141 = r17132 * r17140; double r17142 = r17131 + r17141; double r17143 = r17142 - r17133; return r17143; } double f_of(float a, float b) { float r17144 = a; float r17145 = r17144 * r17144; float r17146 = b; float r17147 = r17146 * r17146; float r17148 = r17145 + r17147; float r17149 = r17148 * r17148; float r17150 = 4.0f; float r17151 = 1.0f; float r17152 = r17151 + r17144; float r17153 = r17145 * r17152; float r17154 = 3.0f; float r17155 = r17154 * r17144; float r17156 = r17151 - r17155; float r17157 = r17147 * r17156; float r17158 = r17153 + r17157; float r17159 = pow(r17158, r17151); float r17160 = r17150 * r17159; float r17161 = r17149 + r17160; float r17162 = r17161 - r17151; return r17162; } double f_od(double a, double b) { double r17163 = a; double r17164 = r17163 * r17163; double r17165 = b; double r17166 = r17165 * r17165; double r17167 = r17164 + r17166; double r17168 = r17167 * r17167; double r17169 = 4.0; double r17170 = 1.0; double r17171 = r17170 + r17163; double r17172 = r17164 * r17171; double r17173 = 3.0; double r17174 = r17173 * r17163; double r17175 = r17170 - r17174; double r17176 = r17166 * r17175; double r17177 = r17172 + r17176; double r17178 = pow(r17177, r17170); double r17179 = r17169 * r17178; double r17180 = r17168 + r17179; double r17181 = r17180 - r17170; return r17181; } 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 r17182, r17183, r17184, r17185, r17186, r17187, r17188, r17189, r17190, r17191, r17192, r17193, r17194, r17195, r17196, r17197, r17198, r17199; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17182); mpfr_init(r17183); mpfr_init(r17184); mpfr_init(r17185); mpfr_init(r17186); mpfr_init(r17187); mpfr_init_set_str(r17188, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17189, "1", 10, MPFR_RNDN); mpfr_init(r17190); mpfr_init(r17191); mpfr_init_set_str(r17192, "3", 10, MPFR_RNDN); mpfr_init(r17193); mpfr_init(r17194); mpfr_init(r17195); mpfr_init(r17196); mpfr_init(r17197); mpfr_init(r17198); mpfr_init(r17199); } double f_im(double a, double b) { mpfr_set_d(r17182, a, MPFR_RNDN); mpfr_sqr(r17183, r17182, MPFR_RNDN); mpfr_set_d(r17184, b, MPFR_RNDN); mpfr_sqr(r17185, r17184, MPFR_RNDN); mpfr_add(r17186, r17183, r17185, MPFR_RNDN); mpfr_sqr(r17187, r17186, MPFR_RNDN); ; ; mpfr_add(r17190, r17189, r17182, MPFR_RNDN); mpfr_mul(r17191, r17183, r17190, MPFR_RNDN); ; mpfr_mul(r17193, r17192, r17182, MPFR_RNDN); mpfr_sub(r17194, r17189, r17193, MPFR_RNDN); mpfr_mul(r17195, r17185, r17194, MPFR_RNDN); mpfr_add(r17196, r17191, r17195, MPFR_RNDN); mpfr_mul(r17197, r17188, r17196, MPFR_RNDN); mpfr_add(r17198, r17187, r17197, MPFR_RNDN); mpfr_sub(r17199, r17198, r17189, MPFR_RNDN); return mpfr_get_d(r17199, MPFR_RNDN); } static mpfr_t r17200, r17201, r17202, r17203, r17204, r17205, r17206, r17207, r17208, r17209, r17210, r17211, r17212, r17213, r17214, r17215, r17216, r17217, r17218; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17200); mpfr_init(r17201); mpfr_init(r17202); mpfr_init(r17203); mpfr_init(r17204); mpfr_init(r17205); mpfr_init_set_str(r17206, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17207, "1", 10, MPFR_RNDN); mpfr_init(r17208); mpfr_init(r17209); mpfr_init_set_str(r17210, "3", 10, MPFR_RNDN); mpfr_init(r17211); mpfr_init(r17212); mpfr_init(r17213); mpfr_init(r17214); mpfr_init(r17215); mpfr_init(r17216); mpfr_init(r17217); mpfr_init(r17218); } double f_fm(double a, double b) { mpfr_set_d(r17200, a, MPFR_RNDN); mpfr_sqr(r17201, r17200, MPFR_RNDN); mpfr_set_d(r17202, b, MPFR_RNDN); mpfr_sqr(r17203, r17202, MPFR_RNDN); mpfr_add(r17204, r17201, r17203, MPFR_RNDN); mpfr_sqr(r17205, r17204, MPFR_RNDN); ; ; mpfr_add(r17208, r17207, r17200, MPFR_RNDN); mpfr_mul(r17209, r17201, r17208, MPFR_RNDN); ; mpfr_mul(r17211, r17210, r17200, MPFR_RNDN); mpfr_sub(r17212, r17207, r17211, MPFR_RNDN); mpfr_mul(r17213, r17203, r17212, MPFR_RNDN); mpfr_add(r17214, r17209, r17213, MPFR_RNDN); mpfr_pow(r17215, r17214, r17207, MPFR_RNDN); mpfr_mul(r17216, r17206, r17215, MPFR_RNDN); mpfr_add(r17217, r17205, r17216, MPFR_RNDN); mpfr_sub(r17218, r17217, r17207, MPFR_RNDN); return mpfr_get_d(r17218, MPFR_RNDN); } static mpfr_t r17219, r17220, r17221, r17222, r17223, r17224, r17225, r17226, r17227, r17228, r17229, r17230, r17231, r17232, r17233, r17234, r17235, r17236, r17237; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17219); mpfr_init(r17220); mpfr_init(r17221); mpfr_init(r17222); mpfr_init(r17223); mpfr_init(r17224); mpfr_init_set_str(r17225, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17226, "1", 10, MPFR_RNDN); mpfr_init(r17227); mpfr_init(r17228); mpfr_init_set_str(r17229, "3", 10, MPFR_RNDN); mpfr_init(r17230); mpfr_init(r17231); mpfr_init(r17232); mpfr_init(r17233); mpfr_init(r17234); mpfr_init(r17235); mpfr_init(r17236); mpfr_init(r17237); } double f_dm(double a, double b) { mpfr_set_d(r17219, a, MPFR_RNDN); mpfr_sqr(r17220, r17219, MPFR_RNDN); mpfr_set_d(r17221, b, MPFR_RNDN); mpfr_sqr(r17222, r17221, MPFR_RNDN); mpfr_add(r17223, r17220, r17222, MPFR_RNDN); mpfr_sqr(r17224, r17223, MPFR_RNDN); ; ; mpfr_add(r17227, r17226, r17219, MPFR_RNDN); mpfr_mul(r17228, r17220, r17227, MPFR_RNDN); ; mpfr_mul(r17230, r17229, r17219, MPFR_RNDN); mpfr_sub(r17231, r17226, r17230, MPFR_RNDN); mpfr_mul(r17232, r17222, r17231, MPFR_RNDN); mpfr_add(r17233, r17228, r17232, MPFR_RNDN); mpfr_pow(r17234, r17233, r17226, MPFR_RNDN); mpfr_mul(r17235, r17225, r17234, MPFR_RNDN); mpfr_add(r17236, r17224, r17235, MPFR_RNDN); mpfr_sub(r17237, r17236, r17226, MPFR_RNDN); return mpfr_get_d(r17237, MPFR_RNDN); }