#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r19215 = 1.0f; float r19216 = 0.1049934947f; float r19217 = x; float r19218 = r19217 * r19217; float r19219 = r19216 * r19218; float r19220 = r19215 + r19219; float r19221 = 0.0424060604f; float r19222 = r19218 * r19218; float r19223 = r19221 * r19222; float r19224 = r19220 + r19223; float r19225 = 0.0072644182f; float r19226 = r19222 * r19218; float r19227 = r19225 * r19226; float r19228 = r19224 + r19227; float r19229 = 0.0005064034f; float r19230 = r19226 * r19218; float r19231 = r19229 * r19230; float r19232 = r19228 + r19231; float r19233 = 0.0001789971f; float r19234 = r19230 * r19218; float r19235 = r19233 * r19234; float r19236 = r19232 + r19235; float r19237 = 0.7715471019f; float r19238 = r19237 * r19218; float r19239 = r19215 + r19238; float r19240 = 0.2909738639f; float r19241 = r19240 * r19222; float r19242 = r19239 + r19241; float r19243 = 0.0694555761f; float r19244 = r19243 * r19226; float r19245 = r19242 + r19244; float r19246 = 0.0140005442f; float r19247 = r19246 * r19230; float r19248 = r19245 + r19247; float r19249 = 0.0008327945f; float r19250 = r19249 * r19234; float r19251 = r19248 + r19250; float r19252 = 2.0f; float r19253 = r19252 * r19233; float r19254 = r19234 * r19218; float r19255 = r19253 * r19254; float r19256 = r19251 + r19255; float r19257 = r19236 / r19256; float r19258 = r19257 * r19217; return r19258; } double f_id(double x) { double r19259 = 1.0; double r19260 = 0.1049934947; double r19261 = x; double r19262 = r19261 * r19261; double r19263 = r19260 * r19262; double r19264 = r19259 + r19263; double r19265 = 0.0424060604; double r19266 = r19262 * r19262; double r19267 = r19265 * r19266; double r19268 = r19264 + r19267; double r19269 = 0.0072644182; double r19270 = r19266 * r19262; double r19271 = r19269 * r19270; double r19272 = r19268 + r19271; double r19273 = 0.0005064034; double r19274 = r19270 * r19262; double r19275 = r19273 * r19274; double r19276 = r19272 + r19275; double r19277 = 0.0001789971; double r19278 = r19274 * r19262; double r19279 = r19277 * r19278; double r19280 = r19276 + r19279; double r19281 = 0.7715471019; double r19282 = r19281 * r19262; double r19283 = r19259 + r19282; double r19284 = 0.2909738639; double r19285 = r19284 * r19266; double r19286 = r19283 + r19285; double r19287 = 0.0694555761; double r19288 = r19287 * r19270; double r19289 = r19286 + r19288; double r19290 = 0.0140005442; double r19291 = r19290 * r19274; double r19292 = r19289 + r19291; double r19293 = 0.0008327945; double r19294 = r19293 * r19278; double r19295 = r19292 + r19294; double r19296 = 2.0; double r19297 = r19296 * r19277; double r19298 = r19278 * r19262; double r19299 = r19297 * r19298; double r19300 = r19295 + r19299; double r19301 = r19280 / r19300; double r19302 = r19301 * r19261; return r19302; } double f_of(float x) { float r19303 = x; float r19304 = r19303 * (r19303 * r19303); float r19305 = r19304 * r19304; float r19306 = 0.0694555761f; float r19307 = 0.0140005442f; float r19308 = r19307 * r19303; float r19309 = r19303 * r19308; float r19310 = r19306 + r19309; float r19311 = r19305 * r19310; float r19312 = r19303 * r19303; float r19313 = r19304 * r19312; float r19314 = r19313 * r19313; float r19315 = 0.0001789971f; float r19316 = r19315 * r19303; float r19317 = 2.0f; float r19318 = r19303 * r19317; float r19319 = r19316 * r19318; float r19320 = 0.0008327945f; float r19321 = r19319 + r19320; float r19322 = r19314 * r19321; float r19323 = r19311 + r19322; float r19324 = 0.7715471019f; float r19325 = r19303 * r19324; float r19326 = r19325 * r19303; float r19327 = 0.2909738639f; float r19328 = r19303 * r19327; float r19329 = r19304 * r19328; float r19330 = r19326 + r19329; float r19331 = 1.0f; float r19332 = r19330 + r19331; float r19333 = r19323 + r19332; float r19334 = r19303 / r19333; float r19335 = 0.1049934947f; float r19336 = r19303 * r19335; float r19337 = r19303 * r19336; float r19338 = 0.0424060604f; float r19339 = r19338 * r19303; float r19340 = r19339 * r19304; float r19341 = r19340 + r19331; float r19342 = r19337 + r19341; float r19343 = r19312 * r19312; float r19344 = r19343 * r19343; float r19345 = r19312 * r19315; float r19346 = 0.0005064034f; float r19347 = r19345 + r19346; float r19348 = r19344 * r19347; float r19349 = 0.0072644182f; float r19350 = r19349 * r19305; float r19351 = r19348 + r19350; float r19352 = r19342 + r19351; float r19353 = r19334 * r19352; return r19353; } double f_od(double x) { double r19354 = x; double r19355 = r19354 * (r19354 * r19354); double r19356 = r19355 * r19355; double r19357 = 0.0694555761; double r19358 = 0.0140005442; double r19359 = r19358 * r19354; double r19360 = r19354 * r19359; double r19361 = r19357 + r19360; double r19362 = r19356 * r19361; double r19363 = r19354 * r19354; double r19364 = r19355 * r19363; double r19365 = r19364 * r19364; double r19366 = 0.0001789971; double r19367 = r19366 * r19354; double r19368 = 2.0; double r19369 = r19354 * r19368; double r19370 = r19367 * r19369; double r19371 = 0.0008327945; double r19372 = r19370 + r19371; double r19373 = r19365 * r19372; double r19374 = r19362 + r19373; double r19375 = 0.7715471019; double r19376 = r19354 * r19375; double r19377 = r19376 * r19354; double r19378 = 0.2909738639; double r19379 = r19354 * r19378; double r19380 = r19355 * r19379; double r19381 = r19377 + r19380; double r19382 = 1.0; double r19383 = r19381 + r19382; double r19384 = r19374 + r19383; double r19385 = r19354 / r19384; double r19386 = 0.1049934947; double r19387 = r19354 * r19386; double r19388 = r19354 * r19387; double r19389 = 0.0424060604; double r19390 = r19389 * r19354; double r19391 = r19390 * r19355; double r19392 = r19391 + r19382; double r19393 = r19388 + r19392; double r19394 = r19363 * r19363; double r19395 = r19394 * r19394; double r19396 = r19363 * r19366; double r19397 = 0.0005064034; double r19398 = r19396 + r19397; double r19399 = r19395 * r19398; double r19400 = 0.0072644182; double r19401 = r19400 * r19356; double r19402 = r19399 + r19401; double r19403 = r19393 + r19402; double r19404 = r19385 * r19403; return r19404; } 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 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, r19432, r19433, r19434, r19435, r19436, r19437, r19438, r19439, r19440, r19441, r19442, r19443, r19444, r19445, r19446, r19447, r19448; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19405, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19406, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19407); mpfr_init(r19408); mpfr_init(r19409); mpfr_init(r19410); mpfr_init_set_str(r19411, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19412); mpfr_init(r19413); mpfr_init(r19414); mpfr_init_set_str(r19415, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19416); mpfr_init(r19417); mpfr_init(r19418); mpfr_init_set_str(r19419, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19420); mpfr_init(r19421); mpfr_init(r19422); mpfr_init_set_str(r19423, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19424); mpfr_init(r19425); mpfr_init(r19426); mpfr_init_set_str(r19427, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19428); mpfr_init(r19429); mpfr_init_set_str(r19430, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19431); mpfr_init(r19432); mpfr_init_set_str(r19433, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19434); mpfr_init(r19435); mpfr_init_set_str(r19436, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19437); mpfr_init(r19438); mpfr_init_set_str(r19439, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19440); mpfr_init(r19441); mpfr_init_set_str(r19442, "2", 10, MPFR_RNDN); mpfr_init(r19443); mpfr_init(r19444); mpfr_init(r19445); mpfr_init(r19446); mpfr_init(r19447); mpfr_init(r19448); } double f_im(double x) { ; ; mpfr_set_d(r19407, x, MPFR_RNDN); mpfr_mul(r19408, r19407, r19407, MPFR_RNDN); mpfr_mul(r19409, r19406, r19408, MPFR_RNDN); mpfr_add(r19410, r19405, r19409, MPFR_RNDN); ; mpfr_mul(r19412, r19408, r19408, MPFR_RNDN); mpfr_mul(r19413, r19411, r19412, MPFR_RNDN); mpfr_add(r19414, r19410, r19413, MPFR_RNDN); ; mpfr_mul(r19416, r19412, r19408, MPFR_RNDN); mpfr_mul(r19417, r19415, r19416, MPFR_RNDN); mpfr_add(r19418, r19414, r19417, MPFR_RNDN); ; mpfr_mul(r19420, r19416, r19408, MPFR_RNDN); mpfr_mul(r19421, r19419, r19420, MPFR_RNDN); mpfr_add(r19422, r19418, r19421, MPFR_RNDN); ; mpfr_mul(r19424, r19420, r19408, MPFR_RNDN); mpfr_mul(r19425, r19423, r19424, MPFR_RNDN); mpfr_add(r19426, r19422, r19425, MPFR_RNDN); ; mpfr_mul(r19428, r19427, r19408, MPFR_RNDN); mpfr_add(r19429, r19405, r19428, MPFR_RNDN); ; mpfr_mul(r19431, r19430, r19412, MPFR_RNDN); mpfr_add(r19432, r19429, r19431, MPFR_RNDN); ; mpfr_mul(r19434, r19433, r19416, MPFR_RNDN); mpfr_add(r19435, r19432, r19434, MPFR_RNDN); ; mpfr_mul(r19437, r19436, r19420, MPFR_RNDN); mpfr_add(r19438, r19435, r19437, MPFR_RNDN); ; mpfr_mul(r19440, r19439, r19424, MPFR_RNDN); mpfr_add(r19441, r19438, r19440, MPFR_RNDN); ; mpfr_mul(r19443, r19442, r19423, MPFR_RNDN); mpfr_mul(r19444, r19424, r19408, MPFR_RNDN); mpfr_mul(r19445, r19443, r19444, MPFR_RNDN); mpfr_add(r19446, r19441, r19445, MPFR_RNDN); mpfr_div(r19447, r19426, r19446, MPFR_RNDN); mpfr_mul(r19448, r19447, r19407, MPFR_RNDN); return mpfr_get_d(r19448, MPFR_RNDN); } static mpfr_t 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, r19483, r19484, r19485, r19486, r19487, r19488, r19489, r19490, r19491, r19492, r19493, r19494, r19495, r19496, r19497, r19498, r19499; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19449); mpfr_init(r19450); mpfr_init(r19451); mpfr_init_set_str(r19452, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19453, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19454); mpfr_init(r19455); mpfr_init(r19456); mpfr_init(r19457); mpfr_init(r19458); mpfr_init(r19459); mpfr_init(r19460); mpfr_init_set_str(r19461, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19462); mpfr_init_set_str(r19463, "2", 10, MPFR_RNDN); mpfr_init(r19464); mpfr_init(r19465); mpfr_init_set_str(r19466, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19467); mpfr_init(r19468); mpfr_init(r19469); mpfr_init_set_str(r19470, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19471); mpfr_init(r19472); mpfr_init_set_str(r19473, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19474); mpfr_init(r19475); mpfr_init(r19476); mpfr_init_set_str(r19477, "1", 10, MPFR_RNDN); mpfr_init(r19478); mpfr_init(r19479); mpfr_init(r19480); mpfr_init_set_str(r19481, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19482); mpfr_init(r19483); mpfr_init_set_str(r19484, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19485); mpfr_init(r19486); mpfr_init(r19487); mpfr_init(r19488); mpfr_init(r19489); mpfr_init(r19490); mpfr_init(r19491); mpfr_init_set_str(r19492, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19493); mpfr_init(r19494); mpfr_init_set_str(r19495, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19496); mpfr_init(r19497); mpfr_init(r19498); mpfr_init(r19499); } double f_fm(double x) { mpfr_set_d(r19449, x, MPFR_RNDN); mpfr_mul(r19450, r19449, r19449, MPFR_RNDN); mpfr_mul(r19450, r19450, r19449, MPFR_RNDN); mpfr_mul(r19451, r19450, r19450, MPFR_RNDN); ; ; mpfr_mul(r19454, r19453, r19449, MPFR_RNDN); mpfr_mul(r19455, r19449, r19454, MPFR_RNDN); mpfr_add(r19456, r19452, r19455, MPFR_RNDN); mpfr_mul(r19457, r19451, r19456, MPFR_RNDN); mpfr_sqr(r19458, r19449, MPFR_RNDN); mpfr_mul(r19459, r19450, r19458, MPFR_RNDN); mpfr_mul(r19460, r19459, r19459, MPFR_RNDN); ; mpfr_mul(r19462, r19461, r19449, MPFR_RNDN); ; mpfr_mul(r19464, r19449, r19463, MPFR_RNDN); mpfr_mul(r19465, r19462, r19464, MPFR_RNDN); ; mpfr_add(r19467, r19465, r19466, MPFR_RNDN); mpfr_mul(r19468, r19460, r19467, MPFR_RNDN); mpfr_add(r19469, r19457, r19468, MPFR_RNDN); ; mpfr_mul(r19471, r19449, r19470, MPFR_RNDN); mpfr_mul(r19472, r19471, r19449, MPFR_RNDN); ; mpfr_mul(r19474, r19449, r19473, MPFR_RNDN); mpfr_mul(r19475, r19450, r19474, MPFR_RNDN); mpfr_add(r19476, r19472, r19475, MPFR_RNDN); ; mpfr_add(r19478, r19476, r19477, MPFR_RNDN); mpfr_add(r19479, r19469, r19478, MPFR_RNDN); mpfr_div(r19480, r19449, r19479, MPFR_RNDN); ; mpfr_mul(r19482, r19449, r19481, MPFR_RNDN); mpfr_mul(r19483, r19449, r19482, MPFR_RNDN); ; mpfr_mul(r19485, r19484, r19449, MPFR_RNDN); mpfr_mul(r19486, r19485, r19450, MPFR_RNDN); mpfr_add(r19487, r19486, r19477, MPFR_RNDN); mpfr_add(r19488, r19483, r19487, MPFR_RNDN); mpfr_mul(r19489, r19458, r19458, MPFR_RNDN); mpfr_sqr(r19490, r19489, MPFR_RNDN); mpfr_mul(r19491, r19458, r19461, MPFR_RNDN); ; mpfr_add(r19493, r19491, r19492, MPFR_RNDN); mpfr_mul(r19494, r19490, r19493, MPFR_RNDN); ; mpfr_mul(r19496, r19495, r19451, MPFR_RNDN); mpfr_add(r19497, r19494, r19496, MPFR_RNDN); mpfr_add(r19498, r19488, r19497, MPFR_RNDN); mpfr_mul(r19499, r19480, r19498, MPFR_RNDN); return mpfr_get_d(r19499, MPFR_RNDN); } static mpfr_t r19500, r19501, r19502, r19503, r19504, r19505, r19506, r19507, r19508, r19509, r19510, r19511, r19512, r19513, r19514, r19515, r19516, r19517, r19518, r19519, r19520, r19521, r19522, r19523, r19524, r19525, r19526, r19527, r19528, r19529, r19530, r19531, r19532, r19533, r19534, r19535, r19536, r19537, r19538, r19539, r19540, r19541, r19542, r19543, r19544, r19545, r19546, r19547, r19548, r19549, r19550; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19500); mpfr_init(r19501); mpfr_init(r19502); mpfr_init_set_str(r19503, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19504, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19505); mpfr_init(r19506); mpfr_init(r19507); mpfr_init(r19508); mpfr_init(r19509); mpfr_init(r19510); mpfr_init(r19511); mpfr_init_set_str(r19512, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19513); mpfr_init_set_str(r19514, "2", 10, MPFR_RNDN); mpfr_init(r19515); mpfr_init(r19516); mpfr_init_set_str(r19517, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19518); mpfr_init(r19519); mpfr_init(r19520); mpfr_init_set_str(r19521, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19522); mpfr_init(r19523); mpfr_init_set_str(r19524, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19525); mpfr_init(r19526); mpfr_init(r19527); mpfr_init_set_str(r19528, "1", 10, MPFR_RNDN); mpfr_init(r19529); mpfr_init(r19530); mpfr_init(r19531); mpfr_init_set_str(r19532, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19533); mpfr_init(r19534); mpfr_init_set_str(r19535, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19536); mpfr_init(r19537); mpfr_init(r19538); mpfr_init(r19539); mpfr_init(r19540); mpfr_init(r19541); mpfr_init(r19542); mpfr_init_set_str(r19543, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19544); mpfr_init(r19545); mpfr_init_set_str(r19546, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19547); mpfr_init(r19548); mpfr_init(r19549); mpfr_init(r19550); } double f_dm(double x) { mpfr_set_d(r19500, x, MPFR_RNDN); mpfr_mul(r19501, r19500, r19500, MPFR_RNDN); mpfr_mul(r19501, r19501, r19500, MPFR_RNDN); mpfr_mul(r19502, r19501, r19501, MPFR_RNDN); ; ; mpfr_mul(r19505, r19504, r19500, MPFR_RNDN); mpfr_mul(r19506, r19500, r19505, MPFR_RNDN); mpfr_add(r19507, r19503, r19506, MPFR_RNDN); mpfr_mul(r19508, r19502, r19507, MPFR_RNDN); mpfr_sqr(r19509, r19500, MPFR_RNDN); mpfr_mul(r19510, r19501, r19509, MPFR_RNDN); mpfr_mul(r19511, r19510, r19510, MPFR_RNDN); ; mpfr_mul(r19513, r19512, r19500, MPFR_RNDN); ; mpfr_mul(r19515, r19500, r19514, MPFR_RNDN); mpfr_mul(r19516, r19513, r19515, MPFR_RNDN); ; mpfr_add(r19518, r19516, r19517, MPFR_RNDN); mpfr_mul(r19519, r19511, r19518, MPFR_RNDN); mpfr_add(r19520, r19508, r19519, MPFR_RNDN); ; mpfr_mul(r19522, r19500, r19521, MPFR_RNDN); mpfr_mul(r19523, r19522, r19500, MPFR_RNDN); ; mpfr_mul(r19525, r19500, r19524, MPFR_RNDN); mpfr_mul(r19526, r19501, r19525, MPFR_RNDN); mpfr_add(r19527, r19523, r19526, MPFR_RNDN); ; mpfr_add(r19529, r19527, r19528, MPFR_RNDN); mpfr_add(r19530, r19520, r19529, MPFR_RNDN); mpfr_div(r19531, r19500, r19530, MPFR_RNDN); ; mpfr_mul(r19533, r19500, r19532, MPFR_RNDN); mpfr_mul(r19534, r19500, r19533, MPFR_RNDN); ; mpfr_mul(r19536, r19535, r19500, MPFR_RNDN); mpfr_mul(r19537, r19536, r19501, MPFR_RNDN); mpfr_add(r19538, r19537, r19528, MPFR_RNDN); mpfr_add(r19539, r19534, r19538, MPFR_RNDN); mpfr_mul(r19540, r19509, r19509, MPFR_RNDN); mpfr_sqr(r19541, r19540, MPFR_RNDN); mpfr_mul(r19542, r19509, r19512, MPFR_RNDN); ; mpfr_add(r19544, r19542, r19543, MPFR_RNDN); mpfr_mul(r19545, r19541, r19544, MPFR_RNDN); ; mpfr_mul(r19547, r19546, r19502, MPFR_RNDN); mpfr_add(r19548, r19545, r19547, MPFR_RNDN); mpfr_add(r19549, r19539, r19548, MPFR_RNDN); mpfr_mul(r19550, r19531, r19549, MPFR_RNDN); return mpfr_get_d(r19550, MPFR_RNDN); }