#include #include #include #include #include char *name = "Jmat.Real.gamma, branch z less than 0.5"; double f_if(float z) { float r9691 = atan2(1.0, 0.0); float r9692 = z; float r9693 = r9691 * r9692; float r9694 = sin(r9693); float r9695 = r9691 / r9694; float r9696 = 2.0f; float r9697 = r9691 * r9696; float r9698 = sqrt(r9697); float r9699 = 1.0f; float r9700 = r9699 - r9692; float r9701 = r9700 - r9699; float r9702 = 7.0f; float r9703 = r9701 + r9702; float r9704 = 0.5f; float r9705 = r9703 + r9704; float r9706 = r9701 + r9704; float r9707 = pow(r9705, r9706); float r9708 = r9698 * r9707; float r9709 = -r9705; float r9710 = exp(r9709); float r9711 = r9708 * r9710; float r9712 = 0.9999999999998099f; float r9713 = 676.5203681218851f; float r9714 = r9701 + r9699; float r9715 = r9713 / r9714; float r9716 = r9712 + r9715; float r9717 = -1259.1392167224028f; float r9718 = r9701 + r9696; float r9719 = r9717 / r9718; float r9720 = r9716 + r9719; float r9721 = 771.3234287776531f; float r9722 = 3.0f; float r9723 = r9701 + r9722; float r9724 = r9721 / r9723; float r9725 = r9720 + r9724; float r9726 = -176.6150291621406f; float r9727 = 4.0f; float r9728 = r9701 + r9727; float r9729 = r9726 / r9728; float r9730 = r9725 + r9729; float r9731 = 12.507343278686905f; float r9732 = 5.0f; float r9733 = r9701 + r9732; float r9734 = r9731 / r9733; float r9735 = r9730 + r9734; float r9736 = -0.13857109526572012f; float r9737 = 6.0f; float r9738 = r9701 + r9737; float r9739 = r9736 / r9738; float r9740 = r9735 + r9739; float r9741 = 9.984369578019572e-06f; float r9742 = r9741 / r9703; float r9743 = r9740 + r9742; float r9744 = 1.5056327351493116e-07f; float r9745 = 8.0f; float r9746 = r9701 + r9745; float r9747 = r9744 / r9746; float r9748 = r9743 + r9747; float r9749 = r9711 * r9748; float r9750 = r9695 * r9749; return r9750; } double f_id(double z) { double r9751 = atan2(1.0, 0.0); double r9752 = z; double r9753 = r9751 * r9752; double r9754 = sin(r9753); double r9755 = r9751 / r9754; double r9756 = 2.0; double r9757 = r9751 * r9756; double r9758 = sqrt(r9757); double r9759 = 1.0; double r9760 = r9759 - r9752; double r9761 = r9760 - r9759; double r9762 = 7.0; double r9763 = r9761 + r9762; double r9764 = 0.5; double r9765 = r9763 + r9764; double r9766 = r9761 + r9764; double r9767 = pow(r9765, r9766); double r9768 = r9758 * r9767; double r9769 = -r9765; double r9770 = exp(r9769); double r9771 = r9768 * r9770; double r9772 = 0.9999999999998099; double r9773 = 676.5203681218851; double r9774 = r9761 + r9759; double r9775 = r9773 / r9774; double r9776 = r9772 + r9775; double r9777 = -1259.1392167224028; double r9778 = r9761 + r9756; double r9779 = r9777 / r9778; double r9780 = r9776 + r9779; double r9781 = 771.3234287776531; double r9782 = 3.0; double r9783 = r9761 + r9782; double r9784 = r9781 / r9783; double r9785 = r9780 + r9784; double r9786 = -176.6150291621406; double r9787 = 4.0; double r9788 = r9761 + r9787; double r9789 = r9786 / r9788; double r9790 = r9785 + r9789; double r9791 = 12.507343278686905; double r9792 = 5.0; double r9793 = r9761 + r9792; double r9794 = r9791 / r9793; double r9795 = r9790 + r9794; double r9796 = -0.13857109526572012; double r9797 = 6.0; double r9798 = r9761 + r9797; double r9799 = r9796 / r9798; double r9800 = r9795 + r9799; double r9801 = 9.984369578019572e-06; double r9802 = r9801 / r9763; double r9803 = r9800 + r9802; double r9804 = 1.5056327351493116e-07; double r9805 = 8.0; double r9806 = r9761 + r9805; double r9807 = r9804 / r9806; double r9808 = r9803 + r9807; double r9809 = r9771 * r9808; double r9810 = r9755 * r9809; return r9810; } double f_of(float z) { float r9811 = 12.507343278686905f; float r9812 = 1.0f; float r9813 = z; float r9814 = r9812 - r9813; float r9815 = 5.0f; float r9816 = r9812 - r9815; float r9817 = r9814 - r9816; float r9818 = r9811 / r9817; float r9819 = -0.13857109526572012f; float r9820 = 6.0f; float r9821 = r9820 + r9812; float r9822 = r9812 + r9813; float r9823 = r9821 - r9822; float r9824 = r9819 / r9823; float r9825 = r9818 + r9824; float r9826 = -176.6150291621406f; float r9827 = 4.0f; float r9828 = r9812 - r9827; float r9829 = r9814 - r9828; float r9830 = r9826 / r9829; float r9831 = 771.3234287776531f; float r9832 = 3.0f; float r9833 = r9812 - r9832; float r9834 = r9814 - r9833; float r9835 = r9831 / r9834; float r9836 = r9830 + r9835; float r9837 = r9825 + r9836; float r9838 = 0.9999999999998099f; float r9839 = -1259.1392167224028f; float r9840 = 2.0f; float r9841 = r9840 + r9812; float r9842 = r9841 - r9822; float r9843 = r9839 / r9842; float r9844 = 676.5203681218851f; float r9845 = r9844 / r9814; float r9846 = r9843 + r9845; float r9847 = r9838 + r9846; float r9848 = r9837 + r9847; float r9849 = 9.984369578019572e-06f; float r9850 = 7.0f; float r9851 = r9812 + r9850; float r9852 = r9851 - r9822; float r9853 = r9849 / r9852; float r9854 = 1.5056327351493116e-07f; float r9855 = 8.0f; float r9856 = r9855 + r9812; float r9857 = r9856 - r9822; float r9858 = r9854 / r9857; float r9859 = r9853 + r9858; float r9860 = r9848 + r9859; float r9861 = atan2(1.0, 0.0); float r9862 = r9840 * r9861; float r9863 = sqrt(r9862); float r9864 = r9861 * r9813; float r9865 = sin(r9864); float r9866 = r9861 / r9865; float r9867 = r9863 * r9866; float r9868 = 0.5f; float r9869 = r9812 + r9868; float r9870 = r9869 - r9822; float r9871 = r9870 + r9850; float r9872 = pow(r9871, r9870); float r9873 = r9869 + r9850; float r9874 = exp(r9873); float r9875 = r9812 / r9874; float r9876 = r9872 * r9875; float r9877 = exp(r9822); float r9878 = r9876 * r9877; float r9879 = r9867 * r9878; float r9880 = r9860 * r9879; return r9880; } double f_od(double z) { double r9881 = 12.507343278686905; double r9882 = 1.0; double r9883 = z; double r9884 = r9882 - r9883; double r9885 = 5.0; double r9886 = r9882 - r9885; double r9887 = r9884 - r9886; double r9888 = r9881 / r9887; double r9889 = -0.13857109526572012; double r9890 = 6.0; double r9891 = r9890 + r9882; double r9892 = r9882 + r9883; double r9893 = r9891 - r9892; double r9894 = r9889 / r9893; double r9895 = r9888 + r9894; double r9896 = -176.6150291621406; double r9897 = 4.0; double r9898 = r9882 - r9897; double r9899 = r9884 - r9898; double r9900 = r9896 / r9899; double r9901 = 771.3234287776531; double r9902 = 3.0; double r9903 = r9882 - r9902; double r9904 = r9884 - r9903; double r9905 = r9901 / r9904; double r9906 = r9900 + r9905; double r9907 = r9895 + r9906; double r9908 = 0.9999999999998099; double r9909 = -1259.1392167224028; double r9910 = 2.0; double r9911 = r9910 + r9882; double r9912 = r9911 - r9892; double r9913 = r9909 / r9912; double r9914 = 676.5203681218851; double r9915 = r9914 / r9884; double r9916 = r9913 + r9915; double r9917 = r9908 + r9916; double r9918 = r9907 + r9917; double r9919 = 9.984369578019572e-06; double r9920 = 7.0; double r9921 = r9882 + r9920; double r9922 = r9921 - r9892; double r9923 = r9919 / r9922; double r9924 = 1.5056327351493116e-07; double r9925 = 8.0; double r9926 = r9925 + r9882; double r9927 = r9926 - r9892; double r9928 = r9924 / r9927; double r9929 = r9923 + r9928; double r9930 = r9918 + r9929; double r9931 = atan2(1.0, 0.0); double r9932 = r9910 * r9931; double r9933 = sqrt(r9932); double r9934 = r9931 * r9883; double r9935 = sin(r9934); double r9936 = r9931 / r9935; double r9937 = r9933 * r9936; double r9938 = 0.5; double r9939 = r9882 + r9938; double r9940 = r9939 - r9892; double r9941 = r9940 + r9920; double r9942 = pow(r9941, r9940); double r9943 = r9939 + r9920; double r9944 = exp(r9943); double r9945 = r9882 / r9944; double r9946 = r9942 * r9945; double r9947 = exp(r9892); double r9948 = r9946 * r9947; double r9949 = r9937 * r9948; double r9950 = r9930 * r9949; return r9950; } 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 r9951, r9952, r9953, r9954, r9955, r9956, r9957, r9958, r9959, r9960, r9961, r9962, r9963, r9964, r9965, r9966, r9967, r9968, r9969, r9970, r9971, r9972, r9973, r9974, r9975, r9976, r9977, r9978, r9979, r9980, r9981, r9982, r9983, r9984, r9985, r9986, r9987, r9988, r9989, r9990, r9991, r9992, r9993, r9994, r9995, r9996, r9997, r9998, r9999, r10000, r10001, r10002, r10003, r10004, r10005, r10006, r10007, r10008, r10009, r10010; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r9951); mpfr_init(r9952); mpfr_init(r9953); mpfr_init(r9954); mpfr_init(r9955); mpfr_init_set_str(r9956, "2", 10, MPFR_RNDN); mpfr_init(r9957); mpfr_init(r9958); mpfr_init_set_str(r9959, "1", 10, MPFR_RNDN); mpfr_init(r9960); mpfr_init(r9961); mpfr_init_set_str(r9962, "7", 10, MPFR_RNDN); mpfr_init(r9963); mpfr_init_set_str(r9964, "0.5", 10, MPFR_RNDN); mpfr_init(r9965); mpfr_init(r9966); mpfr_init(r9967); mpfr_init(r9968); mpfr_init(r9969); mpfr_init(r9970); mpfr_init(r9971); mpfr_init_set_str(r9972, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r9973, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r9974); mpfr_init(r9975); mpfr_init(r9976); mpfr_init_set_str(r9977, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init(r9978); mpfr_init(r9979); mpfr_init(r9980); mpfr_init_set_str(r9981, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r9982, "3", 10, MPFR_RNDN); mpfr_init(r9983); mpfr_init(r9984); mpfr_init(r9985); mpfr_init_set_str(r9986, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r9987, "4", 10, MPFR_RNDN); mpfr_init(r9988); mpfr_init(r9989); mpfr_init(r9990); mpfr_init_set_str(r9991, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r9992, "5", 10, MPFR_RNDN); mpfr_init(r9993); mpfr_init(r9994); mpfr_init(r9995); mpfr_init_set_str(r9996, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r9997, "6", 10, MPFR_RNDN); mpfr_init(r9998); mpfr_init(r9999); mpfr_init(r10000); mpfr_init_set_str(r10001, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init(r10002); mpfr_init(r10003); mpfr_init_set_str(r10004, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r10005, "8", 10, MPFR_RNDN); mpfr_init(r10006); mpfr_init(r10007); mpfr_init(r10008); mpfr_init(r10009); mpfr_init(r10010); } double f_im(double z) { mpfr_const_pi(r9951, MPFR_RNDN); mpfr_set_d(r9952, z, MPFR_RNDN); mpfr_mul(r9953, r9951, r9952, MPFR_RNDN); mpfr_sin(r9954, r9953, MPFR_RNDN); mpfr_div(r9955, r9951, r9954, MPFR_RNDN); ; mpfr_mul(r9957, r9951, r9956, MPFR_RNDN); mpfr_sqrt(r9958, r9957, MPFR_RNDN); ; mpfr_sub(r9960, r9959, r9952, MPFR_RNDN); mpfr_sub(r9961, r9960, r9959, MPFR_RNDN); ; mpfr_add(r9963, r9961, r9962, MPFR_RNDN); ; mpfr_add(r9965, r9963, r9964, MPFR_RNDN); mpfr_add(r9966, r9961, r9964, MPFR_RNDN); mpfr_pow(r9967, r9965, r9966, MPFR_RNDN); mpfr_mul(r9968, r9958, r9967, MPFR_RNDN); mpfr_neg(r9969, r9965, MPFR_RNDN); mpfr_exp(r9970, r9969, MPFR_RNDN); mpfr_mul(r9971, r9968, r9970, MPFR_RNDN); ; ; mpfr_add(r9974, r9961, r9959, MPFR_RNDN); mpfr_div(r9975, r9973, r9974, MPFR_RNDN); mpfr_add(r9976, r9972, r9975, MPFR_RNDN); ; mpfr_add(r9978, r9961, r9956, MPFR_RNDN); mpfr_div(r9979, r9977, r9978, MPFR_RNDN); mpfr_add(r9980, r9976, r9979, MPFR_RNDN); ; ; mpfr_add(r9983, r9961, r9982, MPFR_RNDN); mpfr_div(r9984, r9981, r9983, MPFR_RNDN); mpfr_add(r9985, r9980, r9984, MPFR_RNDN); ; ; mpfr_add(r9988, r9961, r9987, MPFR_RNDN); mpfr_div(r9989, r9986, r9988, MPFR_RNDN); mpfr_add(r9990, r9985, r9989, MPFR_RNDN); ; ; mpfr_add(r9993, r9961, r9992, MPFR_RNDN); mpfr_div(r9994, r9991, r9993, MPFR_RNDN); mpfr_add(r9995, r9990, r9994, MPFR_RNDN); ; ; mpfr_add(r9998, r9961, r9997, MPFR_RNDN); mpfr_div(r9999, r9996, r9998, MPFR_RNDN); mpfr_add(r10000, r9995, r9999, MPFR_RNDN); ; mpfr_div(r10002, r10001, r9963, MPFR_RNDN); mpfr_add(r10003, r10000, r10002, MPFR_RNDN); ; ; mpfr_add(r10006, r9961, r10005, MPFR_RNDN); mpfr_div(r10007, r10004, r10006, MPFR_RNDN); mpfr_add(r10008, r10003, r10007, MPFR_RNDN); mpfr_mul(r10009, r9971, r10008, MPFR_RNDN); mpfr_mul(r10010, r9955, r10009, MPFR_RNDN); return mpfr_get_d(r10010, MPFR_RNDN); } static mpfr_t r10011, r10012, r10013, r10014, r10015, r10016, r10017, r10018, r10019, r10020, r10021, r10022, r10023, r10024, r10025, r10026, r10027, r10028, r10029, r10030, r10031, r10032, r10033, r10034, r10035, r10036, r10037, r10038, r10039, r10040, r10041, r10042, r10043, r10044, r10045, r10046, r10047, r10048, r10049, r10050, r10051, r10052, r10053, r10054, r10055, r10056, r10057, r10058, r10059, r10060, r10061, r10062, r10063, r10064, r10065, r10066, r10067, r10068, r10069, r10070, r10071, r10072, r10073, r10074, r10075, r10076, r10077, r10078, r10079, r10080; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r10011, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r10012, "1", 10, MPFR_RNDN); mpfr_init(r10013); mpfr_init(r10014); mpfr_init_set_str(r10015, "5", 10, MPFR_RNDN); mpfr_init(r10016); mpfr_init(r10017); mpfr_init(r10018); mpfr_init_set_str(r10019, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r10020, "6", 10, MPFR_RNDN); mpfr_init(r10021); mpfr_init(r10022); mpfr_init(r10023); mpfr_init(r10024); mpfr_init(r10025); mpfr_init_set_str(r10026, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r10027, "4", 10, MPFR_RNDN); mpfr_init(r10028); mpfr_init(r10029); mpfr_init(r10030); mpfr_init_set_str(r10031, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r10032, "3", 10, MPFR_RNDN); mpfr_init(r10033); mpfr_init(r10034); mpfr_init(r10035); mpfr_init(r10036); mpfr_init(r10037); mpfr_init_set_str(r10038, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r10039, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r10040, "2", 10, MPFR_RNDN); mpfr_init(r10041); mpfr_init(r10042); mpfr_init(r10043); mpfr_init_set_str(r10044, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r10045); mpfr_init(r10046); mpfr_init(r10047); mpfr_init(r10048); mpfr_init_set_str(r10049, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init_set_str(r10050, "7", 10, MPFR_RNDN); mpfr_init(r10051); mpfr_init(r10052); mpfr_init(r10053); mpfr_init_set_str(r10054, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r10055, "8", 10, MPFR_RNDN); mpfr_init(r10056); mpfr_init(r10057); mpfr_init(r10058); mpfr_init(r10059); mpfr_init(r10060); mpfr_init(r10061); mpfr_init(r10062); mpfr_init(r10063); mpfr_init(r10064); mpfr_init(r10065); mpfr_init(r10066); mpfr_init(r10067); mpfr_init_set_str(r10068, "0.5", 10, MPFR_RNDN); mpfr_init(r10069); mpfr_init(r10070); mpfr_init(r10071); mpfr_init(r10072); mpfr_init(r10073); mpfr_init(r10074); mpfr_init(r10075); mpfr_init(r10076); mpfr_init(r10077); mpfr_init(r10078); mpfr_init(r10079); mpfr_init(r10080); } double f_fm(double z) { ; ; mpfr_set_d(r10013, z, MPFR_RNDN); mpfr_sub(r10014, r10012, r10013, MPFR_RNDN); ; mpfr_sub(r10016, r10012, r10015, MPFR_RNDN); mpfr_sub(r10017, r10014, r10016, MPFR_RNDN); mpfr_div(r10018, r10011, r10017, MPFR_RNDN); ; ; mpfr_add(r10021, r10020, r10012, MPFR_RNDN); mpfr_add(r10022, r10012, r10013, MPFR_RNDN); mpfr_sub(r10023, r10021, r10022, MPFR_RNDN); mpfr_div(r10024, r10019, r10023, MPFR_RNDN); mpfr_add(r10025, r10018, r10024, MPFR_RNDN); ; ; mpfr_sub(r10028, r10012, r10027, MPFR_RNDN); mpfr_sub(r10029, r10014, r10028, MPFR_RNDN); mpfr_div(r10030, r10026, r10029, MPFR_RNDN); ; ; mpfr_sub(r10033, r10012, r10032, MPFR_RNDN); mpfr_sub(r10034, r10014, r10033, MPFR_RNDN); mpfr_div(r10035, r10031, r10034, MPFR_RNDN); mpfr_add(r10036, r10030, r10035, MPFR_RNDN); mpfr_add(r10037, r10025, r10036, MPFR_RNDN); ; ; ; mpfr_add(r10041, r10040, r10012, MPFR_RNDN); mpfr_sub(r10042, r10041, r10022, MPFR_RNDN); mpfr_div(r10043, r10039, r10042, MPFR_RNDN); ; mpfr_div(r10045, r10044, r10014, MPFR_RNDN); mpfr_add(r10046, r10043, r10045, MPFR_RNDN); mpfr_add(r10047, r10038, r10046, MPFR_RNDN); mpfr_add(r10048, r10037, r10047, MPFR_RNDN); ; ; mpfr_add(r10051, r10012, r10050, MPFR_RNDN); mpfr_sub(r10052, r10051, r10022, MPFR_RNDN); mpfr_div(r10053, r10049, r10052, MPFR_RNDN); ; ; mpfr_add(r10056, r10055, r10012, MPFR_RNDN); mpfr_sub(r10057, r10056, r10022, MPFR_RNDN); mpfr_div(r10058, r10054, r10057, MPFR_RNDN); mpfr_add(r10059, r10053, r10058, MPFR_RNDN); mpfr_add(r10060, r10048, r10059, MPFR_RNDN); mpfr_const_pi(r10061, MPFR_RNDN); mpfr_mul(r10062, r10040, r10061, MPFR_RNDN); mpfr_sqrt(r10063, r10062, MPFR_RNDN); mpfr_mul(r10064, r10061, r10013, MPFR_RNDN); mpfr_sin(r10065, r10064, MPFR_RNDN); mpfr_div(r10066, r10061, r10065, MPFR_RNDN); mpfr_mul(r10067, r10063, r10066, MPFR_RNDN); ; mpfr_add(r10069, r10012, r10068, MPFR_RNDN); mpfr_sub(r10070, r10069, r10022, MPFR_RNDN); mpfr_add(r10071, r10070, r10050, MPFR_RNDN); mpfr_pow(r10072, r10071, r10070, MPFR_RNDN); mpfr_add(r10073, r10069, r10050, MPFR_RNDN); mpfr_exp(r10074, r10073, MPFR_RNDN); mpfr_div(r10075, r10012, r10074, MPFR_RNDN); mpfr_mul(r10076, r10072, r10075, MPFR_RNDN); mpfr_exp(r10077, r10022, MPFR_RNDN); mpfr_mul(r10078, r10076, r10077, MPFR_RNDN); mpfr_mul(r10079, r10067, r10078, MPFR_RNDN); mpfr_mul(r10080, r10060, r10079, MPFR_RNDN); return mpfr_get_d(r10080, MPFR_RNDN); } static mpfr_t r10081, r10082, r10083, r10084, r10085, r10086, r10087, r10088, r10089, r10090, r10091, r10092, r10093, r10094, r10095, r10096, r10097, r10098, r10099, r10100, r10101, r10102, r10103, r10104, r10105, r10106, r10107, r10108, r10109, r10110, r10111, r10112, r10113, r10114, r10115, r10116, r10117, r10118, r10119, r10120, r10121, r10122, r10123, r10124, r10125, r10126, r10127, r10128, r10129, r10130, r10131, r10132, r10133, r10134, r10135, r10136, r10137, r10138, r10139, r10140, r10141, r10142, r10143, r10144, r10145, r10146, r10147, r10148, r10149, r10150; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r10081, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r10082, "1", 10, MPFR_RNDN); mpfr_init(r10083); mpfr_init(r10084); mpfr_init_set_str(r10085, "5", 10, MPFR_RNDN); mpfr_init(r10086); mpfr_init(r10087); mpfr_init(r10088); mpfr_init_set_str(r10089, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r10090, "6", 10, MPFR_RNDN); mpfr_init(r10091); mpfr_init(r10092); mpfr_init(r10093); mpfr_init(r10094); mpfr_init(r10095); mpfr_init_set_str(r10096, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r10097, "4", 10, MPFR_RNDN); mpfr_init(r10098); mpfr_init(r10099); mpfr_init(r10100); mpfr_init_set_str(r10101, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r10102, "3", 10, MPFR_RNDN); mpfr_init(r10103); mpfr_init(r10104); mpfr_init(r10105); mpfr_init(r10106); mpfr_init(r10107); mpfr_init_set_str(r10108, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r10109, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r10110, "2", 10, MPFR_RNDN); mpfr_init(r10111); mpfr_init(r10112); mpfr_init(r10113); mpfr_init_set_str(r10114, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r10115); mpfr_init(r10116); mpfr_init(r10117); mpfr_init(r10118); mpfr_init_set_str(r10119, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init_set_str(r10120, "7", 10, MPFR_RNDN); mpfr_init(r10121); mpfr_init(r10122); mpfr_init(r10123); mpfr_init_set_str(r10124, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r10125, "8", 10, MPFR_RNDN); mpfr_init(r10126); mpfr_init(r10127); mpfr_init(r10128); mpfr_init(r10129); mpfr_init(r10130); mpfr_init(r10131); mpfr_init(r10132); mpfr_init(r10133); mpfr_init(r10134); mpfr_init(r10135); mpfr_init(r10136); mpfr_init(r10137); mpfr_init_set_str(r10138, "0.5", 10, MPFR_RNDN); mpfr_init(r10139); mpfr_init(r10140); mpfr_init(r10141); mpfr_init(r10142); mpfr_init(r10143); mpfr_init(r10144); mpfr_init(r10145); mpfr_init(r10146); mpfr_init(r10147); mpfr_init(r10148); mpfr_init(r10149); mpfr_init(r10150); } double f_dm(double z) { ; ; mpfr_set_d(r10083, z, MPFR_RNDN); mpfr_sub(r10084, r10082, r10083, MPFR_RNDN); ; mpfr_sub(r10086, r10082, r10085, MPFR_RNDN); mpfr_sub(r10087, r10084, r10086, MPFR_RNDN); mpfr_div(r10088, r10081, r10087, MPFR_RNDN); ; ; mpfr_add(r10091, r10090, r10082, MPFR_RNDN); mpfr_add(r10092, r10082, r10083, MPFR_RNDN); mpfr_sub(r10093, r10091, r10092, MPFR_RNDN); mpfr_div(r10094, r10089, r10093, MPFR_RNDN); mpfr_add(r10095, r10088, r10094, MPFR_RNDN); ; ; mpfr_sub(r10098, r10082, r10097, MPFR_RNDN); mpfr_sub(r10099, r10084, r10098, MPFR_RNDN); mpfr_div(r10100, r10096, r10099, MPFR_RNDN); ; ; mpfr_sub(r10103, r10082, r10102, MPFR_RNDN); mpfr_sub(r10104, r10084, r10103, MPFR_RNDN); mpfr_div(r10105, r10101, r10104, MPFR_RNDN); mpfr_add(r10106, r10100, r10105, MPFR_RNDN); mpfr_add(r10107, r10095, r10106, MPFR_RNDN); ; ; ; mpfr_add(r10111, r10110, r10082, MPFR_RNDN); mpfr_sub(r10112, r10111, r10092, MPFR_RNDN); mpfr_div(r10113, r10109, r10112, MPFR_RNDN); ; mpfr_div(r10115, r10114, r10084, MPFR_RNDN); mpfr_add(r10116, r10113, r10115, MPFR_RNDN); mpfr_add(r10117, r10108, r10116, MPFR_RNDN); mpfr_add(r10118, r10107, r10117, MPFR_RNDN); ; ; mpfr_add(r10121, r10082, r10120, MPFR_RNDN); mpfr_sub(r10122, r10121, r10092, MPFR_RNDN); mpfr_div(r10123, r10119, r10122, MPFR_RNDN); ; ; mpfr_add(r10126, r10125, r10082, MPFR_RNDN); mpfr_sub(r10127, r10126, r10092, MPFR_RNDN); mpfr_div(r10128, r10124, r10127, MPFR_RNDN); mpfr_add(r10129, r10123, r10128, MPFR_RNDN); mpfr_add(r10130, r10118, r10129, MPFR_RNDN); mpfr_const_pi(r10131, MPFR_RNDN); mpfr_mul(r10132, r10110, r10131, MPFR_RNDN); mpfr_sqrt(r10133, r10132, MPFR_RNDN); mpfr_mul(r10134, r10131, r10083, MPFR_RNDN); mpfr_sin(r10135, r10134, MPFR_RNDN); mpfr_div(r10136, r10131, r10135, MPFR_RNDN); mpfr_mul(r10137, r10133, r10136, MPFR_RNDN); ; mpfr_add(r10139, r10082, r10138, MPFR_RNDN); mpfr_sub(r10140, r10139, r10092, MPFR_RNDN); mpfr_add(r10141, r10140, r10120, MPFR_RNDN); mpfr_pow(r10142, r10141, r10140, MPFR_RNDN); mpfr_add(r10143, r10139, r10120, MPFR_RNDN); mpfr_exp(r10144, r10143, MPFR_RNDN); mpfr_div(r10145, r10082, r10144, MPFR_RNDN); mpfr_mul(r10146, r10142, r10145, MPFR_RNDN); mpfr_exp(r10147, r10092, MPFR_RNDN); mpfr_mul(r10148, r10146, r10147, MPFR_RNDN); mpfr_mul(r10149, r10137, r10148, MPFR_RNDN); mpfr_mul(r10150, r10130, r10149, MPFR_RNDN); return mpfr_get_d(r10150, MPFR_RNDN); }