#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r18175 = 0.5f; float r18176 = re; float r18177 = cos(r18176); float r18178 = r18175 * r18177; float r18179 = 0.0f; float r18180 = im; float r18181 = r18179 - r18180; float r18182 = exp(r18181); float r18183 = exp(r18180); float r18184 = r18182 - r18183; float r18185 = r18178 * r18184; return r18185; } double f_id(double re, double im) { double r18186 = 0.5; double r18187 = re; double r18188 = cos(r18187); double r18189 = r18186 * r18188; double r18190 = 0.0; double r18191 = im; double r18192 = r18190 - r18191; double r18193 = exp(r18192); double r18194 = exp(r18191); double r18195 = r18193 - r18194; double r18196 = r18189 * r18195; return r18196; } double f_of(float re, float im) { float r18197 = 0.01666666753590107f; float r18198 = im; float r18199 = 5.0f; float r18200 = pow(r18198, r18199); float r18201 = r18197 * r18200; float r18202 = 2.0f; float r18203 = r18202 * r18198; float r18204 = 0.3333333432674408f; float r18205 = 3.0f; float r18206 = pow(r18198, r18205); float r18207 = r18204 * r18206; float r18208 = r18203 + r18207; float r18209 = r18201 + r18208; float r18210 = -r18209; float r18211 = re; float r18212 = cos(r18211); float r18213 = 0.5f; float r18214 = r18212 * r18213; float r18215 = r18210 * r18214; return r18215; } double f_od(double re, double im) { double r18216 = 0.01666666753590107; double r18217 = im; double r18218 = 5.0; double r18219 = pow(r18217, r18218); double r18220 = r18216 * r18219; double r18221 = 2.0; double r18222 = r18221 * r18217; double r18223 = 0.3333333432674408; double r18224 = 3.0; double r18225 = pow(r18217, r18224); double r18226 = r18223 * r18225; double r18227 = r18222 + r18226; double r18228 = r18220 + r18227; double r18229 = -r18228; double r18230 = re; double r18231 = cos(r18230); double r18232 = 0.5; double r18233 = r18231 * r18232; double r18234 = r18229 * r18233; return r18234; } 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 r18235, r18236, r18237, r18238, r18239, r18240, r18241, r18242, r18243, r18244, r18245; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18235, "0.5", 10, MPFR_RNDN); mpfr_init(r18236); mpfr_init(r18237); mpfr_init(r18238); mpfr_init_set_str(r18239, "0", 10, MPFR_RNDN); mpfr_init(r18240); mpfr_init(r18241); mpfr_init(r18242); mpfr_init(r18243); mpfr_init(r18244); mpfr_init(r18245); } double f_im(double re, double im) { ; mpfr_set_d(r18236, re, MPFR_RNDN); mpfr_cos(r18237, r18236, MPFR_RNDN); mpfr_mul(r18238, r18235, r18237, MPFR_RNDN); ; mpfr_set_d(r18240, im, MPFR_RNDN); mpfr_sub(r18241, r18239, r18240, MPFR_RNDN); mpfr_exp(r18242, r18241, MPFR_RNDN); mpfr_exp(r18243, r18240, MPFR_RNDN); mpfr_sub(r18244, r18242, r18243, MPFR_RNDN); mpfr_mul(r18245, r18238, r18244, MPFR_RNDN); return mpfr_get_d(r18245, MPFR_RNDN); } static mpfr_t r18246, r18247, r18248, r18249, r18250, r18251, r18252, r18253, r18254, r18255, r18256, r18257, r18258, r18259, r18260, r18261, r18262, r18263, r18264; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18246, "1/60", 10, MPFR_RNDN); mpfr_init(r18247); mpfr_init_set_str(r18248, "5", 10, MPFR_RNDN); mpfr_init(r18249); mpfr_init(r18250); mpfr_init_set_str(r18251, "2", 10, MPFR_RNDN); mpfr_init(r18252); mpfr_init_set_str(r18253, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18254, "3", 10, MPFR_RNDN); mpfr_init(r18255); mpfr_init(r18256); mpfr_init(r18257); mpfr_init(r18258); mpfr_init(r18259); mpfr_init(r18260); mpfr_init(r18261); mpfr_init_set_str(r18262, "0.5", 10, MPFR_RNDN); mpfr_init(r18263); mpfr_init(r18264); } double f_fm(double re, double im) { ; mpfr_set_d(r18247, im, MPFR_RNDN); ; mpfr_pow(r18249, r18247, r18248, MPFR_RNDN); mpfr_mul(r18250, r18246, r18249, MPFR_RNDN); ; mpfr_mul(r18252, r18251, r18247, MPFR_RNDN); ; ; mpfr_pow(r18255, r18247, r18254, MPFR_RNDN); mpfr_mul(r18256, r18253, r18255, MPFR_RNDN); mpfr_add(r18257, r18252, r18256, MPFR_RNDN); mpfr_add(r18258, r18250, r18257, MPFR_RNDN); mpfr_neg(r18259, r18258, MPFR_RNDN); mpfr_set_d(r18260, re, MPFR_RNDN); mpfr_cos(r18261, r18260, MPFR_RNDN); ; mpfr_mul(r18263, r18261, r18262, MPFR_RNDN); mpfr_mul(r18264, r18259, r18263, MPFR_RNDN); return mpfr_get_d(r18264, MPFR_RNDN); } static mpfr_t r18265, r18266, r18267, r18268, r18269, r18270, r18271, r18272, r18273, r18274, r18275, r18276, r18277, r18278, r18279, r18280, r18281, r18282, r18283; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18265, "1/60", 10, MPFR_RNDN); mpfr_init(r18266); mpfr_init_set_str(r18267, "5", 10, MPFR_RNDN); mpfr_init(r18268); mpfr_init(r18269); mpfr_init_set_str(r18270, "2", 10, MPFR_RNDN); mpfr_init(r18271); mpfr_init_set_str(r18272, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18273, "3", 10, MPFR_RNDN); mpfr_init(r18274); mpfr_init(r18275); mpfr_init(r18276); mpfr_init(r18277); mpfr_init(r18278); mpfr_init(r18279); mpfr_init(r18280); mpfr_init_set_str(r18281, "0.5", 10, MPFR_RNDN); mpfr_init(r18282); mpfr_init(r18283); } double f_dm(double re, double im) { ; mpfr_set_d(r18266, im, MPFR_RNDN); ; mpfr_pow(r18268, r18266, r18267, MPFR_RNDN); mpfr_mul(r18269, r18265, r18268, MPFR_RNDN); ; mpfr_mul(r18271, r18270, r18266, MPFR_RNDN); ; ; mpfr_pow(r18274, r18266, r18273, MPFR_RNDN); mpfr_mul(r18275, r18272, r18274, MPFR_RNDN); mpfr_add(r18276, r18271, r18275, MPFR_RNDN); mpfr_add(r18277, r18269, r18276, MPFR_RNDN); mpfr_neg(r18278, r18277, MPFR_RNDN); mpfr_set_d(r18279, re, MPFR_RNDN); mpfr_cos(r18280, r18279, MPFR_RNDN); ; mpfr_mul(r18282, r18280, r18281, MPFR_RNDN); mpfr_mul(r18283, r18278, r18282, MPFR_RNDN); return mpfr_get_d(r18283, MPFR_RNDN); }