#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r19147 = 1.0f; float r19148 = 0.1049934947f; float r19149 = x; float r19150 = r19149 * r19149; float r19151 = r19148 * r19150; float r19152 = r19147 + r19151; float r19153 = 0.0424060604f; float r19154 = r19150 * r19150; float r19155 = r19153 * r19154; float r19156 = r19152 + r19155; float r19157 = 0.0072644182f; float r19158 = r19154 * r19150; float r19159 = r19157 * r19158; float r19160 = r19156 + r19159; float r19161 = 0.0005064034f; float r19162 = r19158 * r19150; float r19163 = r19161 * r19162; float r19164 = r19160 + r19163; float r19165 = 0.0001789971f; float r19166 = r19162 * r19150; float r19167 = r19165 * r19166; float r19168 = r19164 + r19167; float r19169 = 0.7715471019f; float r19170 = r19169 * r19150; float r19171 = r19147 + r19170; float r19172 = 0.2909738639f; float r19173 = r19172 * r19154; float r19174 = r19171 + r19173; float r19175 = 0.0694555761f; float r19176 = r19175 * r19158; float r19177 = r19174 + r19176; float r19178 = 0.0140005442f; float r19179 = r19178 * r19162; float r19180 = r19177 + r19179; float r19181 = 0.0008327945f; float r19182 = r19181 * r19166; float r19183 = r19180 + r19182; float r19184 = 2.0f; float r19185 = r19184 * r19165; float r19186 = r19166 * r19150; float r19187 = r19185 * r19186; float r19188 = r19183 + r19187; float r19189 = r19168 / r19188; float r19190 = r19189 * r19149; return r19190; } double f_id(double x) { double r19191 = 1.0; double r19192 = 0.1049934947; double r19193 = x; double r19194 = r19193 * r19193; double r19195 = r19192 * r19194; double r19196 = r19191 + r19195; double r19197 = 0.0424060604; double r19198 = r19194 * r19194; double r19199 = r19197 * r19198; double r19200 = r19196 + r19199; double r19201 = 0.0072644182; double r19202 = r19198 * r19194; double r19203 = r19201 * r19202; double r19204 = r19200 + r19203; double r19205 = 0.0005064034; double r19206 = r19202 * r19194; double r19207 = r19205 * r19206; double r19208 = r19204 + r19207; double r19209 = 0.0001789971; double r19210 = r19206 * r19194; double r19211 = r19209 * r19210; double r19212 = r19208 + r19211; double r19213 = 0.7715471019; double r19214 = r19213 * r19194; double r19215 = r19191 + r19214; double r19216 = 0.2909738639; double r19217 = r19216 * r19198; double r19218 = r19215 + r19217; double r19219 = 0.0694555761; double r19220 = r19219 * r19202; double r19221 = r19218 + r19220; double r19222 = 0.0140005442; double r19223 = r19222 * r19206; double r19224 = r19221 + r19223; double r19225 = 0.0008327945; double r19226 = r19225 * r19210; double r19227 = r19224 + r19226; double r19228 = 2.0; double r19229 = r19228 * r19209; double r19230 = r19210 * r19194; double r19231 = r19229 * r19230; double r19232 = r19227 + r19231; double r19233 = r19212 / r19232; double r19234 = r19233 * r19193; return r19234; } double f_of(float x) { float r19235 = x; float r19236 = r19235 * (r19235 * r19235); float r19237 = r19236 * r19236; float r19238 = 0.0694555761f; float r19239 = 0.0140005442f; float r19240 = r19239 * r19235; float r19241 = r19235 * r19240; float r19242 = r19238 + r19241; float r19243 = r19237 * r19242; float r19244 = r19235 * r19235; float r19245 = r19236 * r19244; float r19246 = r19245 * r19245; float r19247 = 0.0001789971f; float r19248 = r19247 * r19235; float r19249 = 2.0f; float r19250 = r19235 * r19249; float r19251 = r19248 * r19250; float r19252 = 0.0008327945f; float r19253 = r19251 + r19252; float r19254 = r19246 * r19253; float r19255 = r19243 + r19254; float r19256 = 0.7715471019f; float r19257 = r19235 * r19256; float r19258 = r19257 * r19235; float r19259 = 0.2909738639f; float r19260 = r19235 * r19259; float r19261 = r19236 * r19260; float r19262 = r19258 + r19261; float r19263 = 1.0f; float r19264 = r19262 + r19263; float r19265 = r19255 + r19264; float r19266 = r19235 / r19265; float r19267 = 0.1049934947f; float r19268 = r19235 * r19267; float r19269 = r19235 * r19268; float r19270 = 0.0424060604f; float r19271 = r19270 * r19235; float r19272 = r19271 * r19236; float r19273 = r19272 + r19263; float r19274 = r19269 + r19273; float r19275 = r19244 * r19244; float r19276 = r19275 * r19275; float r19277 = r19244 * r19247; float r19278 = 0.0005064034f; float r19279 = r19277 + r19278; float r19280 = r19276 * r19279; float r19281 = 0.0072644182f; float r19282 = r19281 * r19237; float r19283 = r19280 + r19282; float r19284 = r19274 + r19283; float r19285 = r19266 * r19284; return r19285; } double f_od(double x) { double r19286 = x; double r19287 = r19286 * (r19286 * r19286); double r19288 = r19287 * r19287; double r19289 = 0.0694555761; double r19290 = 0.0140005442; double r19291 = r19290 * r19286; double r19292 = r19286 * r19291; double r19293 = r19289 + r19292; double r19294 = r19288 * r19293; double r19295 = r19286 * r19286; double r19296 = r19287 * r19295; double r19297 = r19296 * r19296; double r19298 = 0.0001789971; double r19299 = r19298 * r19286; double r19300 = 2.0; double r19301 = r19286 * r19300; double r19302 = r19299 * r19301; double r19303 = 0.0008327945; double r19304 = r19302 + r19303; double r19305 = r19297 * r19304; double r19306 = r19294 + r19305; double r19307 = 0.7715471019; double r19308 = r19286 * r19307; double r19309 = r19308 * r19286; double r19310 = 0.2909738639; double r19311 = r19286 * r19310; double r19312 = r19287 * r19311; double r19313 = r19309 + r19312; double r19314 = 1.0; double r19315 = r19313 + r19314; double r19316 = r19306 + r19315; double r19317 = r19286 / r19316; double r19318 = 0.1049934947; double r19319 = r19286 * r19318; double r19320 = r19286 * r19319; double r19321 = 0.0424060604; double r19322 = r19321 * r19286; double r19323 = r19322 * r19287; double r19324 = r19323 + r19314; double r19325 = r19320 + r19324; double r19326 = r19295 * r19295; double r19327 = r19326 * r19326; double r19328 = r19295 * r19298; double r19329 = 0.0005064034; double r19330 = r19328 + r19329; double r19331 = r19327 * r19330; double r19332 = 0.0072644182; double r19333 = r19332 * r19288; double r19334 = r19331 + r19333; double r19335 = r19325 + r19334; double r19336 = r19317 * r19335; return r19336; } 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 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, r19378, r19379, r19380; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19337, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19338, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19339); mpfr_init(r19340); mpfr_init(r19341); mpfr_init(r19342); mpfr_init_set_str(r19343, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19344); mpfr_init(r19345); mpfr_init(r19346); mpfr_init_set_str(r19347, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19348); mpfr_init(r19349); mpfr_init(r19350); mpfr_init_set_str(r19351, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19352); mpfr_init(r19353); mpfr_init(r19354); mpfr_init_set_str(r19355, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19356); mpfr_init(r19357); mpfr_init(r19358); mpfr_init_set_str(r19359, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19360); mpfr_init(r19361); mpfr_init_set_str(r19362, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19363); mpfr_init(r19364); mpfr_init_set_str(r19365, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19366); mpfr_init(r19367); mpfr_init_set_str(r19368, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19369); mpfr_init(r19370); mpfr_init_set_str(r19371, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19372); mpfr_init(r19373); mpfr_init_set_str(r19374, "2", 10, MPFR_RNDN); mpfr_init(r19375); mpfr_init(r19376); mpfr_init(r19377); mpfr_init(r19378); mpfr_init(r19379); mpfr_init(r19380); } double f_im(double x) { ; ; mpfr_set_d(r19339, x, MPFR_RNDN); mpfr_mul(r19340, r19339, r19339, MPFR_RNDN); mpfr_mul(r19341, r19338, r19340, MPFR_RNDN); mpfr_add(r19342, r19337, r19341, MPFR_RNDN); ; mpfr_mul(r19344, r19340, r19340, MPFR_RNDN); mpfr_mul(r19345, r19343, r19344, MPFR_RNDN); mpfr_add(r19346, r19342, r19345, MPFR_RNDN); ; mpfr_mul(r19348, r19344, r19340, MPFR_RNDN); mpfr_mul(r19349, r19347, r19348, MPFR_RNDN); mpfr_add(r19350, r19346, r19349, MPFR_RNDN); ; mpfr_mul(r19352, r19348, r19340, MPFR_RNDN); mpfr_mul(r19353, r19351, r19352, MPFR_RNDN); mpfr_add(r19354, r19350, r19353, MPFR_RNDN); ; mpfr_mul(r19356, r19352, r19340, MPFR_RNDN); mpfr_mul(r19357, r19355, r19356, MPFR_RNDN); mpfr_add(r19358, r19354, r19357, MPFR_RNDN); ; mpfr_mul(r19360, r19359, r19340, MPFR_RNDN); mpfr_add(r19361, r19337, r19360, MPFR_RNDN); ; mpfr_mul(r19363, r19362, r19344, MPFR_RNDN); mpfr_add(r19364, r19361, r19363, MPFR_RNDN); ; mpfr_mul(r19366, r19365, r19348, MPFR_RNDN); mpfr_add(r19367, r19364, r19366, MPFR_RNDN); ; mpfr_mul(r19369, r19368, r19352, MPFR_RNDN); mpfr_add(r19370, r19367, r19369, MPFR_RNDN); ; mpfr_mul(r19372, r19371, r19356, MPFR_RNDN); mpfr_add(r19373, r19370, r19372, MPFR_RNDN); ; mpfr_mul(r19375, r19374, r19355, MPFR_RNDN); mpfr_mul(r19376, r19356, r19340, MPFR_RNDN); mpfr_mul(r19377, r19375, r19376, MPFR_RNDN); mpfr_add(r19378, r19373, r19377, MPFR_RNDN); mpfr_div(r19379, r19358, r19378, MPFR_RNDN); mpfr_mul(r19380, r19379, r19339, MPFR_RNDN); return mpfr_get_d(r19380, MPFR_RNDN); } static mpfr_t 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, r19431; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19381); mpfr_init(r19382); mpfr_init(r19383); mpfr_init_set_str(r19384, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19385, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19386); mpfr_init(r19387); mpfr_init(r19388); mpfr_init(r19389); mpfr_init(r19390); mpfr_init(r19391); mpfr_init(r19392); mpfr_init_set_str(r19393, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19394); mpfr_init_set_str(r19395, "2", 10, MPFR_RNDN); mpfr_init(r19396); mpfr_init(r19397); mpfr_init_set_str(r19398, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19399); mpfr_init(r19400); mpfr_init(r19401); mpfr_init_set_str(r19402, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19403); mpfr_init(r19404); mpfr_init_set_str(r19405, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19406); mpfr_init(r19407); mpfr_init(r19408); mpfr_init_set_str(r19409, "1", 10, MPFR_RNDN); mpfr_init(r19410); mpfr_init(r19411); mpfr_init(r19412); mpfr_init_set_str(r19413, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19414); mpfr_init(r19415); mpfr_init_set_str(r19416, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19417); mpfr_init(r19418); mpfr_init(r19419); mpfr_init(r19420); mpfr_init(r19421); mpfr_init(r19422); mpfr_init(r19423); mpfr_init_set_str(r19424, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19425); mpfr_init(r19426); mpfr_init_set_str(r19427, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19428); mpfr_init(r19429); mpfr_init(r19430); mpfr_init(r19431); } double f_fm(double x) { mpfr_set_d(r19381, x, MPFR_RNDN); mpfr_mul(r19382, r19381, r19381, MPFR_RNDN); mpfr_mul(r19382, r19382, r19381, MPFR_RNDN); mpfr_mul(r19383, r19382, r19382, MPFR_RNDN); ; ; mpfr_mul(r19386, r19385, r19381, MPFR_RNDN); mpfr_mul(r19387, r19381, r19386, MPFR_RNDN); mpfr_add(r19388, r19384, r19387, MPFR_RNDN); mpfr_mul(r19389, r19383, r19388, MPFR_RNDN); mpfr_sqr(r19390, r19381, MPFR_RNDN); mpfr_mul(r19391, r19382, r19390, MPFR_RNDN); mpfr_mul(r19392, r19391, r19391, MPFR_RNDN); ; mpfr_mul(r19394, r19393, r19381, MPFR_RNDN); ; mpfr_mul(r19396, r19381, r19395, MPFR_RNDN); mpfr_mul(r19397, r19394, r19396, MPFR_RNDN); ; mpfr_add(r19399, r19397, r19398, MPFR_RNDN); mpfr_mul(r19400, r19392, r19399, MPFR_RNDN); mpfr_add(r19401, r19389, r19400, MPFR_RNDN); ; mpfr_mul(r19403, r19381, r19402, MPFR_RNDN); mpfr_mul(r19404, r19403, r19381, MPFR_RNDN); ; mpfr_mul(r19406, r19381, r19405, MPFR_RNDN); mpfr_mul(r19407, r19382, r19406, MPFR_RNDN); mpfr_add(r19408, r19404, r19407, MPFR_RNDN); ; mpfr_add(r19410, r19408, r19409, MPFR_RNDN); mpfr_add(r19411, r19401, r19410, MPFR_RNDN); mpfr_div(r19412, r19381, r19411, MPFR_RNDN); ; mpfr_mul(r19414, r19381, r19413, MPFR_RNDN); mpfr_mul(r19415, r19381, r19414, MPFR_RNDN); ; mpfr_mul(r19417, r19416, r19381, MPFR_RNDN); mpfr_mul(r19418, r19417, r19382, MPFR_RNDN); mpfr_add(r19419, r19418, r19409, MPFR_RNDN); mpfr_add(r19420, r19415, r19419, MPFR_RNDN); mpfr_mul(r19421, r19390, r19390, MPFR_RNDN); mpfr_sqr(r19422, r19421, MPFR_RNDN); mpfr_mul(r19423, r19390, r19393, MPFR_RNDN); ; mpfr_add(r19425, r19423, r19424, MPFR_RNDN); mpfr_mul(r19426, r19422, r19425, MPFR_RNDN); ; mpfr_mul(r19428, r19427, r19383, MPFR_RNDN); mpfr_add(r19429, r19426, r19428, MPFR_RNDN); mpfr_add(r19430, r19420, r19429, MPFR_RNDN); mpfr_mul(r19431, r19412, r19430, MPFR_RNDN); return mpfr_get_d(r19431, MPFR_RNDN); } static mpfr_t r19432, r19433, r19434, r19435, r19436, r19437, r19438, r19439, r19440, r19441, r19442, r19443, r19444, r19445, r19446, r19447, r19448, r19449, r19450, r19451, r19452, r19453, r19454, r19455, r19456, r19457, r19458, r19459, r19460, r19461, r19462, r19463, r19464, r19465, r19466, r19467, r19468, r19469, r19470, r19471, r19472, r19473, r19474, r19475, r19476, r19477, r19478, r19479, r19480, r19481, r19482; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19432); mpfr_init(r19433); mpfr_init(r19434); mpfr_init_set_str(r19435, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19436, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19437); mpfr_init(r19438); mpfr_init(r19439); mpfr_init(r19440); mpfr_init(r19441); mpfr_init(r19442); mpfr_init(r19443); mpfr_init_set_str(r19444, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19445); mpfr_init_set_str(r19446, "2", 10, MPFR_RNDN); mpfr_init(r19447); mpfr_init(r19448); mpfr_init_set_str(r19449, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19450); mpfr_init(r19451); mpfr_init(r19452); mpfr_init_set_str(r19453, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19454); mpfr_init(r19455); mpfr_init_set_str(r19456, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19457); mpfr_init(r19458); mpfr_init(r19459); mpfr_init_set_str(r19460, "1", 10, MPFR_RNDN); mpfr_init(r19461); mpfr_init(r19462); mpfr_init(r19463); mpfr_init_set_str(r19464, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19465); mpfr_init(r19466); mpfr_init_set_str(r19467, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19468); mpfr_init(r19469); mpfr_init(r19470); mpfr_init(r19471); mpfr_init(r19472); mpfr_init(r19473); mpfr_init(r19474); mpfr_init_set_str(r19475, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19476); mpfr_init(r19477); mpfr_init_set_str(r19478, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19479); mpfr_init(r19480); mpfr_init(r19481); mpfr_init(r19482); } double f_dm(double x) { mpfr_set_d(r19432, x, MPFR_RNDN); mpfr_mul(r19433, r19432, r19432, MPFR_RNDN); mpfr_mul(r19433, r19433, r19432, MPFR_RNDN); mpfr_mul(r19434, r19433, r19433, MPFR_RNDN); ; ; mpfr_mul(r19437, r19436, r19432, MPFR_RNDN); mpfr_mul(r19438, r19432, r19437, MPFR_RNDN); mpfr_add(r19439, r19435, r19438, MPFR_RNDN); mpfr_mul(r19440, r19434, r19439, MPFR_RNDN); mpfr_sqr(r19441, r19432, MPFR_RNDN); mpfr_mul(r19442, r19433, r19441, MPFR_RNDN); mpfr_mul(r19443, r19442, r19442, MPFR_RNDN); ; mpfr_mul(r19445, r19444, r19432, MPFR_RNDN); ; mpfr_mul(r19447, r19432, r19446, MPFR_RNDN); mpfr_mul(r19448, r19445, r19447, MPFR_RNDN); ; mpfr_add(r19450, r19448, r19449, MPFR_RNDN); mpfr_mul(r19451, r19443, r19450, MPFR_RNDN); mpfr_add(r19452, r19440, r19451, MPFR_RNDN); ; mpfr_mul(r19454, r19432, r19453, MPFR_RNDN); mpfr_mul(r19455, r19454, r19432, MPFR_RNDN); ; mpfr_mul(r19457, r19432, r19456, MPFR_RNDN); mpfr_mul(r19458, r19433, r19457, MPFR_RNDN); mpfr_add(r19459, r19455, r19458, MPFR_RNDN); ; mpfr_add(r19461, r19459, r19460, MPFR_RNDN); mpfr_add(r19462, r19452, r19461, MPFR_RNDN); mpfr_div(r19463, r19432, r19462, MPFR_RNDN); ; mpfr_mul(r19465, r19432, r19464, MPFR_RNDN); mpfr_mul(r19466, r19432, r19465, MPFR_RNDN); ; mpfr_mul(r19468, r19467, r19432, MPFR_RNDN); mpfr_mul(r19469, r19468, r19433, MPFR_RNDN); mpfr_add(r19470, r19469, r19460, MPFR_RNDN); mpfr_add(r19471, r19466, r19470, MPFR_RNDN); mpfr_mul(r19472, r19441, r19441, MPFR_RNDN); mpfr_sqr(r19473, r19472, MPFR_RNDN); mpfr_mul(r19474, r19441, r19444, MPFR_RNDN); ; mpfr_add(r19476, r19474, r19475, MPFR_RNDN); mpfr_mul(r19477, r19473, r19476, MPFR_RNDN); ; mpfr_mul(r19479, r19478, r19434, MPFR_RNDN); mpfr_add(r19480, r19477, r19479, MPFR_RNDN); mpfr_add(r19481, r19471, r19480, MPFR_RNDN); mpfr_mul(r19482, r19463, r19481, MPFR_RNDN); return mpfr_get_d(r19482, MPFR_RNDN); }