#include #include #include #include #include char *name = "math.log/1 on complex, real part"; double f_if(float re, float im) { float r18094 = re; float r18095 = r18094 * r18094; float r18096 = im; float r18097 = r18096 * r18096; float r18098 = r18095 + r18097; float r18099 = sqrt(r18098); float r18100 = log(r18099); return r18100; } double f_id(double re, double im) { double r18101 = re; double r18102 = r18101 * r18101; double r18103 = im; double r18104 = r18103 * r18103; double r18105 = r18102 + r18104; double r18106 = sqrt(r18105); double r18107 = log(r18106); return r18107; } double f_of(float re, float im) { float r18108 = re; float r18109 = -1.1110224014220058e+125f; bool r18110 = r18108 <= r18109; float r18111 = -r18108; float r18112 = log(r18111); float r18113 = -1.219585437824595e-174f; bool r18114 = r18108 <= r18113; float r18115 = r18108 * r18108; float r18116 = im; float r18117 = r18116 * r18116; float r18118 = r18115 + r18117; float r18119 = sqrt(r18118); float r18120 = log(r18119); float r18121 = 1.1913008854416164e-287f; bool r18122 = r18108 <= r18121; float r18123 = log(r18116); float r18124 = 3.2338231936094466e+51f; bool r18125 = r18108 <= r18124; float r18126 = log(r18108); float r18127 = r18125 ? r18120 : r18126; float r18128 = r18122 ? r18123 : r18127; float r18129 = r18114 ? r18120 : r18128; float r18130 = r18110 ? r18112 : r18129; return r18130; } double f_od(double re, double im) { double r18131 = re; double r18132 = -1.1110224014220058e+125; bool r18133 = r18131 <= r18132; double r18134 = -r18131; double r18135 = log(r18134); double r18136 = -1.219585437824595e-174; bool r18137 = r18131 <= r18136; double r18138 = r18131 * r18131; double r18139 = im; double r18140 = r18139 * r18139; double r18141 = r18138 + r18140; double r18142 = sqrt(r18141); double r18143 = log(r18142); double r18144 = 1.1913008854416164e-287; bool r18145 = r18131 <= r18144; double r18146 = log(r18139); double r18147 = 3.2338231936094466e+51; bool r18148 = r18131 <= r18147; double r18149 = log(r18131); double r18150 = r18148 ? r18143 : r18149; double r18151 = r18145 ? r18146 : r18150; double r18152 = r18137 ? r18143 : r18151; double r18153 = r18133 ? r18135 : r18152; return r18153; } 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 r18154, r18155, r18156, r18157, r18158, r18159, r18160; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r18154); mpfr_init(r18155); mpfr_init(r18156); mpfr_init(r18157); mpfr_init(r18158); mpfr_init(r18159); mpfr_init(r18160); } double f_im(double re, double im) { mpfr_set_d(r18154, re, MPFR_RNDN); mpfr_mul(r18155, r18154, r18154, MPFR_RNDN); mpfr_set_d(r18156, im, MPFR_RNDN); mpfr_mul(r18157, r18156, r18156, MPFR_RNDN); mpfr_add(r18158, r18155, r18157, MPFR_RNDN); mpfr_sqrt(r18159, r18158, MPFR_RNDN); mpfr_log(r18160, r18159, MPFR_RNDN); return mpfr_get_d(r18160, MPFR_RNDN); } static mpfr_t r18161, r18162, r18163, r18164, r18165, r18166, r18167, r18168, r18169, r18170, r18171, r18172, r18173, r18174, r18175, r18176, r18177, r18178, r18179, r18180, r18181, r18182, r18183; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r18161); mpfr_init_set_str(r18162, "-1.1110224014220058e+125", 10, MPFR_RNDN); mpfr_init(r18163); mpfr_init(r18164); mpfr_init(r18165); mpfr_init_set_str(r18166, "-1.219585437824595e-174", 10, MPFR_RNDN); mpfr_init(r18167); mpfr_init(r18168); mpfr_init(r18169); mpfr_init(r18170); mpfr_init(r18171); mpfr_init(r18172); mpfr_init(r18173); mpfr_init_set_str(r18174, "1.1913008854416164e-287", 10, MPFR_RNDN); mpfr_init(r18175); mpfr_init(r18176); mpfr_init_set_str(r18177, "3.2338231936094466e+51", 10, MPFR_RNDN); mpfr_init(r18178); mpfr_init(r18179); mpfr_init(r18180); mpfr_init(r18181); mpfr_init(r18182); mpfr_init(r18183); } double f_fm(double re, double im) { mpfr_set_d(r18161, re, MPFR_RNDN); ; mpfr_set_si(r18163, mpfr_cmp(r18161, r18162) <= 0, MPFR_RNDN); mpfr_neg(r18164, r18161, MPFR_RNDN); mpfr_log(r18165, r18164, MPFR_RNDN); ; mpfr_set_si(r18167, mpfr_cmp(r18161, r18166) <= 0, MPFR_RNDN); mpfr_mul(r18168, r18161, r18161, MPFR_RNDN); mpfr_set_d(r18169, im, MPFR_RNDN); mpfr_mul(r18170, r18169, r18169, MPFR_RNDN); mpfr_add(r18171, r18168, r18170, MPFR_RNDN); mpfr_sqrt(r18172, r18171, MPFR_RNDN); mpfr_log(r18173, r18172, MPFR_RNDN); ; mpfr_set_si(r18175, mpfr_cmp(r18161, r18174) <= 0, MPFR_RNDN); mpfr_log(r18176, r18169, MPFR_RNDN); ; mpfr_set_si(r18178, mpfr_cmp(r18161, r18177) <= 0, MPFR_RNDN); mpfr_log(r18179, r18161, MPFR_RNDN); if (mpfr_get_si(r18178, MPFR_RNDN)) { mpfr_set(r18180, r18173, MPFR_RNDN); } else { mpfr_set(r18180, r18179, MPFR_RNDN); }; if (mpfr_get_si(r18175, MPFR_RNDN)) { mpfr_set(r18181, r18176, MPFR_RNDN); } else { mpfr_set(r18181, r18180, MPFR_RNDN); }; if (mpfr_get_si(r18167, MPFR_RNDN)) { mpfr_set(r18182, r18173, MPFR_RNDN); } else { mpfr_set(r18182, r18181, MPFR_RNDN); }; if (mpfr_get_si(r18163, MPFR_RNDN)) { mpfr_set(r18183, r18165, MPFR_RNDN); } else { mpfr_set(r18183, r18182, MPFR_RNDN); }; return mpfr_get_d(r18183, MPFR_RNDN); } static mpfr_t r18184, r18185, r18186, r18187, r18188, r18189, r18190, r18191, r18192, r18193, r18194, r18195, r18196, r18197, r18198, r18199, r18200, r18201, r18202, r18203, r18204, r18205, r18206; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r18184); mpfr_init_set_str(r18185, "-1.1110224014220058e+125", 10, MPFR_RNDN); mpfr_init(r18186); mpfr_init(r18187); mpfr_init(r18188); mpfr_init_set_str(r18189, "-1.219585437824595e-174", 10, MPFR_RNDN); mpfr_init(r18190); mpfr_init(r18191); mpfr_init(r18192); mpfr_init(r18193); mpfr_init(r18194); mpfr_init(r18195); mpfr_init(r18196); mpfr_init_set_str(r18197, "1.1913008854416164e-287", 10, MPFR_RNDN); mpfr_init(r18198); mpfr_init(r18199); mpfr_init_set_str(r18200, "3.2338231936094466e+51", 10, MPFR_RNDN); mpfr_init(r18201); mpfr_init(r18202); mpfr_init(r18203); mpfr_init(r18204); mpfr_init(r18205); mpfr_init(r18206); } double f_dm(double re, double im) { mpfr_set_d(r18184, re, MPFR_RNDN); ; mpfr_set_si(r18186, mpfr_cmp(r18184, r18185) <= 0, MPFR_RNDN); mpfr_neg(r18187, r18184, MPFR_RNDN); mpfr_log(r18188, r18187, MPFR_RNDN); ; mpfr_set_si(r18190, mpfr_cmp(r18184, r18189) <= 0, MPFR_RNDN); mpfr_mul(r18191, r18184, r18184, MPFR_RNDN); mpfr_set_d(r18192, im, MPFR_RNDN); mpfr_mul(r18193, r18192, r18192, MPFR_RNDN); mpfr_add(r18194, r18191, r18193, MPFR_RNDN); mpfr_sqrt(r18195, r18194, MPFR_RNDN); mpfr_log(r18196, r18195, MPFR_RNDN); ; mpfr_set_si(r18198, mpfr_cmp(r18184, r18197) <= 0, MPFR_RNDN); mpfr_log(r18199, r18192, MPFR_RNDN); ; mpfr_set_si(r18201, mpfr_cmp(r18184, r18200) <= 0, MPFR_RNDN); mpfr_log(r18202, r18184, MPFR_RNDN); if (mpfr_get_si(r18201, MPFR_RNDN)) { mpfr_set(r18203, r18196, MPFR_RNDN); } else { mpfr_set(r18203, r18202, MPFR_RNDN); }; if (mpfr_get_si(r18198, MPFR_RNDN)) { mpfr_set(r18204, r18199, MPFR_RNDN); } else { mpfr_set(r18204, r18203, MPFR_RNDN); }; if (mpfr_get_si(r18190, MPFR_RNDN)) { mpfr_set(r18205, r18196, MPFR_RNDN); } else { mpfr_set(r18205, r18204, MPFR_RNDN); }; if (mpfr_get_si(r18186, MPFR_RNDN)) { mpfr_set(r18206, r18188, MPFR_RNDN); } else { mpfr_set(r18206, r18205, MPFR_RNDN); }; return mpfr_get_d(r18206, MPFR_RNDN); }