#include #include #include #include #include char *name = "Jmat.Real.gamma, branch z less than 0.5"; double f_if(float z) { float r26735 = atan2(1.0, 0.0); float r26736 = z; float r26737 = r26735 * r26736; float r26738 = sin(r26737); float r26739 = r26735 / r26738; float r26740 = 2; float r26741 = r26735 * r26740; float r26742 = sqrt(r26741); float r26743 = 1; float r26744 = r26743 - r26736; float r26745 = r26744 - r26743; float r26746 = 7; float r26747 = r26745 + r26746; float r26748 = 0.5; float r26749 = r26747 + r26748; float r26750 = r26745 + r26748; float r26751 = pow(r26749, r26750); float r26752 = r26742 * r26751; float r26753 = -r26749; float r26754 = exp(r26753); float r26755 = r26752 * r26754; float r26756 = 0.9999999999998099; float r26757 = 676.5203681218851; float r26758 = r26745 + r26743; float r26759 = r26757 / r26758; float r26760 = r26756 + r26759; float r26761 = -1259.1392167224028; float r26762 = r26745 + r26740; float r26763 = r26761 / r26762; float r26764 = r26760 + r26763; float r26765 = 771.3234287776531; float r26766 = 3; float r26767 = r26745 + r26766; float r26768 = r26765 / r26767; float r26769 = r26764 + r26768; float r26770 = -176.6150291621406; float r26771 = 4; float r26772 = r26745 + r26771; float r26773 = r26770 / r26772; float r26774 = r26769 + r26773; float r26775 = 12.507343278686905; float r26776 = 5; float r26777 = r26745 + r26776; float r26778 = r26775 / r26777; float r26779 = r26774 + r26778; float r26780 = -0.13857109526572012; float r26781 = 6; float r26782 = r26745 + r26781; float r26783 = r26780 / r26782; float r26784 = r26779 + r26783; float r26785 = 9.984369578019572e-06; float r26786 = r26785 / r26747; float r26787 = r26784 + r26786; float r26788 = 1.5056327351493116e-07; float r26789 = 8; float r26790 = r26745 + r26789; float r26791 = r26788 / r26790; float r26792 = r26787 + r26791; float r26793 = r26755 * r26792; float r26794 = r26739 * r26793; return r26794; } double f_id(double z) { double r26795 = atan2(1.0, 0.0); double r26796 = z; double r26797 = r26795 * r26796; double r26798 = sin(r26797); double r26799 = r26795 / r26798; double r26800 = 2; double r26801 = r26795 * r26800; double r26802 = sqrt(r26801); double r26803 = 1; double r26804 = r26803 - r26796; double r26805 = r26804 - r26803; double r26806 = 7; double r26807 = r26805 + r26806; double r26808 = 0.5; double r26809 = r26807 + r26808; double r26810 = r26805 + r26808; double r26811 = pow(r26809, r26810); double r26812 = r26802 * r26811; double r26813 = -r26809; double r26814 = exp(r26813); double r26815 = r26812 * r26814; double r26816 = 0.9999999999998099; double r26817 = 676.5203681218851; double r26818 = r26805 + r26803; double r26819 = r26817 / r26818; double r26820 = r26816 + r26819; double r26821 = -1259.1392167224028; double r26822 = r26805 + r26800; double r26823 = r26821 / r26822; double r26824 = r26820 + r26823; double r26825 = 771.3234287776531; double r26826 = 3; double r26827 = r26805 + r26826; double r26828 = r26825 / r26827; double r26829 = r26824 + r26828; double r26830 = -176.6150291621406; double r26831 = 4; double r26832 = r26805 + r26831; double r26833 = r26830 / r26832; double r26834 = r26829 + r26833; double r26835 = 12.507343278686905; double r26836 = 5; double r26837 = r26805 + r26836; double r26838 = r26835 / r26837; double r26839 = r26834 + r26838; double r26840 = -0.13857109526572012; double r26841 = 6; double r26842 = r26805 + r26841; double r26843 = r26840 / r26842; double r26844 = r26839 + r26843; double r26845 = 9.984369578019572e-06; double r26846 = r26845 / r26807; double r26847 = r26844 + r26846; double r26848 = 1.5056327351493116e-07; double r26849 = 8; double r26850 = r26805 + r26849; double r26851 = r26848 / r26850; double r26852 = r26847 + r26851; double r26853 = r26815 * r26852; double r26854 = r26799 * r26853; return r26854; } double f_of(float z) { float r26855 = 1.5056327351493116e-07; float r26856 = 8; float r26857 = 1; float r26858 = r26856 + r26857; float r26859 = z; float r26860 = r26857 + r26859; float r26861 = r26858 - r26860; float r26862 = r26855 / r26861; float r26863 = 771.3234287776531; float r26864 = r26857 - r26859; float r26865 = 3; float r26866 = r26857 - r26865; float r26867 = r26864 - r26866; float r26868 = r26863 / r26867; float r26869 = 0.9999999999998099; float r26870 = 676.5203681218851; float r26871 = r26870 / r26864; float r26872 = r26869 + r26871; float r26873 = -1259.1392167224028; float r26874 = 2; float r26875 = r26857 + r26874; float r26876 = r26875 - r26860; float r26877 = r26873 / r26876; float r26878 = r26872 + r26877; float r26879 = r26868 + r26878; float r26880 = 9.984369578019572e-06; float r26881 = 7; float r26882 = r26881 + r26857; float r26883 = r26882 - r26860; float r26884 = r26880 / r26883; float r26885 = -0.13857109526572012; float r26886 = 6; float r26887 = r26886 + r26857; float r26888 = r26887 - r26860; float r26889 = r26885 / r26888; float r26890 = r26884 + r26889; float r26891 = -176.6150291621406; float r26892 = 4; float r26893 = r26857 - r26892; float r26894 = r26864 - r26893; float r26895 = r26891 / r26894; float r26896 = 12.507343278686905; float r26897 = 5; float r26898 = r26857 - r26897; float r26899 = r26864 - r26898; float r26900 = r26896 / r26899; float r26901 = r26895 + r26900; float r26902 = r26890 + r26901; float r26903 = r26879 + r26902; float r26904 = r26862 + r26903; float r26905 = atan2(1.0, 0.0); float r26906 = r26859 * r26905; float r26907 = sin(r26906); float r26908 = r26905 / r26907; float r26909 = 7.5; float r26910 = r26909 - r26859; float r26911 = 0.5; float r26912 = r26911 - r26859; float r26913 = pow(r26910, r26912); float r26914 = sqrt(r26874); float r26915 = r26913 * r26914; float r26916 = r26859 - r26909; float r26917 = exp(r26916); float r26918 = sqrt(r26905); float r26919 = sqrt(r26918); float r26920 = r26919 * r26919; float r26921 = r26917 * r26920; float r26922 = r26915 * r26921; float r26923 = r26908 * r26922; float r26924 = r26904 * r26923; return r26924; } double f_od(double z) { double r26925 = 1.5056327351493116e-07; double r26926 = 8; double r26927 = 1; double r26928 = r26926 + r26927; double r26929 = z; double r26930 = r26927 + r26929; double r26931 = r26928 - r26930; double r26932 = r26925 / r26931; double r26933 = 771.3234287776531; double r26934 = r26927 - r26929; double r26935 = 3; double r26936 = r26927 - r26935; double r26937 = r26934 - r26936; double r26938 = r26933 / r26937; double r26939 = 0.9999999999998099; double r26940 = 676.5203681218851; double r26941 = r26940 / r26934; double r26942 = r26939 + r26941; double r26943 = -1259.1392167224028; double r26944 = 2; double r26945 = r26927 + r26944; double r26946 = r26945 - r26930; double r26947 = r26943 / r26946; double r26948 = r26942 + r26947; double r26949 = r26938 + r26948; double r26950 = 9.984369578019572e-06; double r26951 = 7; double r26952 = r26951 + r26927; double r26953 = r26952 - r26930; double r26954 = r26950 / r26953; double r26955 = -0.13857109526572012; double r26956 = 6; double r26957 = r26956 + r26927; double r26958 = r26957 - r26930; double r26959 = r26955 / r26958; double r26960 = r26954 + r26959; double r26961 = -176.6150291621406; double r26962 = 4; double r26963 = r26927 - r26962; double r26964 = r26934 - r26963; double r26965 = r26961 / r26964; double r26966 = 12.507343278686905; double r26967 = 5; double r26968 = r26927 - r26967; double r26969 = r26934 - r26968; double r26970 = r26966 / r26969; double r26971 = r26965 + r26970; double r26972 = r26960 + r26971; double r26973 = r26949 + r26972; double r26974 = r26932 + r26973; double r26975 = atan2(1.0, 0.0); double r26976 = r26929 * r26975; double r26977 = sin(r26976); double r26978 = r26975 / r26977; double r26979 = 7.5; double r26980 = r26979 - r26929; double r26981 = 0.5; double r26982 = r26981 - r26929; double r26983 = pow(r26980, r26982); double r26984 = sqrt(r26944); double r26985 = r26983 * r26984; double r26986 = r26929 - r26979; double r26987 = exp(r26986); double r26988 = sqrt(r26975); double r26989 = sqrt(r26988); double r26990 = r26989 * r26989; double r26991 = r26987 * r26990; double r26992 = r26985 * r26991; double r26993 = r26978 * r26992; double r26994 = r26974 * r26993; return r26994; } 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 r26995, r26996, r26997, r26998, r26999, r27000, r27001, r27002, r27003, r27004, r27005, r27006, r27007, r27008, r27009, r27010, r27011, r27012, r27013, r27014, r27015, r27016, r27017, r27018, r27019, r27020, r27021, r27022, r27023, r27024, r27025, r27026, r27027, r27028, r27029, r27030, r27031, r27032, r27033, r27034, r27035, r27036, r27037, r27038, r27039, r27040, r27041, r27042, r27043, r27044, r27045, r27046, r27047, r27048, r27049, r27050, r27051, r27052, r27053, r27054; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26995); mpfr_init(r26996); mpfr_init(r26997); mpfr_init(r26998); mpfr_init(r26999); mpfr_init_set_str(r27000, "2", 10, MPFR_RNDN); mpfr_init(r27001); mpfr_init(r27002); mpfr_init_set_str(r27003, "1", 10, MPFR_RNDN); mpfr_init(r27004); mpfr_init(r27005); mpfr_init_set_str(r27006, "7", 10, MPFR_RNDN); mpfr_init(r27007); mpfr_init_set_str(r27008, "0.5", 10, MPFR_RNDN); mpfr_init(r27009); mpfr_init(r27010); mpfr_init(r27011); mpfr_init(r27012); mpfr_init(r27013); mpfr_init(r27014); mpfr_init(r27015); mpfr_init_set_str(r27016, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r27017, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r27018); mpfr_init(r27019); mpfr_init(r27020); mpfr_init_set_str(r27021, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init(r27022); mpfr_init(r27023); mpfr_init(r27024); mpfr_init_set_str(r27025, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r27026, "3", 10, MPFR_RNDN); mpfr_init(r27027); mpfr_init(r27028); mpfr_init(r27029); mpfr_init_set_str(r27030, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r27031, "4", 10, MPFR_RNDN); mpfr_init(r27032); mpfr_init(r27033); mpfr_init(r27034); mpfr_init_set_str(r27035, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r27036, "5", 10, MPFR_RNDN); mpfr_init(r27037); mpfr_init(r27038); mpfr_init(r27039); mpfr_init_set_str(r27040, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r27041, "6", 10, MPFR_RNDN); mpfr_init(r27042); mpfr_init(r27043); mpfr_init(r27044); mpfr_init_set_str(r27045, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init(r27046); mpfr_init(r27047); mpfr_init_set_str(r27048, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r27049, "8", 10, MPFR_RNDN); mpfr_init(r27050); mpfr_init(r27051); mpfr_init(r27052); mpfr_init(r27053); mpfr_init(r27054); } double f_im(double z) { mpfr_const_pi(r26995, MPFR_RNDN); mpfr_set_d(r26996, z, MPFR_RNDN); mpfr_mul(r26997, r26995, r26996, MPFR_RNDN); mpfr_sin(r26998, r26997, MPFR_RNDN); mpfr_div(r26999, r26995, r26998, MPFR_RNDN); ; mpfr_mul(r27001, r26995, r27000, MPFR_RNDN); mpfr_sqrt(r27002, r27001, MPFR_RNDN); ; mpfr_sub(r27004, r27003, r26996, MPFR_RNDN); mpfr_sub(r27005, r27004, r27003, MPFR_RNDN); ; mpfr_add(r27007, r27005, r27006, MPFR_RNDN); ; mpfr_add(r27009, r27007, r27008, MPFR_RNDN); mpfr_add(r27010, r27005, r27008, MPFR_RNDN); mpfr_pow(r27011, r27009, r27010, MPFR_RNDN); mpfr_mul(r27012, r27002, r27011, MPFR_RNDN); mpfr_neg(r27013, r27009, MPFR_RNDN); mpfr_exp(r27014, r27013, MPFR_RNDN); mpfr_mul(r27015, r27012, r27014, MPFR_RNDN); ; ; mpfr_add(r27018, r27005, r27003, MPFR_RNDN); mpfr_div(r27019, r27017, r27018, MPFR_RNDN); mpfr_add(r27020, r27016, r27019, MPFR_RNDN); ; mpfr_add(r27022, r27005, r27000, MPFR_RNDN); mpfr_div(r27023, r27021, r27022, MPFR_RNDN); mpfr_add(r27024, r27020, r27023, MPFR_RNDN); ; ; mpfr_add(r27027, r27005, r27026, MPFR_RNDN); mpfr_div(r27028, r27025, r27027, MPFR_RNDN); mpfr_add(r27029, r27024, r27028, MPFR_RNDN); ; ; mpfr_add(r27032, r27005, r27031, MPFR_RNDN); mpfr_div(r27033, r27030, r27032, MPFR_RNDN); mpfr_add(r27034, r27029, r27033, MPFR_RNDN); ; ; mpfr_add(r27037, r27005, r27036, MPFR_RNDN); mpfr_div(r27038, r27035, r27037, MPFR_RNDN); mpfr_add(r27039, r27034, r27038, MPFR_RNDN); ; ; mpfr_add(r27042, r27005, r27041, MPFR_RNDN); mpfr_div(r27043, r27040, r27042, MPFR_RNDN); mpfr_add(r27044, r27039, r27043, MPFR_RNDN); ; mpfr_div(r27046, r27045, r27007, MPFR_RNDN); mpfr_add(r27047, r27044, r27046, MPFR_RNDN); ; ; mpfr_add(r27050, r27005, r27049, MPFR_RNDN); mpfr_div(r27051, r27048, r27050, MPFR_RNDN); mpfr_add(r27052, r27047, r27051, MPFR_RNDN); mpfr_mul(r27053, r27015, r27052, MPFR_RNDN); mpfr_mul(r27054, r26999, r27053, MPFR_RNDN); return mpfr_get_d(r27054, MPFR_RNDN); } static mpfr_t r27055, r27056, r27057, r27058, r27059, r27060, r27061, r27062, r27063, r27064, r27065, r27066, r27067, r27068, r27069, r27070, r27071, r27072, r27073, r27074, r27075, r27076, r27077, r27078, r27079, r27080, r27081, r27082, r27083, r27084, r27085, r27086, r27087, r27088, r27089, r27090, r27091, r27092, r27093, r27094, r27095, r27096, r27097, r27098, r27099, r27100, r27101, r27102, r27103, r27104, r27105, r27106, r27107, r27108, r27109, r27110, r27111, r27112, r27113, r27114, r27115, r27116, r27117, r27118, r27119, r27120, r27121, r27122, r27123, r27124; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27055, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r27056, "8", 10, MPFR_RNDN); mpfr_init_set_str(r27057, "1", 10, MPFR_RNDN); mpfr_init(r27058); mpfr_init(r27059); mpfr_init(r27060); mpfr_init(r27061); mpfr_init(r27062); mpfr_init_set_str(r27063, "771.3234287776531", 10, MPFR_RNDN); mpfr_init(r27064); mpfr_init_set_str(r27065, "3", 10, MPFR_RNDN); mpfr_init(r27066); mpfr_init(r27067); mpfr_init(r27068); mpfr_init_set_str(r27069, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r27070, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r27071); mpfr_init(r27072); mpfr_init_set_str(r27073, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r27074, "2", 10, MPFR_RNDN); mpfr_init(r27075); mpfr_init(r27076); mpfr_init(r27077); mpfr_init(r27078); mpfr_init(r27079); mpfr_init_set_str(r27080, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init_set_str(r27081, "7", 10, MPFR_RNDN); mpfr_init(r27082); mpfr_init(r27083); mpfr_init(r27084); mpfr_init_set_str(r27085, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r27086, "6", 10, MPFR_RNDN); mpfr_init(r27087); mpfr_init(r27088); mpfr_init(r27089); mpfr_init(r27090); mpfr_init_set_str(r27091, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r27092, "4", 10, MPFR_RNDN); mpfr_init(r27093); mpfr_init(r27094); mpfr_init(r27095); mpfr_init_set_str(r27096, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r27097, "5", 10, MPFR_RNDN); mpfr_init(r27098); mpfr_init(r27099); mpfr_init(r27100); mpfr_init(r27101); mpfr_init(r27102); mpfr_init(r27103); mpfr_init(r27104); mpfr_init(r27105); mpfr_init(r27106); mpfr_init(r27107); mpfr_init(r27108); mpfr_init_set_str(r27109, "7.5", 10, MPFR_RNDN); mpfr_init(r27110); mpfr_init_set_str(r27111, "0.5", 10, MPFR_RNDN); mpfr_init(r27112); mpfr_init(r27113); mpfr_init(r27114); mpfr_init(r27115); mpfr_init(r27116); mpfr_init(r27117); mpfr_init(r27118); mpfr_init(r27119); mpfr_init(r27120); mpfr_init(r27121); mpfr_init(r27122); mpfr_init(r27123); mpfr_init(r27124); } double f_fm(double z) { ; ; ; mpfr_add(r27058, r27056, r27057, MPFR_RNDN); mpfr_set_d(r27059, z, MPFR_RNDN); mpfr_add(r27060, r27057, r27059, MPFR_RNDN); mpfr_sub(r27061, r27058, r27060, MPFR_RNDN); mpfr_div(r27062, r27055, r27061, MPFR_RNDN); ; mpfr_sub(r27064, r27057, r27059, MPFR_RNDN); ; mpfr_sub(r27066, r27057, r27065, MPFR_RNDN); mpfr_sub(r27067, r27064, r27066, MPFR_RNDN); mpfr_div(r27068, r27063, r27067, MPFR_RNDN); ; ; mpfr_div(r27071, r27070, r27064, MPFR_RNDN); mpfr_add(r27072, r27069, r27071, MPFR_RNDN); ; ; mpfr_add(r27075, r27057, r27074, MPFR_RNDN); mpfr_sub(r27076, r27075, r27060, MPFR_RNDN); mpfr_div(r27077, r27073, r27076, MPFR_RNDN); mpfr_add(r27078, r27072, r27077, MPFR_RNDN); mpfr_add(r27079, r27068, r27078, MPFR_RNDN); ; ; mpfr_add(r27082, r27081, r27057, MPFR_RNDN); mpfr_sub(r27083, r27082, r27060, MPFR_RNDN); mpfr_div(r27084, r27080, r27083, MPFR_RNDN); ; ; mpfr_add(r27087, r27086, r27057, MPFR_RNDN); mpfr_sub(r27088, r27087, r27060, MPFR_RNDN); mpfr_div(r27089, r27085, r27088, MPFR_RNDN); mpfr_add(r27090, r27084, r27089, MPFR_RNDN); ; ; mpfr_sub(r27093, r27057, r27092, MPFR_RNDN); mpfr_sub(r27094, r27064, r27093, MPFR_RNDN); mpfr_div(r27095, r27091, r27094, MPFR_RNDN); ; ; mpfr_sub(r27098, r27057, r27097, MPFR_RNDN); mpfr_sub(r27099, r27064, r27098, MPFR_RNDN); mpfr_div(r27100, r27096, r27099, MPFR_RNDN); mpfr_add(r27101, r27095, r27100, MPFR_RNDN); mpfr_add(r27102, r27090, r27101, MPFR_RNDN); mpfr_add(r27103, r27079, r27102, MPFR_RNDN); mpfr_add(r27104, r27062, r27103, MPFR_RNDN); mpfr_const_pi(r27105, MPFR_RNDN); mpfr_mul(r27106, r27059, r27105, MPFR_RNDN); mpfr_sin(r27107, r27106, MPFR_RNDN); mpfr_div(r27108, r27105, r27107, MPFR_RNDN); ; mpfr_sub(r27110, r27109, r27059, MPFR_RNDN); ; mpfr_sub(r27112, r27111, r27059, MPFR_RNDN); mpfr_pow(r27113, r27110, r27112, MPFR_RNDN); mpfr_sqrt(r27114, r27074, MPFR_RNDN); mpfr_mul(r27115, r27113, r27114, MPFR_RNDN); mpfr_sub(r27116, r27059, r27109, MPFR_RNDN); mpfr_exp(r27117, r27116, MPFR_RNDN); mpfr_sqrt(r27118, r27105, MPFR_RNDN); mpfr_sqrt(r27119, r27118, MPFR_RNDN); mpfr_mul(r27120, r27119, r27119, MPFR_RNDN); mpfr_mul(r27121, r27117, r27120, MPFR_RNDN); mpfr_mul(r27122, r27115, r27121, MPFR_RNDN); mpfr_mul(r27123, r27108, r27122, MPFR_RNDN); mpfr_mul(r27124, r27104, r27123, MPFR_RNDN); return mpfr_get_d(r27124, MPFR_RNDN); } static mpfr_t r27125, r27126, r27127, r27128, r27129, r27130, r27131, r27132, r27133, r27134, r27135, r27136, r27137, r27138, r27139, r27140, r27141, r27142, r27143, r27144, r27145, r27146, r27147, r27148, r27149, r27150, r27151, r27152, r27153, r27154, r27155, r27156, r27157, r27158, r27159, r27160, r27161, r27162, r27163, r27164, r27165, r27166, r27167, r27168, r27169, r27170, r27171, r27172, r27173, r27174, r27175, r27176, r27177, r27178, r27179, r27180, r27181, r27182, r27183, r27184, r27185, r27186, r27187, r27188, r27189, r27190, r27191, r27192, r27193, r27194; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27125, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r27126, "8", 10, MPFR_RNDN); mpfr_init_set_str(r27127, "1", 10, MPFR_RNDN); mpfr_init(r27128); mpfr_init(r27129); mpfr_init(r27130); mpfr_init(r27131); mpfr_init(r27132); mpfr_init_set_str(r27133, "771.3234287776531", 10, MPFR_RNDN); mpfr_init(r27134); mpfr_init_set_str(r27135, "3", 10, MPFR_RNDN); mpfr_init(r27136); mpfr_init(r27137); mpfr_init(r27138); mpfr_init_set_str(r27139, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r27140, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r27141); mpfr_init(r27142); mpfr_init_set_str(r27143, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r27144, "2", 10, MPFR_RNDN); mpfr_init(r27145); mpfr_init(r27146); mpfr_init(r27147); mpfr_init(r27148); mpfr_init(r27149); mpfr_init_set_str(r27150, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init_set_str(r27151, "7", 10, MPFR_RNDN); mpfr_init(r27152); mpfr_init(r27153); mpfr_init(r27154); mpfr_init_set_str(r27155, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r27156, "6", 10, MPFR_RNDN); mpfr_init(r27157); mpfr_init(r27158); mpfr_init(r27159); mpfr_init(r27160); mpfr_init_set_str(r27161, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r27162, "4", 10, MPFR_RNDN); mpfr_init(r27163); mpfr_init(r27164); mpfr_init(r27165); mpfr_init_set_str(r27166, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r27167, "5", 10, MPFR_RNDN); mpfr_init(r27168); mpfr_init(r27169); mpfr_init(r27170); mpfr_init(r27171); mpfr_init(r27172); mpfr_init(r27173); mpfr_init(r27174); mpfr_init(r27175); mpfr_init(r27176); mpfr_init(r27177); mpfr_init(r27178); mpfr_init_set_str(r27179, "7.5", 10, MPFR_RNDN); mpfr_init(r27180); mpfr_init_set_str(r27181, "0.5", 10, MPFR_RNDN); mpfr_init(r27182); mpfr_init(r27183); mpfr_init(r27184); mpfr_init(r27185); mpfr_init(r27186); mpfr_init(r27187); mpfr_init(r27188); mpfr_init(r27189); mpfr_init(r27190); mpfr_init(r27191); mpfr_init(r27192); mpfr_init(r27193); mpfr_init(r27194); } double f_dm(double z) { ; ; ; mpfr_add(r27128, r27126, r27127, MPFR_RNDN); mpfr_set_d(r27129, z, MPFR_RNDN); mpfr_add(r27130, r27127, r27129, MPFR_RNDN); mpfr_sub(r27131, r27128, r27130, MPFR_RNDN); mpfr_div(r27132, r27125, r27131, MPFR_RNDN); ; mpfr_sub(r27134, r27127, r27129, MPFR_RNDN); ; mpfr_sub(r27136, r27127, r27135, MPFR_RNDN); mpfr_sub(r27137, r27134, r27136, MPFR_RNDN); mpfr_div(r27138, r27133, r27137, MPFR_RNDN); ; ; mpfr_div(r27141, r27140, r27134, MPFR_RNDN); mpfr_add(r27142, r27139, r27141, MPFR_RNDN); ; ; mpfr_add(r27145, r27127, r27144, MPFR_RNDN); mpfr_sub(r27146, r27145, r27130, MPFR_RNDN); mpfr_div(r27147, r27143, r27146, MPFR_RNDN); mpfr_add(r27148, r27142, r27147, MPFR_RNDN); mpfr_add(r27149, r27138, r27148, MPFR_RNDN); ; ; mpfr_add(r27152, r27151, r27127, MPFR_RNDN); mpfr_sub(r27153, r27152, r27130, MPFR_RNDN); mpfr_div(r27154, r27150, r27153, MPFR_RNDN); ; ; mpfr_add(r27157, r27156, r27127, MPFR_RNDN); mpfr_sub(r27158, r27157, r27130, MPFR_RNDN); mpfr_div(r27159, r27155, r27158, MPFR_RNDN); mpfr_add(r27160, r27154, r27159, MPFR_RNDN); ; ; mpfr_sub(r27163, r27127, r27162, MPFR_RNDN); mpfr_sub(r27164, r27134, r27163, MPFR_RNDN); mpfr_div(r27165, r27161, r27164, MPFR_RNDN); ; ; mpfr_sub(r27168, r27127, r27167, MPFR_RNDN); mpfr_sub(r27169, r27134, r27168, MPFR_RNDN); mpfr_div(r27170, r27166, r27169, MPFR_RNDN); mpfr_add(r27171, r27165, r27170, MPFR_RNDN); mpfr_add(r27172, r27160, r27171, MPFR_RNDN); mpfr_add(r27173, r27149, r27172, MPFR_RNDN); mpfr_add(r27174, r27132, r27173, MPFR_RNDN); mpfr_const_pi(r27175, MPFR_RNDN); mpfr_mul(r27176, r27129, r27175, MPFR_RNDN); mpfr_sin(r27177, r27176, MPFR_RNDN); mpfr_div(r27178, r27175, r27177, MPFR_RNDN); ; mpfr_sub(r27180, r27179, r27129, MPFR_RNDN); ; mpfr_sub(r27182, r27181, r27129, MPFR_RNDN); mpfr_pow(r27183, r27180, r27182, MPFR_RNDN); mpfr_sqrt(r27184, r27144, MPFR_RNDN); mpfr_mul(r27185, r27183, r27184, MPFR_RNDN); mpfr_sub(r27186, r27129, r27179, MPFR_RNDN); mpfr_exp(r27187, r27186, MPFR_RNDN); mpfr_sqrt(r27188, r27175, MPFR_RNDN); mpfr_sqrt(r27189, r27188, MPFR_RNDN); mpfr_mul(r27190, r27189, r27189, MPFR_RNDN); mpfr_mul(r27191, r27187, r27190, MPFR_RNDN); mpfr_mul(r27192, r27185, r27191, MPFR_RNDN); mpfr_mul(r27193, r27178, r27192, MPFR_RNDN); mpfr_mul(r27194, r27174, r27193, MPFR_RNDN); return mpfr_get_d(r27194, MPFR_RNDN); }