#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r20057 = 0.5f; float r20058 = re; float r20059 = sin(r20058); float r20060 = r20057 * r20059; float r20061 = im; float r20062 = -r20061; float r20063 = exp(r20062); float r20064 = exp(r20061); float r20065 = r20063 - r20064; float r20066 = r20060 * r20065; return r20066; } double f_id(double re, double im) { double r20067 = 0.5; double r20068 = re; double r20069 = sin(r20068); double r20070 = r20067 * r20069; double r20071 = im; double r20072 = -r20071; double r20073 = exp(r20072); double r20074 = exp(r20071); double r20075 = r20073 - r20074; double r20076 = r20070 * r20075; return r20076; } double f_of(float re, float im) { float r20077 = 0.5f; float r20078 = re; float r20079 = sin(r20078); float r20080 = r20077 * r20079; float r20081 = 0.01666666753590107f; float r20082 = im; float r20083 = 5.0f; float r20084 = pow(r20082, r20083); float r20085 = r20081 * r20084; float r20086 = 2.0f; float r20087 = r20086 * r20082; float r20088 = 0.3333333432674408f; float r20089 = 3.0f; float r20090 = pow(r20082, r20089); float r20091 = r20088 * r20090; float r20092 = r20087 + r20091; float r20093 = r20085 + r20092; float r20094 = -r20093; float r20095 = r20080 * r20094; return r20095; } double f_od(double re, double im) { double r20096 = 0.5; double r20097 = re; double r20098 = sin(r20097); double r20099 = r20096 * r20098; double r20100 = 0.01666666753590107; double r20101 = im; double r20102 = 5.0; double r20103 = pow(r20101, r20102); double r20104 = r20100 * r20103; double r20105 = 2.0; double r20106 = r20105 * r20101; double r20107 = 0.3333333432674408; double r20108 = 3.0; double r20109 = pow(r20101, r20108); double r20110 = r20107 * r20109; double r20111 = r20106 + r20110; double r20112 = r20104 + r20111; double r20113 = -r20112; double r20114 = r20099 * r20113; return r20114; } 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 r20115, r20116, r20117, r20118, r20119, r20120, r20121, r20122, r20123, r20124; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r20115, "0.5", 10, MPFR_RNDN); mpfr_init(r20116); mpfr_init(r20117); mpfr_init(r20118); mpfr_init(r20119); mpfr_init(r20120); mpfr_init(r20121); mpfr_init(r20122); mpfr_init(r20123); mpfr_init(r20124); } double f_im(double re, double im) { ; mpfr_set_d(r20116, re, MPFR_RNDN); mpfr_sin(r20117, r20116, MPFR_RNDN); mpfr_mul(r20118, r20115, r20117, MPFR_RNDN); mpfr_set_d(r20119, im, MPFR_RNDN); mpfr_neg(r20120, r20119, MPFR_RNDN); mpfr_exp(r20121, r20120, MPFR_RNDN); mpfr_exp(r20122, r20119, MPFR_RNDN); mpfr_sub(r20123, r20121, r20122, MPFR_RNDN); mpfr_mul(r20124, r20118, r20123, MPFR_RNDN); return mpfr_get_d(r20124, MPFR_RNDN); } static mpfr_t r20125, r20126, r20127, r20128, r20129, r20130, r20131, r20132, r20133, r20134, r20135, r20136, r20137, r20138, r20139, r20140, r20141, r20142, r20143; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20125, "0.5", 10, MPFR_RNDN); mpfr_init(r20126); mpfr_init(r20127); mpfr_init(r20128); mpfr_init_set_str(r20129, "1/60", 10, MPFR_RNDN); mpfr_init(r20130); mpfr_init_set_str(r20131, "5", 10, MPFR_RNDN); mpfr_init(r20132); mpfr_init(r20133); mpfr_init_set_str(r20134, "2", 10, MPFR_RNDN); mpfr_init(r20135); mpfr_init_set_str(r20136, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20137, "3", 10, MPFR_RNDN); mpfr_init(r20138); mpfr_init(r20139); mpfr_init(r20140); mpfr_init(r20141); mpfr_init(r20142); mpfr_init(r20143); } double f_fm(double re, double im) { ; mpfr_set_d(r20126, re, MPFR_RNDN); mpfr_sin(r20127, r20126, MPFR_RNDN); mpfr_mul(r20128, r20125, r20127, MPFR_RNDN); ; mpfr_set_d(r20130, im, MPFR_RNDN); ; mpfr_pow(r20132, r20130, r20131, MPFR_RNDN); mpfr_mul(r20133, r20129, r20132, MPFR_RNDN); ; mpfr_mul(r20135, r20134, r20130, MPFR_RNDN); ; ; mpfr_pow(r20138, r20130, r20137, MPFR_RNDN); mpfr_mul(r20139, r20136, r20138, MPFR_RNDN); mpfr_add(r20140, r20135, r20139, MPFR_RNDN); mpfr_add(r20141, r20133, r20140, MPFR_RNDN); mpfr_neg(r20142, r20141, MPFR_RNDN); mpfr_mul(r20143, r20128, r20142, MPFR_RNDN); return mpfr_get_d(r20143, MPFR_RNDN); } static mpfr_t r20144, r20145, r20146, r20147, r20148, r20149, r20150, r20151, r20152, r20153, r20154, r20155, r20156, r20157, r20158, r20159, r20160, r20161, r20162; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20144, "0.5", 10, MPFR_RNDN); mpfr_init(r20145); mpfr_init(r20146); mpfr_init(r20147); mpfr_init_set_str(r20148, "1/60", 10, MPFR_RNDN); mpfr_init(r20149); mpfr_init_set_str(r20150, "5", 10, MPFR_RNDN); mpfr_init(r20151); mpfr_init(r20152); mpfr_init_set_str(r20153, "2", 10, MPFR_RNDN); mpfr_init(r20154); mpfr_init_set_str(r20155, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20156, "3", 10, MPFR_RNDN); mpfr_init(r20157); mpfr_init(r20158); mpfr_init(r20159); mpfr_init(r20160); mpfr_init(r20161); mpfr_init(r20162); } double f_dm(double re, double im) { ; mpfr_set_d(r20145, re, MPFR_RNDN); mpfr_sin(r20146, r20145, MPFR_RNDN); mpfr_mul(r20147, r20144, r20146, MPFR_RNDN); ; mpfr_set_d(r20149, im, MPFR_RNDN); ; mpfr_pow(r20151, r20149, r20150, MPFR_RNDN); mpfr_mul(r20152, r20148, r20151, MPFR_RNDN); ; mpfr_mul(r20154, r20153, r20149, MPFR_RNDN); ; ; mpfr_pow(r20157, r20149, r20156, MPFR_RNDN); mpfr_mul(r20158, r20155, r20157, MPFR_RNDN); mpfr_add(r20159, r20154, r20158, MPFR_RNDN); mpfr_add(r20160, r20152, r20159, MPFR_RNDN); mpfr_neg(r20161, r20160, MPFR_RNDN); mpfr_mul(r20162, r20147, r20161, MPFR_RNDN); return mpfr_get_d(r20162, MPFR_RNDN); }