#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19223 = 0.5f; float r19224 = re; float r19225 = cos(r19224); float r19226 = r19223 * r19225; float r19227 = 0.0f; float r19228 = im; float r19229 = r19227 - r19228; float r19230 = exp(r19229); float r19231 = exp(r19228); float r19232 = r19230 - r19231; float r19233 = r19226 * r19232; return r19233; } double f_id(double re, double im) { double r19234 = 0.5; double r19235 = re; double r19236 = cos(r19235); double r19237 = r19234 * r19236; double r19238 = 0.0; double r19239 = im; double r19240 = r19238 - r19239; double r19241 = exp(r19240); double r19242 = exp(r19239); double r19243 = r19241 - r19242; double r19244 = r19237 * r19243; return r19244; } double f_of(float re, float im) { float r19245 = 0.01666666753590107f; float r19246 = im; float r19247 = 5.0f; float r19248 = pow(r19246, r19247); float r19249 = r19245 * r19248; float r19250 = 2.0f; float r19251 = r19250 * r19246; float r19252 = 0.3333333432674408f; float r19253 = 3.0f; float r19254 = pow(r19246, r19253); float r19255 = r19252 * r19254; float r19256 = r19251 + r19255; float r19257 = r19249 + r19256; float r19258 = -r19257; float r19259 = re; float r19260 = cos(r19259); float r19261 = 0.5f; float r19262 = r19260 * r19261; float r19263 = r19258 * r19262; return r19263; } double f_od(double re, double im) { double r19264 = 0.01666666753590107; double r19265 = im; double r19266 = 5.0; double r19267 = pow(r19265, r19266); double r19268 = r19264 * r19267; double r19269 = 2.0; double r19270 = r19269 * r19265; double r19271 = 0.3333333432674408; double r19272 = 3.0; double r19273 = pow(r19265, r19272); double r19274 = r19271 * r19273; double r19275 = r19270 + r19274; double r19276 = r19268 + r19275; double r19277 = -r19276; double r19278 = re; double r19279 = cos(r19278); double r19280 = 0.5; double r19281 = r19279 * r19280; double r19282 = r19277 * r19281; return r19282; } 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 r19283, r19284, r19285, r19286, r19287, r19288, r19289, r19290, r19291, r19292, r19293; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19283, "0.5", 10, MPFR_RNDN); mpfr_init(r19284); mpfr_init(r19285); mpfr_init(r19286); mpfr_init_set_str(r19287, "0", 10, MPFR_RNDN); mpfr_init(r19288); mpfr_init(r19289); mpfr_init(r19290); mpfr_init(r19291); mpfr_init(r19292); mpfr_init(r19293); } double f_im(double re, double im) { ; mpfr_set_d(r19284, re, MPFR_RNDN); mpfr_cos(r19285, r19284, MPFR_RNDN); mpfr_mul(r19286, r19283, r19285, MPFR_RNDN); ; mpfr_set_d(r19288, im, MPFR_RNDN); mpfr_sub(r19289, r19287, r19288, MPFR_RNDN); mpfr_exp(r19290, r19289, MPFR_RNDN); mpfr_exp(r19291, r19288, MPFR_RNDN); mpfr_sub(r19292, r19290, r19291, MPFR_RNDN); mpfr_mul(r19293, r19286, r19292, MPFR_RNDN); return mpfr_get_d(r19293, MPFR_RNDN); } static mpfr_t r19294, r19295, r19296, r19297, r19298, r19299, r19300, r19301, r19302, r19303, r19304, r19305, r19306, r19307, r19308, r19309, r19310, r19311, r19312; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19294, "1/60", 10, MPFR_RNDN); mpfr_init(r19295); mpfr_init_set_str(r19296, "5", 10, MPFR_RNDN); mpfr_init(r19297); mpfr_init(r19298); mpfr_init_set_str(r19299, "2", 10, MPFR_RNDN); mpfr_init(r19300); mpfr_init_set_str(r19301, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19302, "3", 10, MPFR_RNDN); mpfr_init(r19303); mpfr_init(r19304); mpfr_init(r19305); mpfr_init(r19306); mpfr_init(r19307); mpfr_init(r19308); mpfr_init(r19309); mpfr_init_set_str(r19310, "0.5", 10, MPFR_RNDN); mpfr_init(r19311); mpfr_init(r19312); } double f_fm(double re, double im) { ; mpfr_set_d(r19295, im, MPFR_RNDN); ; mpfr_pow(r19297, r19295, r19296, MPFR_RNDN); mpfr_mul(r19298, r19294, r19297, MPFR_RNDN); ; mpfr_mul(r19300, r19299, r19295, MPFR_RNDN); ; ; mpfr_pow(r19303, r19295, r19302, MPFR_RNDN); mpfr_mul(r19304, r19301, r19303, MPFR_RNDN); mpfr_add(r19305, r19300, r19304, MPFR_RNDN); mpfr_add(r19306, r19298, r19305, MPFR_RNDN); mpfr_neg(r19307, r19306, MPFR_RNDN); mpfr_set_d(r19308, re, MPFR_RNDN); mpfr_cos(r19309, r19308, MPFR_RNDN); ; mpfr_mul(r19311, r19309, r19310, MPFR_RNDN); mpfr_mul(r19312, r19307, r19311, MPFR_RNDN); return mpfr_get_d(r19312, MPFR_RNDN); } static mpfr_t r19313, r19314, r19315, r19316, r19317, r19318, r19319, r19320, r19321, r19322, r19323, r19324, r19325, r19326, r19327, r19328, r19329, r19330, r19331; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19313, "1/60", 10, MPFR_RNDN); mpfr_init(r19314); mpfr_init_set_str(r19315, "5", 10, MPFR_RNDN); mpfr_init(r19316); mpfr_init(r19317); mpfr_init_set_str(r19318, "2", 10, MPFR_RNDN); mpfr_init(r19319); mpfr_init_set_str(r19320, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19321, "3", 10, MPFR_RNDN); mpfr_init(r19322); mpfr_init(r19323); mpfr_init(r19324); mpfr_init(r19325); mpfr_init(r19326); mpfr_init(r19327); mpfr_init(r19328); mpfr_init_set_str(r19329, "0.5", 10, MPFR_RNDN); mpfr_init(r19330); mpfr_init(r19331); } double f_dm(double re, double im) { ; mpfr_set_d(r19314, im, MPFR_RNDN); ; mpfr_pow(r19316, r19314, r19315, MPFR_RNDN); mpfr_mul(r19317, r19313, r19316, MPFR_RNDN); ; mpfr_mul(r19319, r19318, r19314, MPFR_RNDN); ; ; mpfr_pow(r19322, r19314, r19321, MPFR_RNDN); mpfr_mul(r19323, r19320, r19322, MPFR_RNDN); mpfr_add(r19324, r19319, r19323, MPFR_RNDN); mpfr_add(r19325, r19317, r19324, MPFR_RNDN); mpfr_neg(r19326, r19325, MPFR_RNDN); mpfr_set_d(r19327, re, MPFR_RNDN); mpfr_cos(r19328, r19327, MPFR_RNDN); ; mpfr_mul(r19330, r19328, r19329, MPFR_RNDN); mpfr_mul(r19331, r19326, r19330, MPFR_RNDN); return mpfr_get_d(r19331, MPFR_RNDN); }