#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r20181 = 0.5f; float r20182 = re; float r20183 = sin(r20182); float r20184 = r20181 * r20183; float r20185 = im; float r20186 = -r20185; float r20187 = exp(r20186); float r20188 = exp(r20185); float r20189 = r20187 - r20188; float r20190 = r20184 * r20189; return r20190; } double f_id(double re, double im) { double r20191 = 0.5; double r20192 = re; double r20193 = sin(r20192); double r20194 = r20191 * r20193; double r20195 = im; double r20196 = -r20195; double r20197 = exp(r20196); double r20198 = exp(r20195); double r20199 = r20197 - r20198; double r20200 = r20194 * r20199; return r20200; } double f_of(float re, float im) { float r20201 = 0.5f; float r20202 = re; float r20203 = sin(r20202); float r20204 = r20201 * r20203; float r20205 = 0.016666666666666666f; float r20206 = im; float r20207 = 5.0f; float r20208 = pow(r20206, r20207); float r20209 = r20205 * r20208; float r20210 = 2.0f; float r20211 = r20210 * r20206; float r20212 = 0.3333333333333333f; float r20213 = 3.0f; float r20214 = pow(r20206, r20213); float r20215 = r20212 * r20214; float r20216 = r20211 + r20215; float r20217 = r20209 + r20216; float r20218 = -r20217; float r20219 = r20204 * r20218; return r20219; } double f_od(double re, double im) { double r20220 = 0.5; double r20221 = re; double r20222 = sin(r20221); double r20223 = r20220 * r20222; double r20224 = 0.016666666666666666; double r20225 = im; double r20226 = 5.0; double r20227 = pow(r20225, r20226); double r20228 = r20224 * r20227; double r20229 = 2.0; double r20230 = r20229 * r20225; double r20231 = 0.3333333333333333; double r20232 = 3.0; double r20233 = pow(r20225, r20232); double r20234 = r20231 * r20233; double r20235 = r20230 + r20234; double r20236 = r20228 + r20235; double r20237 = -r20236; double r20238 = r20223 * r20237; return r20238; } 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 r20239, r20240, r20241, r20242, r20243, r20244, r20245, r20246, r20247, r20248; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r20239, "0.5", 10, MPFR_RNDN); mpfr_init(r20240); mpfr_init(r20241); mpfr_init(r20242); mpfr_init(r20243); mpfr_init(r20244); mpfr_init(r20245); mpfr_init(r20246); mpfr_init(r20247); mpfr_init(r20248); } double f_im(double re, double im) { ; mpfr_set_d(r20240, re, MPFR_RNDN); mpfr_sin(r20241, r20240, MPFR_RNDN); mpfr_mul(r20242, r20239, r20241, MPFR_RNDN); mpfr_set_d(r20243, im, MPFR_RNDN); mpfr_neg(r20244, r20243, MPFR_RNDN); mpfr_exp(r20245, r20244, MPFR_RNDN); mpfr_exp(r20246, r20243, MPFR_RNDN); mpfr_sub(r20247, r20245, r20246, MPFR_RNDN); mpfr_mul(r20248, r20242, r20247, MPFR_RNDN); return mpfr_get_d(r20248, MPFR_RNDN); } static mpfr_t r20249, r20250, r20251, r20252, r20253, r20254, r20255, r20256, r20257, r20258, r20259, r20260, r20261, r20262, r20263, r20264, r20265, r20266, r20267; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20249, "0.5", 10, MPFR_RNDN); mpfr_init(r20250); mpfr_init(r20251); mpfr_init(r20252); mpfr_init_set_str(r20253, "1/60", 10, MPFR_RNDN); mpfr_init(r20254); mpfr_init_set_str(r20255, "5", 10, MPFR_RNDN); mpfr_init(r20256); mpfr_init(r20257); mpfr_init_set_str(r20258, "2", 10, MPFR_RNDN); mpfr_init(r20259); mpfr_init_set_str(r20260, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20261, "3", 10, MPFR_RNDN); mpfr_init(r20262); mpfr_init(r20263); mpfr_init(r20264); mpfr_init(r20265); mpfr_init(r20266); mpfr_init(r20267); } double f_fm(double re, double im) { ; mpfr_set_d(r20250, re, MPFR_RNDN); mpfr_sin(r20251, r20250, MPFR_RNDN); mpfr_mul(r20252, r20249, r20251, MPFR_RNDN); ; mpfr_set_d(r20254, im, MPFR_RNDN); ; mpfr_pow(r20256, r20254, r20255, MPFR_RNDN); mpfr_mul(r20257, r20253, r20256, MPFR_RNDN); ; mpfr_mul(r20259, r20258, r20254, MPFR_RNDN); ; ; mpfr_pow(r20262, r20254, r20261, MPFR_RNDN); mpfr_mul(r20263, r20260, r20262, MPFR_RNDN); mpfr_add(r20264, r20259, r20263, MPFR_RNDN); mpfr_add(r20265, r20257, r20264, MPFR_RNDN); mpfr_neg(r20266, r20265, MPFR_RNDN); mpfr_mul(r20267, r20252, r20266, MPFR_RNDN); return mpfr_get_d(r20267, MPFR_RNDN); } static mpfr_t r20268, r20269, r20270, r20271, r20272, r20273, r20274, r20275, r20276, r20277, r20278, r20279, r20280, r20281, r20282, r20283, r20284, r20285, r20286; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20268, "0.5", 10, MPFR_RNDN); mpfr_init(r20269); mpfr_init(r20270); mpfr_init(r20271); mpfr_init_set_str(r20272, "1/60", 10, MPFR_RNDN); mpfr_init(r20273); mpfr_init_set_str(r20274, "5", 10, MPFR_RNDN); mpfr_init(r20275); mpfr_init(r20276); mpfr_init_set_str(r20277, "2", 10, MPFR_RNDN); mpfr_init(r20278); mpfr_init_set_str(r20279, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20280, "3", 10, MPFR_RNDN); mpfr_init(r20281); mpfr_init(r20282); mpfr_init(r20283); mpfr_init(r20284); mpfr_init(r20285); mpfr_init(r20286); } double f_dm(double re, double im) { ; mpfr_set_d(r20269, re, MPFR_RNDN); mpfr_sin(r20270, r20269, MPFR_RNDN); mpfr_mul(r20271, r20268, r20270, MPFR_RNDN); ; mpfr_set_d(r20273, im, MPFR_RNDN); ; mpfr_pow(r20275, r20273, r20274, MPFR_RNDN); mpfr_mul(r20276, r20272, r20275, MPFR_RNDN); ; mpfr_mul(r20278, r20277, r20273, MPFR_RNDN); ; ; mpfr_pow(r20281, r20273, r20280, MPFR_RNDN); mpfr_mul(r20282, r20279, r20281, MPFR_RNDN); mpfr_add(r20283, r20278, r20282, MPFR_RNDN); mpfr_add(r20284, r20276, r20283, MPFR_RNDN); mpfr_neg(r20285, r20284, MPFR_RNDN); mpfr_mul(r20286, r20271, r20285, MPFR_RNDN); return mpfr_get_d(r20286, MPFR_RNDN); }