#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r19023 = 1.0f; float r19024 = 0.1049934947f; float r19025 = x; float r19026 = r19025 * r19025; float r19027 = r19024 * r19026; float r19028 = r19023 + r19027; float r19029 = 0.0424060604f; float r19030 = r19026 * r19026; float r19031 = r19029 * r19030; float r19032 = r19028 + r19031; float r19033 = 0.0072644182f; float r19034 = r19030 * r19026; float r19035 = r19033 * r19034; float r19036 = r19032 + r19035; float r19037 = 0.0005064034f; float r19038 = r19034 * r19026; float r19039 = r19037 * r19038; float r19040 = r19036 + r19039; float r19041 = 0.0001789971f; float r19042 = r19038 * r19026; float r19043 = r19041 * r19042; float r19044 = r19040 + r19043; float r19045 = 0.7715471019f; float r19046 = r19045 * r19026; float r19047 = r19023 + r19046; float r19048 = 0.2909738639f; float r19049 = r19048 * r19030; float r19050 = r19047 + r19049; float r19051 = 0.0694555761f; float r19052 = r19051 * r19034; float r19053 = r19050 + r19052; float r19054 = 0.0140005442f; float r19055 = r19054 * r19038; float r19056 = r19053 + r19055; float r19057 = 0.0008327945f; float r19058 = r19057 * r19042; float r19059 = r19056 + r19058; float r19060 = 2.0f; float r19061 = r19060 * r19041; float r19062 = r19042 * r19026; float r19063 = r19061 * r19062; float r19064 = r19059 + r19063; float r19065 = r19044 / r19064; float r19066 = r19065 * r19025; return r19066; } double f_id(double x) { double r19067 = 1.0; double r19068 = 0.1049934947; double r19069 = x; double r19070 = r19069 * r19069; double r19071 = r19068 * r19070; double r19072 = r19067 + r19071; double r19073 = 0.0424060604; double r19074 = r19070 * r19070; double r19075 = r19073 * r19074; double r19076 = r19072 + r19075; double r19077 = 0.0072644182; double r19078 = r19074 * r19070; double r19079 = r19077 * r19078; double r19080 = r19076 + r19079; double r19081 = 0.0005064034; double r19082 = r19078 * r19070; double r19083 = r19081 * r19082; double r19084 = r19080 + r19083; double r19085 = 0.0001789971; double r19086 = r19082 * r19070; double r19087 = r19085 * r19086; double r19088 = r19084 + r19087; double r19089 = 0.7715471019; double r19090 = r19089 * r19070; double r19091 = r19067 + r19090; double r19092 = 0.2909738639; double r19093 = r19092 * r19074; double r19094 = r19091 + r19093; double r19095 = 0.0694555761; double r19096 = r19095 * r19078; double r19097 = r19094 + r19096; double r19098 = 0.0140005442; double r19099 = r19098 * r19082; double r19100 = r19097 + r19099; double r19101 = 0.0008327945; double r19102 = r19101 * r19086; double r19103 = r19100 + r19102; double r19104 = 2.0; double r19105 = r19104 * r19085; double r19106 = r19086 * r19070; double r19107 = r19105 * r19106; double r19108 = r19103 + r19107; double r19109 = r19088 / r19108; double r19110 = r19109 * r19069; return r19110; } double f_of(float x) { float r19111 = x; float r19112 = r19111 * r19111; float r19113 = 0.0001789971f; float r19114 = r19112 * r19113; float r19115 = r19112 * r19112; float r19116 = r19115 * r19115; float r19117 = r19114 * r19116; float r19118 = 0.0424060604f; float r19119 = r19111 * r19118; float r19120 = r19111 * (r19111 * r19111); float r19121 = r19119 * r19120; float r19122 = 1.0f; float r19123 = r19121 + r19122; float r19124 = 0.1049934947f; float r19125 = r19124 * r19112; float r19126 = r19123 + r19125; float r19127 = r19117 + r19126; float r19128 = r19111 * r19127; float r19129 = 0.7715471019f; float r19130 = r19111 * r19129; float r19131 = r19130 * r19111; float r19132 = r19122 + r19131; float r19133 = 0.2909738639f; float r19134 = r19111 * r19133; float r19135 = r19120 * r19134; float r19136 = r19132 + r19135; float r19137 = r19120 * r19120; float r19138 = 0.0694555761f; float r19139 = 0.0140005442f; float r19140 = r19139 * r19111; float r19141 = r19111 * r19140; float r19142 = r19138 + r19141; float r19143 = r19137 * r19142; float r19144 = r19137 * r19115; float r19145 = 2.0f; float r19146 = r19145 * r19113; float r19147 = r19146 * r19112; float r19148 = 0.0008327945f; float r19149 = r19147 + r19148; float r19150 = r19144 * r19149; float r19151 = r19143 + r19150; float r19152 = r19136 + r19151; float r19153 = r19128 / r19152; float r19154 = r19111 * r19137; float r19155 = 0.0005064034f; float r19156 = r19155 * r19112; float r19157 = 0.0072644182f; float r19158 = r19156 + r19157; float r19159 = r19154 * r19158; float r19160 = r19112 * (r19112 * r19112); float r19161 = r19115 * r19160; float r19162 = r19112 * r19146; float r19163 = r19162 + r19148; float r19164 = r19161 * r19163; float r19165 = r19129 * r19111; float r19166 = r19111 * r19165; float r19167 = r19134 * r19120; float r19168 = r19167 + r19122; float r19169 = r19166 + r19168; float r19170 = r19164 + r19169; float r19171 = r19141 + r19138; float r19172 = r19160 * r19171; float r19173 = r19170 + r19172; float r19174 = r19159 / r19173; float r19175 = r19153 + r19174; return r19175; } double f_od(double x) { double r19176 = x; double r19177 = r19176 * r19176; double r19178 = 0.0001789971; double r19179 = r19177 * r19178; double r19180 = r19177 * r19177; double r19181 = r19180 * r19180; double r19182 = r19179 * r19181; double r19183 = 0.0424060604; double r19184 = r19176 * r19183; double r19185 = r19176 * (r19176 * r19176); double r19186 = r19184 * r19185; double r19187 = 1.0; double r19188 = r19186 + r19187; double r19189 = 0.1049934947; double r19190 = r19189 * r19177; double r19191 = r19188 + r19190; double r19192 = r19182 + r19191; double r19193 = r19176 * r19192; double r19194 = 0.7715471019; double r19195 = r19176 * r19194; double r19196 = r19195 * r19176; double r19197 = r19187 + r19196; double r19198 = 0.2909738639; double r19199 = r19176 * r19198; double r19200 = r19185 * r19199; double r19201 = r19197 + r19200; double r19202 = r19185 * r19185; double r19203 = 0.0694555761; double r19204 = 0.0140005442; double r19205 = r19204 * r19176; double r19206 = r19176 * r19205; double r19207 = r19203 + r19206; double r19208 = r19202 * r19207; double r19209 = r19202 * r19180; double r19210 = 2.0; double r19211 = r19210 * r19178; double r19212 = r19211 * r19177; double r19213 = 0.0008327945; double r19214 = r19212 + r19213; double r19215 = r19209 * r19214; double r19216 = r19208 + r19215; double r19217 = r19201 + r19216; double r19218 = r19193 / r19217; double r19219 = r19176 * r19202; double r19220 = 0.0005064034; double r19221 = r19220 * r19177; double r19222 = 0.0072644182; double r19223 = r19221 + r19222; double r19224 = r19219 * r19223; double r19225 = r19177 * (r19177 * r19177); double r19226 = r19180 * r19225; double r19227 = r19177 * r19211; double r19228 = r19227 + r19213; double r19229 = r19226 * r19228; double r19230 = r19194 * r19176; double r19231 = r19176 * r19230; double r19232 = r19199 * r19185; double r19233 = r19232 + r19187; double r19234 = r19231 + r19233; double r19235 = r19229 + r19234; double r19236 = r19206 + r19203; double r19237 = r19225 * r19236; double r19238 = r19235 + r19237; double r19239 = r19224 / r19238; double r19240 = r19218 + r19239; return r19240; } 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 r19241, r19242, r19243, r19244, r19245, r19246, r19247, r19248, r19249, r19250, r19251, r19252, r19253, r19254, r19255, r19256, r19257, r19258, r19259, r19260, r19261, r19262, r19263, r19264, r19265, r19266, r19267, r19268, r19269, r19270, r19271, r19272, r19273, r19274, r19275, r19276, r19277, r19278, r19279, r19280, r19281, r19282, r19283, r19284; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19241, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19242, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19243); mpfr_init(r19244); mpfr_init(r19245); mpfr_init(r19246); mpfr_init_set_str(r19247, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19248); mpfr_init(r19249); mpfr_init(r19250); mpfr_init_set_str(r19251, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19252); mpfr_init(r19253); mpfr_init(r19254); mpfr_init_set_str(r19255, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19256); mpfr_init(r19257); mpfr_init(r19258); mpfr_init_set_str(r19259, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19260); mpfr_init(r19261); mpfr_init(r19262); mpfr_init_set_str(r19263, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19264); mpfr_init(r19265); mpfr_init_set_str(r19266, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19267); mpfr_init(r19268); mpfr_init_set_str(r19269, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19270); mpfr_init(r19271); mpfr_init_set_str(r19272, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19273); mpfr_init(r19274); mpfr_init_set_str(r19275, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19276); mpfr_init(r19277); mpfr_init_set_str(r19278, "2", 10, MPFR_RNDN); mpfr_init(r19279); mpfr_init(r19280); mpfr_init(r19281); mpfr_init(r19282); mpfr_init(r19283); mpfr_init(r19284); } double f_im(double x) { ; ; mpfr_set_d(r19243, x, MPFR_RNDN); mpfr_mul(r19244, r19243, r19243, MPFR_RNDN); mpfr_mul(r19245, r19242, r19244, MPFR_RNDN); mpfr_add(r19246, r19241, r19245, MPFR_RNDN); ; mpfr_mul(r19248, r19244, r19244, MPFR_RNDN); mpfr_mul(r19249, r19247, r19248, MPFR_RNDN); mpfr_add(r19250, r19246, r19249, MPFR_RNDN); ; mpfr_mul(r19252, r19248, r19244, MPFR_RNDN); mpfr_mul(r19253, r19251, r19252, MPFR_RNDN); mpfr_add(r19254, r19250, r19253, MPFR_RNDN); ; mpfr_mul(r19256, r19252, r19244, MPFR_RNDN); mpfr_mul(r19257, r19255, r19256, MPFR_RNDN); mpfr_add(r19258, r19254, r19257, MPFR_RNDN); ; mpfr_mul(r19260, r19256, r19244, MPFR_RNDN); mpfr_mul(r19261, r19259, r19260, MPFR_RNDN); mpfr_add(r19262, r19258, r19261, MPFR_RNDN); ; mpfr_mul(r19264, r19263, r19244, MPFR_RNDN); mpfr_add(r19265, r19241, r19264, MPFR_RNDN); ; mpfr_mul(r19267, r19266, r19248, MPFR_RNDN); mpfr_add(r19268, r19265, r19267, MPFR_RNDN); ; mpfr_mul(r19270, r19269, r19252, MPFR_RNDN); mpfr_add(r19271, r19268, r19270, MPFR_RNDN); ; mpfr_mul(r19273, r19272, r19256, MPFR_RNDN); mpfr_add(r19274, r19271, r19273, MPFR_RNDN); ; mpfr_mul(r19276, r19275, r19260, MPFR_RNDN); mpfr_add(r19277, r19274, r19276, MPFR_RNDN); ; mpfr_mul(r19279, r19278, r19259, MPFR_RNDN); mpfr_mul(r19280, r19260, r19244, MPFR_RNDN); mpfr_mul(r19281, r19279, r19280, MPFR_RNDN); mpfr_add(r19282, r19277, r19281, MPFR_RNDN); mpfr_div(r19283, r19262, r19282, MPFR_RNDN); mpfr_mul(r19284, r19283, r19243, MPFR_RNDN); return mpfr_get_d(r19284, MPFR_RNDN); } static mpfr_t r19285, r19286, r19287, r19288, r19289, r19290, r19291, r19292, r19293, r19294, r19295, r19296, r19297, r19298, r19299, r19300, r19301, r19302, r19303, r19304, r19305, r19306, r19307, r19308, r19309, r19310, r19311, r19312, r19313, r19314, r19315, r19316, r19317, r19318, r19319, r19320, r19321, r19322, r19323, r19324, r19325, r19326, r19327, r19328, r19329, r19330, r19331, r19332, r19333, r19334, r19335, r19336, r19337, r19338, r19339, r19340, r19341, r19342, r19343, r19344, r19345, r19346, r19347, r19348, r19349; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19285); mpfr_init(r19286); mpfr_init_set_str(r19287, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19288); mpfr_init(r19289); mpfr_init(r19290); mpfr_init(r19291); mpfr_init_set_str(r19292, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19293); mpfr_init(r19294); mpfr_init(r19295); mpfr_init_set_str(r19296, "1", 10, MPFR_RNDN); mpfr_init(r19297); mpfr_init_set_str(r19298, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19299); mpfr_init(r19300); mpfr_init(r19301); mpfr_init(r19302); mpfr_init_set_str(r19303, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19304); mpfr_init(r19305); mpfr_init(r19306); mpfr_init_set_str(r19307, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19308); mpfr_init(r19309); mpfr_init(r19310); mpfr_init(r19311); mpfr_init_set_str(r19312, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19313, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19314); mpfr_init(r19315); mpfr_init(r19316); mpfr_init(r19317); mpfr_init(r19318); mpfr_init_set_str(r19319, "2", 10, MPFR_RNDN); mpfr_init(r19320); mpfr_init(r19321); mpfr_init_set_str(r19322, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19323); mpfr_init(r19324); mpfr_init(r19325); mpfr_init(r19326); mpfr_init(r19327); mpfr_init(r19328); mpfr_init_set_str(r19329, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19330); mpfr_init_set_str(r19331, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19332); mpfr_init(r19333); mpfr_init(r19334); mpfr_init(r19335); mpfr_init(r19336); mpfr_init(r19337); mpfr_init(r19338); mpfr_init(r19339); mpfr_init(r19340); mpfr_init(r19341); mpfr_init(r19342); mpfr_init(r19343); mpfr_init(r19344); mpfr_init(r19345); mpfr_init(r19346); mpfr_init(r19347); mpfr_init(r19348); mpfr_init(r19349); } double f_fm(double x) { mpfr_set_d(r19285, x, MPFR_RNDN); mpfr_sqr(r19286, r19285, MPFR_RNDN); ; mpfr_mul(r19288, r19286, r19287, MPFR_RNDN); mpfr_mul(r19289, r19286, r19286, MPFR_RNDN); mpfr_mul(r19290, r19289, r19289, MPFR_RNDN); mpfr_mul(r19291, r19288, r19290, MPFR_RNDN); ; mpfr_mul(r19293, r19285, r19292, MPFR_RNDN); mpfr_mul(r19294, r19285, r19285, MPFR_RNDN); mpfr_mul(r19294, r19294, r19285, MPFR_RNDN); mpfr_mul(r19295, r19293, r19294, MPFR_RNDN); ; mpfr_add(r19297, r19295, r19296, MPFR_RNDN); ; mpfr_mul(r19299, r19298, r19286, MPFR_RNDN); mpfr_add(r19300, r19297, r19299, MPFR_RNDN); mpfr_add(r19301, r19291, r19300, MPFR_RNDN); mpfr_mul(r19302, r19285, r19301, MPFR_RNDN); ; mpfr_mul(r19304, r19285, r19303, MPFR_RNDN); mpfr_mul(r19305, r19304, r19285, MPFR_RNDN); mpfr_add(r19306, r19296, r19305, MPFR_RNDN); ; mpfr_mul(r19308, r19285, r19307, MPFR_RNDN); mpfr_mul(r19309, r19294, r19308, MPFR_RNDN); mpfr_add(r19310, r19306, r19309, MPFR_RNDN); mpfr_mul(r19311, r19294, r19294, MPFR_RNDN); ; ; mpfr_mul(r19314, r19313, r19285, MPFR_RNDN); mpfr_mul(r19315, r19285, r19314, MPFR_RNDN); mpfr_add(r19316, r19312, r19315, MPFR_RNDN); mpfr_mul(r19317, r19311, r19316, MPFR_RNDN); mpfr_mul(r19318, r19311, r19289, MPFR_RNDN); ; mpfr_mul(r19320, r19319, r19287, MPFR_RNDN); mpfr_mul(r19321, r19320, r19286, MPFR_RNDN); ; mpfr_add(r19323, r19321, r19322, MPFR_RNDN); mpfr_mul(r19324, r19318, r19323, MPFR_RNDN); mpfr_add(r19325, r19317, r19324, MPFR_RNDN); mpfr_add(r19326, r19310, r19325, MPFR_RNDN); mpfr_div(r19327, r19302, r19326, MPFR_RNDN); mpfr_mul(r19328, r19285, r19311, MPFR_RNDN); ; mpfr_mul(r19330, r19329, r19286, MPFR_RNDN); ; mpfr_add(r19332, r19330, r19331, MPFR_RNDN); mpfr_mul(r19333, r19328, r19332, MPFR_RNDN); mpfr_mul(r19334, r19286, r19286, MPFR_RNDN); mpfr_mul(r19334, r19334, r19286, MPFR_RNDN); mpfr_mul(r19335, r19289, r19334, MPFR_RNDN); mpfr_mul(r19336, r19286, r19320, MPFR_RNDN); mpfr_add(r19337, r19336, r19322, MPFR_RNDN); mpfr_mul(r19338, r19335, r19337, MPFR_RNDN); mpfr_mul(r19339, r19303, r19285, MPFR_RNDN); mpfr_mul(r19340, r19285, r19339, MPFR_RNDN); mpfr_mul(r19341, r19308, r19294, MPFR_RNDN); mpfr_add(r19342, r19341, r19296, MPFR_RNDN); mpfr_add(r19343, r19340, r19342, MPFR_RNDN); mpfr_add(r19344, r19338, r19343, MPFR_RNDN); mpfr_add(r19345, r19315, r19312, MPFR_RNDN); mpfr_mul(r19346, r19334, r19345, MPFR_RNDN); mpfr_add(r19347, r19344, r19346, MPFR_RNDN); mpfr_div(r19348, r19333, r19347, MPFR_RNDN); mpfr_add(r19349, r19327, r19348, MPFR_RNDN); return mpfr_get_d(r19349, MPFR_RNDN); } static mpfr_t r19350, r19351, r19352, r19353, r19354, r19355, r19356, r19357, r19358, r19359, r19360, r19361, r19362, r19363, r19364, r19365, r19366, r19367, r19368, r19369, r19370, r19371, r19372, r19373, r19374, r19375, r19376, r19377, r19378, r19379, r19380, r19381, r19382, r19383, r19384, r19385, r19386, r19387, r19388, r19389, r19390, r19391, r19392, r19393, r19394, r19395, r19396, r19397, r19398, r19399, r19400, r19401, r19402, r19403, r19404, r19405, r19406, r19407, r19408, r19409, r19410, r19411, r19412, r19413, r19414; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19350); mpfr_init(r19351); mpfr_init_set_str(r19352, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19353); mpfr_init(r19354); mpfr_init(r19355); mpfr_init(r19356); mpfr_init_set_str(r19357, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19358); mpfr_init(r19359); mpfr_init(r19360); mpfr_init_set_str(r19361, "1", 10, MPFR_RNDN); mpfr_init(r19362); mpfr_init_set_str(r19363, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19364); mpfr_init(r19365); mpfr_init(r19366); mpfr_init(r19367); mpfr_init_set_str(r19368, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19369); mpfr_init(r19370); mpfr_init(r19371); mpfr_init_set_str(r19372, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19373); mpfr_init(r19374); mpfr_init(r19375); mpfr_init(r19376); mpfr_init_set_str(r19377, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19378, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19379); mpfr_init(r19380); mpfr_init(r19381); mpfr_init(r19382); mpfr_init(r19383); mpfr_init_set_str(r19384, "2", 10, MPFR_RNDN); mpfr_init(r19385); mpfr_init(r19386); mpfr_init_set_str(r19387, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19388); mpfr_init(r19389); mpfr_init(r19390); mpfr_init(r19391); mpfr_init(r19392); mpfr_init(r19393); mpfr_init_set_str(r19394, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19395); mpfr_init_set_str(r19396, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19397); mpfr_init(r19398); mpfr_init(r19399); mpfr_init(r19400); mpfr_init(r19401); mpfr_init(r19402); mpfr_init(r19403); mpfr_init(r19404); mpfr_init(r19405); mpfr_init(r19406); mpfr_init(r19407); mpfr_init(r19408); mpfr_init(r19409); mpfr_init(r19410); mpfr_init(r19411); mpfr_init(r19412); mpfr_init(r19413); mpfr_init(r19414); } double f_dm(double x) { mpfr_set_d(r19350, x, MPFR_RNDN); mpfr_sqr(r19351, r19350, MPFR_RNDN); ; mpfr_mul(r19353, r19351, r19352, MPFR_RNDN); mpfr_mul(r19354, r19351, r19351, MPFR_RNDN); mpfr_mul(r19355, r19354, r19354, MPFR_RNDN); mpfr_mul(r19356, r19353, r19355, MPFR_RNDN); ; mpfr_mul(r19358, r19350, r19357, MPFR_RNDN); mpfr_mul(r19359, r19350, r19350, MPFR_RNDN); mpfr_mul(r19359, r19359, r19350, MPFR_RNDN); mpfr_mul(r19360, r19358, r19359, MPFR_RNDN); ; mpfr_add(r19362, r19360, r19361, MPFR_RNDN); ; mpfr_mul(r19364, r19363, r19351, MPFR_RNDN); mpfr_add(r19365, r19362, r19364, MPFR_RNDN); mpfr_add(r19366, r19356, r19365, MPFR_RNDN); mpfr_mul(r19367, r19350, r19366, MPFR_RNDN); ; mpfr_mul(r19369, r19350, r19368, MPFR_RNDN); mpfr_mul(r19370, r19369, r19350, MPFR_RNDN); mpfr_add(r19371, r19361, r19370, MPFR_RNDN); ; mpfr_mul(r19373, r19350, r19372, MPFR_RNDN); mpfr_mul(r19374, r19359, r19373, MPFR_RNDN); mpfr_add(r19375, r19371, r19374, MPFR_RNDN); mpfr_mul(r19376, r19359, r19359, MPFR_RNDN); ; ; mpfr_mul(r19379, r19378, r19350, MPFR_RNDN); mpfr_mul(r19380, r19350, r19379, MPFR_RNDN); mpfr_add(r19381, r19377, r19380, MPFR_RNDN); mpfr_mul(r19382, r19376, r19381, MPFR_RNDN); mpfr_mul(r19383, r19376, r19354, MPFR_RNDN); ; mpfr_mul(r19385, r19384, r19352, MPFR_RNDN); mpfr_mul(r19386, r19385, r19351, MPFR_RNDN); ; mpfr_add(r19388, r19386, r19387, MPFR_RNDN); mpfr_mul(r19389, r19383, r19388, MPFR_RNDN); mpfr_add(r19390, r19382, r19389, MPFR_RNDN); mpfr_add(r19391, r19375, r19390, MPFR_RNDN); mpfr_div(r19392, r19367, r19391, MPFR_RNDN); mpfr_mul(r19393, r19350, r19376, MPFR_RNDN); ; mpfr_mul(r19395, r19394, r19351, MPFR_RNDN); ; mpfr_add(r19397, r19395, r19396, MPFR_RNDN); mpfr_mul(r19398, r19393, r19397, MPFR_RNDN); mpfr_mul(r19399, r19351, r19351, MPFR_RNDN); mpfr_mul(r19399, r19399, r19351, MPFR_RNDN); mpfr_mul(r19400, r19354, r19399, MPFR_RNDN); mpfr_mul(r19401, r19351, r19385, MPFR_RNDN); mpfr_add(r19402, r19401, r19387, MPFR_RNDN); mpfr_mul(r19403, r19400, r19402, MPFR_RNDN); mpfr_mul(r19404, r19368, r19350, MPFR_RNDN); mpfr_mul(r19405, r19350, r19404, MPFR_RNDN); mpfr_mul(r19406, r19373, r19359, MPFR_RNDN); mpfr_add(r19407, r19406, r19361, MPFR_RNDN); mpfr_add(r19408, r19405, r19407, MPFR_RNDN); mpfr_add(r19409, r19403, r19408, MPFR_RNDN); mpfr_add(r19410, r19380, r19377, MPFR_RNDN); mpfr_mul(r19411, r19399, r19410, MPFR_RNDN); mpfr_add(r19412, r19409, r19411, MPFR_RNDN); mpfr_div(r19413, r19398, r19412, MPFR_RNDN); mpfr_add(r19414, r19392, r19413, MPFR_RNDN); return mpfr_get_d(r19414, MPFR_RNDN); }