#include #include #include #include #include char *name = "math.log10 on complex, real part"; double f_if(float re, float im) { float r21112 = re; float r21113 = r21112 * r21112; float r21114 = im; float r21115 = r21114 * r21114; float r21116 = r21113 + r21115; float r21117 = sqrt(r21116); float r21118 = log(r21117); float r21119 = 10; float r21120 = log(r21119); float r21121 = r21118 / r21120; return r21121; } double f_id(double re, double im) { double r21122 = re; double r21123 = r21122 * r21122; double r21124 = im; double r21125 = r21124 * r21124; double r21126 = r21123 + r21125; double r21127 = sqrt(r21126); double r21128 = log(r21127); double r21129 = 10; double r21130 = log(r21129); double r21131 = r21128 / r21130; return r21131; } double f_of(float re, float im) { float r21132 = re; float r21133 = -2.4746401519635354e+107; bool r21134 = r21132 <= r21133; float r21135 = 1/2; float r21136 = 10; float r21137 = log(r21136); float r21138 = sqrt(r21137); float r21139 = r21135 / r21138; float r21140 = sqrt(r21139); float r21141 = -2; float r21142 = sqrt(r21135); float r21143 = -1; float r21144 = r21143 / r21132; float r21145 = log(r21144); float r21146 = r21142 * r21145; float r21147 = 1; float r21148 = 3; float r21149 = pow(r21137, r21148); float r21150 = r21147 / r21149; float r21151 = 1/4; float r21152 = pow(r21150, r21151); float r21153 = r21146 * r21152; float r21154 = r21141 * r21153; float r21155 = r21140 * r21154; float r21156 = 2.5859190809973263e+94; bool r21157 = r21132 <= r21156; float r21158 = r21132 * r21132; float r21159 = im; float r21160 = r21159 * r21159; float r21161 = r21158 + r21160; float r21162 = log(r21161); float r21163 = r21162 / r21138; float r21164 = r21140 * r21163; float r21165 = r21140 * r21164; float r21166 = r21147 / r21137; float r21167 = sqrt(r21166); float r21168 = log(r21132); float r21169 = -r21141; float r21170 = r21168 * r21169; float r21171 = r21167 * r21170; float r21172 = r21139 * r21171; float r21173 = r21157 ? r21165 : r21172; float r21174 = r21134 ? r21155 : r21173; return r21174; } double f_od(double re, double im) { double r21175 = re; double r21176 = -2.4746401519635354e+107; bool r21177 = r21175 <= r21176; double r21178 = 1/2; double r21179 = 10; double r21180 = log(r21179); double r21181 = sqrt(r21180); double r21182 = r21178 / r21181; double r21183 = sqrt(r21182); double r21184 = -2; double r21185 = sqrt(r21178); double r21186 = -1; double r21187 = r21186 / r21175; double r21188 = log(r21187); double r21189 = r21185 * r21188; double r21190 = 1; double r21191 = 3; double r21192 = pow(r21180, r21191); double r21193 = r21190 / r21192; double r21194 = 1/4; double r21195 = pow(r21193, r21194); double r21196 = r21189 * r21195; double r21197 = r21184 * r21196; double r21198 = r21183 * r21197; double r21199 = 2.5859190809973263e+94; bool r21200 = r21175 <= r21199; double r21201 = r21175 * r21175; double r21202 = im; double r21203 = r21202 * r21202; double r21204 = r21201 + r21203; double r21205 = log(r21204); double r21206 = r21205 / r21181; double r21207 = r21183 * r21206; double r21208 = r21183 * r21207; double r21209 = r21190 / r21180; double r21210 = sqrt(r21209); double r21211 = log(r21175); double r21212 = -r21184; double r21213 = r21211 * r21212; double r21214 = r21210 * r21213; double r21215 = r21182 * r21214; double r21216 = r21200 ? r21208 : r21215; double r21217 = r21177 ? r21198 : r21216; return r21217; } 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 r21218, r21219, r21220, r21221, r21222, r21223, r21224, r21225, r21226, r21227; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r21218); mpfr_init(r21219); mpfr_init(r21220); mpfr_init(r21221); mpfr_init(r21222); mpfr_init(r21223); mpfr_init(r21224); mpfr_init_set_str(r21225, "10", 10, MPFR_RNDN); mpfr_init(r21226); mpfr_init(r21227); } double f_im(double re, double im) { mpfr_set_d(r21218, re, MPFR_RNDN); mpfr_mul(r21219, r21218, r21218, MPFR_RNDN); mpfr_set_d(r21220, im, MPFR_RNDN); mpfr_mul(r21221, r21220, r21220, MPFR_RNDN); mpfr_add(r21222, r21219, r21221, MPFR_RNDN); mpfr_sqrt(r21223, r21222, MPFR_RNDN); mpfr_log(r21224, r21223, MPFR_RNDN); ; mpfr_log(r21226, r21225, MPFR_RNDN); mpfr_div(r21227, r21224, r21226, MPFR_RNDN); return mpfr_get_d(r21227, MPFR_RNDN); } static mpfr_t r21228, r21229, r21230, r21231, r21232, r21233, r21234, r21235, r21236, r21237, r21238, r21239, r21240, r21241, r21242, r21243, r21244, r21245, r21246, r21247, r21248, r21249, r21250, r21251, r21252, r21253, r21254, r21255, r21256, r21257, r21258, r21259, r21260, r21261, r21262, r21263, r21264, r21265, r21266, r21267, r21268, r21269, r21270; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r21228); mpfr_init_set_str(r21229, "-2.4746401519635354e+107", 10, MPFR_RNDN); mpfr_init(r21230); mpfr_init_set_str(r21231, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r21232, "10", 10, MPFR_RNDN); mpfr_init(r21233); mpfr_init(r21234); mpfr_init(r21235); mpfr_init(r21236); mpfr_init_set_str(r21237, "-2", 10, MPFR_RNDN); mpfr_init(r21238); mpfr_init_set_str(r21239, "-1", 10, MPFR_RNDN); mpfr_init(r21240); mpfr_init(r21241); mpfr_init(r21242); mpfr_init_set_str(r21243, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21244, "3", 10, MPFR_RNDN); mpfr_init(r21245); mpfr_init(r21246); mpfr_init_set_str(r21247, "1/4", 10, MPFR_RNDN); mpfr_init(r21248); mpfr_init(r21249); mpfr_init(r21250); mpfr_init(r21251); mpfr_init_set_str(r21252, "2.5859190809973263e+94", 10, MPFR_RNDN); mpfr_init(r21253); mpfr_init(r21254); mpfr_init(r21255); mpfr_init(r21256); mpfr_init(r21257); mpfr_init(r21258); mpfr_init(r21259); mpfr_init(r21260); mpfr_init(r21261); mpfr_init(r21262); mpfr_init(r21263); mpfr_init(r21264); mpfr_init(r21265); mpfr_init(r21266); mpfr_init(r21267); mpfr_init(r21268); mpfr_init(r21269); mpfr_init(r21270); } double f_fm(double re, double im) { mpfr_set_d(r21228, re, MPFR_RNDN); ; mpfr_set_si(r21230, mpfr_cmp(r21228, r21229) <= 0, MPFR_RNDN); ; ; mpfr_log(r21233, r21232, MPFR_RNDN); mpfr_sqrt(r21234, r21233, MPFR_RNDN); mpfr_div(r21235, r21231, r21234, MPFR_RNDN); mpfr_sqrt(r21236, r21235, MPFR_RNDN); ; mpfr_sqrt(r21238, r21231, MPFR_RNDN); ; mpfr_div(r21240, r21239, r21228, MPFR_RNDN); mpfr_log(r21241, r21240, MPFR_RNDN); mpfr_mul(r21242, r21238, r21241, MPFR_RNDN); ; ; mpfr_pow(r21245, r21233, r21244, MPFR_RNDN); mpfr_div(r21246, r21243, r21245, MPFR_RNDN); ; mpfr_pow(r21248, r21246, r21247, MPFR_RNDN); mpfr_mul(r21249, r21242, r21248, MPFR_RNDN); mpfr_mul(r21250, r21237, r21249, MPFR_RNDN); mpfr_mul(r21251, r21236, r21250, MPFR_RNDN); ; mpfr_set_si(r21253, mpfr_cmp(r21228, r21252) <= 0, MPFR_RNDN); mpfr_mul(r21254, r21228, r21228, MPFR_RNDN); mpfr_set_d(r21255, im, MPFR_RNDN); mpfr_mul(r21256, r21255, r21255, MPFR_RNDN); mpfr_add(r21257, r21254, r21256, MPFR_RNDN); mpfr_log(r21258, r21257, MPFR_RNDN); mpfr_div(r21259, r21258, r21234, MPFR_RNDN); mpfr_mul(r21260, r21236, r21259, MPFR_RNDN); mpfr_mul(r21261, r21236, r21260, MPFR_RNDN); mpfr_div(r21262, r21243, r21233, MPFR_RNDN); mpfr_sqrt(r21263, r21262, MPFR_RNDN); mpfr_log(r21264, r21228, MPFR_RNDN); mpfr_neg(r21265, r21237, MPFR_RNDN); mpfr_mul(r21266, r21264, r21265, MPFR_RNDN); mpfr_mul(r21267, r21263, r21266, MPFR_RNDN); mpfr_mul(r21268, r21235, r21267, MPFR_RNDN); if (mpfr_get_si(r21253, MPFR_RNDN)) { mpfr_set(r21269, r21261, MPFR_RNDN); } else { mpfr_set(r21269, r21268, MPFR_RNDN); }; if (mpfr_get_si(r21230, MPFR_RNDN)) { mpfr_set(r21270, r21251, MPFR_RNDN); } else { mpfr_set(r21270, r21269, MPFR_RNDN); }; return mpfr_get_d(r21270, MPFR_RNDN); } static mpfr_t r21271, r21272, r21273, r21274, r21275, r21276, r21277, r21278, r21279, r21280, r21281, r21282, r21283, r21284, r21285, r21286, r21287, r21288, r21289, r21290, r21291, r21292, r21293, r21294, r21295, r21296, r21297, r21298, r21299, r21300, r21301, r21302, r21303, r21304, r21305, r21306, r21307, r21308, r21309, r21310, r21311, r21312, r21313; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r21271); mpfr_init_set_str(r21272, "-2.4746401519635354e+107", 10, MPFR_RNDN); mpfr_init(r21273); mpfr_init_set_str(r21274, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r21275, "10", 10, MPFR_RNDN); mpfr_init(r21276); mpfr_init(r21277); mpfr_init(r21278); mpfr_init(r21279); mpfr_init_set_str(r21280, "-2", 10, MPFR_RNDN); mpfr_init(r21281); mpfr_init_set_str(r21282, "-1", 10, MPFR_RNDN); mpfr_init(r21283); mpfr_init(r21284); mpfr_init(r21285); mpfr_init_set_str(r21286, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21287, "3", 10, MPFR_RNDN); mpfr_init(r21288); mpfr_init(r21289); mpfr_init_set_str(r21290, "1/4", 10, MPFR_RNDN); mpfr_init(r21291); mpfr_init(r21292); mpfr_init(r21293); mpfr_init(r21294); mpfr_init_set_str(r21295, "2.5859190809973263e+94", 10, MPFR_RNDN); mpfr_init(r21296); mpfr_init(r21297); mpfr_init(r21298); mpfr_init(r21299); mpfr_init(r21300); mpfr_init(r21301); mpfr_init(r21302); mpfr_init(r21303); mpfr_init(r21304); mpfr_init(r21305); mpfr_init(r21306); mpfr_init(r21307); mpfr_init(r21308); mpfr_init(r21309); mpfr_init(r21310); mpfr_init(r21311); mpfr_init(r21312); mpfr_init(r21313); } double f_dm(double re, double im) { mpfr_set_d(r21271, re, MPFR_RNDN); ; mpfr_set_si(r21273, mpfr_cmp(r21271, r21272) <= 0, MPFR_RNDN); ; ; mpfr_log(r21276, r21275, MPFR_RNDN); mpfr_sqrt(r21277, r21276, MPFR_RNDN); mpfr_div(r21278, r21274, r21277, MPFR_RNDN); mpfr_sqrt(r21279, r21278, MPFR_RNDN); ; mpfr_sqrt(r21281, r21274, MPFR_RNDN); ; mpfr_div(r21283, r21282, r21271, MPFR_RNDN); mpfr_log(r21284, r21283, MPFR_RNDN); mpfr_mul(r21285, r21281, r21284, MPFR_RNDN); ; ; mpfr_pow(r21288, r21276, r21287, MPFR_RNDN); mpfr_div(r21289, r21286, r21288, MPFR_RNDN); ; mpfr_pow(r21291, r21289, r21290, MPFR_RNDN); mpfr_mul(r21292, r21285, r21291, MPFR_RNDN); mpfr_mul(r21293, r21280, r21292, MPFR_RNDN); mpfr_mul(r21294, r21279, r21293, MPFR_RNDN); ; mpfr_set_si(r21296, mpfr_cmp(r21271, r21295) <= 0, MPFR_RNDN); mpfr_mul(r21297, r21271, r21271, MPFR_RNDN); mpfr_set_d(r21298, im, MPFR_RNDN); mpfr_mul(r21299, r21298, r21298, MPFR_RNDN); mpfr_add(r21300, r21297, r21299, MPFR_RNDN); mpfr_log(r21301, r21300, MPFR_RNDN); mpfr_div(r21302, r21301, r21277, MPFR_RNDN); mpfr_mul(r21303, r21279, r21302, MPFR_RNDN); mpfr_mul(r21304, r21279, r21303, MPFR_RNDN); mpfr_div(r21305, r21286, r21276, MPFR_RNDN); mpfr_sqrt(r21306, r21305, MPFR_RNDN); mpfr_log(r21307, r21271, MPFR_RNDN); mpfr_neg(r21308, r21280, MPFR_RNDN); mpfr_mul(r21309, r21307, r21308, MPFR_RNDN); mpfr_mul(r21310, r21306, r21309, MPFR_RNDN); mpfr_mul(r21311, r21278, r21310, MPFR_RNDN); if (mpfr_get_si(r21296, MPFR_RNDN)) { mpfr_set(r21312, r21304, MPFR_RNDN); } else { mpfr_set(r21312, r21311, MPFR_RNDN); }; if (mpfr_get_si(r21273, MPFR_RNDN)) { mpfr_set(r21313, r21294, MPFR_RNDN); } else { mpfr_set(r21313, r21312, MPFR_RNDN); }; return mpfr_get_d(r21313, MPFR_RNDN); }