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