#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19351 = 0.5f; float r19352 = re; float r19353 = cos(r19352); float r19354 = r19351 * r19353; float r19355 = 0.0f; float r19356 = im; float r19357 = r19355 - r19356; float r19358 = exp(r19357); float r19359 = exp(r19356); float r19360 = r19358 - r19359; float r19361 = r19354 * r19360; return r19361; } double f_id(double re, double im) { double r19362 = 0.5; double r19363 = re; double r19364 = cos(r19363); double r19365 = r19362 * r19364; double r19366 = 0.0; double r19367 = im; double r19368 = r19366 - r19367; double r19369 = exp(r19368); double r19370 = exp(r19367); double r19371 = r19369 - r19370; double r19372 = r19365 * r19371; return r19372; } double f_of(float re, float im) { float r19373 = 0.016666666666666666f; float r19374 = im; float r19375 = 5.0f; float r19376 = pow(r19374, r19375); float r19377 = r19373 * r19376; float r19378 = 2.0f; float r19379 = r19378 * r19374; float r19380 = 0.3333333333333333f; float r19381 = 3.0f; float r19382 = pow(r19374, r19381); float r19383 = r19380 * r19382; float r19384 = r19379 + r19383; float r19385 = r19377 + r19384; float r19386 = -r19385; float r19387 = re; float r19388 = cos(r19387); float r19389 = 0.5f; float r19390 = r19388 * r19389; float r19391 = r19386 * r19390; return r19391; } double f_od(double re, double im) { double r19392 = 0.016666666666666666; double r19393 = im; double r19394 = 5.0; double r19395 = pow(r19393, r19394); double r19396 = r19392 * r19395; double r19397 = 2.0; double r19398 = r19397 * r19393; double r19399 = 0.3333333333333333; double r19400 = 3.0; double r19401 = pow(r19393, r19400); double r19402 = r19399 * r19401; double r19403 = r19398 + r19402; double r19404 = r19396 + r19403; double r19405 = -r19404; double r19406 = re; double r19407 = cos(r19406); double r19408 = 0.5; double r19409 = r19407 * r19408; double r19410 = r19405 * r19409; return r19410; } 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 r19411, r19412, r19413, r19414, r19415, r19416, r19417, r19418, r19419, r19420, r19421; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19411, "0.5", 10, MPFR_RNDN); mpfr_init(r19412); mpfr_init(r19413); mpfr_init(r19414); mpfr_init_set_str(r19415, "0", 10, MPFR_RNDN); mpfr_init(r19416); mpfr_init(r19417); mpfr_init(r19418); mpfr_init(r19419); mpfr_init(r19420); mpfr_init(r19421); } double f_im(double re, double im) { ; mpfr_set_d(r19412, re, MPFR_RNDN); mpfr_cos(r19413, r19412, MPFR_RNDN); mpfr_mul(r19414, r19411, r19413, MPFR_RNDN); ; mpfr_set_d(r19416, im, MPFR_RNDN); mpfr_sub(r19417, r19415, r19416, MPFR_RNDN); mpfr_exp(r19418, r19417, MPFR_RNDN); mpfr_exp(r19419, r19416, MPFR_RNDN); mpfr_sub(r19420, r19418, r19419, MPFR_RNDN); mpfr_mul(r19421, r19414, r19420, MPFR_RNDN); return mpfr_get_d(r19421, MPFR_RNDN); } static mpfr_t r19422, r19423, r19424, r19425, r19426, r19427, r19428, r19429, r19430, r19431, r19432, r19433, r19434, r19435, r19436, r19437, r19438, r19439, r19440; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19422, "1/60", 10, MPFR_RNDN); mpfr_init(r19423); mpfr_init_set_str(r19424, "5", 10, MPFR_RNDN); mpfr_init(r19425); mpfr_init(r19426); mpfr_init_set_str(r19427, "2", 10, MPFR_RNDN); mpfr_init(r19428); mpfr_init_set_str(r19429, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19430, "3", 10, MPFR_RNDN); mpfr_init(r19431); mpfr_init(r19432); mpfr_init(r19433); mpfr_init(r19434); mpfr_init(r19435); mpfr_init(r19436); mpfr_init(r19437); mpfr_init_set_str(r19438, "0.5", 10, MPFR_RNDN); mpfr_init(r19439); mpfr_init(r19440); } double f_fm(double re, double im) { ; mpfr_set_d(r19423, im, MPFR_RNDN); ; mpfr_pow(r19425, r19423, r19424, MPFR_RNDN); mpfr_mul(r19426, r19422, r19425, MPFR_RNDN); ; mpfr_mul(r19428, r19427, r19423, MPFR_RNDN); ; ; mpfr_pow(r19431, r19423, r19430, MPFR_RNDN); mpfr_mul(r19432, r19429, r19431, MPFR_RNDN); mpfr_add(r19433, r19428, r19432, MPFR_RNDN); mpfr_add(r19434, r19426, r19433, MPFR_RNDN); mpfr_neg(r19435, r19434, MPFR_RNDN); mpfr_set_d(r19436, re, MPFR_RNDN); mpfr_cos(r19437, r19436, MPFR_RNDN); ; mpfr_mul(r19439, r19437, r19438, MPFR_RNDN); mpfr_mul(r19440, r19435, r19439, MPFR_RNDN); return mpfr_get_d(r19440, MPFR_RNDN); } static mpfr_t r19441, r19442, r19443, r19444, r19445, r19446, r19447, r19448, r19449, r19450, r19451, r19452, r19453, r19454, r19455, r19456, r19457, r19458, r19459; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19441, "1/60", 10, MPFR_RNDN); mpfr_init(r19442); mpfr_init_set_str(r19443, "5", 10, MPFR_RNDN); mpfr_init(r19444); mpfr_init(r19445); mpfr_init_set_str(r19446, "2", 10, MPFR_RNDN); mpfr_init(r19447); mpfr_init_set_str(r19448, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19449, "3", 10, MPFR_RNDN); mpfr_init(r19450); mpfr_init(r19451); mpfr_init(r19452); mpfr_init(r19453); mpfr_init(r19454); mpfr_init(r19455); mpfr_init(r19456); mpfr_init_set_str(r19457, "0.5", 10, MPFR_RNDN); mpfr_init(r19458); mpfr_init(r19459); } double f_dm(double re, double im) { ; mpfr_set_d(r19442, im, MPFR_RNDN); ; mpfr_pow(r19444, r19442, r19443, MPFR_RNDN); mpfr_mul(r19445, r19441, r19444, MPFR_RNDN); ; mpfr_mul(r19447, r19446, r19442, MPFR_RNDN); ; ; mpfr_pow(r19450, r19442, r19449, MPFR_RNDN); mpfr_mul(r19451, r19448, r19450, MPFR_RNDN); mpfr_add(r19452, r19447, r19451, MPFR_RNDN); mpfr_add(r19453, r19445, r19452, MPFR_RNDN); mpfr_neg(r19454, r19453, MPFR_RNDN); mpfr_set_d(r19455, re, MPFR_RNDN); mpfr_cos(r19456, r19455, MPFR_RNDN); ; mpfr_mul(r19458, r19456, r19457, MPFR_RNDN); mpfr_mul(r19459, r19454, r19458, MPFR_RNDN); return mpfr_get_d(r19459, MPFR_RNDN); }