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