#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19287 = 0.5f; float r19288 = re; float r19289 = cos(r19288); float r19290 = r19287 * r19289; float r19291 = 0.0f; float r19292 = im; float r19293 = r19291 - r19292; float r19294 = exp(r19293); float r19295 = exp(r19292); float r19296 = r19294 - r19295; float r19297 = r19290 * r19296; return r19297; } double f_id(double re, double im) { double r19298 = 0.5; double r19299 = re; double r19300 = cos(r19299); double r19301 = r19298 * r19300; double r19302 = 0.0; double r19303 = im; double r19304 = r19302 - r19303; double r19305 = exp(r19304); double r19306 = exp(r19303); double r19307 = r19305 - r19306; double r19308 = r19301 * r19307; return r19308; } double f_of(float re, float im) { float r19309 = 0.016666666666666666f; float r19310 = im; float r19311 = 5.0f; float r19312 = pow(r19310, r19311); float r19313 = r19310 * (r19310 * r19310); float r19314 = 0.3333333333333333f; float r19315 = 2.0f; float r19316 = r19315 * r19310; float r19317 = fma(r19313, r19314, r19316); float r19318 = fma(r19309, r19312, r19317); float r19319 = 1.0f; float r19320 = r19318 / r19319; float r19321 = re; float r19322 = cos(r19321); float r19323 = 0.5f; float r19324 = -r19323; float r19325 = r19322 * r19324; float r19326 = r19320 * r19325; return r19326; } double f_od(double re, double im) { double r19327 = 0.016666666666666666; double r19328 = im; double r19329 = 5.0; double r19330 = pow(r19328, r19329); double r19331 = r19328 * (r19328 * r19328); double r19332 = 0.3333333333333333; double r19333 = 2.0; double r19334 = r19333 * r19328; double r19335 = fma(r19331, r19332, r19334); double r19336 = fma(r19327, r19330, r19335); double r19337 = 1.0; double r19338 = r19336 / r19337; double r19339 = re; double r19340 = cos(r19339); double r19341 = 0.5; double r19342 = -r19341; double r19343 = r19340 * r19342; double r19344 = r19338 * r19343; return r19344; } 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 r19345, r19346, r19347, r19348, r19349, r19350, r19351, r19352, r19353, r19354, r19355; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19345, "0.5", 10, MPFR_RNDN); mpfr_init(r19346); mpfr_init(r19347); mpfr_init(r19348); mpfr_init_set_str(r19349, "0", 10, MPFR_RNDN); mpfr_init(r19350); mpfr_init(r19351); mpfr_init(r19352); mpfr_init(r19353); mpfr_init(r19354); mpfr_init(r19355); } double f_im(double re, double im) { ; mpfr_set_d(r19346, re, MPFR_RNDN); mpfr_cos(r19347, r19346, MPFR_RNDN); mpfr_mul(r19348, r19345, r19347, MPFR_RNDN); ; mpfr_set_d(r19350, im, MPFR_RNDN); mpfr_sub(r19351, r19349, r19350, MPFR_RNDN); mpfr_exp(r19352, r19351, MPFR_RNDN); mpfr_exp(r19353, r19350, MPFR_RNDN); mpfr_sub(r19354, r19352, r19353, MPFR_RNDN); mpfr_mul(r19355, r19348, r19354, MPFR_RNDN); return mpfr_get_d(r19355, MPFR_RNDN); } static mpfr_t r19356, r19357, r19358, r19359, r19360, r19361, r19362, r19363, r19364, r19365, r19366, r19367, r19368, r19369, r19370, r19371, r19372, r19373; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19356, "1/60", 10, MPFR_RNDN); mpfr_init(r19357); mpfr_init_set_str(r19358, "5", 10, MPFR_RNDN); mpfr_init(r19359); mpfr_init(r19360); mpfr_init_set_str(r19361, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19362, "2", 10, MPFR_RNDN); mpfr_init(r19363); mpfr_init(r19364); mpfr_init(r19365); mpfr_init_set_str(r19366, "1", 10, MPFR_RNDN); mpfr_init(r19367); mpfr_init(r19368); mpfr_init(r19369); mpfr_init_set_str(r19370, "0.5", 10, MPFR_RNDN); mpfr_init(r19371); mpfr_init(r19372); mpfr_init(r19373); } double f_fm(double re, double im) { ; mpfr_set_d(r19357, im, MPFR_RNDN); ; mpfr_pow(r19359, r19357, r19358, MPFR_RNDN); mpfr_mul(r19360, r19357, r19357, MPFR_RNDN); mpfr_mul(r19360, r19360, r19357, MPFR_RNDN); ; ; mpfr_mul(r19363, r19362, r19357, MPFR_RNDN); mpfr_fma(r19364, r19360, r19361, r19363, MPFR_RNDN); mpfr_fma(r19365, r19356, r19359, r19364, MPFR_RNDN); ; mpfr_div(r19367, r19365, r19366, MPFR_RNDN); mpfr_set_d(r19368, re, MPFR_RNDN); mpfr_cos(r19369, r19368, MPFR_RNDN); ; mpfr_neg(r19371, r19370, MPFR_RNDN); mpfr_mul(r19372, r19369, r19371, MPFR_RNDN); mpfr_mul(r19373, r19367, 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; void setup_mpfr_f_dm() { 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, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19380, "2", 10, MPFR_RNDN); mpfr_init(r19381); mpfr_init(r19382); mpfr_init(r19383); mpfr_init_set_str(r19384, "1", 10, MPFR_RNDN); mpfr_init(r19385); mpfr_init(r19386); mpfr_init(r19387); mpfr_init_set_str(r19388, "0.5", 10, MPFR_RNDN); mpfr_init(r19389); mpfr_init(r19390); mpfr_init(r19391); } double f_dm(double re, double im) { ; mpfr_set_d(r19375, im, MPFR_RNDN); ; mpfr_pow(r19377, r19375, r19376, MPFR_RNDN); mpfr_mul(r19378, r19375, r19375, MPFR_RNDN); mpfr_mul(r19378, r19378, r19375, MPFR_RNDN); ; ; mpfr_mul(r19381, r19380, r19375, MPFR_RNDN); mpfr_fma(r19382, r19378, r19379, r19381, MPFR_RNDN); mpfr_fma(r19383, r19374, r19377, r19382, MPFR_RNDN); ; mpfr_div(r19385, r19383, r19384, MPFR_RNDN); mpfr_set_d(r19386, re, MPFR_RNDN); mpfr_cos(r19387, r19386, MPFR_RNDN); ; mpfr_neg(r19389, r19388, MPFR_RNDN); mpfr_mul(r19390, r19387, r19389, MPFR_RNDN); mpfr_mul(r19391, r19385, r19390, MPFR_RNDN); return mpfr_get_d(r19391, MPFR_RNDN); }