#include #include #include #include #include char *name = "math.log/2 on complex, real part"; double f_if(float re, float im, float base) { float r15618 = re; float r15619 = r15618 * r15618; float r15620 = im; float r15621 = r15620 * r15620; float r15622 = r15619 + r15621; float r15623 = sqrt(r15622); float r15624 = log(r15623); float r15625 = base; float r15626 = log(r15625); float r15627 = r15624 * r15626; float r15628 = atan2(r15620, r15618); float r15629 = 0.0f; float r15630 = r15628 * r15629; float r15631 = r15627 + r15630; float r15632 = r15626 * r15626; float r15633 = r15629 * r15629; float r15634 = r15632 + r15633; float r15635 = r15631 / r15634; return r15635; } double f_id(double re, double im, double base) { double r15636 = re; double r15637 = r15636 * r15636; double r15638 = im; double r15639 = r15638 * r15638; double r15640 = r15637 + r15639; double r15641 = sqrt(r15640); double r15642 = log(r15641); double r15643 = base; double r15644 = log(r15643); double r15645 = r15642 * r15644; double r15646 = atan2(r15638, r15636); double r15647 = 0.0; double r15648 = r15646 * r15647; double r15649 = r15645 + r15648; double r15650 = r15644 * r15644; double r15651 = r15647 * r15647; double r15652 = r15650 + r15651; double r15653 = r15649 / r15652; return r15653; } double f_of(float re, float im, float base) { float r15654 = re; float r15655 = -17576479424512.0f; bool r15656 = r15654 <= r15655; float r15657 = -r15654; float r15658 = log(r15657); float r15659 = base; float r15660 = log(r15659); float r15661 = r15658 / r15660; float r15662 = 32708774.0f; bool r15663 = r15654 <= r15662; float r15664 = r15654 * r15654; float r15665 = im; float r15666 = r15665 * r15665; float r15667 = r15664 + r15666; float r15668 = sqrt(r15667); float r15669 = log(r15668); float r15670 = r15660 * r15669; float r15671 = 0.0f; float r15672 = r15670 + r15671; float r15673 = r15660 * (r15660 * r15660); float r15674 = r15673 * r15673; float r15675 = cbrt(r15674); float r15676 = r15672 / r15675; float r15677 = log(r15654); float r15678 = r15677 / r15660; float r15679 = r15663 ? r15676 : r15678; float r15680 = r15656 ? r15661 : r15679; return r15680; } double f_od(double re, double im, double base) { double r15681 = re; double r15682 = -17576479424512.0; bool r15683 = r15681 <= r15682; double r15684 = -r15681; double r15685 = log(r15684); double r15686 = base; double r15687 = log(r15686); double r15688 = r15685 / r15687; double r15689 = 32708774.0; bool r15690 = r15681 <= r15689; double r15691 = r15681 * r15681; double r15692 = im; double r15693 = r15692 * r15692; double r15694 = r15691 + r15693; double r15695 = sqrt(r15694); double r15696 = log(r15695); double r15697 = r15687 * r15696; double r15698 = 0.0; double r15699 = r15697 + r15698; double r15700 = r15687 * (r15687 * r15687); double r15701 = r15700 * r15700; double r15702 = cbrt(r15701); double r15703 = r15699 / r15702; double r15704 = log(r15681); double r15705 = r15704 / r15687; double r15706 = r15690 ? r15703 : r15705; double r15707 = r15683 ? r15688 : r15706; return r15707; } 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 r15708, r15709, r15710, r15711, r15712, r15713, r15714, r15715, r15716, r15717, r15718, r15719, r15720, r15721, r15722, r15723, r15724, r15725; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r15708); mpfr_init(r15709); mpfr_init(r15710); mpfr_init(r15711); mpfr_init(r15712); mpfr_init(r15713); mpfr_init(r15714); mpfr_init(r15715); mpfr_init(r15716); mpfr_init(r15717); mpfr_init(r15718); mpfr_init_set_str(r15719, "0", 10, MPFR_RNDN); mpfr_init(r15720); mpfr_init(r15721); mpfr_init(r15722); mpfr_init(r15723); mpfr_init(r15724); mpfr_init(r15725); } double f_im(double re, double im, double base) { mpfr_set_d(r15708, re, MPFR_RNDN); mpfr_mul(r15709, r15708, r15708, MPFR_RNDN); mpfr_set_d(r15710, im, MPFR_RNDN); mpfr_mul(r15711, r15710, r15710, MPFR_RNDN); mpfr_add(r15712, r15709, r15711, MPFR_RNDN); mpfr_sqrt(r15713, r15712, MPFR_RNDN); mpfr_log(r15714, r15713, MPFR_RNDN); mpfr_set_d(r15715, base, MPFR_RNDN); mpfr_log(r15716, r15715, MPFR_RNDN); mpfr_mul(r15717, r15714, r15716, MPFR_RNDN); mpfr_atan2(r15718, r15710, r15708, MPFR_RNDN); ; mpfr_mul(r15720, r15718, r15719, MPFR_RNDN); mpfr_add(r15721, r15717, r15720, MPFR_RNDN); mpfr_mul(r15722, r15716, r15716, MPFR_RNDN); mpfr_mul(r15723, r15719, r15719, MPFR_RNDN); mpfr_add(r15724, r15722, r15723, MPFR_RNDN); mpfr_div(r15725, r15721, r15724, MPFR_RNDN); return mpfr_get_d(r15725, MPFR_RNDN); } static mpfr_t r15726, r15727, r15728, r15729, r15730, r15731, r15732, r15733, r15734, r15735, r15736, r15737, r15738, r15739, r15740, r15741, r15742, r15743, r15744, r15745, r15746, r15747, r15748, r15749, r15750, r15751, r15752; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r15726); mpfr_init_set_str(r15727, "-1.7576479f+13", 10, MPFR_RNDN); mpfr_init(r15728); mpfr_init(r15729); mpfr_init(r15730); mpfr_init(r15731); mpfr_init(r15732); mpfr_init(r15733); mpfr_init_set_str(r15734, "32708774.0f0", 10, MPFR_RNDN); mpfr_init(r15735); mpfr_init(r15736); mpfr_init(r15737); mpfr_init(r15738); mpfr_init(r15739); mpfr_init(r15740); mpfr_init(r15741); mpfr_init(r15742); mpfr_init_set_str(r15743, "0", 10, MPFR_RNDN); mpfr_init(r15744); mpfr_init(r15745); mpfr_init(r15746); mpfr_init(r15747); mpfr_init(r15748); mpfr_init(r15749); mpfr_init(r15750); mpfr_init(r15751); mpfr_init(r15752); } double f_fm(double re, double im, double base) { mpfr_set_d(r15726, re, MPFR_RNDN); ; mpfr_set_si(r15728, mpfr_cmp(r15726, r15727) <= 0, MPFR_RNDN); mpfr_neg(r15729, r15726, MPFR_RNDN); mpfr_log(r15730, r15729, MPFR_RNDN); mpfr_set_d(r15731, base, MPFR_RNDN); mpfr_log(r15732, r15731, MPFR_RNDN); mpfr_div(r15733, r15730, r15732, MPFR_RNDN); ; mpfr_set_si(r15735, mpfr_cmp(r15726, r15734) <= 0, MPFR_RNDN); mpfr_sqr(r15736, r15726, MPFR_RNDN); mpfr_set_d(r15737, im, MPFR_RNDN); mpfr_mul(r15738, r15737, r15737, MPFR_RNDN); mpfr_add(r15739, r15736, r15738, MPFR_RNDN); mpfr_sqrt(r15740, r15739, MPFR_RNDN); mpfr_log(r15741, r15740, MPFR_RNDN); mpfr_mul(r15742, r15732, r15741, MPFR_RNDN); ; mpfr_add(r15744, r15742, r15743, MPFR_RNDN); mpfr_mul(r15745, r15732, r15732, MPFR_RNDN); mpfr_mul(r15745, r15745, r15732, MPFR_RNDN); mpfr_mul(r15746, r15745, r15745, MPFR_RNDN); mpfr_cbrt(r15747, r15746, MPFR_RNDN); mpfr_div(r15748, r15744, r15747, MPFR_RNDN); mpfr_log(r15749, r15726, MPFR_RNDN); mpfr_div(r15750, r15749, r15732, MPFR_RNDN); if (mpfr_get_si(r15735, MPFR_RNDN)) { mpfr_set(r15751, r15748, MPFR_RNDN); } else { mpfr_set(r15751, r15750, MPFR_RNDN); }; if (mpfr_get_si(r15728, MPFR_RNDN)) { mpfr_set(r15752, r15733, MPFR_RNDN); } else { mpfr_set(r15752, r15751, MPFR_RNDN); }; return mpfr_get_d(r15752, MPFR_RNDN); } static mpfr_t r15753, r15754, r15755, r15756, r15757, r15758, r15759, r15760, r15761, r15762, r15763, r15764, r15765, r15766, r15767, r15768, r15769, r15770, r15771, r15772, r15773, r15774, r15775, r15776, r15777, r15778, r15779; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15753); mpfr_init_set_str(r15754, "-1.7576479f+13", 10, MPFR_RNDN); mpfr_init(r15755); mpfr_init(r15756); mpfr_init(r15757); mpfr_init(r15758); mpfr_init(r15759); mpfr_init(r15760); mpfr_init_set_str(r15761, "32708774.0f0", 10, MPFR_RNDN); mpfr_init(r15762); mpfr_init(r15763); mpfr_init(r15764); mpfr_init(r15765); mpfr_init(r15766); mpfr_init(r15767); mpfr_init(r15768); mpfr_init(r15769); mpfr_init_set_str(r15770, "0", 10, MPFR_RNDN); mpfr_init(r15771); mpfr_init(r15772); mpfr_init(r15773); mpfr_init(r15774); mpfr_init(r15775); mpfr_init(r15776); mpfr_init(r15777); mpfr_init(r15778); mpfr_init(r15779); } double f_dm(double re, double im, double base) { mpfr_set_d(r15753, re, MPFR_RNDN); ; mpfr_set_si(r15755, mpfr_cmp(r15753, r15754) <= 0, MPFR_RNDN); mpfr_neg(r15756, r15753, MPFR_RNDN); mpfr_log(r15757, r15756, MPFR_RNDN); mpfr_set_d(r15758, base, MPFR_RNDN); mpfr_log(r15759, r15758, MPFR_RNDN); mpfr_div(r15760, r15757, r15759, MPFR_RNDN); ; mpfr_set_si(r15762, mpfr_cmp(r15753, r15761) <= 0, MPFR_RNDN); mpfr_sqr(r15763, r15753, MPFR_RNDN); mpfr_set_d(r15764, im, MPFR_RNDN); mpfr_mul(r15765, r15764, r15764, MPFR_RNDN); mpfr_add(r15766, r15763, r15765, MPFR_RNDN); mpfr_sqrt(r15767, r15766, MPFR_RNDN); mpfr_log(r15768, r15767, MPFR_RNDN); mpfr_mul(r15769, r15759, r15768, MPFR_RNDN); ; mpfr_add(r15771, r15769, r15770, MPFR_RNDN); mpfr_mul(r15772, r15759, r15759, MPFR_RNDN); mpfr_mul(r15772, r15772, r15759, MPFR_RNDN); mpfr_mul(r15773, r15772, r15772, MPFR_RNDN); mpfr_cbrt(r15774, r15773, MPFR_RNDN); mpfr_div(r15775, r15771, r15774, MPFR_RNDN); mpfr_log(r15776, r15753, MPFR_RNDN); mpfr_div(r15777, r15776, r15759, MPFR_RNDN); if (mpfr_get_si(r15762, MPFR_RNDN)) { mpfr_set(r15778, r15775, MPFR_RNDN); } else { mpfr_set(r15778, r15777, MPFR_RNDN); }; if (mpfr_get_si(r15755, MPFR_RNDN)) { mpfr_set(r15779, r15760, MPFR_RNDN); } else { mpfr_set(r15779, r15778, MPFR_RNDN); }; return mpfr_get_d(r15779, MPFR_RNDN); }