#include #include #include #include #include char *name = "Jmat.Real.erfi, branch x less than or equal to 0.5"; double f_if(float x) { float r24128 = 1; float r24129 = atan2(1.0, 0.0); float r24130 = sqrt(r24129); float r24131 = r24128 / r24130; float r24132 = 2; float r24133 = x; float r24134 = fabs(r24133); float r24135 = r24132 * r24134; float r24136 = 3; float r24137 = r24132 / r24136; float r24138 = r24134 * r24134; float r24139 = r24138 * r24134; float r24140 = r24137 * r24139; float r24141 = r24135 + r24140; float r24142 = 5; float r24143 = r24128 / r24142; float r24144 = r24139 * r24134; float r24145 = r24144 * r24134; float r24146 = r24143 * r24145; float r24147 = r24141 + r24146; float r24148 = 21; float r24149 = r24128 / r24148; float r24150 = r24145 * r24134; float r24151 = r24150 * r24134; float r24152 = r24149 * r24151; float r24153 = r24147 + r24152; float r24154 = r24131 * r24153; float r24155 = fabs(r24154); return r24155; } double f_id(double x) { double r24156 = 1; double r24157 = atan2(1.0, 0.0); double r24158 = sqrt(r24157); double r24159 = r24156 / r24158; double r24160 = 2; double r24161 = x; double r24162 = fabs(r24161); double r24163 = r24160 * r24162; double r24164 = 3; double r24165 = r24160 / r24164; double r24166 = r24162 * r24162; double r24167 = r24166 * r24162; double r24168 = r24165 * r24167; double r24169 = r24163 + r24168; double r24170 = 5; double r24171 = r24156 / r24170; double r24172 = r24167 * r24162; double r24173 = r24172 * r24162; double r24174 = r24171 * r24173; double r24175 = r24169 + r24174; double r24176 = 21; double r24177 = r24156 / r24176; double r24178 = r24173 * r24162; double r24179 = r24178 * r24162; double r24180 = r24177 * r24179; double r24181 = r24175 + r24180; double r24182 = r24159 * r24181; double r24183 = fabs(r24182); return r24183; } double f_of(float x) { float r24184 = 1; float r24185 = atan2(1.0, 0.0); float r24186 = sqrt(r24185); float r24187 = r24184 / r24186; float r24188 = 2; float r24189 = x; float r24190 = fabs(r24189); float r24191 = r24188 * r24190; float r24192 = 3; float r24193 = r24188 / r24192; float r24194 = r24190 * r24190; float r24195 = r24194 * r24190; float r24196 = r24193 * r24195; float r24197 = r24191 + r24196; float r24198 = 5; float r24199 = r24184 / r24198; float r24200 = r24195 * r24190; float r24201 = r24200 * r24190; float r24202 = r24199 * r24201; float r24203 = r24197 + r24202; float r24204 = 21; float r24205 = r24190 / r24204; float r24206 = r24192 + r24192; float r24207 = pow(r24190, r24206); float r24208 = r24205 * r24207; float r24209 = pow(r24208, r24184); float r24210 = r24203 + r24209; float r24211 = r24187 * r24210; float r24212 = fabs(r24211); return r24212; } double f_od(double x) { double r24213 = 1; double r24214 = atan2(1.0, 0.0); double r24215 = sqrt(r24214); double r24216 = r24213 / r24215; double r24217 = 2; double r24218 = x; double r24219 = fabs(r24218); double r24220 = r24217 * r24219; double r24221 = 3; double r24222 = r24217 / r24221; double r24223 = r24219 * r24219; double r24224 = r24223 * r24219; double r24225 = r24222 * r24224; double r24226 = r24220 + r24225; double r24227 = 5; double r24228 = r24213 / r24227; double r24229 = r24224 * r24219; double r24230 = r24229 * r24219; double r24231 = r24228 * r24230; double r24232 = r24226 + r24231; double r24233 = 21; double r24234 = r24219 / r24233; double r24235 = r24221 + r24221; double r24236 = pow(r24219, r24235); double r24237 = r24234 * r24236; double r24238 = pow(r24237, r24213); double r24239 = r24232 + r24238; double r24240 = r24216 * r24239; double r24241 = fabs(r24240); return r24241; } 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 r24242, r24243, r24244, 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; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r24242, "1", 10, MPFR_RNDN); mpfr_init(r24243); mpfr_init(r24244); mpfr_init(r24245); mpfr_init_set_str(r24246, "2", 10, MPFR_RNDN); mpfr_init(r24247); mpfr_init(r24248); mpfr_init(r24249); mpfr_init_set_str(r24250, "3", 10, MPFR_RNDN); mpfr_init(r24251); mpfr_init(r24252); mpfr_init(r24253); mpfr_init(r24254); mpfr_init(r24255); mpfr_init_set_str(r24256, "5", 10, MPFR_RNDN); mpfr_init(r24257); mpfr_init(r24258); mpfr_init(r24259); mpfr_init(r24260); mpfr_init(r24261); mpfr_init_set_str(r24262, "21", 10, MPFR_RNDN); mpfr_init(r24263); mpfr_init(r24264); mpfr_init(r24265); mpfr_init(r24266); mpfr_init(r24267); mpfr_init(r24268); mpfr_init(r24269); } double f_im(double x) { ; mpfr_const_pi(r24243, MPFR_RNDN); mpfr_sqrt(r24244, r24243, MPFR_RNDN); mpfr_div(r24245, r24242, r24244, MPFR_RNDN); ; mpfr_set_d(r24247, x, MPFR_RNDN); mpfr_abs(r24248, r24247, MPFR_RNDN); mpfr_mul(r24249, r24246, r24248, MPFR_RNDN); ; mpfr_div(r24251, r24246, r24250, MPFR_RNDN); mpfr_mul(r24252, r24248, r24248, MPFR_RNDN); mpfr_mul(r24253, r24252, r24248, MPFR_RNDN); mpfr_mul(r24254, r24251, r24253, MPFR_RNDN); mpfr_add(r24255, r24249, r24254, MPFR_RNDN); ; mpfr_div(r24257, r24242, r24256, MPFR_RNDN); mpfr_mul(r24258, r24253, r24248, MPFR_RNDN); mpfr_mul(r24259, r24258, r24248, MPFR_RNDN); mpfr_mul(r24260, r24257, r24259, MPFR_RNDN); mpfr_add(r24261, r24255, r24260, MPFR_RNDN); ; mpfr_div(r24263, r24242, r24262, MPFR_RNDN); mpfr_mul(r24264, r24259, r24248, MPFR_RNDN); mpfr_mul(r24265, r24264, r24248, MPFR_RNDN); mpfr_mul(r24266, r24263, r24265, MPFR_RNDN); mpfr_add(r24267, r24261, r24266, MPFR_RNDN); mpfr_mul(r24268, r24245, r24267, MPFR_RNDN); mpfr_abs(r24269, r24268, MPFR_RNDN); return mpfr_get_d(r24269, MPFR_RNDN); } static mpfr_t r24270, r24271, r24272, r24273, r24274, r24275, r24276, r24277, r24278, r24279, r24280, r24281, r24282, r24283, r24284, r24285, r24286, r24287, r24288, r24289, r24290, r24291, r24292, r24293, r24294, r24295, r24296, r24297, r24298; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24270, "1", 10, MPFR_RNDN); mpfr_init(r24271); mpfr_init(r24272); mpfr_init(r24273); mpfr_init_set_str(r24274, "2", 10, MPFR_RNDN); mpfr_init(r24275); mpfr_init(r24276); mpfr_init(r24277); mpfr_init_set_str(r24278, "3", 10, MPFR_RNDN); mpfr_init(r24279); mpfr_init(r24280); mpfr_init(r24281); mpfr_init(r24282); mpfr_init(r24283); mpfr_init_set_str(r24284, "5", 10, MPFR_RNDN); mpfr_init(r24285); mpfr_init(r24286); mpfr_init(r24287); mpfr_init(r24288); mpfr_init(r24289); mpfr_init_set_str(r24290, "21", 10, MPFR_RNDN); mpfr_init(r24291); mpfr_init(r24292); mpfr_init(r24293); mpfr_init(r24294); mpfr_init(r24295); mpfr_init(r24296); mpfr_init(r24297); mpfr_init(r24298); } double f_fm(double x) { ; mpfr_const_pi(r24271, MPFR_RNDN); mpfr_sqrt(r24272, r24271, MPFR_RNDN); mpfr_div(r24273, r24270, r24272, MPFR_RNDN); ; mpfr_set_d(r24275, x, MPFR_RNDN); mpfr_abs(r24276, r24275, MPFR_RNDN); mpfr_mul(r24277, r24274, r24276, MPFR_RNDN); ; mpfr_div(r24279, r24274, r24278, MPFR_RNDN); mpfr_mul(r24280, r24276, r24276, MPFR_RNDN); mpfr_mul(r24281, r24280, r24276, MPFR_RNDN); mpfr_mul(r24282, r24279, r24281, MPFR_RNDN); mpfr_add(r24283, r24277, r24282, MPFR_RNDN); ; mpfr_div(r24285, r24270, r24284, MPFR_RNDN); mpfr_mul(r24286, r24281, r24276, MPFR_RNDN); mpfr_mul(r24287, r24286, r24276, MPFR_RNDN); mpfr_mul(r24288, r24285, r24287, MPFR_RNDN); mpfr_add(r24289, r24283, r24288, MPFR_RNDN); ; mpfr_div(r24291, r24276, r24290, MPFR_RNDN); mpfr_add(r24292, r24278, r24278, MPFR_RNDN); mpfr_pow(r24293, r24276, r24292, MPFR_RNDN); mpfr_mul(r24294, r24291, r24293, MPFR_RNDN); mpfr_pow(r24295, r24294, r24270, MPFR_RNDN); mpfr_add(r24296, r24289, r24295, MPFR_RNDN); mpfr_mul(r24297, r24273, r24296, MPFR_RNDN); mpfr_abs(r24298, r24297, MPFR_RNDN); return mpfr_get_d(r24298, MPFR_RNDN); } static mpfr_t r24299, r24300, r24301, r24302, r24303, r24304, r24305, r24306, r24307, r24308, r24309, r24310, r24311, r24312, r24313, r24314, r24315, r24316, r24317, r24318, r24319, r24320, r24321, r24322, r24323, r24324, r24325, r24326, r24327; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24299, "1", 10, MPFR_RNDN); mpfr_init(r24300); mpfr_init(r24301); mpfr_init(r24302); mpfr_init_set_str(r24303, "2", 10, MPFR_RNDN); mpfr_init(r24304); mpfr_init(r24305); mpfr_init(r24306); mpfr_init_set_str(r24307, "3", 10, MPFR_RNDN); mpfr_init(r24308); mpfr_init(r24309); mpfr_init(r24310); mpfr_init(r24311); mpfr_init(r24312); mpfr_init_set_str(r24313, "5", 10, MPFR_RNDN); mpfr_init(r24314); mpfr_init(r24315); mpfr_init(r24316); mpfr_init(r24317); mpfr_init(r24318); mpfr_init_set_str(r24319, "21", 10, MPFR_RNDN); mpfr_init(r24320); mpfr_init(r24321); mpfr_init(r24322); mpfr_init(r24323); mpfr_init(r24324); mpfr_init(r24325); mpfr_init(r24326); mpfr_init(r24327); } double f_dm(double x) { ; mpfr_const_pi(r24300, MPFR_RNDN); mpfr_sqrt(r24301, r24300, MPFR_RNDN); mpfr_div(r24302, r24299, r24301, MPFR_RNDN); ; mpfr_set_d(r24304, x, MPFR_RNDN); mpfr_abs(r24305, r24304, MPFR_RNDN); mpfr_mul(r24306, r24303, r24305, MPFR_RNDN); ; mpfr_div(r24308, r24303, r24307, MPFR_RNDN); mpfr_mul(r24309, r24305, r24305, MPFR_RNDN); mpfr_mul(r24310, r24309, r24305, MPFR_RNDN); mpfr_mul(r24311, r24308, r24310, MPFR_RNDN); mpfr_add(r24312, r24306, r24311, MPFR_RNDN); ; mpfr_div(r24314, r24299, r24313, MPFR_RNDN); mpfr_mul(r24315, r24310, r24305, MPFR_RNDN); mpfr_mul(r24316, r24315, r24305, MPFR_RNDN); mpfr_mul(r24317, r24314, r24316, MPFR_RNDN); mpfr_add(r24318, r24312, r24317, MPFR_RNDN); ; mpfr_div(r24320, r24305, r24319, MPFR_RNDN); mpfr_add(r24321, r24307, r24307, MPFR_RNDN); mpfr_pow(r24322, r24305, r24321, MPFR_RNDN); mpfr_mul(r24323, r24320, r24322, MPFR_RNDN); mpfr_pow(r24324, r24323, r24299, MPFR_RNDN); mpfr_add(r24325, r24318, r24324, MPFR_RNDN); mpfr_mul(r24326, r24302, r24325, MPFR_RNDN); mpfr_abs(r24327, r24326, MPFR_RNDN); return mpfr_get_d(r24327, MPFR_RNDN); }