#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r18226 = 0.5f; float r18227 = re; float r18228 = cos(r18227); float r18229 = r18226 * r18228; float r18230 = 0.0f; float r18231 = im; float r18232 = r18230 - r18231; float r18233 = exp(r18232); float r18234 = exp(r18231); float r18235 = r18233 - r18234; float r18236 = r18229 * r18235; return r18236; } double f_id(double re, double im) { double r18237 = 0.5; double r18238 = re; double r18239 = cos(r18238); double r18240 = r18237 * r18239; double r18241 = 0.0; double r18242 = im; double r18243 = r18241 - r18242; double r18244 = exp(r18243); double r18245 = exp(r18242); double r18246 = r18244 - r18245; double r18247 = r18240 * r18246; return r18247; } double f_of(float re, float im) { float r18248 = 0.01666666753590107f; float r18249 = im; float r18250 = 5.0f; float r18251 = pow(r18249, r18250); float r18252 = r18248 * r18251; float r18253 = 2.0f; float r18254 = r18253 * r18249; float r18255 = 0.3333333432674408f; float r18256 = 3.0f; float r18257 = pow(r18249, r18256); float r18258 = r18255 * r18257; float r18259 = r18254 + r18258; float r18260 = r18252 + r18259; float r18261 = -r18260; float r18262 = re; float r18263 = cos(r18262); float r18264 = 0.5f; float r18265 = r18263 * r18264; float r18266 = r18261 * r18265; return r18266; } double f_od(double re, double im) { double r18267 = 0.01666666753590107; double r18268 = im; double r18269 = 5.0; double r18270 = pow(r18268, r18269); double r18271 = r18267 * r18270; double r18272 = 2.0; double r18273 = r18272 * r18268; double r18274 = 0.3333333432674408; double r18275 = 3.0; double r18276 = pow(r18268, r18275); double r18277 = r18274 * r18276; double r18278 = r18273 + r18277; double r18279 = r18271 + r18278; double r18280 = -r18279; double r18281 = re; double r18282 = cos(r18281); double r18283 = 0.5; double r18284 = r18282 * r18283; double r18285 = r18280 * r18284; return r18285; } 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 r18286, r18287, r18288, r18289, r18290, r18291, r18292, r18293, r18294, r18295, r18296; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18286, "0.5", 10, MPFR_RNDN); mpfr_init(r18287); mpfr_init(r18288); mpfr_init(r18289); mpfr_init_set_str(r18290, "0", 10, MPFR_RNDN); mpfr_init(r18291); mpfr_init(r18292); mpfr_init(r18293); mpfr_init(r18294); mpfr_init(r18295); mpfr_init(r18296); } double f_im(double re, double im) { ; mpfr_set_d(r18287, re, MPFR_RNDN); mpfr_cos(r18288, r18287, MPFR_RNDN); mpfr_mul(r18289, r18286, r18288, MPFR_RNDN); ; mpfr_set_d(r18291, im, MPFR_RNDN); mpfr_sub(r18292, r18290, r18291, MPFR_RNDN); mpfr_exp(r18293, r18292, MPFR_RNDN); mpfr_exp(r18294, r18291, MPFR_RNDN); mpfr_sub(r18295, r18293, r18294, MPFR_RNDN); mpfr_mul(r18296, r18289, r18295, MPFR_RNDN); return mpfr_get_d(r18296, MPFR_RNDN); } static mpfr_t r18297, r18298, r18299, r18300, r18301, r18302, r18303, r18304, r18305, r18306, r18307, r18308, r18309, r18310, r18311, r18312, r18313, r18314, r18315; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18297, "1/60", 10, MPFR_RNDN); mpfr_init(r18298); mpfr_init_set_str(r18299, "5", 10, MPFR_RNDN); mpfr_init(r18300); mpfr_init(r18301); mpfr_init_set_str(r18302, "2", 10, MPFR_RNDN); mpfr_init(r18303); mpfr_init_set_str(r18304, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18305, "3", 10, MPFR_RNDN); mpfr_init(r18306); mpfr_init(r18307); mpfr_init(r18308); mpfr_init(r18309); mpfr_init(r18310); mpfr_init(r18311); mpfr_init(r18312); mpfr_init_set_str(r18313, "0.5", 10, MPFR_RNDN); mpfr_init(r18314); mpfr_init(r18315); } double f_fm(double re, double im) { ; mpfr_set_d(r18298, im, MPFR_RNDN); ; mpfr_pow(r18300, r18298, r18299, MPFR_RNDN); mpfr_mul(r18301, r18297, r18300, MPFR_RNDN); ; mpfr_mul(r18303, r18302, r18298, MPFR_RNDN); ; ; mpfr_pow(r18306, r18298, r18305, MPFR_RNDN); mpfr_mul(r18307, r18304, r18306, MPFR_RNDN); mpfr_add(r18308, r18303, r18307, MPFR_RNDN); mpfr_add(r18309, r18301, r18308, MPFR_RNDN); mpfr_neg(r18310, r18309, MPFR_RNDN); mpfr_set_d(r18311, re, MPFR_RNDN); mpfr_cos(r18312, r18311, MPFR_RNDN); ; mpfr_mul(r18314, r18312, r18313, MPFR_RNDN); mpfr_mul(r18315, r18310, r18314, MPFR_RNDN); return mpfr_get_d(r18315, MPFR_RNDN); } static mpfr_t r18316, r18317, r18318, r18319, r18320, r18321, r18322, r18323, r18324, r18325, r18326, r18327, r18328, r18329, r18330, r18331, r18332, r18333, r18334; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18316, "1/60", 10, MPFR_RNDN); mpfr_init(r18317); mpfr_init_set_str(r18318, "5", 10, MPFR_RNDN); mpfr_init(r18319); mpfr_init(r18320); mpfr_init_set_str(r18321, "2", 10, MPFR_RNDN); mpfr_init(r18322); mpfr_init_set_str(r18323, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r18324, "3", 10, MPFR_RNDN); mpfr_init(r18325); mpfr_init(r18326); mpfr_init(r18327); mpfr_init(r18328); mpfr_init(r18329); mpfr_init(r18330); mpfr_init(r18331); mpfr_init_set_str(r18332, "0.5", 10, MPFR_RNDN); mpfr_init(r18333); mpfr_init(r18334); } double f_dm(double re, double im) { ; mpfr_set_d(r18317, im, MPFR_RNDN); ; mpfr_pow(r18319, r18317, r18318, MPFR_RNDN); mpfr_mul(r18320, r18316, r18319, MPFR_RNDN); ; mpfr_mul(r18322, r18321, r18317, MPFR_RNDN); ; ; mpfr_pow(r18325, r18317, r18324, MPFR_RNDN); mpfr_mul(r18326, r18323, r18325, MPFR_RNDN); mpfr_add(r18327, r18322, r18326, MPFR_RNDN); mpfr_add(r18328, r18320, r18327, MPFR_RNDN); mpfr_neg(r18329, r18328, MPFR_RNDN); mpfr_set_d(r18330, re, MPFR_RNDN); mpfr_cos(r18331, r18330, MPFR_RNDN); ; mpfr_mul(r18333, r18331, r18332, MPFR_RNDN); mpfr_mul(r18334, r18329, r18333, MPFR_RNDN); return mpfr_get_d(r18334, MPFR_RNDN); }