#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (20)"; double f_if(float f) { float r7734 = 1.0f; float r7735 = atan2(1.0, 0.0); float r7736 = 4.0f; float r7737 = r7735 / r7736; float r7738 = r7734 / r7737; float r7739 = f; float r7740 = r7737 * r7739; float r7741 = exp(r7740); float r7742 = -r7740; float r7743 = exp(r7742); float r7744 = r7741 + r7743; float r7745 = r7741 - r7743; float r7746 = r7744 / r7745; float r7747 = log(r7746); float r7748 = r7738 * r7747; float r7749 = -r7748; return r7749; } double f_id(double f) { double r7750 = 1.0; double r7751 = atan2(1.0, 0.0); double r7752 = 4.0; double r7753 = r7751 / r7752; double r7754 = r7750 / r7753; double r7755 = f; double r7756 = r7753 * r7755; double r7757 = exp(r7756); double r7758 = -r7756; double r7759 = exp(r7758); double r7760 = r7757 + r7759; double r7761 = r7757 - r7759; double r7762 = r7760 / r7761; double r7763 = log(r7762); double r7764 = r7754 * r7763; double r7765 = -r7764; return r7765; } double f_of(float f) { float r7766 = 4.0f; float r7767 = atan2(1.0, 0.0); float r7768 = r7766 / r7767; float r7769 = cbrt(r7768); float r7770 = r7769 * r7769; float r7771 = log(r7770); float r7772 = r7768 * r7771; float r7773 = log(r7769); float r7774 = f; float r7775 = log(r7774); float r7776 = r7773 - r7775; float r7777 = r7768 * r7776; float r7778 = r7772 + r7777; float r7779 = 0.0012152777777777778f; float r7780 = r7767 * r7779; float r7781 = r7767 * r7780; float r7782 = pow(r7774, r7766); float r7783 = r7782 * r7767; float r7784 = r7781 * r7783; float r7785 = 0.08333333333333333f; float r7786 = r7785 * r7767; float r7787 = r7774 * r7774; float r7788 = r7786 * r7787; float r7789 = r7784 - r7788; float r7790 = r7778 - r7789; float r7791 = -r7790; return r7791; } double f_od(double f) { double r7792 = 4.0; double r7793 = atan2(1.0, 0.0); double r7794 = r7792 / r7793; double r7795 = cbrt(r7794); double r7796 = r7795 * r7795; double r7797 = log(r7796); double r7798 = r7794 * r7797; double r7799 = log(r7795); double r7800 = f; double r7801 = log(r7800); double r7802 = r7799 - r7801; double r7803 = r7794 * r7802; double r7804 = r7798 + r7803; double r7805 = 0.0012152777777777778; double r7806 = r7793 * r7805; double r7807 = r7793 * r7806; double r7808 = pow(r7800, r7792); double r7809 = r7808 * r7793; double r7810 = r7807 * r7809; double r7811 = 0.08333333333333333; double r7812 = r7811 * r7793; double r7813 = r7800 * r7800; double r7814 = r7812 * r7813; double r7815 = r7810 - r7814; double r7816 = r7804 - r7815; double r7817 = -r7816; return r7817; } 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 r7818, r7819, r7820, r7821, r7822, r7823, r7824, r7825, r7826, r7827, r7828, r7829, r7830, r7831, r7832, r7833; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r7818, "1", 10, MPFR_RNDN); mpfr_init(r7819); mpfr_init_set_str(r7820, "4", 10, MPFR_RNDN); mpfr_init(r7821); mpfr_init(r7822); mpfr_init(r7823); mpfr_init(r7824); mpfr_init(r7825); mpfr_init(r7826); mpfr_init(r7827); mpfr_init(r7828); mpfr_init(r7829); mpfr_init(r7830); mpfr_init(r7831); mpfr_init(r7832); mpfr_init(r7833); } double f_im(double f) { ; mpfr_const_pi(r7819, MPFR_RNDN); ; mpfr_div(r7821, r7819, r7820, MPFR_RNDN); mpfr_div(r7822, r7818, r7821, MPFR_RNDN); mpfr_set_d(r7823, f, MPFR_RNDN); mpfr_mul(r7824, r7821, r7823, MPFR_RNDN); mpfr_exp(r7825, r7824, MPFR_RNDN); mpfr_neg(r7826, r7824, MPFR_RNDN); mpfr_exp(r7827, r7826, MPFR_RNDN); mpfr_add(r7828, r7825, r7827, MPFR_RNDN); mpfr_sub(r7829, r7825, r7827, MPFR_RNDN); mpfr_div(r7830, r7828, r7829, MPFR_RNDN); mpfr_log(r7831, r7830, MPFR_RNDN); mpfr_mul(r7832, r7822, r7831, MPFR_RNDN); mpfr_neg(r7833, r7832, MPFR_RNDN); return mpfr_get_d(r7833, MPFR_RNDN); } static mpfr_t r7834, r7835, r7836, r7837, r7838, r7839, r7840, r7841, r7842, r7843, r7844, r7845, r7846, r7847, r7848, r7849, r7850, r7851, r7852, r7853, r7854, r7855, r7856, r7857, r7858, r7859; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r7834, "4", 10, MPFR_RNDN); mpfr_init(r7835); mpfr_init(r7836); mpfr_init(r7837); mpfr_init(r7838); mpfr_init(r7839); mpfr_init(r7840); mpfr_init(r7841); mpfr_init(r7842); mpfr_init(r7843); mpfr_init(r7844); mpfr_init(r7845); mpfr_init(r7846); mpfr_init_set_str(r7847, "7/5760", 10, MPFR_RNDN); mpfr_init(r7848); mpfr_init(r7849); mpfr_init(r7850); mpfr_init(r7851); mpfr_init(r7852); mpfr_init_set_str(r7853, "1/12", 10, MPFR_RNDN); mpfr_init(r7854); mpfr_init(r7855); mpfr_init(r7856); mpfr_init(r7857); mpfr_init(r7858); mpfr_init(r7859); } double f_fm(double f) { ; mpfr_const_pi(r7835, MPFR_RNDN); mpfr_div(r7836, r7834, r7835, MPFR_RNDN); mpfr_cbrt(r7837, r7836, MPFR_RNDN); mpfr_mul(r7838, r7837, r7837, MPFR_RNDN); mpfr_log(r7839, r7838, MPFR_RNDN); mpfr_mul(r7840, r7836, r7839, MPFR_RNDN); mpfr_log(r7841, r7837, MPFR_RNDN); mpfr_set_d(r7842, f, MPFR_RNDN); mpfr_log(r7843, r7842, MPFR_RNDN); mpfr_sub(r7844, r7841, r7843, MPFR_RNDN); mpfr_mul(r7845, r7836, r7844, MPFR_RNDN); mpfr_add(r7846, r7840, r7845, MPFR_RNDN); ; mpfr_mul(r7848, r7835, r7847, MPFR_RNDN); mpfr_mul(r7849, r7835, r7848, MPFR_RNDN); mpfr_pow(r7850, r7842, r7834, MPFR_RNDN); mpfr_mul(r7851, r7850, r7835, MPFR_RNDN); mpfr_mul(r7852, r7849, r7851, MPFR_RNDN); ; mpfr_mul(r7854, r7853, r7835, MPFR_RNDN); mpfr_mul(r7855, r7842, r7842, MPFR_RNDN); mpfr_mul(r7856, r7854, r7855, MPFR_RNDN); mpfr_sub(r7857, r7852, r7856, MPFR_RNDN); mpfr_sub(r7858, r7846, r7857, MPFR_RNDN); mpfr_neg(r7859, r7858, MPFR_RNDN); return mpfr_get_d(r7859, MPFR_RNDN); } static mpfr_t r7860, r7861, r7862, r7863, r7864, r7865, r7866, r7867, r7868, r7869, r7870, r7871, r7872, r7873, r7874, r7875, r7876, r7877, r7878, r7879, r7880, r7881, r7882, r7883, r7884, r7885; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r7860, "4", 10, MPFR_RNDN); mpfr_init(r7861); mpfr_init(r7862); mpfr_init(r7863); mpfr_init(r7864); mpfr_init(r7865); mpfr_init(r7866); mpfr_init(r7867); mpfr_init(r7868); mpfr_init(r7869); mpfr_init(r7870); mpfr_init(r7871); mpfr_init(r7872); mpfr_init_set_str(r7873, "7/5760", 10, MPFR_RNDN); mpfr_init(r7874); mpfr_init(r7875); mpfr_init(r7876); mpfr_init(r7877); mpfr_init(r7878); mpfr_init_set_str(r7879, "1/12", 10, MPFR_RNDN); mpfr_init(r7880); mpfr_init(r7881); mpfr_init(r7882); mpfr_init(r7883); mpfr_init(r7884); mpfr_init(r7885); } double f_dm(double f) { ; mpfr_const_pi(r7861, MPFR_RNDN); mpfr_div(r7862, r7860, r7861, MPFR_RNDN); mpfr_cbrt(r7863, r7862, MPFR_RNDN); mpfr_mul(r7864, r7863, r7863, MPFR_RNDN); mpfr_log(r7865, r7864, MPFR_RNDN); mpfr_mul(r7866, r7862, r7865, MPFR_RNDN); mpfr_log(r7867, r7863, MPFR_RNDN); mpfr_set_d(r7868, f, MPFR_RNDN); mpfr_log(r7869, r7868, MPFR_RNDN); mpfr_sub(r7870, r7867, r7869, MPFR_RNDN); mpfr_mul(r7871, r7862, r7870, MPFR_RNDN); mpfr_add(r7872, r7866, r7871, MPFR_RNDN); ; mpfr_mul(r7874, r7861, r7873, MPFR_RNDN); mpfr_mul(r7875, r7861, r7874, MPFR_RNDN); mpfr_pow(r7876, r7868, r7860, MPFR_RNDN); mpfr_mul(r7877, r7876, r7861, MPFR_RNDN); mpfr_mul(r7878, r7875, r7877, MPFR_RNDN); ; mpfr_mul(r7880, r7879, r7861, MPFR_RNDN); mpfr_mul(r7881, r7868, r7868, MPFR_RNDN); mpfr_mul(r7882, r7880, r7881, MPFR_RNDN); mpfr_sub(r7883, r7878, r7882, MPFR_RNDN); mpfr_sub(r7884, r7872, r7883, MPFR_RNDN); mpfr_neg(r7885, r7884, MPFR_RNDN); return mpfr_get_d(r7885, MPFR_RNDN); }