#include #include #include #include #include char *name = "math.log10 on complex, real part"; double f_if(float re, float im) { float r15406 = re; float r15407 = r15406 * r15406; float r15408 = im; float r15409 = r15408 * r15408; float r15410 = r15407 + r15409; float r15411 = sqrt(r15410); float r15412 = log(r15411); float r15413 = 10.0f; float r15414 = log(r15413); float r15415 = r15412 / r15414; return r15415; } double f_id(double re, double im) { double r15416 = re; double r15417 = r15416 * r15416; double r15418 = im; double r15419 = r15418 * r15418; double r15420 = r15417 + r15419; double r15421 = sqrt(r15420); double r15422 = log(r15421); double r15423 = 10.0; double r15424 = log(r15423); double r15425 = r15422 / r15424; return r15425; } double f_of(float re, float im) { float r15426 = re; float r15427 = -694992240640.0f; bool r15428 = r15426 <= r15427; float r15429 = 3.0f; float r15430 = 10.0f; float r15431 = log(r15430); float r15432 = r15429 / r15431; float r15433 = -r15426; float r15434 = cbrt(r15433); float r15435 = log(r15434); float r15436 = r15432 * r15435; float r15437 = 149.7889862060547f; bool r15438 = r15426 <= r15437; float r15439 = r15426 * r15426; float r15440 = im; float r15441 = r15440 * r15440; float r15442 = r15439 + r15441; float r15443 = sqrt(r15442); float r15444 = cbrt(r15443); float r15445 = r15444 * (r15444 * r15444); float r15446 = log(r15445); float r15447 = r15446 / r15431; float r15448 = cbrt(r15426); float r15449 = log(r15448); float r15450 = r15431 / r15429; float r15451 = r15449 / r15450; float r15452 = r15438 ? r15447 : r15451; float r15453 = r15428 ? r15436 : r15452; return r15453; } double f_od(double re, double im) { double r15454 = re; double r15455 = -694992240640.0; bool r15456 = r15454 <= r15455; double r15457 = 3.0; double r15458 = 10.0; double r15459 = log(r15458); double r15460 = r15457 / r15459; double r15461 = -r15454; double r15462 = cbrt(r15461); double r15463 = log(r15462); double r15464 = r15460 * r15463; double r15465 = 149.7889862060547; bool r15466 = r15454 <= r15465; double r15467 = r15454 * r15454; double r15468 = im; double r15469 = r15468 * r15468; double r15470 = r15467 + r15469; double r15471 = sqrt(r15470); double r15472 = cbrt(r15471); double r15473 = r15472 * (r15472 * r15472); double r15474 = log(r15473); double r15475 = r15474 / r15459; double r15476 = cbrt(r15454); double r15477 = log(r15476); double r15478 = r15459 / r15457; double r15479 = r15477 / r15478; double r15480 = r15466 ? r15475 : r15479; double r15481 = r15456 ? r15464 : r15480; return r15481; } 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 r15482, r15483, r15484, r15485, r15486, r15487, r15488, r15489, r15490, r15491; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r15482); mpfr_init(r15483); mpfr_init(r15484); mpfr_init(r15485); mpfr_init(r15486); mpfr_init(r15487); mpfr_init(r15488); mpfr_init_set_str(r15489, "10", 10, MPFR_RNDN); mpfr_init(r15490); mpfr_init(r15491); } double f_im(double re, double im) { mpfr_set_d(r15482, re, MPFR_RNDN); mpfr_mul(r15483, r15482, r15482, MPFR_RNDN); mpfr_set_d(r15484, im, MPFR_RNDN); mpfr_mul(r15485, r15484, r15484, MPFR_RNDN); mpfr_add(r15486, r15483, r15485, MPFR_RNDN); mpfr_sqrt(r15487, r15486, MPFR_RNDN); mpfr_log(r15488, r15487, MPFR_RNDN); ; mpfr_log(r15490, r15489, MPFR_RNDN); mpfr_div(r15491, r15488, r15490, MPFR_RNDN); return mpfr_get_d(r15491, MPFR_RNDN); } static mpfr_t r15492, r15493, r15494, r15495, r15496, r15497, r15498, r15499, r15500, r15501, r15502, r15503, r15504, r15505, r15506, r15507, r15508, r15509, r15510, r15511, r15512, r15513, r15514, r15515, r15516, r15517, r15518, r15519; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r15492); mpfr_init_set_str(r15493, "-6.9499224f+11", 10, MPFR_RNDN); mpfr_init(r15494); mpfr_init_set_str(r15495, "3", 10, MPFR_RNDN); mpfr_init_set_str(r15496, "10", 10, MPFR_RNDN); mpfr_init(r15497); mpfr_init(r15498); mpfr_init(r15499); mpfr_init(r15500); mpfr_init(r15501); mpfr_init(r15502); mpfr_init_set_str(r15503, "149.78899f0", 10, MPFR_RNDN); mpfr_init(r15504); mpfr_init(r15505); mpfr_init(r15506); mpfr_init(r15507); mpfr_init(r15508); mpfr_init(r15509); mpfr_init(r15510); mpfr_init(r15511); mpfr_init(r15512); mpfr_init(r15513); mpfr_init(r15514); mpfr_init(r15515); mpfr_init(r15516); mpfr_init(r15517); mpfr_init(r15518); mpfr_init(r15519); } double f_fm(double re, double im) { mpfr_set_d(r15492, re, MPFR_RNDN); ; mpfr_set_si(r15494, mpfr_cmp(r15492, r15493) <= 0, MPFR_RNDN); ; ; mpfr_log(r15497, r15496, MPFR_RNDN); mpfr_div(r15498, r15495, r15497, MPFR_RNDN); mpfr_neg(r15499, r15492, MPFR_RNDN); mpfr_cbrt(r15500, r15499, MPFR_RNDN); mpfr_log(r15501, r15500, MPFR_RNDN); mpfr_mul(r15502, r15498, r15501, MPFR_RNDN); ; mpfr_set_si(r15504, mpfr_cmp(r15492, r15503) <= 0, MPFR_RNDN); mpfr_sqr(r15505, r15492, MPFR_RNDN); mpfr_set_d(r15506, im, MPFR_RNDN); mpfr_mul(r15507, r15506, r15506, MPFR_RNDN); mpfr_add(r15508, r15505, r15507, MPFR_RNDN); mpfr_sqrt(r15509, r15508, MPFR_RNDN); mpfr_cbrt(r15510, r15509, MPFR_RNDN); mpfr_mul(r15511, r15510, r15510, MPFR_RNDN); mpfr_mul(r15511, r15511, r15510, MPFR_RNDN); mpfr_log(r15512, r15511, MPFR_RNDN); mpfr_div(r15513, r15512, r15497, MPFR_RNDN); mpfr_cbrt(r15514, r15492, MPFR_RNDN); mpfr_log(r15515, r15514, MPFR_RNDN); mpfr_div(r15516, r15497, r15495, MPFR_RNDN); mpfr_div(r15517, r15515, r15516, MPFR_RNDN); if (mpfr_get_si(r15504, MPFR_RNDN)) { mpfr_set(r15518, r15513, MPFR_RNDN); } else { mpfr_set(r15518, r15517, MPFR_RNDN); }; if (mpfr_get_si(r15494, MPFR_RNDN)) { mpfr_set(r15519, r15502, MPFR_RNDN); } else { mpfr_set(r15519, r15518, MPFR_RNDN); }; return mpfr_get_d(r15519, MPFR_RNDN); } static mpfr_t r15520, r15521, r15522, r15523, r15524, r15525, r15526, r15527, r15528, r15529, r15530, r15531, r15532, r15533, r15534, r15535, r15536, r15537, r15538, r15539, r15540, r15541, r15542, r15543, r15544, r15545, r15546, r15547; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15520); mpfr_init_set_str(r15521, "-6.9499224f+11", 10, MPFR_RNDN); mpfr_init(r15522); mpfr_init_set_str(r15523, "3", 10, MPFR_RNDN); mpfr_init_set_str(r15524, "10", 10, MPFR_RNDN); mpfr_init(r15525); mpfr_init(r15526); mpfr_init(r15527); mpfr_init(r15528); mpfr_init(r15529); mpfr_init(r15530); mpfr_init_set_str(r15531, "149.78899f0", 10, MPFR_RNDN); mpfr_init(r15532); mpfr_init(r15533); mpfr_init(r15534); mpfr_init(r15535); mpfr_init(r15536); mpfr_init(r15537); mpfr_init(r15538); mpfr_init(r15539); mpfr_init(r15540); mpfr_init(r15541); mpfr_init(r15542); mpfr_init(r15543); mpfr_init(r15544); mpfr_init(r15545); mpfr_init(r15546); mpfr_init(r15547); } double f_dm(double re, double im) { mpfr_set_d(r15520, re, MPFR_RNDN); ; mpfr_set_si(r15522, mpfr_cmp(r15520, r15521) <= 0, MPFR_RNDN); ; ; mpfr_log(r15525, r15524, MPFR_RNDN); mpfr_div(r15526, r15523, r15525, MPFR_RNDN); mpfr_neg(r15527, r15520, MPFR_RNDN); mpfr_cbrt(r15528, r15527, MPFR_RNDN); mpfr_log(r15529, r15528, MPFR_RNDN); mpfr_mul(r15530, r15526, r15529, MPFR_RNDN); ; mpfr_set_si(r15532, mpfr_cmp(r15520, r15531) <= 0, MPFR_RNDN); mpfr_sqr(r15533, r15520, MPFR_RNDN); mpfr_set_d(r15534, im, MPFR_RNDN); mpfr_mul(r15535, r15534, r15534, MPFR_RNDN); mpfr_add(r15536, r15533, r15535, MPFR_RNDN); mpfr_sqrt(r15537, r15536, MPFR_RNDN); mpfr_cbrt(r15538, r15537, MPFR_RNDN); mpfr_mul(r15539, r15538, r15538, MPFR_RNDN); mpfr_mul(r15539, r15539, r15538, MPFR_RNDN); mpfr_log(r15540, r15539, MPFR_RNDN); mpfr_div(r15541, r15540, r15525, MPFR_RNDN); mpfr_cbrt(r15542, r15520, MPFR_RNDN); mpfr_log(r15543, r15542, MPFR_RNDN); mpfr_div(r15544, r15525, r15523, MPFR_RNDN); mpfr_div(r15545, r15543, r15544, MPFR_RNDN); if (mpfr_get_si(r15532, MPFR_RNDN)) { mpfr_set(r15546, r15541, MPFR_RNDN); } else { mpfr_set(r15546, r15545, MPFR_RNDN); }; if (mpfr_get_si(r15522, MPFR_RNDN)) { mpfr_set(r15547, r15530, MPFR_RNDN); } else { mpfr_set(r15547, r15546, MPFR_RNDN); }; return mpfr_get_d(r15547, MPFR_RNDN); }