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