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