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