#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r20403 = 0.5f; float r20404 = re; float r20405 = cos(r20404); float r20406 = r20403 * r20405; float r20407 = 0.0f; float r20408 = im; float r20409 = r20407 - r20408; float r20410 = exp(r20409); float r20411 = exp(r20408); float r20412 = r20410 - r20411; float r20413 = r20406 * r20412; return r20413; } double f_id(double re, double im) { double r20414 = 0.5; double r20415 = re; double r20416 = cos(r20415); double r20417 = r20414 * r20416; double r20418 = 0.0; double r20419 = im; double r20420 = r20418 - r20419; double r20421 = exp(r20420); double r20422 = exp(r20419); double r20423 = r20421 - r20422; double r20424 = r20417 * r20423; return r20424; } double f_of(float re, float im) { float r20425 = im; float r20426 = r20425 * (r20425 * r20425); float r20427 = 0.3333333333333333f; float r20428 = r20426 * r20427; float r20429 = 0.016666666666666666f; float r20430 = 5.0f; float r20431 = pow(r20425, r20430); float r20432 = r20429 * r20431; float r20433 = 2.0f; float r20434 = r20425 * r20433; float r20435 = r20432 + r20434; float r20436 = r20428 + r20435; float r20437 = 1.0f; float r20438 = r20436 / r20437; float r20439 = 0.5f; float r20440 = -r20439; float r20441 = re; float r20442 = cos(r20441); float r20443 = r20440 * r20442; float r20444 = r20438 * r20443; return r20444; } double f_od(double re, double im) { double r20445 = im; double r20446 = r20445 * (r20445 * r20445); double r20447 = 0.3333333333333333; double r20448 = r20446 * r20447; double r20449 = 0.016666666666666666; double r20450 = 5.0; double r20451 = pow(r20445, r20450); double r20452 = r20449 * r20451; double r20453 = 2.0; double r20454 = r20445 * r20453; double r20455 = r20452 + r20454; double r20456 = r20448 + r20455; double r20457 = 1.0; double r20458 = r20456 / r20457; double r20459 = 0.5; double r20460 = -r20459; double r20461 = re; double r20462 = cos(r20461); double r20463 = r20460 * r20462; double r20464 = r20458 * r20463; return r20464; } 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 r20465, r20466, r20467, r20468, r20469, r20470, r20471, r20472, r20473, r20474, r20475; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r20465, "0.5", 10, MPFR_RNDN); mpfr_init(r20466); mpfr_init(r20467); mpfr_init(r20468); mpfr_init_set_str(r20469, "0", 10, MPFR_RNDN); mpfr_init(r20470); mpfr_init(r20471); mpfr_init(r20472); mpfr_init(r20473); mpfr_init(r20474); mpfr_init(r20475); } double f_im(double re, double im) { ; mpfr_set_d(r20466, re, MPFR_RNDN); mpfr_cos(r20467, r20466, MPFR_RNDN); mpfr_mul(r20468, r20465, r20467, MPFR_RNDN); ; mpfr_set_d(r20470, im, MPFR_RNDN); mpfr_sub(r20471, r20469, r20470, MPFR_RNDN); mpfr_exp(r20472, r20471, MPFR_RNDN); mpfr_exp(r20473, r20470, MPFR_RNDN); mpfr_sub(r20474, r20472, r20473, MPFR_RNDN); mpfr_mul(r20475, r20468, r20474, MPFR_RNDN); return mpfr_get_d(r20475, MPFR_RNDN); } static mpfr_t r20476, r20477, r20478, r20479, r20480, r20481, r20482, r20483, r20484, r20485, r20486, r20487, r20488, r20489, r20490, r20491, r20492, r20493, r20494, r20495; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r20476); mpfr_init(r20477); mpfr_init_set_str(r20478, "1/3", 10, MPFR_RNDN); mpfr_init(r20479); mpfr_init_set_str(r20480, "1/60", 10, MPFR_RNDN); mpfr_init_set_str(r20481, "5", 10, MPFR_RNDN); mpfr_init(r20482); mpfr_init(r20483); mpfr_init_set_str(r20484, "2", 10, MPFR_RNDN); mpfr_init(r20485); mpfr_init(r20486); mpfr_init(r20487); mpfr_init_set_str(r20488, "1", 10, MPFR_RNDN); mpfr_init(r20489); mpfr_init_set_str(r20490, "0.5", 10, MPFR_RNDN); mpfr_init(r20491); mpfr_init(r20492); mpfr_init(r20493); mpfr_init(r20494); mpfr_init(r20495); } double f_fm(double re, double im) { mpfr_set_d(r20476, im, MPFR_RNDN); mpfr_mul(r20477, r20476, r20476, MPFR_RNDN); mpfr_mul(r20477, r20477, r20476, MPFR_RNDN); ; mpfr_mul(r20479, r20477, r20478, MPFR_RNDN); ; ; mpfr_pow(r20482, r20476, r20481, MPFR_RNDN); mpfr_mul(r20483, r20480, r20482, MPFR_RNDN); ; mpfr_mul(r20485, r20476, r20484, MPFR_RNDN); mpfr_add(r20486, r20483, r20485, MPFR_RNDN); mpfr_add(r20487, r20479, r20486, MPFR_RNDN); ; mpfr_div(r20489, r20487, r20488, MPFR_RNDN); ; mpfr_neg(r20491, r20490, MPFR_RNDN); mpfr_set_d(r20492, re, MPFR_RNDN); mpfr_cos(r20493, r20492, MPFR_RNDN); mpfr_mul(r20494, r20491, r20493, MPFR_RNDN); mpfr_mul(r20495, r20489, r20494, MPFR_RNDN); return mpfr_get_d(r20495, MPFR_RNDN); } static mpfr_t r20496, r20497, r20498, r20499, r20500, r20501, r20502, r20503, r20504, r20505, r20506, r20507, r20508, r20509, r20510, r20511, r20512, r20513, r20514, r20515; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r20496); mpfr_init(r20497); mpfr_init_set_str(r20498, "1/3", 10, MPFR_RNDN); mpfr_init(r20499); mpfr_init_set_str(r20500, "1/60", 10, MPFR_RNDN); mpfr_init_set_str(r20501, "5", 10, MPFR_RNDN); mpfr_init(r20502); mpfr_init(r20503); mpfr_init_set_str(r20504, "2", 10, MPFR_RNDN); mpfr_init(r20505); mpfr_init(r20506); mpfr_init(r20507); mpfr_init_set_str(r20508, "1", 10, MPFR_RNDN); mpfr_init(r20509); mpfr_init_set_str(r20510, "0.5", 10, MPFR_RNDN); mpfr_init(r20511); mpfr_init(r20512); mpfr_init(r20513); mpfr_init(r20514); mpfr_init(r20515); } double f_dm(double re, double im) { mpfr_set_d(r20496, im, MPFR_RNDN); mpfr_mul(r20497, r20496, r20496, MPFR_RNDN); mpfr_mul(r20497, r20497, r20496, MPFR_RNDN); ; mpfr_mul(r20499, r20497, r20498, MPFR_RNDN); ; ; mpfr_pow(r20502, r20496, r20501, MPFR_RNDN); mpfr_mul(r20503, r20500, r20502, MPFR_RNDN); ; mpfr_mul(r20505, r20496, r20504, MPFR_RNDN); mpfr_add(r20506, r20503, r20505, MPFR_RNDN); mpfr_add(r20507, r20499, r20506, MPFR_RNDN); ; mpfr_div(r20509, r20507, r20508, MPFR_RNDN); ; mpfr_neg(r20511, r20510, MPFR_RNDN); mpfr_set_d(r20512, re, MPFR_RNDN); mpfr_cos(r20513, r20512, MPFR_RNDN); mpfr_mul(r20514, r20511, r20513, MPFR_RNDN); mpfr_mul(r20515, r20509, r20514, MPFR_RNDN); return mpfr_get_d(r20515, MPFR_RNDN); }