#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (20)"; double f_if(float f) { float r24119 = 1; float r24120 = atan2(1.0, 0.0); float r24121 = 4; float r24122 = r24120 / r24121; float r24123 = r24119 / r24122; float r24124 = f; float r24125 = r24122 * r24124; float r24126 = exp(r24125); float r24127 = -r24125; float r24128 = exp(r24127); float r24129 = r24126 + r24128; float r24130 = r24126 - r24128; float r24131 = r24129 / r24130; float r24132 = log(r24131); float r24133 = r24123 * r24132; float r24134 = -r24133; return r24134; } double f_id(double f) { double r24135 = 1; double r24136 = atan2(1.0, 0.0); double r24137 = 4; double r24138 = r24136 / r24137; double r24139 = r24135 / r24138; double r24140 = f; double r24141 = r24138 * r24140; double r24142 = exp(r24141); double r24143 = -r24141; double r24144 = exp(r24143); double r24145 = r24142 + r24144; double r24146 = r24142 - r24144; double r24147 = r24145 / r24146; double r24148 = log(r24147); double r24149 = r24139 * r24148; double r24150 = -r24149; return r24150; } double f_of(float f) { float r24151 = 7/5760; float r24152 = atan2(1.0, 0.0); float r24153 = r24152 * r24152; float r24154 = r24151 * r24153; float r24155 = -r24152; float r24156 = f; float r24157 = 4; float r24158 = pow(r24156, r24157); float r24159 = r24155 * r24158; float r24160 = 1/12; float r24161 = r24156 * r24160; float r24162 = r24156 * r24152; float r24163 = r24161 * r24162; float r24164 = fma(r24154, r24159, r24163); float r24165 = r24157 / r24152; float r24166 = log(r24156); float r24167 = log(r24157); float r24168 = r24166 - r24167; float r24169 = r24165 * r24168; float r24170 = log(r24152); float r24171 = r24165 * r24170; float r24172 = r24169 + r24171; float r24173 = r24164 - r24172; float r24174 = 0; float r24175 = r24173 + r24174; float r24176 = -r24175; return r24176; } double f_od(double f) { double r24177 = 7/5760; double r24178 = atan2(1.0, 0.0); double r24179 = r24178 * r24178; double r24180 = r24177 * r24179; double r24181 = -r24178; double r24182 = f; double r24183 = 4; double r24184 = pow(r24182, r24183); double r24185 = r24181 * r24184; double r24186 = 1/12; double r24187 = r24182 * r24186; double r24188 = r24182 * r24178; double r24189 = r24187 * r24188; double r24190 = fma(r24180, r24185, r24189); double r24191 = r24183 / r24178; double r24192 = log(r24182); double r24193 = log(r24183); double r24194 = r24192 - r24193; double r24195 = r24191 * r24194; double r24196 = log(r24178); double r24197 = r24191 * r24196; double r24198 = r24195 + r24197; double r24199 = r24190 - r24198; double r24200 = 0; double r24201 = r24199 + r24200; double r24202 = -r24201; return r24202; } 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 r24203, r24204, r24205, r24206, r24207, r24208, r24209, r24210, r24211, r24212, r24213, r24214, r24215, r24216, r24217, r24218; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r24203, "1", 10, MPFR_RNDN); mpfr_init(r24204); mpfr_init_set_str(r24205, "4", 10, MPFR_RNDN); mpfr_init(r24206); mpfr_init(r24207); mpfr_init(r24208); mpfr_init(r24209); mpfr_init(r24210); mpfr_init(r24211); mpfr_init(r24212); mpfr_init(r24213); mpfr_init(r24214); mpfr_init(r24215); mpfr_init(r24216); mpfr_init(r24217); mpfr_init(r24218); } double f_im(double f) { ; mpfr_const_pi(r24204, MPFR_RNDN); ; mpfr_div(r24206, r24204, r24205, MPFR_RNDN); mpfr_div(r24207, r24203, r24206, MPFR_RNDN); mpfr_set_d(r24208, f, MPFR_RNDN); mpfr_mul(r24209, r24206, r24208, MPFR_RNDN); mpfr_exp(r24210, r24209, MPFR_RNDN); mpfr_neg(r24211, r24209, MPFR_RNDN); mpfr_exp(r24212, r24211, MPFR_RNDN); mpfr_add(r24213, r24210, r24212, MPFR_RNDN); mpfr_sub(r24214, r24210, r24212, MPFR_RNDN); mpfr_div(r24215, r24213, r24214, MPFR_RNDN); mpfr_log(r24216, r24215, MPFR_RNDN); mpfr_mul(r24217, r24207, r24216, MPFR_RNDN); mpfr_neg(r24218, r24217, MPFR_RNDN); return mpfr_get_d(r24218, MPFR_RNDN); } static mpfr_t r24219, r24220, r24221, r24222, r24223, r24224, r24225, r24226, r24227, r24228, r24229, r24230, r24231, r24232, r24233, r24234, r24235, r24236, r24237, r24238, r24239, r24240, r24241, r24242, r24243, r24244; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r24219, "7/5760", 10, MPFR_RNDN); mpfr_init(r24220); mpfr_init(r24221); mpfr_init(r24222); mpfr_init(r24223); mpfr_init(r24224); mpfr_init_set_str(r24225, "4", 10, MPFR_RNDN); mpfr_init(r24226); mpfr_init(r24227); mpfr_init_set_str(r24228, "1/12", 10, MPFR_RNDN); mpfr_init(r24229); mpfr_init(r24230); mpfr_init(r24231); mpfr_init(r24232); mpfr_init(r24233); mpfr_init(r24234); mpfr_init(r24235); mpfr_init(r24236); mpfr_init(r24237); mpfr_init(r24238); mpfr_init(r24239); mpfr_init(r24240); mpfr_init(r24241); mpfr_init_set_str(r24242, "0", 10, MPFR_RNDN); mpfr_init(r24243); mpfr_init(r24244); } double f_fm(double f) { ; mpfr_const_pi(r24220, MPFR_RNDN); mpfr_mul(r24221, r24220, r24220, MPFR_RNDN); mpfr_mul(r24222, r24219, r24221, MPFR_RNDN); mpfr_neg(r24223, r24220, MPFR_RNDN); mpfr_set_d(r24224, f, MPFR_RNDN); ; mpfr_pow(r24226, r24224, r24225, MPFR_RNDN); mpfr_mul(r24227, r24223, r24226, MPFR_RNDN); ; mpfr_mul(r24229, r24224, r24228, MPFR_RNDN); mpfr_mul(r24230, r24224, r24220, MPFR_RNDN); mpfr_mul(r24231, r24229, r24230, MPFR_RNDN); mpfr_fma(r24232, r24222, r24227, r24231, MPFR_RNDN); mpfr_div(r24233, r24225, r24220, MPFR_RNDN); mpfr_log(r24234, r24224, MPFR_RNDN); mpfr_log(r24235, r24225, MPFR_RNDN); mpfr_sub(r24236, r24234, r24235, MPFR_RNDN); mpfr_mul(r24237, r24233, r24236, MPFR_RNDN); mpfr_log(r24238, r24220, MPFR_RNDN); mpfr_mul(r24239, r24233, r24238, MPFR_RNDN); mpfr_add(r24240, r24237, r24239, MPFR_RNDN); mpfr_sub(r24241, r24232, r24240, MPFR_RNDN); ; mpfr_add(r24243, r24241, r24242, MPFR_RNDN); mpfr_neg(r24244, r24243, MPFR_RNDN); return mpfr_get_d(r24244, MPFR_RNDN); } static mpfr_t r24245, r24246, r24247, r24248, r24249, r24250, r24251, r24252, r24253, r24254, r24255, r24256, r24257, r24258, r24259, r24260, r24261, r24262, r24263, r24264, r24265, r24266, r24267, r24268, r24269, r24270; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r24245, "7/5760", 10, MPFR_RNDN); mpfr_init(r24246); mpfr_init(r24247); mpfr_init(r24248); mpfr_init(r24249); mpfr_init(r24250); mpfr_init_set_str(r24251, "4", 10, MPFR_RNDN); mpfr_init(r24252); mpfr_init(r24253); mpfr_init_set_str(r24254, "1/12", 10, MPFR_RNDN); mpfr_init(r24255); mpfr_init(r24256); mpfr_init(r24257); mpfr_init(r24258); mpfr_init(r24259); mpfr_init(r24260); mpfr_init(r24261); mpfr_init(r24262); mpfr_init(r24263); mpfr_init(r24264); mpfr_init(r24265); mpfr_init(r24266); mpfr_init(r24267); mpfr_init_set_str(r24268, "0", 10, MPFR_RNDN); mpfr_init(r24269); mpfr_init(r24270); } double f_dm(double f) { ; mpfr_const_pi(r24246, MPFR_RNDN); mpfr_mul(r24247, r24246, r24246, MPFR_RNDN); mpfr_mul(r24248, r24245, r24247, MPFR_RNDN); mpfr_neg(r24249, r24246, MPFR_RNDN); mpfr_set_d(r24250, f, MPFR_RNDN); ; mpfr_pow(r24252, r24250, r24251, MPFR_RNDN); mpfr_mul(r24253, r24249, r24252, MPFR_RNDN); ; mpfr_mul(r24255, r24250, r24254, MPFR_RNDN); mpfr_mul(r24256, r24250, r24246, MPFR_RNDN); mpfr_mul(r24257, r24255, r24256, MPFR_RNDN); mpfr_fma(r24258, r24248, r24253, r24257, MPFR_RNDN); mpfr_div(r24259, r24251, r24246, MPFR_RNDN); mpfr_log(r24260, r24250, MPFR_RNDN); mpfr_log(r24261, r24251, MPFR_RNDN); mpfr_sub(r24262, r24260, r24261, MPFR_RNDN); mpfr_mul(r24263, r24259, r24262, MPFR_RNDN); mpfr_log(r24264, r24246, MPFR_RNDN); mpfr_mul(r24265, r24259, r24264, MPFR_RNDN); mpfr_add(r24266, r24263, r24265, MPFR_RNDN); mpfr_sub(r24267, r24258, r24266, MPFR_RNDN); ; mpfr_add(r24269, r24267, r24268, MPFR_RNDN); mpfr_neg(r24270, r24269, MPFR_RNDN); return mpfr_get_d(r24270, MPFR_RNDN); }