#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19195 = 0.5f; float r19196 = re; float r19197 = cos(r19196); float r19198 = r19195 * r19197; float r19199 = 0.0f; float r19200 = im; float r19201 = r19199 - r19200; float r19202 = exp(r19201); float r19203 = exp(r19200); float r19204 = r19202 - r19203; float r19205 = r19198 * r19204; return r19205; } double f_id(double re, double im) { double r19206 = 0.5; double r19207 = re; double r19208 = cos(r19207); double r19209 = r19206 * r19208; double r19210 = 0.0; double r19211 = im; double r19212 = r19210 - r19211; double r19213 = exp(r19212); double r19214 = exp(r19211); double r19215 = r19213 - r19214; double r19216 = r19209 * r19215; return r19216; } double f_of(float re, float im) { float r19217 = 0.016666666666666666f; float r19218 = im; float r19219 = 5.0f; float r19220 = pow(r19218, r19219); float r19221 = r19217 * r19220; float r19222 = 2.0f; float r19223 = r19222 * r19218; float r19224 = 0.3333333333333333f; float r19225 = 3.0f; float r19226 = pow(r19218, r19225); float r19227 = r19224 * r19226; float r19228 = r19223 + r19227; float r19229 = r19221 + r19228; float r19230 = -r19229; float r19231 = re; float r19232 = cos(r19231); float r19233 = 0.5f; float r19234 = r19232 * r19233; float r19235 = r19230 * r19234; return r19235; } double f_od(double re, double im) { double r19236 = 0.016666666666666666; double r19237 = im; double r19238 = 5.0; double r19239 = pow(r19237, r19238); double r19240 = r19236 * r19239; double r19241 = 2.0; double r19242 = r19241 * r19237; double r19243 = 0.3333333333333333; double r19244 = 3.0; double r19245 = pow(r19237, r19244); double r19246 = r19243 * r19245; double r19247 = r19242 + r19246; double r19248 = r19240 + r19247; double r19249 = -r19248; double r19250 = re; double r19251 = cos(r19250); double r19252 = 0.5; double r19253 = r19251 * r19252; double r19254 = r19249 * r19253; return r19254; } 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 r19255, r19256, r19257, r19258, r19259, r19260, r19261, r19262, r19263, r19264, r19265; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19255, "0.5", 10, MPFR_RNDN); mpfr_init(r19256); mpfr_init(r19257); mpfr_init(r19258); mpfr_init_set_str(r19259, "0", 10, MPFR_RNDN); mpfr_init(r19260); mpfr_init(r19261); mpfr_init(r19262); mpfr_init(r19263); mpfr_init(r19264); mpfr_init(r19265); } double f_im(double re, double im) { ; mpfr_set_d(r19256, re, MPFR_RNDN); mpfr_cos(r19257, r19256, MPFR_RNDN); mpfr_mul(r19258, r19255, r19257, MPFR_RNDN); ; mpfr_set_d(r19260, im, MPFR_RNDN); mpfr_sub(r19261, r19259, r19260, MPFR_RNDN); mpfr_exp(r19262, r19261, MPFR_RNDN); mpfr_exp(r19263, r19260, MPFR_RNDN); mpfr_sub(r19264, r19262, r19263, MPFR_RNDN); mpfr_mul(r19265, r19258, r19264, MPFR_RNDN); return mpfr_get_d(r19265, MPFR_RNDN); } static mpfr_t r19266, r19267, r19268, r19269, r19270, r19271, r19272, r19273, r19274, r19275, r19276, r19277, r19278, r19279, r19280, r19281, r19282, r19283, r19284; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19266, "1/60", 10, MPFR_RNDN); mpfr_init(r19267); mpfr_init_set_str(r19268, "5", 10, MPFR_RNDN); mpfr_init(r19269); mpfr_init(r19270); mpfr_init_set_str(r19271, "2", 10, MPFR_RNDN); mpfr_init(r19272); mpfr_init_set_str(r19273, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19274, "3", 10, MPFR_RNDN); mpfr_init(r19275); mpfr_init(r19276); mpfr_init(r19277); mpfr_init(r19278); mpfr_init(r19279); mpfr_init(r19280); mpfr_init(r19281); mpfr_init_set_str(r19282, "0.5", 10, MPFR_RNDN); mpfr_init(r19283); mpfr_init(r19284); } double f_fm(double re, double im) { ; mpfr_set_d(r19267, im, MPFR_RNDN); ; mpfr_pow(r19269, r19267, r19268, MPFR_RNDN); mpfr_mul(r19270, r19266, r19269, MPFR_RNDN); ; mpfr_mul(r19272, r19271, r19267, MPFR_RNDN); ; ; mpfr_pow(r19275, r19267, r19274, MPFR_RNDN); mpfr_mul(r19276, r19273, r19275, MPFR_RNDN); mpfr_add(r19277, r19272, r19276, MPFR_RNDN); mpfr_add(r19278, r19270, r19277, MPFR_RNDN); mpfr_neg(r19279, r19278, MPFR_RNDN); mpfr_set_d(r19280, re, MPFR_RNDN); mpfr_cos(r19281, r19280, MPFR_RNDN); ; mpfr_mul(r19283, r19281, r19282, MPFR_RNDN); mpfr_mul(r19284, r19279, r19283, MPFR_RNDN); return mpfr_get_d(r19284, MPFR_RNDN); } static mpfr_t r19285, r19286, r19287, r19288, r19289, r19290, r19291, r19292, r19293, r19294, r19295, r19296, r19297, r19298, r19299, r19300, r19301, r19302, r19303; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19285, "1/60", 10, MPFR_RNDN); mpfr_init(r19286); mpfr_init_set_str(r19287, "5", 10, MPFR_RNDN); mpfr_init(r19288); mpfr_init(r19289); mpfr_init_set_str(r19290, "2", 10, MPFR_RNDN); mpfr_init(r19291); mpfr_init_set_str(r19292, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19293, "3", 10, MPFR_RNDN); mpfr_init(r19294); mpfr_init(r19295); mpfr_init(r19296); mpfr_init(r19297); mpfr_init(r19298); mpfr_init(r19299); mpfr_init(r19300); mpfr_init_set_str(r19301, "0.5", 10, MPFR_RNDN); mpfr_init(r19302); mpfr_init(r19303); } double f_dm(double re, double im) { ; mpfr_set_d(r19286, im, MPFR_RNDN); ; mpfr_pow(r19288, r19286, r19287, MPFR_RNDN); mpfr_mul(r19289, r19285, r19288, MPFR_RNDN); ; mpfr_mul(r19291, r19290, r19286, MPFR_RNDN); ; ; mpfr_pow(r19294, r19286, r19293, MPFR_RNDN); mpfr_mul(r19295, r19292, r19294, MPFR_RNDN); mpfr_add(r19296, r19291, r19295, MPFR_RNDN); mpfr_add(r19297, r19289, r19296, MPFR_RNDN); mpfr_neg(r19298, r19297, MPFR_RNDN); mpfr_set_d(r19299, re, MPFR_RNDN); mpfr_cos(r19300, r19299, MPFR_RNDN); ; mpfr_mul(r19302, r19300, r19301, MPFR_RNDN); mpfr_mul(r19303, r19298, r19302, MPFR_RNDN); return mpfr_get_d(r19303, MPFR_RNDN); }