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