#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r20108 = 0.5f; float r20109 = re; float r20110 = cos(r20109); float r20111 = r20108 * r20110; float r20112 = 0.0f; float r20113 = im; float r20114 = r20112 - r20113; float r20115 = exp(r20114); float r20116 = exp(r20113); float r20117 = r20115 - r20116; float r20118 = r20111 * r20117; return r20118; } double f_id(double re, double im) { double r20119 = 0.5; double r20120 = re; double r20121 = cos(r20120); double r20122 = r20119 * r20121; double r20123 = 0.0; double r20124 = im; double r20125 = r20123 - r20124; double r20126 = exp(r20125); double r20127 = exp(r20124); double r20128 = r20126 - r20127; double r20129 = r20122 * r20128; return r20129; } double f_of(float re, float im) { float r20130 = 0.016666666666666666f; float r20131 = im; float r20132 = 5.0f; float r20133 = pow(r20131, r20132); float r20134 = r20130 * r20133; float r20135 = 2.0f; float r20136 = r20135 * r20131; float r20137 = 0.3333333333333333f; float r20138 = 3.0f; float r20139 = pow(r20131, r20138); float r20140 = r20137 * r20139; float r20141 = r20136 + r20140; float r20142 = r20134 + r20141; float r20143 = -r20142; float r20144 = re; float r20145 = cos(r20144); float r20146 = 0.5f; float r20147 = r20145 * r20146; float r20148 = r20143 * r20147; return r20148; } double f_od(double re, double im) { double r20149 = 0.016666666666666666; double r20150 = im; double r20151 = 5.0; double r20152 = pow(r20150, r20151); double r20153 = r20149 * r20152; double r20154 = 2.0; double r20155 = r20154 * r20150; double r20156 = 0.3333333333333333; double r20157 = 3.0; double r20158 = pow(r20150, r20157); double r20159 = r20156 * r20158; double r20160 = r20155 + r20159; double r20161 = r20153 + r20160; double r20162 = -r20161; double r20163 = re; double r20164 = cos(r20163); double r20165 = 0.5; double r20166 = r20164 * r20165; double r20167 = r20162 * r20166; return r20167; } 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 r20168, r20169, r20170, r20171, r20172, r20173, r20174, r20175, r20176, r20177, r20178; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r20168, "0.5", 10, MPFR_RNDN); mpfr_init(r20169); mpfr_init(r20170); mpfr_init(r20171); mpfr_init_set_str(r20172, "0", 10, MPFR_RNDN); mpfr_init(r20173); mpfr_init(r20174); mpfr_init(r20175); mpfr_init(r20176); mpfr_init(r20177); mpfr_init(r20178); } double f_im(double re, double im) { ; mpfr_set_d(r20169, re, MPFR_RNDN); mpfr_cos(r20170, r20169, MPFR_RNDN); mpfr_mul(r20171, r20168, r20170, MPFR_RNDN); ; mpfr_set_d(r20173, im, MPFR_RNDN); mpfr_sub(r20174, r20172, r20173, MPFR_RNDN); mpfr_exp(r20175, r20174, MPFR_RNDN); mpfr_exp(r20176, r20173, MPFR_RNDN); mpfr_sub(r20177, r20175, r20176, MPFR_RNDN); mpfr_mul(r20178, r20171, r20177, MPFR_RNDN); return mpfr_get_d(r20178, MPFR_RNDN); } static mpfr_t r20179, r20180, r20181, r20182, r20183, r20184, r20185, r20186, r20187, r20188, r20189, r20190, r20191, r20192, r20193, r20194, r20195, r20196, r20197; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20179, "1/60", 10, MPFR_RNDN); mpfr_init(r20180); mpfr_init_set_str(r20181, "5", 10, MPFR_RNDN); mpfr_init(r20182); mpfr_init(r20183); mpfr_init_set_str(r20184, "2", 10, MPFR_RNDN); mpfr_init(r20185); mpfr_init_set_str(r20186, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20187, "3", 10, MPFR_RNDN); mpfr_init(r20188); mpfr_init(r20189); mpfr_init(r20190); mpfr_init(r20191); mpfr_init(r20192); mpfr_init(r20193); mpfr_init(r20194); mpfr_init_set_str(r20195, "0.5", 10, MPFR_RNDN); mpfr_init(r20196); mpfr_init(r20197); } double f_fm(double re, double im) { ; mpfr_set_d(r20180, im, MPFR_RNDN); ; mpfr_pow(r20182, r20180, r20181, MPFR_RNDN); mpfr_mul(r20183, r20179, r20182, MPFR_RNDN); ; mpfr_mul(r20185, r20184, r20180, MPFR_RNDN); ; ; mpfr_pow(r20188, r20180, r20187, MPFR_RNDN); mpfr_mul(r20189, r20186, r20188, MPFR_RNDN); mpfr_add(r20190, r20185, r20189, MPFR_RNDN); mpfr_add(r20191, r20183, r20190, MPFR_RNDN); mpfr_neg(r20192, r20191, MPFR_RNDN); mpfr_set_d(r20193, re, MPFR_RNDN); mpfr_cos(r20194, r20193, MPFR_RNDN); ; mpfr_mul(r20196, r20194, r20195, MPFR_RNDN); mpfr_mul(r20197, r20192, r20196, MPFR_RNDN); return mpfr_get_d(r20197, MPFR_RNDN); } static mpfr_t r20198, r20199, r20200, r20201, r20202, r20203, r20204, r20205, r20206, r20207, r20208, r20209, r20210, r20211, r20212, r20213, r20214, r20215, r20216; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20198, "1/60", 10, MPFR_RNDN); mpfr_init(r20199); mpfr_init_set_str(r20200, "5", 10, MPFR_RNDN); mpfr_init(r20201); mpfr_init(r20202); mpfr_init_set_str(r20203, "2", 10, MPFR_RNDN); mpfr_init(r20204); mpfr_init_set_str(r20205, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20206, "3", 10, MPFR_RNDN); mpfr_init(r20207); mpfr_init(r20208); mpfr_init(r20209); mpfr_init(r20210); mpfr_init(r20211); mpfr_init(r20212); mpfr_init(r20213); mpfr_init_set_str(r20214, "0.5", 10, MPFR_RNDN); mpfr_init(r20215); mpfr_init(r20216); } double f_dm(double re, double im) { ; mpfr_set_d(r20199, im, MPFR_RNDN); ; mpfr_pow(r20201, r20199, r20200, MPFR_RNDN); mpfr_mul(r20202, r20198, r20201, MPFR_RNDN); ; mpfr_mul(r20204, r20203, r20199, MPFR_RNDN); ; ; mpfr_pow(r20207, r20199, r20206, MPFR_RNDN); mpfr_mul(r20208, r20205, r20207, MPFR_RNDN); mpfr_add(r20209, r20204, r20208, MPFR_RNDN); mpfr_add(r20210, r20202, r20209, MPFR_RNDN); mpfr_neg(r20211, r20210, MPFR_RNDN); mpfr_set_d(r20212, re, MPFR_RNDN); mpfr_cos(r20213, r20212, MPFR_RNDN); ; mpfr_mul(r20215, r20213, r20214, MPFR_RNDN); mpfr_mul(r20216, r20211, r20215, MPFR_RNDN); return mpfr_get_d(r20216, MPFR_RNDN); }