#include #include #include #include #include char *name = "math.log10 on complex, real part"; double f_if(float re, float im) { float r8701 = re; float r8702 = r8701 * r8701; float r8703 = im; float r8704 = r8703 * r8703; float r8705 = r8702 + r8704; float r8706 = sqrt(r8705); float r8707 = log(r8706); float r8708 = 10; float r8709 = log(r8708); float r8710 = r8707 / r8709; return r8710; } double f_id(double re, double im) { double r8711 = re; double r8712 = r8711 * r8711; double r8713 = im; double r8714 = r8713 * r8713; double r8715 = r8712 + r8714; double r8716 = sqrt(r8715); double r8717 = log(r8716); double r8718 = 10; double r8719 = log(r8718); double r8720 = r8717 / r8719; return r8720; } double f_of(float re, float im) { float r8721 = re; float r8722 = -r8721; float r8723 = -7.347556294497716e+55; bool r8724 = r8722 <= r8723; float r8725 = 1/2; float r8726 = 10; float r8727 = log(r8726); float r8728 = sqrt(r8727); float r8729 = r8725 / r8728; float r8730 = log(r8721); float r8731 = -2; float r8732 = -r8731; float r8733 = r8730 * r8732; float r8734 = 1; float r8735 = r8734 / r8727; float r8736 = sqrt(r8735); float r8737 = r8733 * r8736; float r8738 = r8729 * r8737; float r8739 = -2.390514674954853e-254; bool r8740 = r8722 <= r8739; float r8741 = r8721 * r8721; float r8742 = im; float r8743 = r8742 * r8742; float r8744 = r8741 + r8743; float r8745 = log(r8744); float r8746 = r8728 / r8745; float r8747 = r8729 / r8746; float r8748 = 4.888860994581046e-300; bool r8749 = r8722 <= r8748; float r8750 = 2; float r8751 = log(r8742); float r8752 = r8751 * r8736; float r8753 = r8750 * r8752; float r8754 = r8729 * r8753; float r8755 = 3.7386430575018087e+48; bool r8756 = r8722 <= r8755; float r8757 = -1; float r8758 = r8757 / r8721; float r8759 = log(r8758); float r8760 = r8759 * r8736; float r8761 = r8731 * r8760; float r8762 = r8729 * r8761; float r8763 = r8756 ? r8747 : r8762; float r8764 = r8749 ? r8754 : r8763; float r8765 = r8740 ? r8747 : r8764; float r8766 = r8724 ? r8738 : r8765; return r8766; } double f_od(double re, double im) { double r8767 = re; double r8768 = -r8767; double r8769 = -7.347556294497716e+55; bool r8770 = r8768 <= r8769; double r8771 = 1/2; double r8772 = 10; double r8773 = log(r8772); double r8774 = sqrt(r8773); double r8775 = r8771 / r8774; double r8776 = log(r8767); double r8777 = -2; double r8778 = -r8777; double r8779 = r8776 * r8778; double r8780 = 1; double r8781 = r8780 / r8773; double r8782 = sqrt(r8781); double r8783 = r8779 * r8782; double r8784 = r8775 * r8783; double r8785 = -2.390514674954853e-254; bool r8786 = r8768 <= r8785; double r8787 = r8767 * r8767; double r8788 = im; double r8789 = r8788 * r8788; double r8790 = r8787 + r8789; double r8791 = log(r8790); double r8792 = r8774 / r8791; double r8793 = r8775 / r8792; double r8794 = 4.888860994581046e-300; bool r8795 = r8768 <= r8794; double r8796 = 2; double r8797 = log(r8788); double r8798 = r8797 * r8782; double r8799 = r8796 * r8798; double r8800 = r8775 * r8799; double r8801 = 3.7386430575018087e+48; bool r8802 = r8768 <= r8801; double r8803 = -1; double r8804 = r8803 / r8767; double r8805 = log(r8804); double r8806 = r8805 * r8782; double r8807 = r8777 * r8806; double r8808 = r8775 * r8807; double r8809 = r8802 ? r8793 : r8808; double r8810 = r8795 ? r8800 : r8809; double r8811 = r8786 ? r8793 : r8810; double r8812 = r8770 ? r8784 : r8811; return r8812; } 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 r8813, r8814, r8815, r8816, r8817, r8818, r8819, r8820, r8821, r8822; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r8813); mpfr_init(r8814); mpfr_init(r8815); mpfr_init(r8816); mpfr_init(r8817); mpfr_init(r8818); mpfr_init(r8819); mpfr_init_set_str(r8820, "10", 10, MPFR_RNDN); mpfr_init(r8821); mpfr_init(r8822); } double f_im(double re, double im) { mpfr_set_d(r8813, re, MPFR_RNDN); mpfr_mul(r8814, r8813, r8813, MPFR_RNDN); mpfr_set_d(r8815, im, MPFR_RNDN); mpfr_mul(r8816, r8815, r8815, MPFR_RNDN); mpfr_add(r8817, r8814, r8816, MPFR_RNDN); mpfr_sqrt(r8818, r8817, MPFR_RNDN); mpfr_log(r8819, r8818, MPFR_RNDN); ; mpfr_log(r8821, r8820, MPFR_RNDN); mpfr_div(r8822, r8819, r8821, MPFR_RNDN); return mpfr_get_d(r8822, MPFR_RNDN); } static mpfr_t r8823, r8824, r8825, r8826, r8827, r8828, r8829, r8830, r8831, r8832, r8833, r8834, r8835, r8836, r8837, r8838, r8839, r8840, r8841, r8842, r8843, r8844, r8845, r8846, r8847, r8848, r8849, r8850, r8851, r8852, r8853, r8854, r8855, r8856, r8857, r8858, r8859, r8860, r8861, r8862, r8863, r8864, r8865, r8866, r8867, r8868; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r8823); mpfr_init(r8824); mpfr_init_set_str(r8825, "-7.347556294497716e+55", 10, MPFR_RNDN); mpfr_init(r8826); mpfr_init_set_str(r8827, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r8828, "10", 10, MPFR_RNDN); mpfr_init(r8829); mpfr_init(r8830); mpfr_init(r8831); mpfr_init(r8832); mpfr_init_set_str(r8833, "-2", 10, MPFR_RNDN); mpfr_init(r8834); mpfr_init(r8835); mpfr_init_set_str(r8836, "1", 10, MPFR_RNDN); mpfr_init(r8837); mpfr_init(r8838); mpfr_init(r8839); mpfr_init(r8840); mpfr_init_set_str(r8841, "-2.390514674954853e-254", 10, MPFR_RNDN); mpfr_init(r8842); mpfr_init(r8843); mpfr_init(r8844); mpfr_init(r8845); mpfr_init(r8846); mpfr_init(r8847); mpfr_init(r8848); mpfr_init(r8849); mpfr_init_set_str(r8850, "4.888860994581046e-300", 10, MPFR_RNDN); mpfr_init(r8851); mpfr_init_set_str(r8852, "2", 10, MPFR_RNDN); mpfr_init(r8853); mpfr_init(r8854); mpfr_init(r8855); mpfr_init(r8856); mpfr_init_set_str(r8857, "3.7386430575018087e+48", 10, MPFR_RNDN); mpfr_init(r8858); mpfr_init_set_str(r8859, "-1", 10, MPFR_RNDN); mpfr_init(r8860); mpfr_init(r8861); mpfr_init(r8862); mpfr_init(r8863); mpfr_init(r8864); mpfr_init(r8865); mpfr_init(r8866); mpfr_init(r8867); mpfr_init(r8868); } double f_fm(double re, double im) { mpfr_set_d(r8823, re, MPFR_RNDN); mpfr_neg(r8824, r8823, MPFR_RNDN); ; mpfr_set_si(r8826, mpfr_cmp(r8824, r8825) <= 0, MPFR_RNDN); ; ; mpfr_log(r8829, r8828, MPFR_RNDN); mpfr_sqrt(r8830, r8829, MPFR_RNDN); mpfr_div(r8831, r8827, r8830, MPFR_RNDN); mpfr_log(r8832, r8823, MPFR_RNDN); ; mpfr_neg(r8834, r8833, MPFR_RNDN); mpfr_mul(r8835, r8832, r8834, MPFR_RNDN); ; mpfr_div(r8837, r8836, r8829, MPFR_RNDN); mpfr_sqrt(r8838, r8837, MPFR_RNDN); mpfr_mul(r8839, r8835, r8838, MPFR_RNDN); mpfr_mul(r8840, r8831, r8839, MPFR_RNDN); ; mpfr_set_si(r8842, mpfr_cmp(r8824, r8841) <= 0, MPFR_RNDN); mpfr_mul(r8843, r8823, r8823, MPFR_RNDN); mpfr_set_d(r8844, im, MPFR_RNDN); mpfr_mul(r8845, r8844, r8844, MPFR_RNDN); mpfr_add(r8846, r8843, r8845, MPFR_RNDN); mpfr_log(r8847, r8846, MPFR_RNDN); mpfr_div(r8848, r8830, r8847, MPFR_RNDN); mpfr_div(r8849, r8831, r8848, MPFR_RNDN); ; mpfr_set_si(r8851, mpfr_cmp(r8824, r8850) <= 0, MPFR_RNDN); ; mpfr_log(r8853, r8844, MPFR_RNDN); mpfr_mul(r8854, r8853, r8838, MPFR_RNDN); mpfr_mul(r8855, r8852, r8854, MPFR_RNDN); mpfr_mul(r8856, r8831, r8855, MPFR_RNDN); ; mpfr_set_si(r8858, mpfr_cmp(r8824, r8857) <= 0, MPFR_RNDN); ; mpfr_div(r8860, r8859, r8823, MPFR_RNDN); mpfr_log(r8861, r8860, MPFR_RNDN); mpfr_mul(r8862, r8861, r8838, MPFR_RNDN); mpfr_mul(r8863, r8833, r8862, MPFR_RNDN); mpfr_mul(r8864, r8831, r8863, MPFR_RNDN); if (mpfr_get_si(r8858, MPFR_RNDN)) { mpfr_set(r8865, r8849, MPFR_RNDN); } else { mpfr_set(r8865, r8864, MPFR_RNDN); }; if (mpfr_get_si(r8851, MPFR_RNDN)) { mpfr_set(r8866, r8856, MPFR_RNDN); } else { mpfr_set(r8866, r8865, MPFR_RNDN); }; if (mpfr_get_si(r8842, MPFR_RNDN)) { mpfr_set(r8867, r8849, MPFR_RNDN); } else { mpfr_set(r8867, r8866, MPFR_RNDN); }; if (mpfr_get_si(r8826, MPFR_RNDN)) { mpfr_set(r8868, r8840, MPFR_RNDN); } else { mpfr_set(r8868, r8867, MPFR_RNDN); }; return mpfr_get_d(r8868, MPFR_RNDN); } static mpfr_t r8869, r8870, r8871, r8872, r8873, r8874, r8875, r8876, r8877, r8878, r8879, r8880, r8881, r8882, r8883, r8884, r8885, r8886, r8887, r8888, r8889, r8890, r8891, r8892, r8893, r8894, r8895, r8896, r8897, r8898, r8899, r8900, r8901, r8902, r8903, r8904, r8905, r8906, r8907, r8908, r8909, r8910, r8911, r8912, r8913, r8914; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r8869); mpfr_init(r8870); mpfr_init_set_str(r8871, "-7.347556294497716e+55", 10, MPFR_RNDN); mpfr_init(r8872); mpfr_init_set_str(r8873, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r8874, "10", 10, MPFR_RNDN); mpfr_init(r8875); mpfr_init(r8876); mpfr_init(r8877); mpfr_init(r8878); mpfr_init_set_str(r8879, "-2", 10, MPFR_RNDN); mpfr_init(r8880); mpfr_init(r8881); mpfr_init_set_str(r8882, "1", 10, MPFR_RNDN); mpfr_init(r8883); mpfr_init(r8884); mpfr_init(r8885); mpfr_init(r8886); mpfr_init_set_str(r8887, "-2.390514674954853e-254", 10, MPFR_RNDN); mpfr_init(r8888); mpfr_init(r8889); mpfr_init(r8890); mpfr_init(r8891); mpfr_init(r8892); mpfr_init(r8893); mpfr_init(r8894); mpfr_init(r8895); mpfr_init_set_str(r8896, "4.888860994581046e-300", 10, MPFR_RNDN); mpfr_init(r8897); mpfr_init_set_str(r8898, "2", 10, MPFR_RNDN); mpfr_init(r8899); mpfr_init(r8900); mpfr_init(r8901); mpfr_init(r8902); mpfr_init_set_str(r8903, "3.7386430575018087e+48", 10, MPFR_RNDN); mpfr_init(r8904); mpfr_init_set_str(r8905, "-1", 10, MPFR_RNDN); mpfr_init(r8906); mpfr_init(r8907); mpfr_init(r8908); mpfr_init(r8909); mpfr_init(r8910); mpfr_init(r8911); mpfr_init(r8912); mpfr_init(r8913); mpfr_init(r8914); } double f_dm(double re, double im) { mpfr_set_d(r8869, re, MPFR_RNDN); mpfr_neg(r8870, r8869, MPFR_RNDN); ; mpfr_set_si(r8872, mpfr_cmp(r8870, r8871) <= 0, MPFR_RNDN); ; ; mpfr_log(r8875, r8874, MPFR_RNDN); mpfr_sqrt(r8876, r8875, MPFR_RNDN); mpfr_div(r8877, r8873, r8876, MPFR_RNDN); mpfr_log(r8878, r8869, MPFR_RNDN); ; mpfr_neg(r8880, r8879, MPFR_RNDN); mpfr_mul(r8881, r8878, r8880, MPFR_RNDN); ; mpfr_div(r8883, r8882, r8875, MPFR_RNDN); mpfr_sqrt(r8884, r8883, MPFR_RNDN); mpfr_mul(r8885, r8881, r8884, MPFR_RNDN); mpfr_mul(r8886, r8877, r8885, MPFR_RNDN); ; mpfr_set_si(r8888, mpfr_cmp(r8870, r8887) <= 0, MPFR_RNDN); mpfr_mul(r8889, r8869, r8869, MPFR_RNDN); mpfr_set_d(r8890, im, MPFR_RNDN); mpfr_mul(r8891, r8890, r8890, MPFR_RNDN); mpfr_add(r8892, r8889, r8891, MPFR_RNDN); mpfr_log(r8893, r8892, MPFR_RNDN); mpfr_div(r8894, r8876, r8893, MPFR_RNDN); mpfr_div(r8895, r8877, r8894, MPFR_RNDN); ; mpfr_set_si(r8897, mpfr_cmp(r8870, r8896) <= 0, MPFR_RNDN); ; mpfr_log(r8899, r8890, MPFR_RNDN); mpfr_mul(r8900, r8899, r8884, MPFR_RNDN); mpfr_mul(r8901, r8898, r8900, MPFR_RNDN); mpfr_mul(r8902, r8877, r8901, MPFR_RNDN); ; mpfr_set_si(r8904, mpfr_cmp(r8870, r8903) <= 0, MPFR_RNDN); ; mpfr_div(r8906, r8905, r8869, MPFR_RNDN); mpfr_log(r8907, r8906, MPFR_RNDN); mpfr_mul(r8908, r8907, r8884, MPFR_RNDN); mpfr_mul(r8909, r8879, r8908, MPFR_RNDN); mpfr_mul(r8910, r8877, r8909, MPFR_RNDN); if (mpfr_get_si(r8904, MPFR_RNDN)) { mpfr_set(r8911, r8895, MPFR_RNDN); } else { mpfr_set(r8911, r8910, MPFR_RNDN); }; if (mpfr_get_si(r8897, MPFR_RNDN)) { mpfr_set(r8912, r8902, MPFR_RNDN); } else { mpfr_set(r8912, r8911, MPFR_RNDN); }; if (mpfr_get_si(r8888, MPFR_RNDN)) { mpfr_set(r8913, r8895, MPFR_RNDN); } else { mpfr_set(r8913, r8912, MPFR_RNDN); }; if (mpfr_get_si(r8872, MPFR_RNDN)) { mpfr_set(r8914, r8886, MPFR_RNDN); } else { mpfr_set(r8914, r8913, MPFR_RNDN); }; return mpfr_get_d(r8914, MPFR_RNDN); }