#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 r26163 = 1; float r26164 = atan2(1.0, 0.0); float r26165 = sqrt(r26164); float r26166 = r26163 / r26165; float r26167 = 2; float r26168 = x; float r26169 = fabs(r26168); float r26170 = r26167 * r26169; float r26171 = 3; float r26172 = r26167 / r26171; float r26173 = r26169 * r26169; float r26174 = r26173 * r26169; float r26175 = r26172 * r26174; float r26176 = r26170 + r26175; float r26177 = 5; float r26178 = r26163 / r26177; float r26179 = r26174 * r26169; float r26180 = r26179 * r26169; float r26181 = r26178 * r26180; float r26182 = r26176 + r26181; float r26183 = 21; float r26184 = r26163 / r26183; float r26185 = r26180 * r26169; float r26186 = r26185 * r26169; float r26187 = r26184 * r26186; float r26188 = r26182 + r26187; float r26189 = r26166 * r26188; float r26190 = fabs(r26189); return r26190; } double f_id(double x) { double r26191 = 1; double r26192 = atan2(1.0, 0.0); double r26193 = sqrt(r26192); double r26194 = r26191 / r26193; double r26195 = 2; double r26196 = x; double r26197 = fabs(r26196); double r26198 = r26195 * r26197; double r26199 = 3; double r26200 = r26195 / r26199; double r26201 = r26197 * r26197; double r26202 = r26201 * r26197; double r26203 = r26200 * r26202; double r26204 = r26198 + r26203; double r26205 = 5; double r26206 = r26191 / r26205; double r26207 = r26202 * r26197; double r26208 = r26207 * r26197; double r26209 = r26206 * r26208; double r26210 = r26204 + r26209; double r26211 = 21; double r26212 = r26191 / r26211; double r26213 = r26208 * r26197; double r26214 = r26213 * r26197; double r26215 = r26212 * r26214; double r26216 = r26210 + r26215; double r26217 = r26194 * r26216; double r26218 = fabs(r26217); return r26218; } double f_of(float x) { float r26219 = 1; float r26220 = atan2(1.0, 0.0); float r26221 = sqrt(r26220); float r26222 = r26219 / r26221; float r26223 = 2; float r26224 = x; float r26225 = fabs(r26224); float r26226 = r26223 * r26225; float r26227 = 3; float r26228 = r26223 / r26227; float r26229 = r26225 * r26225; float r26230 = r26229 * r26225; float r26231 = r26228 * r26230; float r26232 = r26226 + r26231; float r26233 = 5; float r26234 = r26219 / r26233; float r26235 = r26230 * r26225; float r26236 = r26235 * r26225; float r26237 = r26234 * r26236; float r26238 = r26232 + r26237; float r26239 = 21; float r26240 = r26219 / r26239; float r26241 = r26236 * r26225; float r26242 = r26241 * r26225; float r26243 = r26240 * r26242; float r26244 = r26238 + r26243; float r26245 = r26222 * r26244; float r26246 = fabs(r26245); return r26246; } double f_od(double x) { double r26247 = 1; double r26248 = atan2(1.0, 0.0); double r26249 = sqrt(r26248); double r26250 = r26247 / r26249; double r26251 = 2; double r26252 = x; double r26253 = fabs(r26252); double r26254 = r26251 * r26253; double r26255 = 3; double r26256 = r26251 / r26255; double r26257 = r26253 * r26253; double r26258 = r26257 * r26253; double r26259 = r26256 * r26258; double r26260 = r26254 + r26259; double r26261 = 5; double r26262 = r26247 / r26261; double r26263 = r26258 * r26253; double r26264 = r26263 * r26253; double r26265 = r26262 * r26264; double r26266 = r26260 + r26265; double r26267 = 21; double r26268 = r26247 / r26267; double r26269 = r26264 * r26253; double r26270 = r26269 * r26253; double r26271 = r26268 * r26270; double r26272 = r26266 + r26271; double r26273 = r26250 * r26272; double r26274 = fabs(r26273); return r26274; } 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 r26275, r26276, r26277, r26278, r26279, r26280, r26281, r26282, r26283, r26284, r26285, r26286, r26287, r26288, r26289, r26290, r26291, r26292, r26293, r26294, r26295, r26296, r26297, r26298, r26299, r26300, r26301, r26302; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r26275, "1", 10, MPFR_RNDN); mpfr_init(r26276); mpfr_init(r26277); mpfr_init(r26278); mpfr_init_set_str(r26279, "2", 10, MPFR_RNDN); mpfr_init(r26280); mpfr_init(r26281); mpfr_init(r26282); mpfr_init_set_str(r26283, "3", 10, MPFR_RNDN); mpfr_init(r26284); mpfr_init(r26285); mpfr_init(r26286); mpfr_init(r26287); mpfr_init(r26288); mpfr_init_set_str(r26289, "5", 10, MPFR_RNDN); mpfr_init(r26290); mpfr_init(r26291); mpfr_init(r26292); mpfr_init(r26293); mpfr_init(r26294); mpfr_init_set_str(r26295, "21", 10, MPFR_RNDN); mpfr_init(r26296); mpfr_init(r26297); mpfr_init(r26298); mpfr_init(r26299); mpfr_init(r26300); mpfr_init(r26301); mpfr_init(r26302); } double f_im(double x) { ; mpfr_const_pi(r26276, MPFR_RNDN); mpfr_sqrt(r26277, r26276, MPFR_RNDN); mpfr_div(r26278, r26275, r26277, MPFR_RNDN); ; mpfr_set_d(r26280, x, MPFR_RNDN); mpfr_abs(r26281, r26280, MPFR_RNDN); mpfr_mul(r26282, r26279, r26281, MPFR_RNDN); ; mpfr_div(r26284, r26279, r26283, MPFR_RNDN); mpfr_mul(r26285, r26281, r26281, MPFR_RNDN); mpfr_mul(r26286, r26285, r26281, MPFR_RNDN); mpfr_mul(r26287, r26284, r26286, MPFR_RNDN); mpfr_add(r26288, r26282, r26287, MPFR_RNDN); ; mpfr_div(r26290, r26275, r26289, MPFR_RNDN); mpfr_mul(r26291, r26286, r26281, MPFR_RNDN); mpfr_mul(r26292, r26291, r26281, MPFR_RNDN); mpfr_mul(r26293, r26290, r26292, MPFR_RNDN); mpfr_add(r26294, r26288, r26293, MPFR_RNDN); ; mpfr_div(r26296, r26275, r26295, MPFR_RNDN); mpfr_mul(r26297, r26292, r26281, MPFR_RNDN); mpfr_mul(r26298, r26297, r26281, MPFR_RNDN); mpfr_mul(r26299, r26296, r26298, MPFR_RNDN); mpfr_add(r26300, r26294, r26299, MPFR_RNDN); mpfr_mul(r26301, r26278, r26300, MPFR_RNDN); mpfr_abs(r26302, r26301, MPFR_RNDN); return mpfr_get_d(r26302, MPFR_RNDN); } static mpfr_t r26303, r26304, r26305, r26306, r26307, r26308, r26309, r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318, r26319, r26320, r26321, r26322, r26323, r26324, r26325, r26326, r26327, r26328, r26329, r26330; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26303, "1", 10, MPFR_RNDN); mpfr_init(r26304); mpfr_init(r26305); mpfr_init(r26306); mpfr_init_set_str(r26307, "2", 10, MPFR_RNDN); mpfr_init(r26308); mpfr_init(r26309); mpfr_init(r26310); mpfr_init_set_str(r26311, "3", 10, MPFR_RNDN); mpfr_init(r26312); mpfr_init(r26313); mpfr_init(r26314); mpfr_init(r26315); mpfr_init(r26316); mpfr_init_set_str(r26317, "5", 10, MPFR_RNDN); mpfr_init(r26318); mpfr_init(r26319); mpfr_init(r26320); mpfr_init(r26321); mpfr_init(r26322); mpfr_init_set_str(r26323, "21", 10, MPFR_RNDN); mpfr_init(r26324); mpfr_init(r26325); mpfr_init(r26326); mpfr_init(r26327); mpfr_init(r26328); mpfr_init(r26329); mpfr_init(r26330); } double f_fm(double x) { ; mpfr_const_pi(r26304, MPFR_RNDN); mpfr_sqrt(r26305, r26304, MPFR_RNDN); mpfr_div(r26306, r26303, r26305, MPFR_RNDN); ; mpfr_set_d(r26308, x, MPFR_RNDN); mpfr_abs(r26309, r26308, MPFR_RNDN); mpfr_mul(r26310, r26307, r26309, MPFR_RNDN); ; mpfr_div(r26312, r26307, r26311, MPFR_RNDN); mpfr_mul(r26313, r26309, r26309, MPFR_RNDN); mpfr_mul(r26314, r26313, r26309, MPFR_RNDN); mpfr_mul(r26315, r26312, r26314, MPFR_RNDN); mpfr_add(r26316, r26310, r26315, MPFR_RNDN); ; mpfr_div(r26318, r26303, r26317, MPFR_RNDN); mpfr_mul(r26319, r26314, r26309, MPFR_RNDN); mpfr_mul(r26320, r26319, r26309, MPFR_RNDN); mpfr_mul(r26321, r26318, r26320, MPFR_RNDN); mpfr_add(r26322, r26316, r26321, MPFR_RNDN); ; mpfr_div(r26324, r26303, r26323, MPFR_RNDN); mpfr_mul(r26325, r26320, r26309, MPFR_RNDN); mpfr_mul(r26326, r26325, r26309, MPFR_RNDN); mpfr_mul(r26327, r26324, r26326, MPFR_RNDN); mpfr_add(r26328, r26322, r26327, MPFR_RNDN); mpfr_mul(r26329, r26306, r26328, MPFR_RNDN); mpfr_abs(r26330, r26329, MPFR_RNDN); return mpfr_get_d(r26330, MPFR_RNDN); } static mpfr_t r26331, r26332, r26333, r26334, r26335, r26336, r26337, r26338, r26339, r26340, r26341, r26342, r26343, r26344, r26345, r26346, r26347, r26348, r26349, r26350, r26351, r26352, r26353, r26354, r26355, r26356, r26357, r26358; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26331, "1", 10, MPFR_RNDN); mpfr_init(r26332); mpfr_init(r26333); mpfr_init(r26334); mpfr_init_set_str(r26335, "2", 10, MPFR_RNDN); mpfr_init(r26336); mpfr_init(r26337); mpfr_init(r26338); mpfr_init_set_str(r26339, "3", 10, MPFR_RNDN); mpfr_init(r26340); mpfr_init(r26341); mpfr_init(r26342); mpfr_init(r26343); mpfr_init(r26344); mpfr_init_set_str(r26345, "5", 10, MPFR_RNDN); mpfr_init(r26346); mpfr_init(r26347); mpfr_init(r26348); mpfr_init(r26349); mpfr_init(r26350); mpfr_init_set_str(r26351, "21", 10, MPFR_RNDN); mpfr_init(r26352); mpfr_init(r26353); mpfr_init(r26354); mpfr_init(r26355); mpfr_init(r26356); mpfr_init(r26357); mpfr_init(r26358); } double f_dm(double x) { ; mpfr_const_pi(r26332, MPFR_RNDN); mpfr_sqrt(r26333, r26332, MPFR_RNDN); mpfr_div(r26334, r26331, r26333, MPFR_RNDN); ; mpfr_set_d(r26336, x, MPFR_RNDN); mpfr_abs(r26337, r26336, MPFR_RNDN); mpfr_mul(r26338, r26335, r26337, MPFR_RNDN); ; mpfr_div(r26340, r26335, r26339, MPFR_RNDN); mpfr_mul(r26341, r26337, r26337, MPFR_RNDN); mpfr_mul(r26342, r26341, r26337, MPFR_RNDN); mpfr_mul(r26343, r26340, r26342, MPFR_RNDN); mpfr_add(r26344, r26338, r26343, MPFR_RNDN); ; mpfr_div(r26346, r26331, r26345, MPFR_RNDN); mpfr_mul(r26347, r26342, r26337, MPFR_RNDN); mpfr_mul(r26348, r26347, r26337, MPFR_RNDN); mpfr_mul(r26349, r26346, r26348, MPFR_RNDN); mpfr_add(r26350, r26344, r26349, MPFR_RNDN); ; mpfr_div(r26352, r26331, r26351, MPFR_RNDN); mpfr_mul(r26353, r26348, r26337, MPFR_RNDN); mpfr_mul(r26354, r26353, r26337, MPFR_RNDN); mpfr_mul(r26355, r26352, r26354, MPFR_RNDN); mpfr_add(r26356, r26350, r26355, MPFR_RNDN); mpfr_mul(r26357, r26334, r26356, MPFR_RNDN); mpfr_abs(r26358, r26357, MPFR_RNDN); return mpfr_get_d(r26358, MPFR_RNDN); }