#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r26987 = 1; float r26988 = 0.1049934947; float r26989 = x; float r26990 = r26989 * r26989; float r26991 = r26988 * r26990; float r26992 = r26987 + r26991; float r26993 = 0.0424060604; float r26994 = r26990 * r26990; float r26995 = r26993 * r26994; float r26996 = r26992 + r26995; float r26997 = 0.0072644182; float r26998 = r26994 * r26990; float r26999 = r26997 * r26998; float r27000 = r26996 + r26999; float r27001 = 0.0005064034; float r27002 = r26998 * r26990; float r27003 = r27001 * r27002; float r27004 = r27000 + r27003; float r27005 = 0.0001789971; float r27006 = r27002 * r26990; float r27007 = r27005 * r27006; float r27008 = r27004 + r27007; float r27009 = 0.7715471019; float r27010 = r27009 * r26990; float r27011 = r26987 + r27010; float r27012 = 0.2909738639; float r27013 = r27012 * r26994; float r27014 = r27011 + r27013; float r27015 = 0.0694555761; float r27016 = r27015 * r26998; float r27017 = r27014 + r27016; float r27018 = 0.0140005442; float r27019 = r27018 * r27002; float r27020 = r27017 + r27019; float r27021 = 0.0008327945; float r27022 = r27021 * r27006; float r27023 = r27020 + r27022; float r27024 = 2; float r27025 = r27024 * r27005; float r27026 = r27006 * r26990; float r27027 = r27025 * r27026; float r27028 = r27023 + r27027; float r27029 = r27008 / r27028; float r27030 = r27029 * r26989; return r27030; } double f_id(double x) { double r27031 = 1; double r27032 = 0.1049934947; double r27033 = x; double r27034 = r27033 * r27033; double r27035 = r27032 * r27034; double r27036 = r27031 + r27035; double r27037 = 0.0424060604; double r27038 = r27034 * r27034; double r27039 = r27037 * r27038; double r27040 = r27036 + r27039; double r27041 = 0.0072644182; double r27042 = r27038 * r27034; double r27043 = r27041 * r27042; double r27044 = r27040 + r27043; double r27045 = 0.0005064034; double r27046 = r27042 * r27034; double r27047 = r27045 * r27046; double r27048 = r27044 + r27047; double r27049 = 0.0001789971; double r27050 = r27046 * r27034; double r27051 = r27049 * r27050; double r27052 = r27048 + r27051; double r27053 = 0.7715471019; double r27054 = r27053 * r27034; double r27055 = r27031 + r27054; double r27056 = 0.2909738639; double r27057 = r27056 * r27038; double r27058 = r27055 + r27057; double r27059 = 0.0694555761; double r27060 = r27059 * r27042; double r27061 = r27058 + r27060; double r27062 = 0.0140005442; double r27063 = r27062 * r27046; double r27064 = r27061 + r27063; double r27065 = 0.0008327945; double r27066 = r27065 * r27050; double r27067 = r27064 + r27066; double r27068 = 2; double r27069 = r27068 * r27049; double r27070 = r27050 * r27034; double r27071 = r27069 * r27070; double r27072 = r27067 + r27071; double r27073 = r27052 / r27072; double r27074 = r27073 * r27033; return r27074; } double f_of(float x) { float r27075 = x; float r27076 = -512314.0678876898; bool r27077 = r27075 <= r27076; float r27078 = 723.4643878500376; bool r27079 = r27075 <= r27078; float r27080 = !r27079; bool r27081 = r27077 || r27080; float r27082 = 0.2514179000665375; float r27083 = 4; float r27084 = pow(r27075, r27083); float r27085 = r27082 / r27084; float r27086 = 0.15298196345929327; float r27087 = 6; float r27088 = pow(r27075, r27087); float r27089 = r27086 / r27088; float r27090 = r27085 + r27089; float r27091 = 0.5; float r27092 = r27091 / r27075; float r27093 = fma(r27090, r27075, r27092); float r27094 = r27075 * r27075; float r27095 = r27094 * r27094; float r27096 = 0.0072644182; float r27097 = r27075 * r27096; float r27098 = r27097 * r27075; float r27099 = 0.0424060604; float r27100 = 0.1049934947; float r27101 = 1; float r27102 = fma(r27094, r27100, r27101); float r27103 = fma(r27099, r27095, r27102); float r27104 = fma(r27095, r27098, r27103); float r27105 = 0.0001789971; float r27106 = r27105 * r27094; float r27107 = 0.0005064034; float r27108 = r27106 + r27107; float r27109 = 3; float r27110 = r27109 + r27101; float r27111 = pow(r27094, r27110); float r27112 = r27108 * r27111; float r27113 = r27104 + r27112; float r27114 = sqrt(r27113); float r27115 = 0.0694555761; float r27116 = r27115 * r27075; float r27117 = r27116 * r27095; float r27118 = 0.2909738639; float r27119 = 0.7715471019; float r27120 = r27075 * r27119; float r27121 = fma(r27075, r27120, r27101); float r27122 = fma(r27095, r27118, r27121); float r27123 = fma(r27075, r27117, r27122); float r27124 = 2; float r27125 = r27124 * r27105; float r27126 = 0.0008327945; float r27127 = r27126 * r27075; float r27128 = r27127 * r27075; float r27129 = fma(r27095, r27125, r27128); float r27130 = 0.0140005442; float r27131 = r27130 * r27111; float r27132 = fma(r27111, r27129, r27131); float r27133 = r27123 + r27132; float r27134 = r27114 / r27133; float r27135 = r27094 * r27096; float r27136 = fma(r27095, r27135, r27103); float r27137 = r27112 + r27136; float r27138 = sqrt(r27137); float r27139 = r27134 * r27138; float r27140 = r27075 * r27139; float r27141 = r27081 ? r27093 : r27140; return r27141; } double f_od(double x) { double r27142 = x; double r27143 = -512314.0678876898; bool r27144 = r27142 <= r27143; double r27145 = 723.4643878500376; bool r27146 = r27142 <= r27145; double r27147 = !r27146; bool r27148 = r27144 || r27147; double r27149 = 0.2514179000665375; double r27150 = 4; double r27151 = pow(r27142, r27150); double r27152 = r27149 / r27151; double r27153 = 0.15298196345929327; double r27154 = 6; double r27155 = pow(r27142, r27154); double r27156 = r27153 / r27155; double r27157 = r27152 + r27156; double r27158 = 0.5; double r27159 = r27158 / r27142; double r27160 = fma(r27157, r27142, r27159); double r27161 = r27142 * r27142; double r27162 = r27161 * r27161; double r27163 = 0.0072644182; double r27164 = r27142 * r27163; double r27165 = r27164 * r27142; double r27166 = 0.0424060604; double r27167 = 0.1049934947; double r27168 = 1; double r27169 = fma(r27161, r27167, r27168); double r27170 = fma(r27166, r27162, r27169); double r27171 = fma(r27162, r27165, r27170); double r27172 = 0.0001789971; double r27173 = r27172 * r27161; double r27174 = 0.0005064034; double r27175 = r27173 + r27174; double r27176 = 3; double r27177 = r27176 + r27168; double r27178 = pow(r27161, r27177); double r27179 = r27175 * r27178; double r27180 = r27171 + r27179; double r27181 = sqrt(r27180); double r27182 = 0.0694555761; double r27183 = r27182 * r27142; double r27184 = r27183 * r27162; double r27185 = 0.2909738639; double r27186 = 0.7715471019; double r27187 = r27142 * r27186; double r27188 = fma(r27142, r27187, r27168); double r27189 = fma(r27162, r27185, r27188); double r27190 = fma(r27142, r27184, r27189); double r27191 = 2; double r27192 = r27191 * r27172; double r27193 = 0.0008327945; double r27194 = r27193 * r27142; double r27195 = r27194 * r27142; double r27196 = fma(r27162, r27192, r27195); double r27197 = 0.0140005442; double r27198 = r27197 * r27178; double r27199 = fma(r27178, r27196, r27198); double r27200 = r27190 + r27199; double r27201 = r27181 / r27200; double r27202 = r27161 * r27163; double r27203 = fma(r27162, r27202, r27170); double r27204 = r27179 + r27203; double r27205 = sqrt(r27204); double r27206 = r27201 * r27205; double r27207 = r27142 * r27206; double r27208 = r27148 ? r27160 : r27207; return r27208; } 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 r27209, r27210, r27211, r27212, r27213, r27214, r27215, r27216, r27217, r27218, r27219, r27220, 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; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r27209, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27210, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r27211); mpfr_init(r27212); mpfr_init(r27213); mpfr_init(r27214); mpfr_init_set_str(r27215, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r27216); mpfr_init(r27217); mpfr_init(r27218); mpfr_init_set_str(r27219, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r27220); mpfr_init(r27221); mpfr_init(r27222); mpfr_init_set_str(r27223, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r27224); mpfr_init(r27225); mpfr_init(r27226); mpfr_init_set_str(r27227, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r27228); mpfr_init(r27229); mpfr_init(r27230); mpfr_init_set_str(r27231, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r27232); mpfr_init(r27233); mpfr_init_set_str(r27234, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r27235); mpfr_init(r27236); mpfr_init_set_str(r27237, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r27238); mpfr_init(r27239); mpfr_init_set_str(r27240, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r27241); mpfr_init(r27242); mpfr_init_set_str(r27243, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r27244); mpfr_init(r27245); mpfr_init_set_str(r27246, "2", 10, MPFR_RNDN); mpfr_init(r27247); mpfr_init(r27248); mpfr_init(r27249); mpfr_init(r27250); mpfr_init(r27251); mpfr_init(r27252); } double f_im(double x) { ; ; mpfr_set_d(r27211, x, MPFR_RNDN); mpfr_mul(r27212, r27211, r27211, MPFR_RNDN); mpfr_mul(r27213, r27210, r27212, MPFR_RNDN); mpfr_add(r27214, r27209, r27213, MPFR_RNDN); ; mpfr_mul(r27216, r27212, r27212, MPFR_RNDN); mpfr_mul(r27217, r27215, r27216, MPFR_RNDN); mpfr_add(r27218, r27214, r27217, MPFR_RNDN); ; mpfr_mul(r27220, r27216, r27212, MPFR_RNDN); mpfr_mul(r27221, r27219, r27220, MPFR_RNDN); mpfr_add(r27222, r27218, r27221, MPFR_RNDN); ; mpfr_mul(r27224, r27220, r27212, MPFR_RNDN); mpfr_mul(r27225, r27223, r27224, MPFR_RNDN); mpfr_add(r27226, r27222, r27225, MPFR_RNDN); ; mpfr_mul(r27228, r27224, r27212, MPFR_RNDN); mpfr_mul(r27229, r27227, r27228, MPFR_RNDN); mpfr_add(r27230, r27226, r27229, MPFR_RNDN); ; mpfr_mul(r27232, r27231, r27212, MPFR_RNDN); mpfr_add(r27233, r27209, r27232, MPFR_RNDN); ; mpfr_mul(r27235, r27234, r27216, MPFR_RNDN); mpfr_add(r27236, r27233, r27235, MPFR_RNDN); ; mpfr_mul(r27238, r27237, r27220, MPFR_RNDN); mpfr_add(r27239, r27236, r27238, MPFR_RNDN); ; mpfr_mul(r27241, r27240, r27224, MPFR_RNDN); mpfr_add(r27242, r27239, r27241, MPFR_RNDN); ; mpfr_mul(r27244, r27243, r27228, MPFR_RNDN); mpfr_add(r27245, r27242, r27244, MPFR_RNDN); ; mpfr_mul(r27247, r27246, r27227, MPFR_RNDN); mpfr_mul(r27248, r27228, r27212, MPFR_RNDN); mpfr_mul(r27249, r27247, r27248, MPFR_RNDN); mpfr_add(r27250, r27245, r27249, MPFR_RNDN); mpfr_div(r27251, r27230, r27250, MPFR_RNDN); mpfr_mul(r27252, r27251, r27211, MPFR_RNDN); return mpfr_get_d(r27252, MPFR_RNDN); } static mpfr_t 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, 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r27253); mpfr_init_set_str(r27254, "-512314.0678876898", 10, MPFR_RNDN); mpfr_init(r27255); mpfr_init_set_str(r27256, "723.4643878500376", 10, MPFR_RNDN); mpfr_init(r27257); mpfr_init(r27258); mpfr_init(r27259); mpfr_init_set_str(r27260, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init_set_str(r27261, "4", 10, MPFR_RNDN); mpfr_init(r27262); mpfr_init(r27263); mpfr_init_set_str(r27264, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r27265, "6", 10, MPFR_RNDN); mpfr_init(r27266); mpfr_init(r27267); mpfr_init(r27268); mpfr_init_set_str(r27269, "0.5", 10, MPFR_RNDN); mpfr_init(r27270); mpfr_init(r27271); mpfr_init(r27272); mpfr_init(r27273); mpfr_init_set_str(r27274, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r27275); mpfr_init(r27276); mpfr_init_set_str(r27277, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r27278, "0.1049934947", 10, MPFR_RNDN); mpfr_init_set_str(r27279, "1", 10, MPFR_RNDN); mpfr_init(r27280); mpfr_init(r27281); mpfr_init(r27282); mpfr_init_set_str(r27283, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r27284); mpfr_init_set_str(r27285, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r27286); mpfr_init_set_str(r27287, "3", 10, MPFR_RNDN); mpfr_init(r27288); mpfr_init(r27289); mpfr_init(r27290); mpfr_init(r27291); mpfr_init(r27292); mpfr_init_set_str(r27293, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r27294); mpfr_init(r27295); mpfr_init_set_str(r27296, "0.2909738639", 10, MPFR_RNDN); mpfr_init_set_str(r27297, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r27298); mpfr_init(r27299); mpfr_init(r27300); mpfr_init(r27301); mpfr_init_set_str(r27302, "2", 10, MPFR_RNDN); mpfr_init(r27303); mpfr_init_set_str(r27304, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r27305); mpfr_init(r27306); mpfr_init(r27307); mpfr_init_set_str(r27308, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r27309); mpfr_init(r27310); mpfr_init(r27311); mpfr_init(r27312); mpfr_init(r27313); mpfr_init(r27314); mpfr_init(r27315); mpfr_init(r27316); mpfr_init(r27317); mpfr_init(r27318); mpfr_init(r27319); } double f_fm(double x) { mpfr_set_d(r27253, x, MPFR_RNDN); ; mpfr_set_si(r27255, mpfr_cmp(r27253, r27254) <= 0, MPFR_RNDN); ; mpfr_set_si(r27257, mpfr_cmp(r27253, r27256) <= 0, MPFR_RNDN); mpfr_set_si(r27258, !mpfr_get_si(r27257, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r27259, mpfr_get_si(r27255, MPFR_RNDN) || mpfr_get_si(r27258, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r27262, r27253, r27261, MPFR_RNDN); mpfr_div(r27263, r27260, r27262, MPFR_RNDN); ; ; mpfr_pow(r27266, r27253, r27265, MPFR_RNDN); mpfr_div(r27267, r27264, r27266, MPFR_RNDN); mpfr_add(r27268, r27263, r27267, MPFR_RNDN); ; mpfr_div(r27270, r27269, r27253, MPFR_RNDN); mpfr_fma(r27271, r27268, r27253, r27270, MPFR_RNDN); mpfr_mul(r27272, r27253, r27253, MPFR_RNDN); mpfr_mul(r27273, r27272, r27272, MPFR_RNDN); ; mpfr_mul(r27275, r27253, r27274, MPFR_RNDN); mpfr_mul(r27276, r27275, r27253, MPFR_RNDN); ; ; ; mpfr_fma(r27280, r27272, r27278, r27279, MPFR_RNDN); mpfr_fma(r27281, r27277, r27273, r27280, MPFR_RNDN); mpfr_fma(r27282, r27273, r27276, r27281, MPFR_RNDN); ; mpfr_mul(r27284, r27283, r27272, MPFR_RNDN); ; mpfr_add(r27286, r27284, r27285, MPFR_RNDN); ; mpfr_add(r27288, r27287, r27279, MPFR_RNDN); mpfr_pow(r27289, r27272, r27288, MPFR_RNDN); mpfr_mul(r27290, r27286, r27289, MPFR_RNDN); mpfr_add(r27291, r27282, r27290, MPFR_RNDN); mpfr_sqrt(r27292, r27291, MPFR_RNDN); ; mpfr_mul(r27294, r27293, r27253, MPFR_RNDN); mpfr_mul(r27295, r27294, r27273, MPFR_RNDN); ; ; mpfr_mul(r27298, r27253, r27297, MPFR_RNDN); mpfr_fma(r27299, r27253, r27298, r27279, MPFR_RNDN); mpfr_fma(r27300, r27273, r27296, r27299, MPFR_RNDN); mpfr_fma(r27301, r27253, r27295, r27300, MPFR_RNDN); ; mpfr_mul(r27303, r27302, r27283, MPFR_RNDN); ; mpfr_mul(r27305, r27304, r27253, MPFR_RNDN); mpfr_mul(r27306, r27305, r27253, MPFR_RNDN); mpfr_fma(r27307, r27273, r27303, r27306, MPFR_RNDN); ; mpfr_mul(r27309, r27308, r27289, MPFR_RNDN); mpfr_fma(r27310, r27289, r27307, r27309, MPFR_RNDN); mpfr_add(r27311, r27301, r27310, MPFR_RNDN); mpfr_div(r27312, r27292, r27311, MPFR_RNDN); mpfr_mul(r27313, r27272, r27274, MPFR_RNDN); mpfr_fma(r27314, r27273, r27313, r27281, MPFR_RNDN); mpfr_add(r27315, r27290, r27314, MPFR_RNDN); mpfr_sqrt(r27316, r27315, MPFR_RNDN); mpfr_mul(r27317, r27312, r27316, MPFR_RNDN); mpfr_mul(r27318, r27253, r27317, MPFR_RNDN); if (mpfr_get_si(r27259, MPFR_RNDN)) { mpfr_set(r27319, r27271, MPFR_RNDN); } else { mpfr_set(r27319, r27318, MPFR_RNDN); }; return mpfr_get_d(r27319, MPFR_RNDN); } static mpfr_t 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, r27361, r27362, r27363, r27364, r27365, r27366, r27367, r27368, r27369, r27370, r27371, r27372, r27373, r27374, r27375, r27376, r27377, r27378, r27379, r27380, r27381, r27382, r27383, r27384, r27385, r27386; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r27320); mpfr_init_set_str(r27321, "-512314.0678876898", 10, MPFR_RNDN); mpfr_init(r27322); mpfr_init_set_str(r27323, "723.4643878500376", 10, MPFR_RNDN); mpfr_init(r27324); mpfr_init(r27325); mpfr_init(r27326); mpfr_init_set_str(r27327, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init_set_str(r27328, "4", 10, MPFR_RNDN); mpfr_init(r27329); mpfr_init(r27330); mpfr_init_set_str(r27331, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r27332, "6", 10, MPFR_RNDN); mpfr_init(r27333); mpfr_init(r27334); mpfr_init(r27335); mpfr_init_set_str(r27336, "0.5", 10, MPFR_RNDN); mpfr_init(r27337); mpfr_init(r27338); mpfr_init(r27339); mpfr_init(r27340); mpfr_init_set_str(r27341, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r27342); mpfr_init(r27343); mpfr_init_set_str(r27344, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r27345, "0.1049934947", 10, MPFR_RNDN); mpfr_init_set_str(r27346, "1", 10, MPFR_RNDN); mpfr_init(r27347); mpfr_init(r27348); mpfr_init(r27349); mpfr_init_set_str(r27350, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r27351); mpfr_init_set_str(r27352, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r27353); mpfr_init_set_str(r27354, "3", 10, MPFR_RNDN); mpfr_init(r27355); mpfr_init(r27356); mpfr_init(r27357); mpfr_init(r27358); mpfr_init(r27359); mpfr_init_set_str(r27360, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r27361); mpfr_init(r27362); mpfr_init_set_str(r27363, "0.2909738639", 10, MPFR_RNDN); mpfr_init_set_str(r27364, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r27365); mpfr_init(r27366); mpfr_init(r27367); mpfr_init(r27368); mpfr_init_set_str(r27369, "2", 10, MPFR_RNDN); mpfr_init(r27370); mpfr_init_set_str(r27371, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r27372); mpfr_init(r27373); mpfr_init(r27374); mpfr_init_set_str(r27375, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r27376); mpfr_init(r27377); mpfr_init(r27378); mpfr_init(r27379); mpfr_init(r27380); mpfr_init(r27381); mpfr_init(r27382); mpfr_init(r27383); mpfr_init(r27384); mpfr_init(r27385); mpfr_init(r27386); } double f_dm(double x) { mpfr_set_d(r27320, x, MPFR_RNDN); ; mpfr_set_si(r27322, mpfr_cmp(r27320, r27321) <= 0, MPFR_RNDN); ; mpfr_set_si(r27324, mpfr_cmp(r27320, r27323) <= 0, MPFR_RNDN); mpfr_set_si(r27325, !mpfr_get_si(r27324, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r27326, mpfr_get_si(r27322, MPFR_RNDN) || mpfr_get_si(r27325, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r27329, r27320, r27328, MPFR_RNDN); mpfr_div(r27330, r27327, r27329, MPFR_RNDN); ; ; mpfr_pow(r27333, r27320, r27332, MPFR_RNDN); mpfr_div(r27334, r27331, r27333, MPFR_RNDN); mpfr_add(r27335, r27330, r27334, MPFR_RNDN); ; mpfr_div(r27337, r27336, r27320, MPFR_RNDN); mpfr_fma(r27338, r27335, r27320, r27337, MPFR_RNDN); mpfr_mul(r27339, r27320, r27320, MPFR_RNDN); mpfr_mul(r27340, r27339, r27339, MPFR_RNDN); ; mpfr_mul(r27342, r27320, r27341, MPFR_RNDN); mpfr_mul(r27343, r27342, r27320, MPFR_RNDN); ; ; ; mpfr_fma(r27347, r27339, r27345, r27346, MPFR_RNDN); mpfr_fma(r27348, r27344, r27340, r27347, MPFR_RNDN); mpfr_fma(r27349, r27340, r27343, r27348, MPFR_RNDN); ; mpfr_mul(r27351, r27350, r27339, MPFR_RNDN); ; mpfr_add(r27353, r27351, r27352, MPFR_RNDN); ; mpfr_add(r27355, r27354, r27346, MPFR_RNDN); mpfr_pow(r27356, r27339, r27355, MPFR_RNDN); mpfr_mul(r27357, r27353, r27356, MPFR_RNDN); mpfr_add(r27358, r27349, r27357, MPFR_RNDN); mpfr_sqrt(r27359, r27358, MPFR_RNDN); ; mpfr_mul(r27361, r27360, r27320, MPFR_RNDN); mpfr_mul(r27362, r27361, r27340, MPFR_RNDN); ; ; mpfr_mul(r27365, r27320, r27364, MPFR_RNDN); mpfr_fma(r27366, r27320, r27365, r27346, MPFR_RNDN); mpfr_fma(r27367, r27340, r27363, r27366, MPFR_RNDN); mpfr_fma(r27368, r27320, r27362, r27367, MPFR_RNDN); ; mpfr_mul(r27370, r27369, r27350, MPFR_RNDN); ; mpfr_mul(r27372, r27371, r27320, MPFR_RNDN); mpfr_mul(r27373, r27372, r27320, MPFR_RNDN); mpfr_fma(r27374, r27340, r27370, r27373, MPFR_RNDN); ; mpfr_mul(r27376, r27375, r27356, MPFR_RNDN); mpfr_fma(r27377, r27356, r27374, r27376, MPFR_RNDN); mpfr_add(r27378, r27368, r27377, MPFR_RNDN); mpfr_div(r27379, r27359, r27378, MPFR_RNDN); mpfr_mul(r27380, r27339, r27341, MPFR_RNDN); mpfr_fma(r27381, r27340, r27380, r27348, MPFR_RNDN); mpfr_add(r27382, r27357, r27381, MPFR_RNDN); mpfr_sqrt(r27383, r27382, MPFR_RNDN); mpfr_mul(r27384, r27379, r27383, MPFR_RNDN); mpfr_mul(r27385, r27320, r27384, MPFR_RNDN); if (mpfr_get_si(r27326, MPFR_RNDN)) { mpfr_set(r27386, r27338, MPFR_RNDN); } else { mpfr_set(r27386, r27385, MPFR_RNDN); }; return mpfr_get_d(r27386, MPFR_RNDN); }