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