#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 r25064 = 1; float r25065 = atan2(1.0, 0.0); float r25066 = sqrt(r25065); float r25067 = r25064 / r25066; float r25068 = 2; float r25069 = x; float r25070 = fabs(r25069); float r25071 = r25068 * r25070; float r25072 = 3; float r25073 = r25068 / r25072; float r25074 = r25070 * r25070; float r25075 = r25074 * r25070; float r25076 = r25073 * r25075; float r25077 = r25071 + r25076; float r25078 = 5; float r25079 = r25064 / r25078; float r25080 = r25075 * r25070; float r25081 = r25080 * r25070; float r25082 = r25079 * r25081; float r25083 = r25077 + r25082; float r25084 = 21; float r25085 = r25064 / r25084; float r25086 = r25081 * r25070; float r25087 = r25086 * r25070; float r25088 = r25085 * r25087; float r25089 = r25083 + r25088; float r25090 = r25067 * r25089; float r25091 = fabs(r25090); return r25091; } double f_id(double x) { double r25092 = 1; double r25093 = atan2(1.0, 0.0); double r25094 = sqrt(r25093); double r25095 = r25092 / r25094; double r25096 = 2; double r25097 = x; double r25098 = fabs(r25097); double r25099 = r25096 * r25098; double r25100 = 3; double r25101 = r25096 / r25100; double r25102 = r25098 * r25098; double r25103 = r25102 * r25098; double r25104 = r25101 * r25103; double r25105 = r25099 + r25104; double r25106 = 5; double r25107 = r25092 / r25106; double r25108 = r25103 * r25098; double r25109 = r25108 * r25098; double r25110 = r25107 * r25109; double r25111 = r25105 + r25110; double r25112 = 21; double r25113 = r25092 / r25112; double r25114 = r25109 * r25098; double r25115 = r25114 * r25098; double r25116 = r25113 * r25115; double r25117 = r25111 + r25116; double r25118 = r25095 * r25117; double r25119 = fabs(r25118); return r25119; } double f_of(float x) { float r25120 = 1; float r25121 = atan2(1.0, 0.0); float r25122 = sqrt(r25121); float r25123 = r25120 / r25122; float r25124 = 2; float r25125 = x; float r25126 = fabs(r25125); float r25127 = r25124 * r25126; float r25128 = r25126 * r25126; float r25129 = 3/2; float r25130 = r25129 / r25126; float r25131 = r25128 / r25130; float r25132 = pow(r25131, r25120); float r25133 = r25127 + r25132; float r25134 = 5; float r25135 = r25120 / r25134; float r25136 = r25128 * r25126; float r25137 = r25136 * r25126; float r25138 = r25137 * r25126; float r25139 = r25135 * r25138; float r25140 = r25133 + r25139; float r25141 = 21; float r25142 = r25120 / r25141; float r25143 = r25138 * r25126; float r25144 = r25143 * r25126; float r25145 = r25142 * r25144; float r25146 = r25140 + r25145; float r25147 = r25123 * r25146; float r25148 = fabs(r25147); return r25148; } double f_od(double x) { double r25149 = 1; double r25150 = atan2(1.0, 0.0); double r25151 = sqrt(r25150); double r25152 = r25149 / r25151; double r25153 = 2; double r25154 = x; double r25155 = fabs(r25154); double r25156 = r25153 * r25155; double r25157 = r25155 * r25155; double r25158 = 3/2; double r25159 = r25158 / r25155; double r25160 = r25157 / r25159; double r25161 = pow(r25160, r25149); double r25162 = r25156 + r25161; double r25163 = 5; double r25164 = r25149 / r25163; double r25165 = r25157 * r25155; double r25166 = r25165 * r25155; double r25167 = r25166 * r25155; double r25168 = r25164 * r25167; double r25169 = r25162 + r25168; double r25170 = 21; double r25171 = r25149 / r25170; double r25172 = r25167 * r25155; double r25173 = r25172 * r25155; double r25174 = r25171 * r25173; double r25175 = r25169 + r25174; double r25176 = r25152 * r25175; double r25177 = fabs(r25176); return r25177; } 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 r25178, r25179, r25180, r25181, r25182, r25183, r25184, r25185, r25186, r25187, r25188, r25189, r25190, r25191, r25192, r25193, r25194, r25195, r25196, r25197, r25198, r25199, r25200, r25201, r25202, r25203, r25204, r25205; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r25178, "1", 10, MPFR_RNDN); mpfr_init(r25179); mpfr_init(r25180); mpfr_init(r25181); mpfr_init_set_str(r25182, "2", 10, MPFR_RNDN); mpfr_init(r25183); mpfr_init(r25184); mpfr_init(r25185); mpfr_init_set_str(r25186, "3", 10, MPFR_RNDN); mpfr_init(r25187); mpfr_init(r25188); mpfr_init(r25189); mpfr_init(r25190); mpfr_init(r25191); mpfr_init_set_str(r25192, "5", 10, MPFR_RNDN); mpfr_init(r25193); mpfr_init(r25194); mpfr_init(r25195); mpfr_init(r25196); mpfr_init(r25197); mpfr_init_set_str(r25198, "21", 10, MPFR_RNDN); mpfr_init(r25199); mpfr_init(r25200); mpfr_init(r25201); mpfr_init(r25202); mpfr_init(r25203); mpfr_init(r25204); mpfr_init(r25205); } double f_im(double x) { ; mpfr_const_pi(r25179, MPFR_RNDN); mpfr_sqrt(r25180, r25179, MPFR_RNDN); mpfr_div(r25181, r25178, r25180, MPFR_RNDN); ; mpfr_set_d(r25183, x, MPFR_RNDN); mpfr_abs(r25184, r25183, MPFR_RNDN); mpfr_mul(r25185, r25182, r25184, MPFR_RNDN); ; mpfr_div(r25187, r25182, r25186, MPFR_RNDN); mpfr_mul(r25188, r25184, r25184, MPFR_RNDN); mpfr_mul(r25189, r25188, r25184, MPFR_RNDN); mpfr_mul(r25190, r25187, r25189, MPFR_RNDN); mpfr_add(r25191, r25185, r25190, MPFR_RNDN); ; mpfr_div(r25193, r25178, r25192, MPFR_RNDN); mpfr_mul(r25194, r25189, r25184, MPFR_RNDN); mpfr_mul(r25195, r25194, r25184, MPFR_RNDN); mpfr_mul(r25196, r25193, r25195, MPFR_RNDN); mpfr_add(r25197, r25191, r25196, MPFR_RNDN); ; mpfr_div(r25199, r25178, r25198, MPFR_RNDN); mpfr_mul(r25200, r25195, r25184, MPFR_RNDN); mpfr_mul(r25201, r25200, r25184, MPFR_RNDN); mpfr_mul(r25202, r25199, r25201, MPFR_RNDN); mpfr_add(r25203, r25197, r25202, MPFR_RNDN); mpfr_mul(r25204, r25181, r25203, MPFR_RNDN); mpfr_abs(r25205, r25204, MPFR_RNDN); return mpfr_get_d(r25205, MPFR_RNDN); } static mpfr_t r25206, r25207, r25208, r25209, r25210, r25211, r25212, r25213, r25214, r25215, r25216, r25217, r25218, r25219, r25220, r25221, r25222, r25223, r25224, r25225, r25226, r25227, r25228, r25229, r25230, r25231, r25232, r25233, r25234; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r25206, "1", 10, MPFR_RNDN); mpfr_init(r25207); mpfr_init(r25208); mpfr_init(r25209); mpfr_init_set_str(r25210, "2", 10, MPFR_RNDN); mpfr_init(r25211); mpfr_init(r25212); mpfr_init(r25213); mpfr_init(r25214); mpfr_init_set_str(r25215, "3/2", 10, MPFR_RNDN); mpfr_init(r25216); mpfr_init(r25217); mpfr_init(r25218); mpfr_init(r25219); mpfr_init_set_str(r25220, "5", 10, MPFR_RNDN); mpfr_init(r25221); mpfr_init(r25222); mpfr_init(r25223); mpfr_init(r25224); mpfr_init(r25225); mpfr_init(r25226); mpfr_init_set_str(r25227, "21", 10, MPFR_RNDN); mpfr_init(r25228); mpfr_init(r25229); mpfr_init(r25230); mpfr_init(r25231); mpfr_init(r25232); mpfr_init(r25233); mpfr_init(r25234); } double f_fm(double x) { ; mpfr_const_pi(r25207, MPFR_RNDN); mpfr_sqrt(r25208, r25207, MPFR_RNDN); mpfr_div(r25209, r25206, r25208, MPFR_RNDN); ; mpfr_set_d(r25211, x, MPFR_RNDN); mpfr_abs(r25212, r25211, MPFR_RNDN); mpfr_mul(r25213, r25210, r25212, MPFR_RNDN); mpfr_mul(r25214, r25212, r25212, MPFR_RNDN); ; mpfr_div(r25216, r25215, r25212, MPFR_RNDN); mpfr_div(r25217, r25214, r25216, MPFR_RNDN); mpfr_pow(r25218, r25217, r25206, MPFR_RNDN); mpfr_add(r25219, r25213, r25218, MPFR_RNDN); ; mpfr_div(r25221, r25206, r25220, MPFR_RNDN); mpfr_mul(r25222, r25214, r25212, MPFR_RNDN); mpfr_mul(r25223, r25222, r25212, MPFR_RNDN); mpfr_mul(r25224, r25223, r25212, MPFR_RNDN); mpfr_mul(r25225, r25221, r25224, MPFR_RNDN); mpfr_add(r25226, r25219, r25225, MPFR_RNDN); ; mpfr_div(r25228, r25206, r25227, MPFR_RNDN); mpfr_mul(r25229, r25224, r25212, MPFR_RNDN); mpfr_mul(r25230, r25229, r25212, MPFR_RNDN); mpfr_mul(r25231, r25228, r25230, MPFR_RNDN); mpfr_add(r25232, r25226, r25231, MPFR_RNDN); mpfr_mul(r25233, r25209, r25232, MPFR_RNDN); mpfr_abs(r25234, r25233, MPFR_RNDN); return mpfr_get_d(r25234, MPFR_RNDN); } static mpfr_t r25235, r25236, r25237, r25238, r25239, r25240, r25241, r25242, r25243, r25244, r25245, r25246, r25247, r25248, r25249, r25250, r25251, r25252, r25253, r25254, r25255, r25256, r25257, r25258, r25259, r25260, r25261, r25262, r25263; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r25235, "1", 10, MPFR_RNDN); mpfr_init(r25236); mpfr_init(r25237); mpfr_init(r25238); mpfr_init_set_str(r25239, "2", 10, MPFR_RNDN); mpfr_init(r25240); mpfr_init(r25241); mpfr_init(r25242); mpfr_init(r25243); mpfr_init_set_str(r25244, "3/2", 10, MPFR_RNDN); mpfr_init(r25245); mpfr_init(r25246); mpfr_init(r25247); mpfr_init(r25248); mpfr_init_set_str(r25249, "5", 10, MPFR_RNDN); mpfr_init(r25250); mpfr_init(r25251); mpfr_init(r25252); mpfr_init(r25253); mpfr_init(r25254); mpfr_init(r25255); mpfr_init_set_str(r25256, "21", 10, MPFR_RNDN); mpfr_init(r25257); mpfr_init(r25258); mpfr_init(r25259); mpfr_init(r25260); mpfr_init(r25261); mpfr_init(r25262); mpfr_init(r25263); } double f_dm(double x) { ; mpfr_const_pi(r25236, MPFR_RNDN); mpfr_sqrt(r25237, r25236, MPFR_RNDN); mpfr_div(r25238, r25235, r25237, MPFR_RNDN); ; mpfr_set_d(r25240, x, MPFR_RNDN); mpfr_abs(r25241, r25240, MPFR_RNDN); mpfr_mul(r25242, r25239, r25241, MPFR_RNDN); mpfr_mul(r25243, r25241, r25241, MPFR_RNDN); ; mpfr_div(r25245, r25244, r25241, MPFR_RNDN); mpfr_div(r25246, r25243, r25245, MPFR_RNDN); mpfr_pow(r25247, r25246, r25235, MPFR_RNDN); mpfr_add(r25248, r25242, r25247, MPFR_RNDN); ; mpfr_div(r25250, r25235, r25249, MPFR_RNDN); mpfr_mul(r25251, r25243, r25241, MPFR_RNDN); mpfr_mul(r25252, r25251, r25241, MPFR_RNDN); mpfr_mul(r25253, r25252, r25241, MPFR_RNDN); mpfr_mul(r25254, r25250, r25253, MPFR_RNDN); mpfr_add(r25255, r25248, r25254, MPFR_RNDN); ; mpfr_div(r25257, r25235, r25256, MPFR_RNDN); mpfr_mul(r25258, r25253, r25241, MPFR_RNDN); mpfr_mul(r25259, r25258, r25241, MPFR_RNDN); mpfr_mul(r25260, r25257, r25259, MPFR_RNDN); mpfr_add(r25261, r25255, r25260, MPFR_RNDN); mpfr_mul(r25262, r25238, r25261, MPFR_RNDN); mpfr_abs(r25263, r25262, MPFR_RNDN); return mpfr_get_d(r25263, MPFR_RNDN); }