#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r25167 = 0.5; float r25168 = re; float r25169 = cos(r25168); float r25170 = r25167 * r25169; float r25171 = 0; float r25172 = im; float r25173 = r25171 - r25172; float r25174 = exp(r25173); float r25175 = exp(r25172); float r25176 = r25174 - r25175; float r25177 = r25170 * r25176; return r25177; } double f_id(double re, double im) { double r25178 = 0.5; double r25179 = re; double r25180 = cos(r25179); double r25181 = r25178 * r25180; double r25182 = 0; double r25183 = im; double r25184 = r25182 - r25183; double r25185 = exp(r25184); double r25186 = exp(r25183); double r25187 = r25185 - r25186; double r25188 = r25181 * r25187; return r25188; } double f_of(float re, float im) { float r25189 = im; float r25190 = 1/9; float r25191 = 2; float r25192 = pow(r25189, r25191); float r25193 = r25190 * r25192; float r25194 = 3; float r25195 = log(r25194); float r25196 = r25193 + r25195; float r25197 = 1/162; float r25198 = 4; float r25199 = pow(r25189, r25198); float r25200 = r25197 * r25199; float r25201 = r25196 - r25200; float r25202 = expm1(r25201); float r25203 = 5; float r25204 = pow(r25189, r25203); float r25205 = 1/60; float r25206 = r25204 * r25205; float r25207 = fma(r25189, r25202, r25206); float r25208 = re; float r25209 = cos(r25208); float r25210 = 0.5; float r25211 = -r25210; float r25212 = r25209 * r25211; float r25213 = r25207 * r25212; return r25213; } double f_od(double re, double im) { double r25214 = im; double r25215 = 1/9; double r25216 = 2; double r25217 = pow(r25214, r25216); double r25218 = r25215 * r25217; double r25219 = 3; double r25220 = log(r25219); double r25221 = r25218 + r25220; double r25222 = 1/162; double r25223 = 4; double r25224 = pow(r25214, r25223); double r25225 = r25222 * r25224; double r25226 = r25221 - r25225; double r25227 = expm1(r25226); double r25228 = 5; double r25229 = pow(r25214, r25228); double r25230 = 1/60; double r25231 = r25229 * r25230; double r25232 = fma(r25214, r25227, r25231); double r25233 = re; double r25234 = cos(r25233); double r25235 = 0.5; double r25236 = -r25235; double r25237 = r25234 * r25236; double r25238 = r25232 * r25237; return r25238; } 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 r25239, r25240, r25241, r25242, r25243, r25244, r25245, r25246, r25247, r25248, r25249; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r25239, "0.5", 10, MPFR_RNDN); mpfr_init(r25240); mpfr_init(r25241); mpfr_init(r25242); mpfr_init_set_str(r25243, "0", 10, MPFR_RNDN); mpfr_init(r25244); mpfr_init(r25245); mpfr_init(r25246); mpfr_init(r25247); mpfr_init(r25248); mpfr_init(r25249); } double f_im(double re, double im) { ; mpfr_set_d(r25240, re, MPFR_RNDN); mpfr_cos(r25241, r25240, MPFR_RNDN); mpfr_mul(r25242, r25239, r25241, MPFR_RNDN); ; mpfr_set_d(r25244, im, MPFR_RNDN); mpfr_sub(r25245, r25243, r25244, MPFR_RNDN); mpfr_exp(r25246, r25245, MPFR_RNDN); mpfr_exp(r25247, r25244, MPFR_RNDN); mpfr_sub(r25248, r25246, r25247, MPFR_RNDN); mpfr_mul(r25249, r25242, r25248, MPFR_RNDN); return mpfr_get_d(r25249, MPFR_RNDN); } static mpfr_t r25250, r25251, r25252, r25253, r25254, r25255, r25256, r25257, r25258, r25259, r25260, r25261, r25262, r25263, r25264, r25265, r25266, r25267, r25268, r25269, r25270, r25271, r25272, r25273, r25274; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25250); mpfr_init_set_str(r25251, "1/9", 10, MPFR_RNDN); mpfr_init_set_str(r25252, "2", 10, MPFR_RNDN); mpfr_init(r25253); mpfr_init(r25254); mpfr_init_set_str(r25255, "3", 10, MPFR_RNDN); mpfr_init(r25256); mpfr_init(r25257); mpfr_init_set_str(r25258, "1/162", 10, MPFR_RNDN); mpfr_init_set_str(r25259, "4", 10, MPFR_RNDN); mpfr_init(r25260); mpfr_init(r25261); mpfr_init(r25262); mpfr_init(r25263); mpfr_init_set_str(r25264, "5", 10, MPFR_RNDN); mpfr_init(r25265); mpfr_init_set_str(r25266, "1/60", 10, MPFR_RNDN); mpfr_init(r25267); mpfr_init(r25268); mpfr_init(r25269); mpfr_init(r25270); mpfr_init_set_str(r25271, "0.5", 10, MPFR_RNDN); mpfr_init(r25272); mpfr_init(r25273); mpfr_init(r25274); } double f_fm(double re, double im) { mpfr_set_d(r25250, im, MPFR_RNDN); ; ; mpfr_pow(r25253, r25250, r25252, MPFR_RNDN); mpfr_mul(r25254, r25251, r25253, MPFR_RNDN); ; mpfr_log(r25256, r25255, MPFR_RNDN); mpfr_add(r25257, r25254, r25256, MPFR_RNDN); ; ; mpfr_pow(r25260, r25250, r25259, MPFR_RNDN); mpfr_mul(r25261, r25258, r25260, MPFR_RNDN); mpfr_sub(r25262, r25257, r25261, MPFR_RNDN); mpfr_expm1(r25263, r25262, MPFR_RNDN); ; mpfr_pow(r25265, r25250, r25264, MPFR_RNDN); ; mpfr_mul(r25267, r25265, r25266, MPFR_RNDN); mpfr_fma(r25268, r25250, r25263, r25267, MPFR_RNDN); mpfr_set_d(r25269, re, MPFR_RNDN); mpfr_cos(r25270, r25269, MPFR_RNDN); ; mpfr_neg(r25272, r25271, MPFR_RNDN); mpfr_mul(r25273, r25270, r25272, MPFR_RNDN); mpfr_mul(r25274, r25268, r25273, MPFR_RNDN); return mpfr_get_d(r25274, MPFR_RNDN); } static mpfr_t r25275, r25276, r25277, r25278, r25279, r25280, r25281, r25282, r25283, r25284, r25285, r25286, r25287, r25288, r25289, r25290, r25291, r25292, r25293, r25294, r25295, r25296, r25297, r25298, r25299; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r25275); mpfr_init_set_str(r25276, "1/9", 10, MPFR_RNDN); mpfr_init_set_str(r25277, "2", 10, MPFR_RNDN); mpfr_init(r25278); mpfr_init(r25279); mpfr_init_set_str(r25280, "3", 10, MPFR_RNDN); mpfr_init(r25281); mpfr_init(r25282); mpfr_init_set_str(r25283, "1/162", 10, MPFR_RNDN); mpfr_init_set_str(r25284, "4", 10, MPFR_RNDN); mpfr_init(r25285); mpfr_init(r25286); mpfr_init(r25287); mpfr_init(r25288); mpfr_init_set_str(r25289, "5", 10, MPFR_RNDN); mpfr_init(r25290); mpfr_init_set_str(r25291, "1/60", 10, MPFR_RNDN); mpfr_init(r25292); mpfr_init(r25293); mpfr_init(r25294); mpfr_init(r25295); mpfr_init_set_str(r25296, "0.5", 10, MPFR_RNDN); mpfr_init(r25297); mpfr_init(r25298); mpfr_init(r25299); } double f_dm(double re, double im) { mpfr_set_d(r25275, im, MPFR_RNDN); ; ; mpfr_pow(r25278, r25275, r25277, MPFR_RNDN); mpfr_mul(r25279, r25276, r25278, MPFR_RNDN); ; mpfr_log(r25281, r25280, MPFR_RNDN); mpfr_add(r25282, r25279, r25281, MPFR_RNDN); ; ; mpfr_pow(r25285, r25275, r25284, MPFR_RNDN); mpfr_mul(r25286, r25283, r25285, MPFR_RNDN); mpfr_sub(r25287, r25282, r25286, MPFR_RNDN); mpfr_expm1(r25288, r25287, MPFR_RNDN); ; mpfr_pow(r25290, r25275, r25289, MPFR_RNDN); ; mpfr_mul(r25292, r25290, r25291, MPFR_RNDN); mpfr_fma(r25293, r25275, r25288, r25292, MPFR_RNDN); mpfr_set_d(r25294, re, MPFR_RNDN); mpfr_cos(r25295, r25294, MPFR_RNDN); ; mpfr_neg(r25297, r25296, MPFR_RNDN); mpfr_mul(r25298, r25295, r25297, MPFR_RNDN); mpfr_mul(r25299, r25293, r25298, MPFR_RNDN); return mpfr_get_d(r25299, MPFR_RNDN); }