#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19315 = 0.5f; float r19316 = re; float r19317 = cos(r19316); float r19318 = r19315 * r19317; float r19319 = 0.0f; float r19320 = im; float r19321 = r19319 - r19320; float r19322 = exp(r19321); float r19323 = exp(r19320); float r19324 = r19322 - r19323; float r19325 = r19318 * r19324; return r19325; } double f_id(double re, double im) { double r19326 = 0.5; double r19327 = re; double r19328 = cos(r19327); double r19329 = r19326 * r19328; double r19330 = 0.0; double r19331 = im; double r19332 = r19330 - r19331; double r19333 = exp(r19332); double r19334 = exp(r19331); double r19335 = r19333 - r19334; double r19336 = r19329 * r19335; return r19336; } double f_of(float re, float im) { float r19337 = 0.01666666753590107f; float r19338 = im; float r19339 = 5.0f; float r19340 = pow(r19338, r19339); float r19341 = r19337 * r19340; float r19342 = 2.0f; float r19343 = r19342 * r19338; float r19344 = 0.3333333432674408f; float r19345 = 3.0f; float r19346 = pow(r19338, r19345); float r19347 = r19344 * r19346; float r19348 = r19343 + r19347; float r19349 = r19341 + r19348; float r19350 = -r19349; float r19351 = re; float r19352 = cos(r19351); float r19353 = 0.5f; float r19354 = r19352 * r19353; float r19355 = r19350 * r19354; return r19355; } double f_od(double re, double im) { double r19356 = 0.01666666753590107; double r19357 = im; double r19358 = 5.0; double r19359 = pow(r19357, r19358); double r19360 = r19356 * r19359; double r19361 = 2.0; double r19362 = r19361 * r19357; double r19363 = 0.3333333432674408; double r19364 = 3.0; double r19365 = pow(r19357, r19364); double r19366 = r19363 * r19365; double r19367 = r19362 + r19366; double r19368 = r19360 + r19367; double r19369 = -r19368; double r19370 = re; double r19371 = cos(r19370); double r19372 = 0.5; double r19373 = r19371 * r19372; double r19374 = r19369 * r19373; return r19374; } 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 r19375, r19376, r19377, r19378, r19379, r19380, r19381, r19382, r19383, r19384, r19385; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19375, "0.5", 10, MPFR_RNDN); mpfr_init(r19376); mpfr_init(r19377); mpfr_init(r19378); mpfr_init_set_str(r19379, "0", 10, MPFR_RNDN); mpfr_init(r19380); mpfr_init(r19381); mpfr_init(r19382); mpfr_init(r19383); mpfr_init(r19384); mpfr_init(r19385); } double f_im(double re, double im) { ; mpfr_set_d(r19376, re, MPFR_RNDN); mpfr_cos(r19377, r19376, MPFR_RNDN); mpfr_mul(r19378, r19375, r19377, MPFR_RNDN); ; mpfr_set_d(r19380, im, MPFR_RNDN); mpfr_sub(r19381, r19379, r19380, MPFR_RNDN); mpfr_exp(r19382, r19381, MPFR_RNDN); mpfr_exp(r19383, r19380, MPFR_RNDN); mpfr_sub(r19384, r19382, r19383, MPFR_RNDN); mpfr_mul(r19385, r19378, r19384, MPFR_RNDN); return mpfr_get_d(r19385, MPFR_RNDN); } static mpfr_t r19386, r19387, r19388, r19389, r19390, r19391, r19392, r19393, r19394, r19395, r19396, r19397, r19398, r19399, r19400, r19401, r19402, r19403, r19404; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19386, "1/60", 10, MPFR_RNDN); mpfr_init(r19387); mpfr_init_set_str(r19388, "5", 10, MPFR_RNDN); mpfr_init(r19389); mpfr_init(r19390); mpfr_init_set_str(r19391, "2", 10, MPFR_RNDN); mpfr_init(r19392); mpfr_init_set_str(r19393, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19394, "3", 10, MPFR_RNDN); mpfr_init(r19395); mpfr_init(r19396); mpfr_init(r19397); mpfr_init(r19398); mpfr_init(r19399); mpfr_init(r19400); mpfr_init(r19401); mpfr_init_set_str(r19402, "0.5", 10, MPFR_RNDN); mpfr_init(r19403); mpfr_init(r19404); } double f_fm(double re, double im) { ; mpfr_set_d(r19387, im, MPFR_RNDN); ; mpfr_pow(r19389, r19387, r19388, MPFR_RNDN); mpfr_mul(r19390, r19386, r19389, MPFR_RNDN); ; mpfr_mul(r19392, r19391, r19387, MPFR_RNDN); ; ; mpfr_pow(r19395, r19387, r19394, MPFR_RNDN); mpfr_mul(r19396, r19393, r19395, MPFR_RNDN); mpfr_add(r19397, r19392, r19396, MPFR_RNDN); mpfr_add(r19398, r19390, r19397, MPFR_RNDN); mpfr_neg(r19399, r19398, MPFR_RNDN); mpfr_set_d(r19400, re, MPFR_RNDN); mpfr_cos(r19401, r19400, MPFR_RNDN); ; mpfr_mul(r19403, r19401, r19402, MPFR_RNDN); mpfr_mul(r19404, r19399, r19403, MPFR_RNDN); return mpfr_get_d(r19404, MPFR_RNDN); } static mpfr_t r19405, r19406, r19407, r19408, r19409, r19410, r19411, r19412, r19413, r19414, r19415, r19416, r19417, r19418, r19419, r19420, r19421, r19422, r19423; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19405, "1/60", 10, MPFR_RNDN); mpfr_init(r19406); mpfr_init_set_str(r19407, "5", 10, MPFR_RNDN); mpfr_init(r19408); mpfr_init(r19409); mpfr_init_set_str(r19410, "2", 10, MPFR_RNDN); mpfr_init(r19411); mpfr_init_set_str(r19412, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19413, "3", 10, MPFR_RNDN); mpfr_init(r19414); mpfr_init(r19415); mpfr_init(r19416); mpfr_init(r19417); mpfr_init(r19418); mpfr_init(r19419); mpfr_init(r19420); mpfr_init_set_str(r19421, "0.5", 10, MPFR_RNDN); mpfr_init(r19422); mpfr_init(r19423); } double f_dm(double re, double im) { ; mpfr_set_d(r19406, im, MPFR_RNDN); ; mpfr_pow(r19408, r19406, r19407, MPFR_RNDN); mpfr_mul(r19409, r19405, r19408, MPFR_RNDN); ; mpfr_mul(r19411, r19410, r19406, MPFR_RNDN); ; ; mpfr_pow(r19414, r19406, r19413, MPFR_RNDN); mpfr_mul(r19415, r19412, r19414, MPFR_RNDN); mpfr_add(r19416, r19411, r19415, MPFR_RNDN); mpfr_add(r19417, r19409, r19416, MPFR_RNDN); mpfr_neg(r19418, r19417, MPFR_RNDN); mpfr_set_d(r19419, re, MPFR_RNDN); mpfr_cos(r19420, r19419, MPFR_RNDN); ; mpfr_mul(r19422, r19420, r19421, MPFR_RNDN); mpfr_mul(r19423, r19418, r19422, MPFR_RNDN); return mpfr_get_d(r19423, MPFR_RNDN); }