#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r20037 = 0.5f; float r20038 = re; float r20039 = sin(r20038); float r20040 = r20037 * r20039; float r20041 = im; float r20042 = -r20041; float r20043 = exp(r20042); float r20044 = exp(r20041); float r20045 = r20043 - r20044; float r20046 = r20040 * r20045; return r20046; } double f_id(double re, double im) { double r20047 = 0.5; double r20048 = re; double r20049 = sin(r20048); double r20050 = r20047 * r20049; double r20051 = im; double r20052 = -r20051; double r20053 = exp(r20052); double r20054 = exp(r20051); double r20055 = r20053 - r20054; double r20056 = r20050 * r20055; return r20056; } double f_of(float re, float im) { float r20057 = 0.5f; float r20058 = re; float r20059 = sin(r20058); float r20060 = r20057 * r20059; float r20061 = 0.016666666666666666f; float r20062 = im; float r20063 = 5.0f; float r20064 = pow(r20062, r20063); float r20065 = r20061 * r20064; float r20066 = 2.0f; float r20067 = r20066 * r20062; float r20068 = 0.3333333333333333f; float r20069 = 3.0f; float r20070 = pow(r20062, r20069); float r20071 = r20068 * r20070; float r20072 = r20067 + r20071; float r20073 = r20065 + r20072; float r20074 = -r20073; float r20075 = r20060 * r20074; return r20075; } double f_od(double re, double im) { double r20076 = 0.5; double r20077 = re; double r20078 = sin(r20077); double r20079 = r20076 * r20078; double r20080 = 0.016666666666666666; double r20081 = im; double r20082 = 5.0; double r20083 = pow(r20081, r20082); double r20084 = r20080 * r20083; double r20085 = 2.0; double r20086 = r20085 * r20081; double r20087 = 0.3333333333333333; double r20088 = 3.0; double r20089 = pow(r20081, r20088); double r20090 = r20087 * r20089; double r20091 = r20086 + r20090; double r20092 = r20084 + r20091; double r20093 = -r20092; double r20094 = r20079 * r20093; return r20094; } 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 r20095, r20096, r20097, r20098, r20099, r20100, r20101, r20102, r20103, r20104; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r20095, "0.5", 10, MPFR_RNDN); mpfr_init(r20096); mpfr_init(r20097); mpfr_init(r20098); mpfr_init(r20099); mpfr_init(r20100); mpfr_init(r20101); mpfr_init(r20102); mpfr_init(r20103); mpfr_init(r20104); } double f_im(double re, double im) { ; mpfr_set_d(r20096, re, MPFR_RNDN); mpfr_sin(r20097, r20096, MPFR_RNDN); mpfr_mul(r20098, r20095, r20097, MPFR_RNDN); mpfr_set_d(r20099, im, MPFR_RNDN); mpfr_neg(r20100, r20099, MPFR_RNDN); mpfr_exp(r20101, r20100, MPFR_RNDN); mpfr_exp(r20102, r20099, MPFR_RNDN); mpfr_sub(r20103, r20101, r20102, MPFR_RNDN); mpfr_mul(r20104, r20098, r20103, MPFR_RNDN); return mpfr_get_d(r20104, MPFR_RNDN); } static mpfr_t r20105, r20106, r20107, r20108, r20109, r20110, r20111, r20112, r20113, r20114, r20115, r20116, r20117, r20118, r20119, r20120, r20121, r20122, r20123; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20105, "0.5", 10, MPFR_RNDN); mpfr_init(r20106); mpfr_init(r20107); mpfr_init(r20108); mpfr_init_set_str(r20109, "1/60", 10, MPFR_RNDN); mpfr_init(r20110); mpfr_init_set_str(r20111, "5", 10, MPFR_RNDN); mpfr_init(r20112); mpfr_init(r20113); mpfr_init_set_str(r20114, "2", 10, MPFR_RNDN); mpfr_init(r20115); mpfr_init_set_str(r20116, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20117, "3", 10, MPFR_RNDN); mpfr_init(r20118); mpfr_init(r20119); mpfr_init(r20120); mpfr_init(r20121); mpfr_init(r20122); mpfr_init(r20123); } double f_fm(double re, double im) { ; mpfr_set_d(r20106, re, MPFR_RNDN); mpfr_sin(r20107, r20106, MPFR_RNDN); mpfr_mul(r20108, r20105, r20107, MPFR_RNDN); ; mpfr_set_d(r20110, im, MPFR_RNDN); ; mpfr_pow(r20112, r20110, r20111, MPFR_RNDN); mpfr_mul(r20113, r20109, r20112, MPFR_RNDN); ; mpfr_mul(r20115, r20114, r20110, MPFR_RNDN); ; ; mpfr_pow(r20118, r20110, r20117, MPFR_RNDN); mpfr_mul(r20119, r20116, r20118, MPFR_RNDN); mpfr_add(r20120, r20115, r20119, MPFR_RNDN); mpfr_add(r20121, r20113, r20120, MPFR_RNDN); mpfr_neg(r20122, r20121, MPFR_RNDN); mpfr_mul(r20123, r20108, r20122, MPFR_RNDN); return mpfr_get_d(r20123, MPFR_RNDN); } static mpfr_t r20124, r20125, r20126, r20127, r20128, r20129, r20130, r20131, r20132, r20133, r20134, r20135, r20136, r20137, r20138, r20139, r20140, r20141, r20142; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20124, "0.5", 10, MPFR_RNDN); mpfr_init(r20125); mpfr_init(r20126); mpfr_init(r20127); mpfr_init_set_str(r20128, "1/60", 10, MPFR_RNDN); mpfr_init(r20129); mpfr_init_set_str(r20130, "5", 10, MPFR_RNDN); mpfr_init(r20131); mpfr_init(r20132); mpfr_init_set_str(r20133, "2", 10, MPFR_RNDN); mpfr_init(r20134); mpfr_init_set_str(r20135, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20136, "3", 10, MPFR_RNDN); mpfr_init(r20137); mpfr_init(r20138); mpfr_init(r20139); mpfr_init(r20140); mpfr_init(r20141); mpfr_init(r20142); } double f_dm(double re, double im) { ; mpfr_set_d(r20125, re, MPFR_RNDN); mpfr_sin(r20126, r20125, MPFR_RNDN); mpfr_mul(r20127, r20124, r20126, MPFR_RNDN); ; mpfr_set_d(r20129, im, MPFR_RNDN); ; mpfr_pow(r20131, r20129, r20130, MPFR_RNDN); mpfr_mul(r20132, r20128, r20131, MPFR_RNDN); ; mpfr_mul(r20134, r20133, r20129, MPFR_RNDN); ; ; mpfr_pow(r20137, r20129, r20136, MPFR_RNDN); mpfr_mul(r20138, r20135, r20137, MPFR_RNDN); mpfr_add(r20139, r20134, r20138, MPFR_RNDN); mpfr_add(r20140, r20132, r20139, MPFR_RNDN); mpfr_neg(r20141, r20140, MPFR_RNDN); mpfr_mul(r20142, r20127, r20141, MPFR_RNDN); return mpfr_get_d(r20142, MPFR_RNDN); }