#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r23094 = 0.5f; float r23095 = re; float r23096 = cos(r23095); float r23097 = r23094 * r23096; float r23098 = 0.0f; float r23099 = im; float r23100 = r23098 - r23099; float r23101 = exp(r23100); float r23102 = exp(r23099); float r23103 = r23101 - r23102; float r23104 = r23097 * r23103; return r23104; } double f_id(double re, double im) { double r23105 = 0.5; double r23106 = re; double r23107 = cos(r23106); double r23108 = r23105 * r23107; double r23109 = 0.0; double r23110 = im; double r23111 = r23109 - r23110; double r23112 = exp(r23111); double r23113 = exp(r23110); double r23114 = r23112 - r23113; double r23115 = r23108 * r23114; return r23115; } double f_of(float re, float im) { float r23116 = im; float r23117 = 1.9183430624335127f; bool r23118 = r23116 <= r23117; float r23119 = 0.016666666666666666f; float r23120 = 5.0f; float r23121 = pow(r23116, r23120); float r23122 = r23119 * r23121; float r23123 = 2.0f; float r23124 = r23123 * r23116; float r23125 = 0.3333333333333333f; float r23126 = 3.0f; float r23127 = pow(r23116, r23126); float r23128 = r23125 * r23127; float r23129 = r23124 + r23128; float r23130 = r23122 + r23129; float r23131 = -r23130; float r23132 = 0.5f; float r23133 = re; float r23134 = cos(r23133); float r23135 = r23132 * r23134; float r23136 = r23131 * r23135; float r23137 = -r23116; float r23138 = exp(r23137); float r23139 = exp(r23116); float r23140 = r23138 - r23139; float r23141 = r23140 * r23135; float r23142 = r23118 ? r23136 : r23141; return r23142; } double f_od(double re, double im) { double r23143 = im; double r23144 = 1.9183430624335127; bool r23145 = r23143 <= r23144; double r23146 = 0.016666666666666666; double r23147 = 5.0; double r23148 = pow(r23143, r23147); double r23149 = r23146 * r23148; double r23150 = 2.0; double r23151 = r23150 * r23143; double r23152 = 0.3333333333333333; double r23153 = 3.0; double r23154 = pow(r23143, r23153); double r23155 = r23152 * r23154; double r23156 = r23151 + r23155; double r23157 = r23149 + r23156; double r23158 = -r23157; double r23159 = 0.5; double r23160 = re; double r23161 = cos(r23160); double r23162 = r23159 * r23161; double r23163 = r23158 * r23162; double r23164 = -r23143; double r23165 = exp(r23164); double r23166 = exp(r23143); double r23167 = r23165 - r23166; double r23168 = r23167 * r23162; double r23169 = r23145 ? r23163 : r23168; return r23169; } 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 r23170, r23171, r23172, r23173, r23174, r23175, r23176, r23177, r23178, r23179, r23180; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r23170, "0.5", 10, MPFR_RNDN); mpfr_init(r23171); mpfr_init(r23172); mpfr_init(r23173); mpfr_init_set_str(r23174, "0", 10, MPFR_RNDN); mpfr_init(r23175); mpfr_init(r23176); mpfr_init(r23177); mpfr_init(r23178); mpfr_init(r23179); mpfr_init(r23180); } double f_im(double re, double im) { ; mpfr_set_d(r23171, re, MPFR_RNDN); mpfr_cos(r23172, r23171, MPFR_RNDN); mpfr_mul(r23173, r23170, r23172, MPFR_RNDN); ; mpfr_set_d(r23175, im, MPFR_RNDN); mpfr_sub(r23176, r23174, r23175, MPFR_RNDN); mpfr_exp(r23177, r23176, MPFR_RNDN); mpfr_exp(r23178, r23175, MPFR_RNDN); mpfr_sub(r23179, r23177, r23178, MPFR_RNDN); mpfr_mul(r23180, r23173, r23179, MPFR_RNDN); return mpfr_get_d(r23180, MPFR_RNDN); } static mpfr_t r23181, r23182, r23183, r23184, r23185, r23186, r23187, r23188, r23189, r23190, r23191, r23192, r23193, r23194, r23195, r23196, r23197, r23198, r23199, r23200, r23201, r23202, r23203, r23204, r23205, r23206, r23207; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r23181); mpfr_init_set_str(r23182, "1.9183430624335127", 10, MPFR_RNDN); mpfr_init(r23183); mpfr_init_set_str(r23184, "1/60", 10, MPFR_RNDN); mpfr_init_set_str(r23185, "5", 10, MPFR_RNDN); mpfr_init(r23186); mpfr_init(r23187); mpfr_init_set_str(r23188, "2", 10, MPFR_RNDN); mpfr_init(r23189); mpfr_init_set_str(r23190, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r23191, "3", 10, MPFR_RNDN); mpfr_init(r23192); mpfr_init(r23193); mpfr_init(r23194); mpfr_init(r23195); mpfr_init(r23196); mpfr_init_set_str(r23197, "0.5", 10, MPFR_RNDN); mpfr_init(r23198); mpfr_init(r23199); mpfr_init(r23200); mpfr_init(r23201); mpfr_init(r23202); mpfr_init(r23203); mpfr_init(r23204); mpfr_init(r23205); mpfr_init(r23206); mpfr_init(r23207); } double f_fm(double re, double im) { mpfr_set_d(r23181, im, MPFR_RNDN); ; mpfr_set_si(r23183, mpfr_cmp(r23181, r23182) <= 0, MPFR_RNDN); ; ; mpfr_pow(r23186, r23181, r23185, MPFR_RNDN); mpfr_mul(r23187, r23184, r23186, MPFR_RNDN); ; mpfr_mul(r23189, r23188, r23181, MPFR_RNDN); ; ; mpfr_pow(r23192, r23181, r23191, MPFR_RNDN); mpfr_mul(r23193, r23190, r23192, MPFR_RNDN); mpfr_add(r23194, r23189, r23193, MPFR_RNDN); mpfr_add(r23195, r23187, r23194, MPFR_RNDN); mpfr_neg(r23196, r23195, MPFR_RNDN); ; mpfr_set_d(r23198, re, MPFR_RNDN); mpfr_cos(r23199, r23198, MPFR_RNDN); mpfr_mul(r23200, r23197, r23199, MPFR_RNDN); mpfr_mul(r23201, r23196, r23200, MPFR_RNDN); mpfr_neg(r23202, r23181, MPFR_RNDN); mpfr_exp(r23203, r23202, MPFR_RNDN); mpfr_exp(r23204, r23181, MPFR_RNDN); mpfr_sub(r23205, r23203, r23204, MPFR_RNDN); mpfr_mul(r23206, r23205, r23200, MPFR_RNDN); if (mpfr_get_si(r23183, MPFR_RNDN)) { mpfr_set(r23207, r23201, MPFR_RNDN); } else { mpfr_set(r23207, r23206, MPFR_RNDN); }; return mpfr_get_d(r23207, MPFR_RNDN); } static mpfr_t r23208, r23209, r23210, r23211, r23212, r23213, r23214, r23215, r23216, r23217, r23218, r23219, r23220, r23221, r23222, r23223, r23224, r23225, r23226, r23227, r23228, r23229, r23230, r23231, r23232, r23233, r23234; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r23208); mpfr_init_set_str(r23209, "1.9183430624335127", 10, MPFR_RNDN); mpfr_init(r23210); mpfr_init_set_str(r23211, "1/60", 10, MPFR_RNDN); mpfr_init_set_str(r23212, "5", 10, MPFR_RNDN); mpfr_init(r23213); mpfr_init(r23214); mpfr_init_set_str(r23215, "2", 10, MPFR_RNDN); mpfr_init(r23216); mpfr_init_set_str(r23217, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r23218, "3", 10, MPFR_RNDN); mpfr_init(r23219); mpfr_init(r23220); mpfr_init(r23221); mpfr_init(r23222); mpfr_init(r23223); mpfr_init_set_str(r23224, "0.5", 10, MPFR_RNDN); mpfr_init(r23225); mpfr_init(r23226); mpfr_init(r23227); mpfr_init(r23228); mpfr_init(r23229); mpfr_init(r23230); mpfr_init(r23231); mpfr_init(r23232); mpfr_init(r23233); mpfr_init(r23234); } double f_dm(double re, double im) { mpfr_set_d(r23208, im, MPFR_RNDN); ; mpfr_set_si(r23210, mpfr_cmp(r23208, r23209) <= 0, MPFR_RNDN); ; ; mpfr_pow(r23213, r23208, r23212, MPFR_RNDN); mpfr_mul(r23214, r23211, r23213, MPFR_RNDN); ; mpfr_mul(r23216, r23215, r23208, MPFR_RNDN); ; ; mpfr_pow(r23219, r23208, r23218, MPFR_RNDN); mpfr_mul(r23220, r23217, r23219, MPFR_RNDN); mpfr_add(r23221, r23216, r23220, MPFR_RNDN); mpfr_add(r23222, r23214, r23221, MPFR_RNDN); mpfr_neg(r23223, r23222, MPFR_RNDN); ; mpfr_set_d(r23225, re, MPFR_RNDN); mpfr_cos(r23226, r23225, MPFR_RNDN); mpfr_mul(r23227, r23224, r23226, MPFR_RNDN); mpfr_mul(r23228, r23223, r23227, MPFR_RNDN); mpfr_neg(r23229, r23208, MPFR_RNDN); mpfr_exp(r23230, r23229, MPFR_RNDN); mpfr_exp(r23231, r23208, MPFR_RNDN); mpfr_sub(r23232, r23230, r23231, MPFR_RNDN); mpfr_mul(r23233, r23232, r23227, MPFR_RNDN); if (mpfr_get_si(r23210, MPFR_RNDN)) { mpfr_set(r23234, r23228, MPFR_RNDN); } else { mpfr_set(r23234, r23233, MPFR_RNDN); }; return mpfr_get_d(r23234, MPFR_RNDN); }