#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 r27195 = 1; float r27196 = atan2(1.0, 0.0); float r27197 = sqrt(r27196); float r27198 = r27195 / r27197; float r27199 = 2; float r27200 = x; float r27201 = fabs(r27200); float r27202 = r27199 * r27201; float r27203 = 3; float r27204 = r27199 / r27203; float r27205 = r27201 * r27201; float r27206 = r27205 * r27201; float r27207 = r27204 * r27206; float r27208 = r27202 + r27207; float r27209 = 5; float r27210 = r27195 / r27209; float r27211 = r27206 * r27201; float r27212 = r27211 * r27201; float r27213 = r27210 * r27212; float r27214 = r27208 + r27213; float r27215 = 21; float r27216 = r27195 / r27215; float r27217 = r27212 * r27201; float r27218 = r27217 * r27201; float r27219 = r27216 * r27218; float r27220 = r27214 + r27219; float r27221 = r27198 * r27220; float r27222 = fabs(r27221); return r27222; } double f_id(double x) { double r27223 = 1; double r27224 = atan2(1.0, 0.0); double r27225 = sqrt(r27224); double r27226 = r27223 / r27225; double r27227 = 2; double r27228 = x; double r27229 = fabs(r27228); double r27230 = r27227 * r27229; double r27231 = 3; double r27232 = r27227 / r27231; double r27233 = r27229 * r27229; double r27234 = r27233 * r27229; double r27235 = r27232 * r27234; double r27236 = r27230 + r27235; double r27237 = 5; double r27238 = r27223 / r27237; double r27239 = r27234 * r27229; double r27240 = r27239 * r27229; double r27241 = r27238 * r27240; double r27242 = r27236 + r27241; double r27243 = 21; double r27244 = r27223 / r27243; double r27245 = r27240 * r27229; double r27246 = r27245 * r27229; double r27247 = r27244 * r27246; double r27248 = r27242 + r27247; double r27249 = r27226 * r27248; double r27250 = fabs(r27249); return r27250; } double f_of(float x) { float r27251 = 1; float r27252 = atan2(1.0, 0.0); float r27253 = sqrt(r27252); float r27254 = r27251 / r27253; float r27255 = 2; float r27256 = x; float r27257 = fabs(r27256); float r27258 = r27255 * r27257; float r27259 = 3; float r27260 = r27255 / r27259; float r27261 = r27257 * r27257; float r27262 = r27261 * r27257; float r27263 = r27260 * r27262; float r27264 = r27258 + r27263; float r27265 = 5; float r27266 = r27251 / r27265; float r27267 = r27262 * r27257; float r27268 = r27267 * r27257; float r27269 = r27266 * r27268; float r27270 = r27264 + r27269; float r27271 = 21; float r27272 = r27257 / r27271; float r27273 = pow(r27257, r27259); float r27274 = r27273 * r27273; float r27275 = r27272 * r27274; float r27276 = pow(r27275, r27251); float r27277 = r27270 + r27276; float r27278 = r27254 * r27277; float r27279 = fabs(r27278); return r27279; } double f_od(double x) { double r27280 = 1; double r27281 = atan2(1.0, 0.0); double r27282 = sqrt(r27281); double r27283 = r27280 / r27282; double r27284 = 2; double r27285 = x; double r27286 = fabs(r27285); double r27287 = r27284 * r27286; double r27288 = 3; double r27289 = r27284 / r27288; double r27290 = r27286 * r27286; double r27291 = r27290 * r27286; double r27292 = r27289 * r27291; double r27293 = r27287 + r27292; double r27294 = 5; double r27295 = r27280 / r27294; double r27296 = r27291 * r27286; double r27297 = r27296 * r27286; double r27298 = r27295 * r27297; double r27299 = r27293 + r27298; double r27300 = 21; double r27301 = r27286 / r27300; double r27302 = pow(r27286, r27288); double r27303 = r27302 * r27302; double r27304 = r27301 * r27303; double r27305 = pow(r27304, r27280); double r27306 = r27299 + r27305; double r27307 = r27283 * r27306; double r27308 = fabs(r27307); return r27308; } 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 r27309, r27310, r27311, r27312, r27313, r27314, r27315, r27316, r27317, r27318, r27319, r27320, r27321, r27322, r27323, r27324, r27325, r27326, r27327, r27328, r27329, r27330, r27331, r27332, r27333, r27334, r27335, r27336; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r27309, "1", 10, MPFR_RNDN); mpfr_init(r27310); mpfr_init(r27311); mpfr_init(r27312); mpfr_init_set_str(r27313, "2", 10, MPFR_RNDN); mpfr_init(r27314); mpfr_init(r27315); mpfr_init(r27316); mpfr_init_set_str(r27317, "3", 10, MPFR_RNDN); mpfr_init(r27318); mpfr_init(r27319); mpfr_init(r27320); mpfr_init(r27321); mpfr_init(r27322); mpfr_init_set_str(r27323, "5", 10, MPFR_RNDN); mpfr_init(r27324); mpfr_init(r27325); mpfr_init(r27326); mpfr_init(r27327); mpfr_init(r27328); mpfr_init_set_str(r27329, "21", 10, MPFR_RNDN); mpfr_init(r27330); mpfr_init(r27331); mpfr_init(r27332); mpfr_init(r27333); mpfr_init(r27334); mpfr_init(r27335); mpfr_init(r27336); } double f_im(double x) { ; mpfr_const_pi(r27310, MPFR_RNDN); mpfr_sqrt(r27311, r27310, MPFR_RNDN); mpfr_div(r27312, r27309, r27311, MPFR_RNDN); ; mpfr_set_d(r27314, x, MPFR_RNDN); mpfr_abs(r27315, r27314, MPFR_RNDN); mpfr_mul(r27316, r27313, r27315, MPFR_RNDN); ; mpfr_div(r27318, r27313, r27317, MPFR_RNDN); mpfr_mul(r27319, r27315, r27315, MPFR_RNDN); mpfr_mul(r27320, r27319, r27315, MPFR_RNDN); mpfr_mul(r27321, r27318, r27320, MPFR_RNDN); mpfr_add(r27322, r27316, r27321, MPFR_RNDN); ; mpfr_div(r27324, r27309, r27323, MPFR_RNDN); mpfr_mul(r27325, r27320, r27315, MPFR_RNDN); mpfr_mul(r27326, r27325, r27315, MPFR_RNDN); mpfr_mul(r27327, r27324, r27326, MPFR_RNDN); mpfr_add(r27328, r27322, r27327, MPFR_RNDN); ; mpfr_div(r27330, r27309, r27329, MPFR_RNDN); mpfr_mul(r27331, r27326, r27315, MPFR_RNDN); mpfr_mul(r27332, r27331, r27315, MPFR_RNDN); mpfr_mul(r27333, r27330, r27332, MPFR_RNDN); mpfr_add(r27334, r27328, r27333, MPFR_RNDN); mpfr_mul(r27335, r27312, r27334, MPFR_RNDN); mpfr_abs(r27336, r27335, MPFR_RNDN); return mpfr_get_d(r27336, MPFR_RNDN); } static mpfr_t r27337, r27338, r27339, r27340, r27341, r27342, r27343, r27344, r27345, r27346, r27347, r27348, r27349, r27350, r27351, r27352, r27353, r27354, r27355, r27356, r27357, r27358, r27359, r27360, r27361, r27362, r27363, r27364, r27365; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27337, "1", 10, MPFR_RNDN); mpfr_init(r27338); mpfr_init(r27339); mpfr_init(r27340); mpfr_init_set_str(r27341, "2", 10, MPFR_RNDN); mpfr_init(r27342); mpfr_init(r27343); mpfr_init(r27344); mpfr_init_set_str(r27345, "3", 10, MPFR_RNDN); mpfr_init(r27346); mpfr_init(r27347); mpfr_init(r27348); mpfr_init(r27349); mpfr_init(r27350); mpfr_init_set_str(r27351, "5", 10, MPFR_RNDN); mpfr_init(r27352); mpfr_init(r27353); mpfr_init(r27354); mpfr_init(r27355); mpfr_init(r27356); mpfr_init_set_str(r27357, "21", 10, MPFR_RNDN); mpfr_init(r27358); mpfr_init(r27359); mpfr_init(r27360); mpfr_init(r27361); mpfr_init(r27362); mpfr_init(r27363); mpfr_init(r27364); mpfr_init(r27365); } double f_fm(double x) { ; mpfr_const_pi(r27338, MPFR_RNDN); mpfr_sqrt(r27339, r27338, MPFR_RNDN); mpfr_div(r27340, r27337, r27339, MPFR_RNDN); ; mpfr_set_d(r27342, x, MPFR_RNDN); mpfr_abs(r27343, r27342, MPFR_RNDN); mpfr_mul(r27344, r27341, r27343, MPFR_RNDN); ; mpfr_div(r27346, r27341, r27345, MPFR_RNDN); mpfr_mul(r27347, r27343, r27343, MPFR_RNDN); mpfr_mul(r27348, r27347, r27343, MPFR_RNDN); mpfr_mul(r27349, r27346, r27348, MPFR_RNDN); mpfr_add(r27350, r27344, r27349, MPFR_RNDN); ; mpfr_div(r27352, r27337, r27351, MPFR_RNDN); mpfr_mul(r27353, r27348, r27343, MPFR_RNDN); mpfr_mul(r27354, r27353, r27343, MPFR_RNDN); mpfr_mul(r27355, r27352, r27354, MPFR_RNDN); mpfr_add(r27356, r27350, r27355, MPFR_RNDN); ; mpfr_div(r27358, r27343, r27357, MPFR_RNDN); mpfr_pow(r27359, r27343, r27345, MPFR_RNDN); mpfr_mul(r27360, r27359, r27359, MPFR_RNDN); mpfr_mul(r27361, r27358, r27360, MPFR_RNDN); mpfr_pow(r27362, r27361, r27337, MPFR_RNDN); mpfr_add(r27363, r27356, r27362, MPFR_RNDN); mpfr_mul(r27364, r27340, r27363, MPFR_RNDN); mpfr_abs(r27365, r27364, MPFR_RNDN); return mpfr_get_d(r27365, MPFR_RNDN); } static mpfr_t r27366, r27367, r27368, r27369, r27370, r27371, r27372, r27373, r27374, r27375, r27376, r27377, r27378, r27379, r27380, r27381, r27382, r27383, r27384, r27385, r27386, r27387, r27388, r27389, r27390, r27391, r27392, r27393, r27394; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27366, "1", 10, MPFR_RNDN); mpfr_init(r27367); mpfr_init(r27368); mpfr_init(r27369); mpfr_init_set_str(r27370, "2", 10, MPFR_RNDN); mpfr_init(r27371); mpfr_init(r27372); mpfr_init(r27373); mpfr_init_set_str(r27374, "3", 10, MPFR_RNDN); mpfr_init(r27375); mpfr_init(r27376); mpfr_init(r27377); mpfr_init(r27378); mpfr_init(r27379); mpfr_init_set_str(r27380, "5", 10, MPFR_RNDN); mpfr_init(r27381); mpfr_init(r27382); mpfr_init(r27383); mpfr_init(r27384); mpfr_init(r27385); mpfr_init_set_str(r27386, "21", 10, MPFR_RNDN); mpfr_init(r27387); mpfr_init(r27388); mpfr_init(r27389); mpfr_init(r27390); mpfr_init(r27391); mpfr_init(r27392); mpfr_init(r27393); mpfr_init(r27394); } double f_dm(double x) { ; mpfr_const_pi(r27367, MPFR_RNDN); mpfr_sqrt(r27368, r27367, MPFR_RNDN); mpfr_div(r27369, r27366, r27368, MPFR_RNDN); ; mpfr_set_d(r27371, x, MPFR_RNDN); mpfr_abs(r27372, r27371, MPFR_RNDN); mpfr_mul(r27373, r27370, r27372, MPFR_RNDN); ; mpfr_div(r27375, r27370, r27374, MPFR_RNDN); mpfr_mul(r27376, r27372, r27372, MPFR_RNDN); mpfr_mul(r27377, r27376, r27372, MPFR_RNDN); mpfr_mul(r27378, r27375, r27377, MPFR_RNDN); mpfr_add(r27379, r27373, r27378, MPFR_RNDN); ; mpfr_div(r27381, r27366, r27380, MPFR_RNDN); mpfr_mul(r27382, r27377, r27372, MPFR_RNDN); mpfr_mul(r27383, r27382, r27372, MPFR_RNDN); mpfr_mul(r27384, r27381, r27383, MPFR_RNDN); mpfr_add(r27385, r27379, r27384, MPFR_RNDN); ; mpfr_div(r27387, r27372, r27386, MPFR_RNDN); mpfr_pow(r27388, r27372, r27374, MPFR_RNDN); mpfr_mul(r27389, r27388, r27388, MPFR_RNDN); mpfr_mul(r27390, r27387, r27389, MPFR_RNDN); mpfr_pow(r27391, r27390, r27366, MPFR_RNDN); mpfr_add(r27392, r27385, r27391, MPFR_RNDN); mpfr_mul(r27393, r27369, r27392, MPFR_RNDN); mpfr_abs(r27394, r27393, MPFR_RNDN); return mpfr_get_d(r27394, MPFR_RNDN); }