#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r25100 = 0.5; float r25101 = re; float r25102 = cos(r25101); float r25103 = r25100 * r25102; float r25104 = 0; float r25105 = im; float r25106 = r25104 - r25105; float r25107 = exp(r25106); float r25108 = exp(r25105); float r25109 = r25107 - r25108; float r25110 = r25103 * r25109; return r25110; } double f_id(double re, double im) { double r25111 = 0.5; double r25112 = re; double r25113 = cos(r25112); double r25114 = r25111 * r25113; double r25115 = 0; double r25116 = im; double r25117 = r25115 - r25116; double r25118 = exp(r25117); double r25119 = exp(r25116); double r25120 = r25118 - r25119; double r25121 = r25114 * r25120; return r25121; } double f_of(float re, float im) { float r25122 = 0.5; float r25123 = re; float r25124 = cos(r25123); float r25125 = r25122 * r25124; float r25126 = 1/60; float r25127 = im; float r25128 = 5; float r25129 = pow(r25127, r25128); float r25130 = r25126 * r25129; float r25131 = 2; float r25132 = r25131 * r25127; float r25133 = 1/3; float r25134 = 3; float r25135 = pow(r25127, r25134); float r25136 = r25133 * r25135; float r25137 = r25132 + r25136; float r25138 = r25130 + r25137; float r25139 = -r25138; float r25140 = r25125 * r25139; return r25140; } double f_od(double re, double im) { double r25141 = 0.5; double r25142 = re; double r25143 = cos(r25142); double r25144 = r25141 * r25143; double r25145 = 1/60; double r25146 = im; double r25147 = 5; double r25148 = pow(r25146, r25147); double r25149 = r25145 * r25148; double r25150 = 2; double r25151 = r25150 * r25146; double r25152 = 1/3; double r25153 = 3; double r25154 = pow(r25146, r25153); double r25155 = r25152 * r25154; double r25156 = r25151 + r25155; double r25157 = r25149 + r25156; double r25158 = -r25157; double r25159 = r25144 * r25158; return r25159; } 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 r25160, r25161, r25162, r25163, r25164, r25165, r25166, r25167, r25168, r25169, r25170; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r25160, "0.5", 10, MPFR_RNDN); mpfr_init(r25161); mpfr_init(r25162); mpfr_init(r25163); mpfr_init_set_str(r25164, "0", 10, MPFR_RNDN); mpfr_init(r25165); mpfr_init(r25166); mpfr_init(r25167); mpfr_init(r25168); mpfr_init(r25169); mpfr_init(r25170); } double f_im(double re, double im) { ; mpfr_set_d(r25161, re, MPFR_RNDN); mpfr_cos(r25162, r25161, MPFR_RNDN); mpfr_mul(r25163, r25160, r25162, MPFR_RNDN); ; mpfr_set_d(r25165, im, MPFR_RNDN); mpfr_sub(r25166, r25164, r25165, MPFR_RNDN); mpfr_exp(r25167, r25166, MPFR_RNDN); mpfr_exp(r25168, r25165, MPFR_RNDN); mpfr_sub(r25169, r25167, r25168, MPFR_RNDN); mpfr_mul(r25170, r25163, r25169, MPFR_RNDN); return mpfr_get_d(r25170, MPFR_RNDN); } static mpfr_t r25171, r25172, r25173, r25174, r25175, r25176, r25177, r25178, r25179, r25180, r25181, r25182, r25183, r25184, r25185, r25186, r25187, r25188, r25189; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r25171, "0.5", 10, MPFR_RNDN); mpfr_init(r25172); mpfr_init(r25173); mpfr_init(r25174); mpfr_init_set_str(r25175, "1/60", 10, MPFR_RNDN); mpfr_init(r25176); mpfr_init_set_str(r25177, "5", 10, MPFR_RNDN); mpfr_init(r25178); mpfr_init(r25179); mpfr_init_set_str(r25180, "2", 10, MPFR_RNDN); mpfr_init(r25181); mpfr_init_set_str(r25182, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r25183, "3", 10, MPFR_RNDN); mpfr_init(r25184); mpfr_init(r25185); mpfr_init(r25186); mpfr_init(r25187); mpfr_init(r25188); mpfr_init(r25189); } double f_fm(double re, double im) { ; mpfr_set_d(r25172, re, MPFR_RNDN); mpfr_cos(r25173, r25172, MPFR_RNDN); mpfr_mul(r25174, r25171, r25173, MPFR_RNDN); ; mpfr_set_d(r25176, im, MPFR_RNDN); ; mpfr_pow(r25178, r25176, r25177, MPFR_RNDN); mpfr_mul(r25179, r25175, r25178, MPFR_RNDN); ; mpfr_mul(r25181, r25180, r25176, MPFR_RNDN); ; ; mpfr_pow(r25184, r25176, r25183, MPFR_RNDN); mpfr_mul(r25185, r25182, r25184, MPFR_RNDN); mpfr_add(r25186, r25181, r25185, MPFR_RNDN); mpfr_add(r25187, r25179, r25186, MPFR_RNDN); mpfr_neg(r25188, r25187, MPFR_RNDN); mpfr_mul(r25189, r25174, r25188, MPFR_RNDN); return mpfr_get_d(r25189, MPFR_RNDN); } static mpfr_t r25190, r25191, r25192, r25193, r25194, r25195, r25196, r25197, r25198, r25199, r25200, r25201, r25202, r25203, r25204, r25205, r25206, r25207, r25208; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r25190, "0.5", 10, MPFR_RNDN); mpfr_init(r25191); mpfr_init(r25192); mpfr_init(r25193); mpfr_init_set_str(r25194, "1/60", 10, MPFR_RNDN); mpfr_init(r25195); mpfr_init_set_str(r25196, "5", 10, MPFR_RNDN); mpfr_init(r25197); mpfr_init(r25198); mpfr_init_set_str(r25199, "2", 10, MPFR_RNDN); mpfr_init(r25200); mpfr_init_set_str(r25201, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r25202, "3", 10, MPFR_RNDN); mpfr_init(r25203); mpfr_init(r25204); mpfr_init(r25205); mpfr_init(r25206); mpfr_init(r25207); mpfr_init(r25208); } double f_dm(double re, double im) { ; mpfr_set_d(r25191, re, MPFR_RNDN); mpfr_cos(r25192, r25191, MPFR_RNDN); mpfr_mul(r25193, r25190, r25192, MPFR_RNDN); ; mpfr_set_d(r25195, im, MPFR_RNDN); ; mpfr_pow(r25197, r25195, r25196, MPFR_RNDN); mpfr_mul(r25198, r25194, r25197, MPFR_RNDN); ; mpfr_mul(r25200, r25199, r25195, MPFR_RNDN); ; ; mpfr_pow(r25203, r25195, r25202, MPFR_RNDN); mpfr_mul(r25204, r25201, r25203, MPFR_RNDN); mpfr_add(r25205, r25200, r25204, MPFR_RNDN); mpfr_add(r25206, r25198, r25205, MPFR_RNDN); mpfr_neg(r25207, r25206, MPFR_RNDN); mpfr_mul(r25208, r25193, r25207, MPFR_RNDN); return mpfr_get_d(r25208, MPFR_RNDN); }