#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r18203 = 0.5f; float r18204 = re; float r18205 = cos(r18204); float r18206 = r18203 * r18205; float r18207 = 0.0f; float r18208 = im; float r18209 = r18207 - r18208; float r18210 = exp(r18209); float r18211 = exp(r18208); float r18212 = r18210 - r18211; float r18213 = r18206 * r18212; return r18213; } double f_id(double re, double im) { double r18214 = 0.5; double r18215 = re; double r18216 = cos(r18215); double r18217 = r18214 * r18216; double r18218 = 0.0; double r18219 = im; double r18220 = r18218 - r18219; double r18221 = exp(r18220); double r18222 = exp(r18219); double r18223 = r18221 - r18222; double r18224 = r18217 * r18223; return r18224; } double f_of(float re, float im) { float r18225 = 0.01666666753590107f; float r18226 = im; float r18227 = 5.0f; float r18228 = pow(r18226, r18227); float r18229 = r18225 * r18228; float r18230 = 2.0f; float r18231 = r18230 * r18226; float r18232 = 0.3333333432674408f; float r18233 = 3.0f; float r18234 = pow(r18226, r18233); float r18235 = r18232 * r18234; float r18236 = r18231 + r18235; float r18237 = r18229 + r18236; float r18238 = -r18237; float r18239 = re; float r18240 = cos(r18239); float r18241 = 0.5f; float r18242 = r18240 * r18241; float r18243 = r18238 * r18242; return r18243; } double f_od(double re, double im) { double r18244 = 0.01666666753590107; double r18245 = im; double r18246 = 5.0; double r18247 = pow(r18245, r18246); double r18248 = r18244 * r18247; double r18249 = 2.0; double r18250 = r18249 * r18245; double r18251 = 0.3333333432674408; double r18252 = 3.0; double r18253 = pow(r18245, r18252); double r18254 = r18251 * r18253; double r18255 = r18250 + r18254; double r18256 = r18248 + r18255; double r18257 = -r18256; double r18258 = re; double r18259 = cos(r18258); double r18260 = 0.5; double r18261 = r18259 * r18260; double r18262 = r18257 * r18261; return r18262; } 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 r18263, r18264, r18265, r18266, r18267, r18268, r18269, r18270, r18271, r18272, r18273; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18263, "0.5", 10, MPFR_RNDN); mpfr_init(r18264); mpfr_init(r18265); mpfr_init(r18266); mpfr_init_set_str(r18267, "0", 10, MPFR_RNDN); mpfr_init(r18268); mpfr_init(r18269); mpfr_init(r18270); mpfr_init(r18271); mpfr_init(r18272); mpfr_init(r18273); } double f_im(double re, double im) { ; mpfr_set_d(r18264, re, MPFR_RNDN); mpfr_cos(r18265, r18264, MPFR_RNDN); mpfr_mul(r18266, r18263, r18265, MPFR_RNDN); ; mpfr_set_d(r18268, im, MPFR_RNDN); mpfr_sub(r18269, r18267, r18268, MPFR_RNDN); mpfr_exp(r18270, r18269, MPFR_RNDN); mpfr_exp(r18271, r18268, MPFR_RNDN); mpfr_sub(r18272, r18270, r18271, MPFR_RNDN); mpfr_mul(r18273, r18266, r18272, MPFR_RNDN); return mpfr_get_d(r18273, MPFR_RNDN); } static mpfr_t r18274, r18275, r18276, r18277, r18278, r18279, r18280, r18281, r18282, r18283, r18284, r18285, r18286, r18287, r18288, r18289, r18290, r18291, r18292; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18274, "1/60", 10, MPFR_RNDN); mpfr_init(r18275); mpfr_init_set_str(r18276, "5", 10, MPFR_RNDN); mpfr_init(r18277); mpfr_init(r18278); mpfr_init_set_str(r18279, "2", 10, MPFR_RNDN); mpfr_init(r18280); mpfr_init_set_str(r18281, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18282, "3", 10, MPFR_RNDN); mpfr_init(r18283); mpfr_init(r18284); mpfr_init(r18285); mpfr_init(r18286); mpfr_init(r18287); mpfr_init(r18288); mpfr_init(r18289); mpfr_init_set_str(r18290, "0.5", 10, MPFR_RNDN); mpfr_init(r18291); mpfr_init(r18292); } double f_fm(double re, double im) { ; mpfr_set_d(r18275, im, MPFR_RNDN); ; mpfr_pow(r18277, r18275, r18276, MPFR_RNDN); mpfr_mul(r18278, r18274, r18277, MPFR_RNDN); ; mpfr_mul(r18280, r18279, r18275, MPFR_RNDN); ; ; mpfr_pow(r18283, r18275, r18282, MPFR_RNDN); mpfr_mul(r18284, r18281, r18283, MPFR_RNDN); mpfr_add(r18285, r18280, r18284, MPFR_RNDN); mpfr_add(r18286, r18278, r18285, MPFR_RNDN); mpfr_neg(r18287, r18286, MPFR_RNDN); mpfr_set_d(r18288, re, MPFR_RNDN); mpfr_cos(r18289, r18288, MPFR_RNDN); ; mpfr_mul(r18291, r18289, r18290, MPFR_RNDN); mpfr_mul(r18292, r18287, r18291, MPFR_RNDN); return mpfr_get_d(r18292, MPFR_RNDN); } static mpfr_t r18293, r18294, r18295, r18296, r18297, r18298, r18299, r18300, r18301, r18302, r18303, r18304, r18305, r18306, r18307, r18308, r18309, r18310, r18311; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18293, "1/60", 10, MPFR_RNDN); mpfr_init(r18294); mpfr_init_set_str(r18295, "5", 10, MPFR_RNDN); mpfr_init(r18296); mpfr_init(r18297); mpfr_init_set_str(r18298, "2", 10, MPFR_RNDN); mpfr_init(r18299); mpfr_init_set_str(r18300, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18301, "3", 10, MPFR_RNDN); mpfr_init(r18302); mpfr_init(r18303); mpfr_init(r18304); mpfr_init(r18305); mpfr_init(r18306); mpfr_init(r18307); mpfr_init(r18308); mpfr_init_set_str(r18309, "0.5", 10, MPFR_RNDN); mpfr_init(r18310); mpfr_init(r18311); } double f_dm(double re, double im) { ; mpfr_set_d(r18294, im, MPFR_RNDN); ; mpfr_pow(r18296, r18294, r18295, MPFR_RNDN); mpfr_mul(r18297, r18293, r18296, MPFR_RNDN); ; mpfr_mul(r18299, r18298, r18294, MPFR_RNDN); ; ; mpfr_pow(r18302, r18294, r18301, MPFR_RNDN); mpfr_mul(r18303, r18300, r18302, MPFR_RNDN); mpfr_add(r18304, r18299, r18303, MPFR_RNDN); mpfr_add(r18305, r18297, r18304, MPFR_RNDN); mpfr_neg(r18306, r18305, MPFR_RNDN); mpfr_set_d(r18307, re, MPFR_RNDN); mpfr_cos(r18308, r18307, MPFR_RNDN); ; mpfr_mul(r18310, r18308, r18309, MPFR_RNDN); mpfr_mul(r18311, r18306, r18310, MPFR_RNDN); return mpfr_get_d(r18311, MPFR_RNDN); }