#include #include #include #include #include char *name = "Jmat.Real.gamma, branch z less than 0.5"; double f_if(float z) { float r26901 = atan2(1.0, 0.0); float r26902 = z; float r26903 = r26901 * r26902; float r26904 = sin(r26903); float r26905 = r26901 / r26904; float r26906 = 2; float r26907 = r26901 * r26906; float r26908 = sqrt(r26907); float r26909 = 1; float r26910 = r26909 - r26902; float r26911 = r26910 - r26909; float r26912 = 7; float r26913 = r26911 + r26912; float r26914 = 0.5; float r26915 = r26913 + r26914; float r26916 = r26911 + r26914; float r26917 = pow(r26915, r26916); float r26918 = r26908 * r26917; float r26919 = -r26915; float r26920 = exp(r26919); float r26921 = r26918 * r26920; float r26922 = 0.9999999999998099; float r26923 = 676.5203681218851; float r26924 = r26911 + r26909; float r26925 = r26923 / r26924; float r26926 = r26922 + r26925; float r26927 = -1259.1392167224028; float r26928 = r26911 + r26906; float r26929 = r26927 / r26928; float r26930 = r26926 + r26929; float r26931 = 771.3234287776531; float r26932 = 3; float r26933 = r26911 + r26932; float r26934 = r26931 / r26933; float r26935 = r26930 + r26934; float r26936 = -176.6150291621406; float r26937 = 4; float r26938 = r26911 + r26937; float r26939 = r26936 / r26938; float r26940 = r26935 + r26939; float r26941 = 12.507343278686905; float r26942 = 5; float r26943 = r26911 + r26942; float r26944 = r26941 / r26943; float r26945 = r26940 + r26944; float r26946 = -0.13857109526572012; float r26947 = 6; float r26948 = r26911 + r26947; float r26949 = r26946 / r26948; float r26950 = r26945 + r26949; float r26951 = 9.984369578019572e-06; float r26952 = r26951 / r26913; float r26953 = r26950 + r26952; float r26954 = 1.5056327351493116e-07; float r26955 = 8; float r26956 = r26911 + r26955; float r26957 = r26954 / r26956; float r26958 = r26953 + r26957; float r26959 = r26921 * r26958; float r26960 = r26905 * r26959; return r26960; } double f_id(double z) { double r26961 = atan2(1.0, 0.0); double r26962 = z; double r26963 = r26961 * r26962; double r26964 = sin(r26963); double r26965 = r26961 / r26964; double r26966 = 2; double r26967 = r26961 * r26966; double r26968 = sqrt(r26967); double r26969 = 1; double r26970 = r26969 - r26962; double r26971 = r26970 - r26969; double r26972 = 7; double r26973 = r26971 + r26972; double r26974 = 0.5; double r26975 = r26973 + r26974; double r26976 = r26971 + r26974; double r26977 = pow(r26975, r26976); double r26978 = r26968 * r26977; double r26979 = -r26975; double r26980 = exp(r26979); double r26981 = r26978 * r26980; double r26982 = 0.9999999999998099; double r26983 = 676.5203681218851; double r26984 = r26971 + r26969; double r26985 = r26983 / r26984; double r26986 = r26982 + r26985; double r26987 = -1259.1392167224028; double r26988 = r26971 + r26966; double r26989 = r26987 / r26988; double r26990 = r26986 + r26989; double r26991 = 771.3234287776531; double r26992 = 3; double r26993 = r26971 + r26992; double r26994 = r26991 / r26993; double r26995 = r26990 + r26994; double r26996 = -176.6150291621406; double r26997 = 4; double r26998 = r26971 + r26997; double r26999 = r26996 / r26998; double r27000 = r26995 + r26999; double r27001 = 12.507343278686905; double r27002 = 5; double r27003 = r26971 + r27002; double r27004 = r27001 / r27003; double r27005 = r27000 + r27004; double r27006 = -0.13857109526572012; double r27007 = 6; double r27008 = r26971 + r27007; double r27009 = r27006 / r27008; double r27010 = r27005 + r27009; double r27011 = 9.984369578019572e-06; double r27012 = r27011 / r26973; double r27013 = r27010 + r27012; double r27014 = 1.5056327351493116e-07; double r27015 = 8; double r27016 = r26971 + r27015; double r27017 = r27014 / r27016; double r27018 = r27013 + r27017; double r27019 = r26981 * r27018; double r27020 = r26965 * r27019; return r27020; } double f_of(float z) { float r27021 = 1.5056327351493116e-07; float r27022 = 8; float r27023 = 1; float r27024 = r27022 + r27023; float r27025 = z; float r27026 = r27023 + r27025; float r27027 = r27024 - r27026; float r27028 = r27021 / r27027; float r27029 = 771.3234287776531; float r27030 = r27023 - r27025; float r27031 = 3; float r27032 = r27023 - r27031; float r27033 = r27030 - r27032; float r27034 = r27029 / r27033; float r27035 = 0.9999999999998099; float r27036 = 676.5203681218851; float r27037 = r27036 / r27030; float r27038 = r27035 + r27037; float r27039 = -1259.1392167224028; float r27040 = 2; float r27041 = -r27025; float r27042 = r27040 + r27041; float r27043 = r27039 / r27042; float r27044 = r27038 + r27043; float r27045 = pow(r27044, r27031); float r27046 = cbrt(r27045); float r27047 = r27034 + r27046; float r27048 = 9.984369578019572e-06; float r27049 = 7; float r27050 = r27049 + r27023; float r27051 = r27050 - r27026; float r27052 = r27048 / r27051; float r27053 = -0.13857109526572012; float r27054 = 6; float r27055 = r27054 + r27023; float r27056 = r27055 - r27026; float r27057 = r27053 / r27056; float r27058 = r27052 + r27057; float r27059 = -176.6150291621406; float r27060 = 4; float r27061 = r27023 - r27060; float r27062 = r27030 - r27061; float r27063 = r27059 / r27062; float r27064 = 12.507343278686905; float r27065 = 5; float r27066 = r27023 - r27065; float r27067 = r27030 - r27066; float r27068 = r27064 / r27067; float r27069 = r27063 + r27068; float r27070 = r27058 + r27069; float r27071 = r27047 + r27070; float r27072 = r27028 + r27071; float r27073 = atan2(1.0, 0.0); float r27074 = r27025 * r27073; float r27075 = sin(r27074); float r27076 = r27073 / r27075; float r27077 = 7.5; float r27078 = r27077 - r27025; float r27079 = 0.5; float r27080 = r27079 - r27025; float r27081 = pow(r27078, r27080); float r27082 = sqrt(r27040); float r27083 = r27081 * r27082; float r27084 = r27025 - r27077; float r27085 = exp(r27084); float r27086 = sqrt(r27073); float r27087 = r27085 * r27086; float r27088 = r27083 * r27087; float r27089 = r27076 * r27088; float r27090 = r27072 * r27089; return r27090; } double f_od(double z) { double r27091 = 1.5056327351493116e-07; double r27092 = 8; double r27093 = 1; double r27094 = r27092 + r27093; double r27095 = z; double r27096 = r27093 + r27095; double r27097 = r27094 - r27096; double r27098 = r27091 / r27097; double r27099 = 771.3234287776531; double r27100 = r27093 - r27095; double r27101 = 3; double r27102 = r27093 - r27101; double r27103 = r27100 - r27102; double r27104 = r27099 / r27103; double r27105 = 0.9999999999998099; double r27106 = 676.5203681218851; double r27107 = r27106 / r27100; double r27108 = r27105 + r27107; double r27109 = -1259.1392167224028; double r27110 = 2; double r27111 = -r27095; double r27112 = r27110 + r27111; double r27113 = r27109 / r27112; double r27114 = r27108 + r27113; double r27115 = pow(r27114, r27101); double r27116 = cbrt(r27115); double r27117 = r27104 + r27116; double r27118 = 9.984369578019572e-06; double r27119 = 7; double r27120 = r27119 + r27093; double r27121 = r27120 - r27096; double r27122 = r27118 / r27121; double r27123 = -0.13857109526572012; double r27124 = 6; double r27125 = r27124 + r27093; double r27126 = r27125 - r27096; double r27127 = r27123 / r27126; double r27128 = r27122 + r27127; double r27129 = -176.6150291621406; double r27130 = 4; double r27131 = r27093 - r27130; double r27132 = r27100 - r27131; double r27133 = r27129 / r27132; double r27134 = 12.507343278686905; double r27135 = 5; double r27136 = r27093 - r27135; double r27137 = r27100 - r27136; double r27138 = r27134 / r27137; double r27139 = r27133 + r27138; double r27140 = r27128 + r27139; double r27141 = r27117 + r27140; double r27142 = r27098 + r27141; double r27143 = atan2(1.0, 0.0); double r27144 = r27095 * r27143; double r27145 = sin(r27144); double r27146 = r27143 / r27145; double r27147 = 7.5; double r27148 = r27147 - r27095; double r27149 = 0.5; double r27150 = r27149 - r27095; double r27151 = pow(r27148, r27150); double r27152 = sqrt(r27110); double r27153 = r27151 * r27152; double r27154 = r27095 - r27147; double r27155 = exp(r27154); double r27156 = sqrt(r27143); double r27157 = r27155 * r27156; double r27158 = r27153 * r27157; double r27159 = r27146 * r27158; double r27160 = r27142 * r27159; return r27160; } 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 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, r27195, r27196, r27197, r27198, r27199, r27200, r27201, r27202, r27203, r27204, r27205, r27206, r27207, r27208, r27209, r27210, r27211, r27212, r27213, r27214, r27215, r27216, r27217, r27218, r27219, r27220; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r27161); mpfr_init(r27162); mpfr_init(r27163); mpfr_init(r27164); mpfr_init(r27165); mpfr_init_set_str(r27166, "2", 10, MPFR_RNDN); mpfr_init(r27167); mpfr_init(r27168); mpfr_init_set_str(r27169, "1", 10, MPFR_RNDN); mpfr_init(r27170); mpfr_init(r27171); mpfr_init_set_str(r27172, "7", 10, MPFR_RNDN); mpfr_init(r27173); mpfr_init_set_str(r27174, "0.5", 10, MPFR_RNDN); mpfr_init(r27175); mpfr_init(r27176); mpfr_init(r27177); mpfr_init(r27178); mpfr_init(r27179); mpfr_init(r27180); mpfr_init(r27181); mpfr_init_set_str(r27182, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r27183, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r27184); mpfr_init(r27185); mpfr_init(r27186); mpfr_init_set_str(r27187, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init(r27188); mpfr_init(r27189); mpfr_init(r27190); mpfr_init_set_str(r27191, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r27192, "3", 10, MPFR_RNDN); mpfr_init(r27193); mpfr_init(r27194); mpfr_init(r27195); mpfr_init_set_str(r27196, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r27197, "4", 10, MPFR_RNDN); mpfr_init(r27198); mpfr_init(r27199); mpfr_init(r27200); mpfr_init_set_str(r27201, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r27202, "5", 10, MPFR_RNDN); mpfr_init(r27203); mpfr_init(r27204); mpfr_init(r27205); mpfr_init_set_str(r27206, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r27207, "6", 10, MPFR_RNDN); mpfr_init(r27208); mpfr_init(r27209); mpfr_init(r27210); mpfr_init_set_str(r27211, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init(r27212); mpfr_init(r27213); mpfr_init_set_str(r27214, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r27215, "8", 10, MPFR_RNDN); mpfr_init(r27216); mpfr_init(r27217); mpfr_init(r27218); mpfr_init(r27219); mpfr_init(r27220); } double f_im(double z) { mpfr_const_pi(r27161, MPFR_RNDN); mpfr_set_d(r27162, z, MPFR_RNDN); mpfr_mul(r27163, r27161, r27162, MPFR_RNDN); mpfr_sin(r27164, r27163, MPFR_RNDN); mpfr_div(r27165, r27161, r27164, MPFR_RNDN); ; mpfr_mul(r27167, r27161, r27166, MPFR_RNDN); mpfr_sqrt(r27168, r27167, MPFR_RNDN); ; mpfr_sub(r27170, r27169, r27162, MPFR_RNDN); mpfr_sub(r27171, r27170, r27169, MPFR_RNDN); ; mpfr_add(r27173, r27171, r27172, MPFR_RNDN); ; mpfr_add(r27175, r27173, r27174, MPFR_RNDN); mpfr_add(r27176, r27171, r27174, MPFR_RNDN); mpfr_pow(r27177, r27175, r27176, MPFR_RNDN); mpfr_mul(r27178, r27168, r27177, MPFR_RNDN); mpfr_neg(r27179, r27175, MPFR_RNDN); mpfr_exp(r27180, r27179, MPFR_RNDN); mpfr_mul(r27181, r27178, r27180, MPFR_RNDN); ; ; mpfr_add(r27184, r27171, r27169, MPFR_RNDN); mpfr_div(r27185, r27183, r27184, MPFR_RNDN); mpfr_add(r27186, r27182, r27185, MPFR_RNDN); ; mpfr_add(r27188, r27171, r27166, MPFR_RNDN); mpfr_div(r27189, r27187, r27188, MPFR_RNDN); mpfr_add(r27190, r27186, r27189, MPFR_RNDN); ; ; mpfr_add(r27193, r27171, r27192, MPFR_RNDN); mpfr_div(r27194, r27191, r27193, MPFR_RNDN); mpfr_add(r27195, r27190, r27194, MPFR_RNDN); ; ; mpfr_add(r27198, r27171, r27197, MPFR_RNDN); mpfr_div(r27199, r27196, r27198, MPFR_RNDN); mpfr_add(r27200, r27195, r27199, MPFR_RNDN); ; ; mpfr_add(r27203, r27171, r27202, MPFR_RNDN); mpfr_div(r27204, r27201, r27203, MPFR_RNDN); mpfr_add(r27205, r27200, r27204, MPFR_RNDN); ; ; mpfr_add(r27208, r27171, r27207, MPFR_RNDN); mpfr_div(r27209, r27206, r27208, MPFR_RNDN); mpfr_add(r27210, r27205, r27209, MPFR_RNDN); ; mpfr_div(r27212, r27211, r27173, MPFR_RNDN); mpfr_add(r27213, r27210, r27212, MPFR_RNDN); ; ; mpfr_add(r27216, r27171, r27215, MPFR_RNDN); mpfr_div(r27217, r27214, r27216, MPFR_RNDN); mpfr_add(r27218, r27213, r27217, MPFR_RNDN); mpfr_mul(r27219, r27181, r27218, MPFR_RNDN); mpfr_mul(r27220, r27165, r27219, MPFR_RNDN); return mpfr_get_d(r27220, MPFR_RNDN); } static mpfr_t r27221, r27222, r27223, r27224, r27225, r27226, r27227, r27228, r27229, r27230, r27231, r27232, r27233, r27234, r27235, r27236, r27237, r27238, r27239, r27240, r27241, r27242, r27243, r27244, r27245, r27246, r27247, r27248, r27249, r27250, r27251, r27252, r27253, r27254, r27255, r27256, r27257, r27258, r27259, r27260, r27261, r27262, r27263, r27264, r27265, r27266, r27267, r27268, r27269, r27270, r27271, r27272, r27273, r27274, r27275, r27276, r27277, r27278, r27279, r27280, r27281, r27282, r27283, r27284, r27285, r27286, r27287, r27288, r27289, r27290; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27221, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r27222, "8", 10, MPFR_RNDN); mpfr_init_set_str(r27223, "1", 10, MPFR_RNDN); mpfr_init(r27224); mpfr_init(r27225); mpfr_init(r27226); mpfr_init(r27227); mpfr_init(r27228); mpfr_init_set_str(r27229, "771.3234287776531", 10, MPFR_RNDN); mpfr_init(r27230); mpfr_init_set_str(r27231, "3", 10, MPFR_RNDN); mpfr_init(r27232); mpfr_init(r27233); mpfr_init(r27234); mpfr_init_set_str(r27235, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r27236, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r27237); mpfr_init(r27238); mpfr_init_set_str(r27239, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r27240, "2", 10, MPFR_RNDN); mpfr_init(r27241); mpfr_init(r27242); mpfr_init(r27243); mpfr_init(r27244); mpfr_init(r27245); mpfr_init(r27246); mpfr_init(r27247); mpfr_init_set_str(r27248, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init_set_str(r27249, "7", 10, MPFR_RNDN); mpfr_init(r27250); mpfr_init(r27251); mpfr_init(r27252); mpfr_init_set_str(r27253, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r27254, "6", 10, MPFR_RNDN); mpfr_init(r27255); mpfr_init(r27256); mpfr_init(r27257); mpfr_init(r27258); mpfr_init_set_str(r27259, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r27260, "4", 10, MPFR_RNDN); mpfr_init(r27261); mpfr_init(r27262); mpfr_init(r27263); mpfr_init_set_str(r27264, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r27265, "5", 10, MPFR_RNDN); mpfr_init(r27266); mpfr_init(r27267); mpfr_init(r27268); mpfr_init(r27269); mpfr_init(r27270); mpfr_init(r27271); mpfr_init(r27272); mpfr_init(r27273); mpfr_init(r27274); mpfr_init(r27275); mpfr_init(r27276); mpfr_init_set_str(r27277, "7.5", 10, MPFR_RNDN); mpfr_init(r27278); mpfr_init_set_str(r27279, "0.5", 10, MPFR_RNDN); mpfr_init(r27280); mpfr_init(r27281); mpfr_init(r27282); mpfr_init(r27283); mpfr_init(r27284); mpfr_init(r27285); mpfr_init(r27286); mpfr_init(r27287); mpfr_init(r27288); mpfr_init(r27289); mpfr_init(r27290); } double f_fm(double z) { ; ; ; mpfr_add(r27224, r27222, r27223, MPFR_RNDN); mpfr_set_d(r27225, z, MPFR_RNDN); mpfr_add(r27226, r27223, r27225, MPFR_RNDN); mpfr_sub(r27227, r27224, r27226, MPFR_RNDN); mpfr_div(r27228, r27221, r27227, MPFR_RNDN); ; mpfr_sub(r27230, r27223, r27225, MPFR_RNDN); ; mpfr_sub(r27232, r27223, r27231, MPFR_RNDN); mpfr_sub(r27233, r27230, r27232, MPFR_RNDN); mpfr_div(r27234, r27229, r27233, MPFR_RNDN); ; ; mpfr_div(r27237, r27236, r27230, MPFR_RNDN); mpfr_add(r27238, r27235, r27237, MPFR_RNDN); ; ; mpfr_neg(r27241, r27225, MPFR_RNDN); mpfr_add(r27242, r27240, r27241, MPFR_RNDN); mpfr_div(r27243, r27239, r27242, MPFR_RNDN); mpfr_add(r27244, r27238, r27243, MPFR_RNDN); mpfr_pow(r27245, r27244, r27231, MPFR_RNDN); mpfr_cbrt(r27246, r27245, MPFR_RNDN); mpfr_add(r27247, r27234, r27246, MPFR_RNDN); ; ; mpfr_add(r27250, r27249, r27223, MPFR_RNDN); mpfr_sub(r27251, r27250, r27226, MPFR_RNDN); mpfr_div(r27252, r27248, r27251, MPFR_RNDN); ; ; mpfr_add(r27255, r27254, r27223, MPFR_RNDN); mpfr_sub(r27256, r27255, r27226, MPFR_RNDN); mpfr_div(r27257, r27253, r27256, MPFR_RNDN); mpfr_add(r27258, r27252, r27257, MPFR_RNDN); ; ; mpfr_sub(r27261, r27223, r27260, MPFR_RNDN); mpfr_sub(r27262, r27230, r27261, MPFR_RNDN); mpfr_div(r27263, r27259, r27262, MPFR_RNDN); ; ; mpfr_sub(r27266, r27223, r27265, MPFR_RNDN); mpfr_sub(r27267, r27230, r27266, MPFR_RNDN); mpfr_div(r27268, r27264, r27267, MPFR_RNDN); mpfr_add(r27269, r27263, r27268, MPFR_RNDN); mpfr_add(r27270, r27258, r27269, MPFR_RNDN); mpfr_add(r27271, r27247, r27270, MPFR_RNDN); mpfr_add(r27272, r27228, r27271, MPFR_RNDN); mpfr_const_pi(r27273, MPFR_RNDN); mpfr_mul(r27274, r27225, r27273, MPFR_RNDN); mpfr_sin(r27275, r27274, MPFR_RNDN); mpfr_div(r27276, r27273, r27275, MPFR_RNDN); ; mpfr_sub(r27278, r27277, r27225, MPFR_RNDN); ; mpfr_sub(r27280, r27279, r27225, MPFR_RNDN); mpfr_pow(r27281, r27278, r27280, MPFR_RNDN); mpfr_sqrt(r27282, r27240, MPFR_RNDN); mpfr_mul(r27283, r27281, r27282, MPFR_RNDN); mpfr_sub(r27284, r27225, r27277, MPFR_RNDN); mpfr_exp(r27285, r27284, MPFR_RNDN); mpfr_sqrt(r27286, r27273, MPFR_RNDN); mpfr_mul(r27287, r27285, r27286, MPFR_RNDN); mpfr_mul(r27288, r27283, r27287, MPFR_RNDN); mpfr_mul(r27289, r27276, r27288, MPFR_RNDN); mpfr_mul(r27290, r27272, r27289, MPFR_RNDN); return mpfr_get_d(r27290, MPFR_RNDN); } static mpfr_t r27291, r27292, r27293, r27294, r27295, r27296, r27297, r27298, r27299, r27300, r27301, r27302, r27303, r27304, r27305, r27306, r27307, r27308, r27309, r27310, r27311, r27312, r27313, r27314, r27315, r27316, r27317, r27318, r27319, r27320, r27321, r27322, r27323, r27324, r27325, r27326, r27327, r27328, r27329, r27330, r27331, r27332, r27333, r27334, r27335, r27336, r27337, r27338, r27339, r27340, r27341, r27342, r27343, r27344, r27345, r27346, r27347, r27348, r27349, r27350, r27351, r27352, r27353, r27354, r27355, r27356, r27357, r27358, r27359, r27360; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27291, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r27292, "8", 10, MPFR_RNDN); mpfr_init_set_str(r27293, "1", 10, MPFR_RNDN); mpfr_init(r27294); mpfr_init(r27295); mpfr_init(r27296); mpfr_init(r27297); mpfr_init(r27298); mpfr_init_set_str(r27299, "771.3234287776531", 10, MPFR_RNDN); mpfr_init(r27300); mpfr_init_set_str(r27301, "3", 10, MPFR_RNDN); mpfr_init(r27302); mpfr_init(r27303); mpfr_init(r27304); mpfr_init_set_str(r27305, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r27306, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r27307); mpfr_init(r27308); mpfr_init_set_str(r27309, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r27310, "2", 10, MPFR_RNDN); mpfr_init(r27311); mpfr_init(r27312); mpfr_init(r27313); mpfr_init(r27314); mpfr_init(r27315); mpfr_init(r27316); mpfr_init(r27317); mpfr_init_set_str(r27318, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init_set_str(r27319, "7", 10, MPFR_RNDN); mpfr_init(r27320); mpfr_init(r27321); mpfr_init(r27322); mpfr_init_set_str(r27323, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r27324, "6", 10, MPFR_RNDN); mpfr_init(r27325); mpfr_init(r27326); mpfr_init(r27327); mpfr_init(r27328); mpfr_init_set_str(r27329, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r27330, "4", 10, MPFR_RNDN); mpfr_init(r27331); mpfr_init(r27332); mpfr_init(r27333); mpfr_init_set_str(r27334, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r27335, "5", 10, MPFR_RNDN); mpfr_init(r27336); mpfr_init(r27337); mpfr_init(r27338); mpfr_init(r27339); mpfr_init(r27340); mpfr_init(r27341); mpfr_init(r27342); mpfr_init(r27343); mpfr_init(r27344); mpfr_init(r27345); mpfr_init(r27346); mpfr_init_set_str(r27347, "7.5", 10, MPFR_RNDN); mpfr_init(r27348); mpfr_init_set_str(r27349, "0.5", 10, MPFR_RNDN); mpfr_init(r27350); mpfr_init(r27351); mpfr_init(r27352); mpfr_init(r27353); mpfr_init(r27354); mpfr_init(r27355); mpfr_init(r27356); mpfr_init(r27357); mpfr_init(r27358); mpfr_init(r27359); mpfr_init(r27360); } double f_dm(double z) { ; ; ; mpfr_add(r27294, r27292, r27293, MPFR_RNDN); mpfr_set_d(r27295, z, MPFR_RNDN); mpfr_add(r27296, r27293, r27295, MPFR_RNDN); mpfr_sub(r27297, r27294, r27296, MPFR_RNDN); mpfr_div(r27298, r27291, r27297, MPFR_RNDN); ; mpfr_sub(r27300, r27293, r27295, MPFR_RNDN); ; mpfr_sub(r27302, r27293, r27301, MPFR_RNDN); mpfr_sub(r27303, r27300, r27302, MPFR_RNDN); mpfr_div(r27304, r27299, r27303, MPFR_RNDN); ; ; mpfr_div(r27307, r27306, r27300, MPFR_RNDN); mpfr_add(r27308, r27305, r27307, MPFR_RNDN); ; ; mpfr_neg(r27311, r27295, MPFR_RNDN); mpfr_add(r27312, r27310, r27311, MPFR_RNDN); mpfr_div(r27313, r27309, r27312, MPFR_RNDN); mpfr_add(r27314, r27308, r27313, MPFR_RNDN); mpfr_pow(r27315, r27314, r27301, MPFR_RNDN); mpfr_cbrt(r27316, r27315, MPFR_RNDN); mpfr_add(r27317, r27304, r27316, MPFR_RNDN); ; ; mpfr_add(r27320, r27319, r27293, MPFR_RNDN); mpfr_sub(r27321, r27320, r27296, MPFR_RNDN); mpfr_div(r27322, r27318, r27321, MPFR_RNDN); ; ; mpfr_add(r27325, r27324, r27293, MPFR_RNDN); mpfr_sub(r27326, r27325, r27296, MPFR_RNDN); mpfr_div(r27327, r27323, r27326, MPFR_RNDN); mpfr_add(r27328, r27322, r27327, MPFR_RNDN); ; ; mpfr_sub(r27331, r27293, r27330, MPFR_RNDN); mpfr_sub(r27332, r27300, r27331, MPFR_RNDN); mpfr_div(r27333, r27329, r27332, MPFR_RNDN); ; ; mpfr_sub(r27336, r27293, r27335, MPFR_RNDN); mpfr_sub(r27337, r27300, r27336, MPFR_RNDN); mpfr_div(r27338, r27334, r27337, MPFR_RNDN); mpfr_add(r27339, r27333, r27338, MPFR_RNDN); mpfr_add(r27340, r27328, r27339, MPFR_RNDN); mpfr_add(r27341, r27317, r27340, MPFR_RNDN); mpfr_add(r27342, r27298, r27341, MPFR_RNDN); mpfr_const_pi(r27343, MPFR_RNDN); mpfr_mul(r27344, r27295, r27343, MPFR_RNDN); mpfr_sin(r27345, r27344, MPFR_RNDN); mpfr_div(r27346, r27343, r27345, MPFR_RNDN); ; mpfr_sub(r27348, r27347, r27295, MPFR_RNDN); ; mpfr_sub(r27350, r27349, r27295, MPFR_RNDN); mpfr_pow(r27351, r27348, r27350, MPFR_RNDN); mpfr_sqrt(r27352, r27310, MPFR_RNDN); mpfr_mul(r27353, r27351, r27352, MPFR_RNDN); mpfr_sub(r27354, r27295, r27347, MPFR_RNDN); mpfr_exp(r27355, r27354, MPFR_RNDN); mpfr_sqrt(r27356, r27343, MPFR_RNDN); mpfr_mul(r27357, r27355, r27356, MPFR_RNDN); mpfr_mul(r27358, r27353, r27357, MPFR_RNDN); mpfr_mul(r27359, r27346, r27358, MPFR_RNDN); mpfr_mul(r27360, r27342, r27359, MPFR_RNDN); return mpfr_get_d(r27360, MPFR_RNDN); }