#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r25115 = 1; float r25116 = 0.1049934947; float r25117 = x; float r25118 = r25117 * r25117; float r25119 = r25116 * r25118; float r25120 = r25115 + r25119; float r25121 = 0.0424060604; float r25122 = r25118 * r25118; float r25123 = r25121 * r25122; float r25124 = r25120 + r25123; float r25125 = 0.0072644182; float r25126 = r25122 * r25118; float r25127 = r25125 * r25126; float r25128 = r25124 + r25127; float r25129 = 0.0005064034; float r25130 = r25126 * r25118; float r25131 = r25129 * r25130; float r25132 = r25128 + r25131; float r25133 = 0.0001789971; float r25134 = r25130 * r25118; float r25135 = r25133 * r25134; float r25136 = r25132 + r25135; float r25137 = 0.7715471019; float r25138 = r25137 * r25118; float r25139 = r25115 + r25138; float r25140 = 0.2909738639; float r25141 = r25140 * r25122; float r25142 = r25139 + r25141; float r25143 = 0.0694555761; float r25144 = r25143 * r25126; float r25145 = r25142 + r25144; float r25146 = 0.0140005442; float r25147 = r25146 * r25130; float r25148 = r25145 + r25147; float r25149 = 0.0008327945; float r25150 = r25149 * r25134; float r25151 = r25148 + r25150; float r25152 = 2; float r25153 = r25152 * r25133; float r25154 = r25134 * r25118; float r25155 = r25153 * r25154; float r25156 = r25151 + r25155; float r25157 = r25136 / r25156; float r25158 = r25157 * r25117; return r25158; } double f_id(double x) { double r25159 = 1; double r25160 = 0.1049934947; double r25161 = x; double r25162 = r25161 * r25161; double r25163 = r25160 * r25162; double r25164 = r25159 + r25163; double r25165 = 0.0424060604; double r25166 = r25162 * r25162; double r25167 = r25165 * r25166; double r25168 = r25164 + r25167; double r25169 = 0.0072644182; double r25170 = r25166 * r25162; double r25171 = r25169 * r25170; double r25172 = r25168 + r25171; double r25173 = 0.0005064034; double r25174 = r25170 * r25162; double r25175 = r25173 * r25174; double r25176 = r25172 + r25175; double r25177 = 0.0001789971; double r25178 = r25174 * r25162; double r25179 = r25177 * r25178; double r25180 = r25176 + r25179; double r25181 = 0.7715471019; double r25182 = r25181 * r25162; double r25183 = r25159 + r25182; double r25184 = 0.2909738639; double r25185 = r25184 * r25166; double r25186 = r25183 + r25185; double r25187 = 0.0694555761; double r25188 = r25187 * r25170; double r25189 = r25186 + r25188; double r25190 = 0.0140005442; double r25191 = r25190 * r25174; double r25192 = r25189 + r25191; double r25193 = 0.0008327945; double r25194 = r25193 * r25178; double r25195 = r25192 + r25194; double r25196 = 2; double r25197 = r25196 * r25177; double r25198 = r25178 * r25162; double r25199 = r25197 * r25198; double r25200 = r25195 + r25199; double r25201 = r25180 / r25200; double r25202 = r25201 * r25161; return r25202; } double f_of(float x) { float r25203 = x; float r25204 = -1006989.1225998843; bool r25205 = r25203 <= r25204; float r25206 = 0.2514179000665375; float r25207 = 4; float r25208 = pow(r25203, r25207); float r25209 = r25206 / r25208; float r25210 = 0.15298196345929327; float r25211 = 6; float r25212 = pow(r25203, r25211); float r25213 = r25210 / r25212; float r25214 = r25209 + r25213; float r25215 = 0.5; float r25216 = r25215 / r25203; float r25217 = fma(r25214, r25203, r25216); float r25218 = 60103344.098054335; bool r25219 = r25203 <= r25218; float r25220 = r25203 * r25203; float r25221 = 3; float r25222 = 1; float r25223 = r25221 + r25222; float r25224 = pow(r25220, r25223); float r25225 = 0.0001789971; float r25226 = r25224 * r25225; float r25227 = 0.0005064034; float r25228 = r25224 * r25227; float r25229 = fma(r25226, r25220, r25228); float r25230 = r25220 * r25220; float r25231 = 0.0072644182; float r25232 = r25203 * r25231; float r25233 = r25232 * r25203; float r25234 = 0.0424060604; float r25235 = 0.1049934947; float r25236 = fma(r25220, r25235, r25222); float r25237 = fma(r25234, r25230, r25236); float r25238 = fma(r25230, r25233, r25237); float r25239 = r25229 + r25238; float r25240 = 2; float r25241 = r25225 * r25240; float r25242 = 0.0008327945; float r25243 = r25242 * r25220; float r25244 = fma(r25241, r25230, r25243); float r25245 = 0.0140005442; float r25246 = r25224 * r25245; float r25247 = fma(r25224, r25244, r25246); float r25248 = 0.0694555761; float r25249 = r25203 * r25248; float r25250 = r25203 * r25249; float r25251 = 0.2909738639; float r25252 = 0.7715471019; float r25253 = fma(r25252, r25220, r25222); float r25254 = fma(r25230, r25251, r25253); float r25255 = fma(r25250, r25230, r25254); float r25256 = r25247 + r25255; float r25257 = r25239 / r25256; float r25258 = sqrt(r25257); float r25259 = r25258 * r25258; float r25260 = r25259 * r25203; float r25261 = r25219 ? r25260 : r25217; float r25262 = r25205 ? r25217 : r25261; return r25262; } double f_od(double x) { double r25263 = x; double r25264 = -1006989.1225998843; bool r25265 = r25263 <= r25264; double r25266 = 0.2514179000665375; double r25267 = 4; double r25268 = pow(r25263, r25267); double r25269 = r25266 / r25268; double r25270 = 0.15298196345929327; double r25271 = 6; double r25272 = pow(r25263, r25271); double r25273 = r25270 / r25272; double r25274 = r25269 + r25273; double r25275 = 0.5; double r25276 = r25275 / r25263; double r25277 = fma(r25274, r25263, r25276); double r25278 = 60103344.098054335; bool r25279 = r25263 <= r25278; double r25280 = r25263 * r25263; double r25281 = 3; double r25282 = 1; double r25283 = r25281 + r25282; double r25284 = pow(r25280, r25283); double r25285 = 0.0001789971; double r25286 = r25284 * r25285; double r25287 = 0.0005064034; double r25288 = r25284 * r25287; double r25289 = fma(r25286, r25280, r25288); double r25290 = r25280 * r25280; double r25291 = 0.0072644182; double r25292 = r25263 * r25291; double r25293 = r25292 * r25263; double r25294 = 0.0424060604; double r25295 = 0.1049934947; double r25296 = fma(r25280, r25295, r25282); double r25297 = fma(r25294, r25290, r25296); double r25298 = fma(r25290, r25293, r25297); double r25299 = r25289 + r25298; double r25300 = 2; double r25301 = r25285 * r25300; double r25302 = 0.0008327945; double r25303 = r25302 * r25280; double r25304 = fma(r25301, r25290, r25303); double r25305 = 0.0140005442; double r25306 = r25284 * r25305; double r25307 = fma(r25284, r25304, r25306); double r25308 = 0.0694555761; double r25309 = r25263 * r25308; double r25310 = r25263 * r25309; double r25311 = 0.2909738639; double r25312 = 0.7715471019; double r25313 = fma(r25312, r25280, r25282); double r25314 = fma(r25290, r25311, r25313); double r25315 = fma(r25310, r25290, r25314); double r25316 = r25307 + r25315; double r25317 = r25299 / r25316; double r25318 = sqrt(r25317); double r25319 = r25318 * r25318; double r25320 = r25319 * r25263; double r25321 = r25279 ? r25320 : r25277; double r25322 = r25265 ? r25277 : r25321; return r25322; } 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 r25323, r25324, r25325, r25326, r25327, r25328, r25329, r25330, r25331, r25332, r25333, r25334, r25335, r25336, r25337, r25338, r25339, r25340, r25341, r25342, r25343, r25344, r25345, r25346, r25347, r25348, r25349, r25350, r25351, r25352, r25353, r25354, r25355, r25356, r25357, r25358, r25359, r25360, r25361, r25362, r25363, r25364, r25365, r25366; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r25323, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25324, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r25325); mpfr_init(r25326); mpfr_init(r25327); mpfr_init(r25328); mpfr_init_set_str(r25329, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r25330); mpfr_init(r25331); mpfr_init(r25332); mpfr_init_set_str(r25333, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r25334); mpfr_init(r25335); mpfr_init(r25336); mpfr_init_set_str(r25337, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r25338); mpfr_init(r25339); mpfr_init(r25340); mpfr_init_set_str(r25341, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r25342); mpfr_init(r25343); mpfr_init(r25344); mpfr_init_set_str(r25345, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r25346); mpfr_init(r25347); mpfr_init_set_str(r25348, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r25349); mpfr_init(r25350); mpfr_init_set_str(r25351, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r25352); mpfr_init(r25353); mpfr_init_set_str(r25354, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r25355); mpfr_init(r25356); mpfr_init_set_str(r25357, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r25358); mpfr_init(r25359); mpfr_init_set_str(r25360, "2", 10, MPFR_RNDN); mpfr_init(r25361); mpfr_init(r25362); mpfr_init(r25363); mpfr_init(r25364); mpfr_init(r25365); mpfr_init(r25366); } double f_im(double x) { ; ; mpfr_set_d(r25325, x, MPFR_RNDN); mpfr_mul(r25326, r25325, r25325, MPFR_RNDN); mpfr_mul(r25327, r25324, r25326, MPFR_RNDN); mpfr_add(r25328, r25323, r25327, MPFR_RNDN); ; mpfr_mul(r25330, r25326, r25326, MPFR_RNDN); mpfr_mul(r25331, r25329, r25330, MPFR_RNDN); mpfr_add(r25332, r25328, r25331, MPFR_RNDN); ; mpfr_mul(r25334, r25330, r25326, MPFR_RNDN); mpfr_mul(r25335, r25333, r25334, MPFR_RNDN); mpfr_add(r25336, r25332, r25335, MPFR_RNDN); ; mpfr_mul(r25338, r25334, r25326, MPFR_RNDN); mpfr_mul(r25339, r25337, r25338, MPFR_RNDN); mpfr_add(r25340, r25336, r25339, MPFR_RNDN); ; mpfr_mul(r25342, r25338, r25326, MPFR_RNDN); mpfr_mul(r25343, r25341, r25342, MPFR_RNDN); mpfr_add(r25344, r25340, r25343, MPFR_RNDN); ; mpfr_mul(r25346, r25345, r25326, MPFR_RNDN); mpfr_add(r25347, r25323, r25346, MPFR_RNDN); ; mpfr_mul(r25349, r25348, r25330, MPFR_RNDN); mpfr_add(r25350, r25347, r25349, MPFR_RNDN); ; mpfr_mul(r25352, r25351, r25334, MPFR_RNDN); mpfr_add(r25353, r25350, r25352, MPFR_RNDN); ; mpfr_mul(r25355, r25354, r25338, MPFR_RNDN); mpfr_add(r25356, r25353, r25355, MPFR_RNDN); ; mpfr_mul(r25358, r25357, r25342, MPFR_RNDN); mpfr_add(r25359, r25356, r25358, MPFR_RNDN); ; mpfr_mul(r25361, r25360, r25341, MPFR_RNDN); mpfr_mul(r25362, r25342, r25326, MPFR_RNDN); mpfr_mul(r25363, r25361, r25362, MPFR_RNDN); mpfr_add(r25364, r25359, r25363, MPFR_RNDN); mpfr_div(r25365, r25344, r25364, MPFR_RNDN); mpfr_mul(r25366, r25365, r25325, MPFR_RNDN); return mpfr_get_d(r25366, MPFR_RNDN); } static mpfr_t r25367, r25368, r25369, r25370, r25371, r25372, r25373, r25374, r25375, r25376, r25377, r25378, r25379, r25380, r25381, r25382, r25383, r25384, r25385, r25386, r25387, r25388, r25389, r25390, r25391, r25392, r25393, r25394, r25395, r25396, r25397, r25398, r25399, r25400, r25401, r25402, r25403, r25404, r25405, r25406, r25407, r25408, r25409, r25410, r25411, r25412, r25413, r25414, r25415, r25416, r25417, r25418, r25419, r25420, r25421, r25422, r25423, r25424, r25425, r25426; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r25367); mpfr_init_set_str(r25368, "-1006989.1225998843", 10, MPFR_RNDN); mpfr_init(r25369); mpfr_init_set_str(r25370, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init_set_str(r25371, "4", 10, MPFR_RNDN); mpfr_init(r25372); mpfr_init(r25373); mpfr_init_set_str(r25374, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r25375, "6", 10, MPFR_RNDN); mpfr_init(r25376); mpfr_init(r25377); mpfr_init(r25378); mpfr_init_set_str(r25379, "0.5", 10, MPFR_RNDN); mpfr_init(r25380); mpfr_init(r25381); mpfr_init_set_str(r25382, "60103344.098054335", 10, MPFR_RNDN); mpfr_init(r25383); mpfr_init(r25384); mpfr_init_set_str(r25385, "3", 10, MPFR_RNDN); mpfr_init_set_str(r25386, "1", 10, MPFR_RNDN); mpfr_init(r25387); mpfr_init(r25388); mpfr_init_set_str(r25389, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r25390); mpfr_init_set_str(r25391, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r25392); mpfr_init(r25393); mpfr_init(r25394); mpfr_init_set_str(r25395, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r25396); mpfr_init(r25397); mpfr_init_set_str(r25398, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r25399, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r25400); mpfr_init(r25401); mpfr_init(r25402); mpfr_init(r25403); mpfr_init_set_str(r25404, "2", 10, MPFR_RNDN); mpfr_init(r25405); mpfr_init_set_str(r25406, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r25407); mpfr_init(r25408); mpfr_init_set_str(r25409, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r25410); mpfr_init(r25411); mpfr_init_set_str(r25412, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r25413); mpfr_init(r25414); mpfr_init_set_str(r25415, "0.2909738639", 10, MPFR_RNDN); mpfr_init_set_str(r25416, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r25417); mpfr_init(r25418); mpfr_init(r25419); mpfr_init(r25420); mpfr_init(r25421); mpfr_init(r25422); mpfr_init(r25423); mpfr_init(r25424); mpfr_init(r25425); mpfr_init(r25426); } double f_fm(double x) { mpfr_set_d(r25367, x, MPFR_RNDN); ; mpfr_set_si(r25369, mpfr_cmp(r25367, r25368) <= 0, MPFR_RNDN); ; ; mpfr_pow(r25372, r25367, r25371, MPFR_RNDN); mpfr_div(r25373, r25370, r25372, MPFR_RNDN); ; ; mpfr_pow(r25376, r25367, r25375, MPFR_RNDN); mpfr_div(r25377, r25374, r25376, MPFR_RNDN); mpfr_add(r25378, r25373, r25377, MPFR_RNDN); ; mpfr_div(r25380, r25379, r25367, MPFR_RNDN); mpfr_fma(r25381, r25378, r25367, r25380, MPFR_RNDN); ; mpfr_set_si(r25383, mpfr_cmp(r25367, r25382) <= 0, MPFR_RNDN); mpfr_mul(r25384, r25367, r25367, MPFR_RNDN); ; ; mpfr_add(r25387, r25385, r25386, MPFR_RNDN); mpfr_pow(r25388, r25384, r25387, MPFR_RNDN); ; mpfr_mul(r25390, r25388, r25389, MPFR_RNDN); ; mpfr_mul(r25392, r25388, r25391, MPFR_RNDN); mpfr_fma(r25393, r25390, r25384, r25392, MPFR_RNDN); mpfr_mul(r25394, r25384, r25384, MPFR_RNDN); ; mpfr_mul(r25396, r25367, r25395, MPFR_RNDN); mpfr_mul(r25397, r25396, r25367, MPFR_RNDN); ; ; mpfr_fma(r25400, r25384, r25399, r25386, MPFR_RNDN); mpfr_fma(r25401, r25398, r25394, r25400, MPFR_RNDN); mpfr_fma(r25402, r25394, r25397, r25401, MPFR_RNDN); mpfr_add(r25403, r25393, r25402, MPFR_RNDN); ; mpfr_mul(r25405, r25389, r25404, MPFR_RNDN); ; mpfr_mul(r25407, r25406, r25384, MPFR_RNDN); mpfr_fma(r25408, r25405, r25394, r25407, MPFR_RNDN); ; mpfr_mul(r25410, r25388, r25409, MPFR_RNDN); mpfr_fma(r25411, r25388, r25408, r25410, MPFR_RNDN); ; mpfr_mul(r25413, r25367, r25412, MPFR_RNDN); mpfr_mul(r25414, r25367, r25413, MPFR_RNDN); ; ; mpfr_fma(r25417, r25416, r25384, r25386, MPFR_RNDN); mpfr_fma(r25418, r25394, r25415, r25417, MPFR_RNDN); mpfr_fma(r25419, r25414, r25394, r25418, MPFR_RNDN); mpfr_add(r25420, r25411, r25419, MPFR_RNDN); mpfr_div(r25421, r25403, r25420, MPFR_RNDN); mpfr_sqrt(r25422, r25421, MPFR_RNDN); mpfr_mul(r25423, r25422, r25422, MPFR_RNDN); mpfr_mul(r25424, r25423, r25367, MPFR_RNDN); if (mpfr_get_si(r25383, MPFR_RNDN)) { mpfr_set(r25425, r25424, MPFR_RNDN); } else { mpfr_set(r25425, r25381, MPFR_RNDN); }; if (mpfr_get_si(r25369, MPFR_RNDN)) { mpfr_set(r25426, r25381, MPFR_RNDN); } else { mpfr_set(r25426, r25425, MPFR_RNDN); }; return mpfr_get_d(r25426, MPFR_RNDN); } static mpfr_t r25427, r25428, r25429, r25430, r25431, r25432, r25433, r25434, r25435, r25436, r25437, r25438, r25439, r25440, r25441, r25442, r25443, r25444, r25445, r25446, r25447, r25448, r25449, r25450, r25451, r25452, r25453, r25454, r25455, r25456, r25457, r25458, r25459, r25460, r25461, r25462, r25463, r25464, r25465, r25466, r25467, r25468, r25469, r25470, r25471, r25472, r25473, r25474, r25475, r25476, r25477, r25478, r25479, r25480, r25481, r25482, r25483, r25484, r25485, r25486; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r25427); mpfr_init_set_str(r25428, "-1006989.1225998843", 10, MPFR_RNDN); mpfr_init(r25429); mpfr_init_set_str(r25430, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init_set_str(r25431, "4", 10, MPFR_RNDN); mpfr_init(r25432); mpfr_init(r25433); mpfr_init_set_str(r25434, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r25435, "6", 10, MPFR_RNDN); mpfr_init(r25436); mpfr_init(r25437); mpfr_init(r25438); mpfr_init_set_str(r25439, "0.5", 10, MPFR_RNDN); mpfr_init(r25440); mpfr_init(r25441); mpfr_init_set_str(r25442, "60103344.098054335", 10, MPFR_RNDN); mpfr_init(r25443); mpfr_init(r25444); mpfr_init_set_str(r25445, "3", 10, MPFR_RNDN); mpfr_init_set_str(r25446, "1", 10, MPFR_RNDN); mpfr_init(r25447); mpfr_init(r25448); mpfr_init_set_str(r25449, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r25450); mpfr_init_set_str(r25451, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r25452); mpfr_init(r25453); mpfr_init(r25454); mpfr_init_set_str(r25455, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r25456); mpfr_init(r25457); mpfr_init_set_str(r25458, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r25459, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r25460); mpfr_init(r25461); mpfr_init(r25462); mpfr_init(r25463); mpfr_init_set_str(r25464, "2", 10, MPFR_RNDN); mpfr_init(r25465); mpfr_init_set_str(r25466, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r25467); mpfr_init(r25468); mpfr_init_set_str(r25469, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r25470); mpfr_init(r25471); mpfr_init_set_str(r25472, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r25473); mpfr_init(r25474); mpfr_init_set_str(r25475, "0.2909738639", 10, MPFR_RNDN); mpfr_init_set_str(r25476, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r25477); mpfr_init(r25478); mpfr_init(r25479); mpfr_init(r25480); mpfr_init(r25481); mpfr_init(r25482); mpfr_init(r25483); mpfr_init(r25484); mpfr_init(r25485); mpfr_init(r25486); } double f_dm(double x) { mpfr_set_d(r25427, x, MPFR_RNDN); ; mpfr_set_si(r25429, mpfr_cmp(r25427, r25428) <= 0, MPFR_RNDN); ; ; mpfr_pow(r25432, r25427, r25431, MPFR_RNDN); mpfr_div(r25433, r25430, r25432, MPFR_RNDN); ; ; mpfr_pow(r25436, r25427, r25435, MPFR_RNDN); mpfr_div(r25437, r25434, r25436, MPFR_RNDN); mpfr_add(r25438, r25433, r25437, MPFR_RNDN); ; mpfr_div(r25440, r25439, r25427, MPFR_RNDN); mpfr_fma(r25441, r25438, r25427, r25440, MPFR_RNDN); ; mpfr_set_si(r25443, mpfr_cmp(r25427, r25442) <= 0, MPFR_RNDN); mpfr_mul(r25444, r25427, r25427, MPFR_RNDN); ; ; mpfr_add(r25447, r25445, r25446, MPFR_RNDN); mpfr_pow(r25448, r25444, r25447, MPFR_RNDN); ; mpfr_mul(r25450, r25448, r25449, MPFR_RNDN); ; mpfr_mul(r25452, r25448, r25451, MPFR_RNDN); mpfr_fma(r25453, r25450, r25444, r25452, MPFR_RNDN); mpfr_mul(r25454, r25444, r25444, MPFR_RNDN); ; mpfr_mul(r25456, r25427, r25455, MPFR_RNDN); mpfr_mul(r25457, r25456, r25427, MPFR_RNDN); ; ; mpfr_fma(r25460, r25444, r25459, r25446, MPFR_RNDN); mpfr_fma(r25461, r25458, r25454, r25460, MPFR_RNDN); mpfr_fma(r25462, r25454, r25457, r25461, MPFR_RNDN); mpfr_add(r25463, r25453, r25462, MPFR_RNDN); ; mpfr_mul(r25465, r25449, r25464, MPFR_RNDN); ; mpfr_mul(r25467, r25466, r25444, MPFR_RNDN); mpfr_fma(r25468, r25465, r25454, r25467, MPFR_RNDN); ; mpfr_mul(r25470, r25448, r25469, MPFR_RNDN); mpfr_fma(r25471, r25448, r25468, r25470, MPFR_RNDN); ; mpfr_mul(r25473, r25427, r25472, MPFR_RNDN); mpfr_mul(r25474, r25427, r25473, MPFR_RNDN); ; ; mpfr_fma(r25477, r25476, r25444, r25446, MPFR_RNDN); mpfr_fma(r25478, r25454, r25475, r25477, MPFR_RNDN); mpfr_fma(r25479, r25474, r25454, r25478, MPFR_RNDN); mpfr_add(r25480, r25471, r25479, MPFR_RNDN); mpfr_div(r25481, r25463, r25480, MPFR_RNDN); mpfr_sqrt(r25482, r25481, MPFR_RNDN); mpfr_mul(r25483, r25482, r25482, MPFR_RNDN); mpfr_mul(r25484, r25483, r25427, MPFR_RNDN); if (mpfr_get_si(r25443, MPFR_RNDN)) { mpfr_set(r25485, r25484, MPFR_RNDN); } else { mpfr_set(r25485, r25441, MPFR_RNDN); }; if (mpfr_get_si(r25429, MPFR_RNDN)) { mpfr_set(r25486, r25441, MPFR_RNDN); } else { mpfr_set(r25486, r25485, MPFR_RNDN); }; return mpfr_get_d(r25486, MPFR_RNDN); }