#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19228 = 0.5f; float r19229 = re; float r19230 = cos(r19229); float r19231 = r19228 * r19230; float r19232 = 0.0f; float r19233 = im; float r19234 = r19232 - r19233; float r19235 = exp(r19234); float r19236 = exp(r19233); float r19237 = r19235 - r19236; float r19238 = r19231 * r19237; return r19238; } double f_id(double re, double im) { double r19239 = 0.5; double r19240 = re; double r19241 = cos(r19240); double r19242 = r19239 * r19241; double r19243 = 0.0; double r19244 = im; double r19245 = r19243 - r19244; double r19246 = exp(r19245); double r19247 = exp(r19244); double r19248 = r19246 - r19247; double r19249 = r19242 * r19248; return r19249; } double f_of(float re, float im) { float r19250 = 0.01666666753590107f; float r19251 = im; float r19252 = 5.0f; float r19253 = pow(r19251, r19252); float r19254 = r19250 * r19253; float r19255 = 2.0f; float r19256 = r19255 * r19251; float r19257 = 0.3333333432674408f; float r19258 = 3.0f; float r19259 = pow(r19251, r19258); float r19260 = r19257 * r19259; float r19261 = r19256 + r19260; float r19262 = r19254 + r19261; float r19263 = -r19262; float r19264 = re; float r19265 = cos(r19264); float r19266 = 0.5f; float r19267 = r19265 * r19266; float r19268 = r19263 * r19267; return r19268; } double f_od(double re, double im) { double r19269 = 0.01666666753590107; double r19270 = im; double r19271 = 5.0; double r19272 = pow(r19270, r19271); double r19273 = r19269 * r19272; double r19274 = 2.0; double r19275 = r19274 * r19270; double r19276 = 0.3333333432674408; double r19277 = 3.0; double r19278 = pow(r19270, r19277); double r19279 = r19276 * r19278; double r19280 = r19275 + r19279; double r19281 = r19273 + r19280; double r19282 = -r19281; double r19283 = re; double r19284 = cos(r19283); double r19285 = 0.5; double r19286 = r19284 * r19285; double r19287 = r19282 * r19286; return r19287; } 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 r19288, r19289, r19290, r19291, r19292, r19293, r19294, r19295, r19296, r19297, r19298; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19288, "0.5", 10, MPFR_RNDN); mpfr_init(r19289); mpfr_init(r19290); mpfr_init(r19291); mpfr_init_set_str(r19292, "0", 10, MPFR_RNDN); mpfr_init(r19293); mpfr_init(r19294); mpfr_init(r19295); mpfr_init(r19296); mpfr_init(r19297); mpfr_init(r19298); } double f_im(double re, double im) { ; mpfr_set_d(r19289, re, MPFR_RNDN); mpfr_cos(r19290, r19289, MPFR_RNDN); mpfr_mul(r19291, r19288, r19290, MPFR_RNDN); ; mpfr_set_d(r19293, im, MPFR_RNDN); mpfr_sub(r19294, r19292, r19293, MPFR_RNDN); mpfr_exp(r19295, r19294, MPFR_RNDN); mpfr_exp(r19296, r19293, MPFR_RNDN); mpfr_sub(r19297, r19295, r19296, MPFR_RNDN); mpfr_mul(r19298, r19291, r19297, MPFR_RNDN); return mpfr_get_d(r19298, MPFR_RNDN); } static mpfr_t r19299, r19300, r19301, r19302, r19303, r19304, r19305, r19306, r19307, r19308, r19309, r19310, r19311, r19312, r19313, r19314, r19315, r19316, r19317; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19299, "1/60", 10, MPFR_RNDN); mpfr_init(r19300); mpfr_init_set_str(r19301, "5", 10, MPFR_RNDN); mpfr_init(r19302); mpfr_init(r19303); mpfr_init_set_str(r19304, "2", 10, MPFR_RNDN); mpfr_init(r19305); mpfr_init_set_str(r19306, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19307, "3", 10, MPFR_RNDN); mpfr_init(r19308); mpfr_init(r19309); mpfr_init(r19310); mpfr_init(r19311); mpfr_init(r19312); mpfr_init(r19313); mpfr_init(r19314); mpfr_init_set_str(r19315, "0.5", 10, MPFR_RNDN); mpfr_init(r19316); mpfr_init(r19317); } double f_fm(double re, double im) { ; mpfr_set_d(r19300, im, MPFR_RNDN); ; mpfr_pow(r19302, r19300, r19301, MPFR_RNDN); mpfr_mul(r19303, r19299, r19302, MPFR_RNDN); ; mpfr_mul(r19305, r19304, r19300, MPFR_RNDN); ; ; mpfr_pow(r19308, r19300, r19307, MPFR_RNDN); mpfr_mul(r19309, r19306, r19308, MPFR_RNDN); mpfr_add(r19310, r19305, r19309, MPFR_RNDN); mpfr_add(r19311, r19303, r19310, MPFR_RNDN); mpfr_neg(r19312, r19311, MPFR_RNDN); mpfr_set_d(r19313, re, MPFR_RNDN); mpfr_cos(r19314, r19313, MPFR_RNDN); ; mpfr_mul(r19316, r19314, r19315, MPFR_RNDN); mpfr_mul(r19317, r19312, r19316, MPFR_RNDN); return mpfr_get_d(r19317, MPFR_RNDN); } static mpfr_t r19318, r19319, r19320, r19321, r19322, r19323, r19324, r19325, r19326, r19327, r19328, r19329, r19330, r19331, r19332, r19333, r19334, r19335, r19336; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19318, "1/60", 10, MPFR_RNDN); mpfr_init(r19319); mpfr_init_set_str(r19320, "5", 10, MPFR_RNDN); mpfr_init(r19321); mpfr_init(r19322); mpfr_init_set_str(r19323, "2", 10, MPFR_RNDN); mpfr_init(r19324); mpfr_init_set_str(r19325, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19326, "3", 10, MPFR_RNDN); mpfr_init(r19327); mpfr_init(r19328); mpfr_init(r19329); mpfr_init(r19330); mpfr_init(r19331); mpfr_init(r19332); mpfr_init(r19333); mpfr_init_set_str(r19334, "0.5", 10, MPFR_RNDN); mpfr_init(r19335); mpfr_init(r19336); } double f_dm(double re, double im) { ; mpfr_set_d(r19319, im, MPFR_RNDN); ; mpfr_pow(r19321, r19319, r19320, MPFR_RNDN); mpfr_mul(r19322, r19318, r19321, MPFR_RNDN); ; mpfr_mul(r19324, r19323, r19319, MPFR_RNDN); ; ; mpfr_pow(r19327, r19319, r19326, MPFR_RNDN); mpfr_mul(r19328, r19325, r19327, MPFR_RNDN); mpfr_add(r19329, r19324, r19328, MPFR_RNDN); mpfr_add(r19330, r19322, r19329, MPFR_RNDN); mpfr_neg(r19331, r19330, MPFR_RNDN); mpfr_set_d(r19332, re, MPFR_RNDN); mpfr_cos(r19333, r19332, MPFR_RNDN); ; mpfr_mul(r19335, r19333, r19334, MPFR_RNDN); mpfr_mul(r19336, r19331, r19335, MPFR_RNDN); return mpfr_get_d(r19336, MPFR_RNDN); }