#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19303 = 0.5f; float r19304 = re; float r19305 = cos(r19304); float r19306 = r19303 * r19305; float r19307 = 0.0f; float r19308 = im; float r19309 = r19307 - r19308; float r19310 = exp(r19309); float r19311 = exp(r19308); float r19312 = r19310 - r19311; float r19313 = r19306 * r19312; return r19313; } double f_id(double re, double im) { double r19314 = 0.5; double r19315 = re; double r19316 = cos(r19315); double r19317 = r19314 * r19316; double r19318 = 0.0; double r19319 = im; double r19320 = r19318 - r19319; double r19321 = exp(r19320); double r19322 = exp(r19319); double r19323 = r19321 - r19322; double r19324 = r19317 * r19323; return r19324; } double f_of(float re, float im) { float r19325 = 0.016666666666666666f; float r19326 = im; float r19327 = 5.0f; float r19328 = pow(r19326, r19327); float r19329 = r19325 * r19328; float r19330 = 2.0f; float r19331 = r19330 * r19326; float r19332 = 0.3333333333333333f; float r19333 = 3.0f; float r19334 = pow(r19326, r19333); float r19335 = r19332 * r19334; float r19336 = r19331 + r19335; float r19337 = r19329 + r19336; float r19338 = -r19337; float r19339 = re; float r19340 = cos(r19339); float r19341 = 0.5f; float r19342 = r19340 * r19341; float r19343 = r19338 * r19342; return r19343; } double f_od(double re, double im) { double r19344 = 0.016666666666666666; double r19345 = im; double r19346 = 5.0; double r19347 = pow(r19345, r19346); double r19348 = r19344 * r19347; double r19349 = 2.0; double r19350 = r19349 * r19345; double r19351 = 0.3333333333333333; double r19352 = 3.0; double r19353 = pow(r19345, r19352); double r19354 = r19351 * r19353; double r19355 = r19350 + r19354; double r19356 = r19348 + r19355; double r19357 = -r19356; double r19358 = re; double r19359 = cos(r19358); double r19360 = 0.5; double r19361 = r19359 * r19360; double r19362 = r19357 * r19361; return r19362; } 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 r19363, r19364, r19365, r19366, r19367, r19368, r19369, r19370, r19371, r19372, r19373; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19363, "0.5", 10, MPFR_RNDN); mpfr_init(r19364); mpfr_init(r19365); mpfr_init(r19366); mpfr_init_set_str(r19367, "0", 10, MPFR_RNDN); mpfr_init(r19368); mpfr_init(r19369); mpfr_init(r19370); mpfr_init(r19371); mpfr_init(r19372); mpfr_init(r19373); } double f_im(double re, double im) { ; mpfr_set_d(r19364, re, MPFR_RNDN); mpfr_cos(r19365, r19364, MPFR_RNDN); mpfr_mul(r19366, r19363, r19365, MPFR_RNDN); ; mpfr_set_d(r19368, im, MPFR_RNDN); mpfr_sub(r19369, r19367, r19368, MPFR_RNDN); mpfr_exp(r19370, r19369, MPFR_RNDN); mpfr_exp(r19371, r19368, MPFR_RNDN); mpfr_sub(r19372, r19370, r19371, MPFR_RNDN); mpfr_mul(r19373, r19366, r19372, MPFR_RNDN); return mpfr_get_d(r19373, MPFR_RNDN); } static mpfr_t r19374, r19375, r19376, r19377, r19378, r19379, r19380, r19381, r19382, r19383, r19384, r19385, r19386, r19387, r19388, r19389, r19390, r19391, r19392; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19374, "1/60", 10, MPFR_RNDN); mpfr_init(r19375); mpfr_init_set_str(r19376, "5", 10, MPFR_RNDN); mpfr_init(r19377); mpfr_init(r19378); mpfr_init_set_str(r19379, "2", 10, MPFR_RNDN); mpfr_init(r19380); mpfr_init_set_str(r19381, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19382, "3", 10, MPFR_RNDN); mpfr_init(r19383); mpfr_init(r19384); mpfr_init(r19385); mpfr_init(r19386); mpfr_init(r19387); mpfr_init(r19388); mpfr_init(r19389); mpfr_init_set_str(r19390, "0.5", 10, MPFR_RNDN); mpfr_init(r19391); mpfr_init(r19392); } double f_fm(double re, double im) { ; mpfr_set_d(r19375, im, MPFR_RNDN); ; mpfr_pow(r19377, r19375, r19376, MPFR_RNDN); mpfr_mul(r19378, r19374, r19377, MPFR_RNDN); ; mpfr_mul(r19380, r19379, r19375, MPFR_RNDN); ; ; mpfr_pow(r19383, r19375, r19382, MPFR_RNDN); mpfr_mul(r19384, r19381, r19383, MPFR_RNDN); mpfr_add(r19385, r19380, r19384, MPFR_RNDN); mpfr_add(r19386, r19378, r19385, MPFR_RNDN); mpfr_neg(r19387, r19386, MPFR_RNDN); mpfr_set_d(r19388, re, MPFR_RNDN); mpfr_cos(r19389, r19388, MPFR_RNDN); ; mpfr_mul(r19391, r19389, r19390, MPFR_RNDN); mpfr_mul(r19392, r19387, r19391, MPFR_RNDN); return mpfr_get_d(r19392, MPFR_RNDN); } static mpfr_t r19393, r19394, r19395, r19396, r19397, r19398, r19399, r19400, r19401, r19402, r19403, r19404, r19405, r19406, r19407, r19408, r19409, r19410, r19411; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19393, "1/60", 10, MPFR_RNDN); mpfr_init(r19394); mpfr_init_set_str(r19395, "5", 10, MPFR_RNDN); mpfr_init(r19396); mpfr_init(r19397); mpfr_init_set_str(r19398, "2", 10, MPFR_RNDN); mpfr_init(r19399); mpfr_init_set_str(r19400, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19401, "3", 10, MPFR_RNDN); mpfr_init(r19402); mpfr_init(r19403); mpfr_init(r19404); mpfr_init(r19405); mpfr_init(r19406); mpfr_init(r19407); mpfr_init(r19408); mpfr_init_set_str(r19409, "0.5", 10, MPFR_RNDN); mpfr_init(r19410); mpfr_init(r19411); } double f_dm(double re, double im) { ; mpfr_set_d(r19394, im, MPFR_RNDN); ; mpfr_pow(r19396, r19394, r19395, MPFR_RNDN); mpfr_mul(r19397, r19393, r19396, MPFR_RNDN); ; mpfr_mul(r19399, r19398, r19394, MPFR_RNDN); ; ; mpfr_pow(r19402, r19394, r19401, MPFR_RNDN); mpfr_mul(r19403, r19400, r19402, MPFR_RNDN); mpfr_add(r19404, r19399, r19403, MPFR_RNDN); mpfr_add(r19405, r19397, r19404, MPFR_RNDN); mpfr_neg(r19406, r19405, MPFR_RNDN); mpfr_set_d(r19407, re, MPFR_RNDN); mpfr_cos(r19408, r19407, MPFR_RNDN); ; mpfr_mul(r19410, r19408, r19409, MPFR_RNDN); mpfr_mul(r19411, r19406, r19410, MPFR_RNDN); return mpfr_get_d(r19411, MPFR_RNDN); }