#include #include #include #include #include char *name = "Toniolo and Linder, Equation (13)"; double f_if(float n, float U, float t, float l, float Om, float U_) { float r19344 = 2.0f; float r19345 = n; float r19346 = r19344 * r19345; float r19347 = U; float r19348 = r19346 * r19347; float r19349 = t; float r19350 = l; float r19351 = r19350 * r19350; float r19352 = Om; float r19353 = r19351 / r19352; float r19354 = r19344 * r19353; float r19355 = r19349 - r19354; float r19356 = r19350 / r19352; float r19357 = r19356 * r19356; float r19358 = r19345 * r19357; float r19359 = U_; float r19360 = r19347 - r19359; float r19361 = r19358 * r19360; float r19362 = r19355 - r19361; float r19363 = r19348 * r19362; float r19364 = sqrt(r19363); return r19364; } double f_id(double n, double U, double t, double l, double Om, double U_) { double r19365 = 2.0; double r19366 = n; double r19367 = r19365 * r19366; double r19368 = U; double r19369 = r19367 * r19368; double r19370 = t; double r19371 = l; double r19372 = r19371 * r19371; double r19373 = Om; double r19374 = r19372 / r19373; double r19375 = r19365 * r19374; double r19376 = r19370 - r19375; double r19377 = r19371 / r19373; double r19378 = r19377 * r19377; double r19379 = r19366 * r19378; double r19380 = U_; double r19381 = r19368 - r19380; double r19382 = r19379 * r19381; double r19383 = r19376 - r19382; double r19384 = r19369 * r19383; double r19385 = sqrt(r19384); return r19385; } double f_of(float n, float U, float t, float l, float Om, float U_) { float r19386 = n; float r19387 = 4.5602004537906414e-307f; bool r19388 = r19386 <= r19387; float r19389 = 2.0f; float r19390 = r19389 * r19386; float r19391 = U; float r19392 = t; float r19393 = l; float r19394 = Om; float r19395 = r19394 / r19393; float r19396 = r19393 / r19395; float r19397 = r19389 * r19396; float r19398 = r19392 - r19397; float r19399 = r19393 / r19394; float r19400 = r19399 * r19399; float r19401 = r19386 * r19400; float r19402 = U_; float r19403 = r19391 - r19402; float r19404 = r19401 * r19403; float r19405 = r19398 - r19404; float r19406 = r19391 * r19405; float r19407 = r19390 * r19406; float r19408 = sqrt(r19407); float r19409 = cbrt(r19408); float r19410 = r19409 * (r19409 * r19409); float r19411 = sqrt(r19390); float r19412 = sqrt(r19406); float r19413 = r19411 * r19412; float r19414 = r19388 ? r19410 : r19413; return r19414; } double f_od(double n, double U, double t, double l, double Om, double U_) { double r19415 = n; double r19416 = 4.5602004537906414e-307; bool r19417 = r19415 <= r19416; double r19418 = 2.0; double r19419 = r19418 * r19415; double r19420 = U; double r19421 = t; double r19422 = l; double r19423 = Om; double r19424 = r19423 / r19422; double r19425 = r19422 / r19424; double r19426 = r19418 * r19425; double r19427 = r19421 - r19426; double r19428 = r19422 / r19423; double r19429 = r19428 * r19428; double r19430 = r19415 * r19429; double r19431 = U_; double r19432 = r19420 - r19431; double r19433 = r19430 * r19432; double r19434 = r19427 - r19433; double r19435 = r19420 * r19434; double r19436 = r19419 * r19435; double r19437 = sqrt(r19436); double r19438 = cbrt(r19437); double r19439 = r19438 * (r19438 * r19438); double r19440 = sqrt(r19419); double r19441 = sqrt(r19435); double r19442 = r19440 * r19441; double r19443 = r19417 ? r19439 : r19442; return r19443; } 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 r19444, r19445, r19446, r19447, r19448, r19449, r19450, r19451, r19452, r19453, r19454, r19455, r19456, r19457, r19458, r19459, r19460, r19461, r19462, r19463, r19464; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r19444, "2", 10, MPFR_RNDN); mpfr_init(r19445); mpfr_init(r19446); mpfr_init(r19447); mpfr_init(r19448); mpfr_init(r19449); mpfr_init(r19450); mpfr_init(r19451); mpfr_init(r19452); mpfr_init(r19453); mpfr_init(r19454); mpfr_init(r19455); mpfr_init(r19456); mpfr_init(r19457); mpfr_init(r19458); mpfr_init(r19459); mpfr_init(r19460); mpfr_init(r19461); mpfr_init(r19462); mpfr_init(r19463); mpfr_init(r19464); } double f_im(double n, double U, double t, double l, double Om, double U_) { ; mpfr_set_d(r19445, n, MPFR_RNDN); mpfr_mul(r19446, r19444, r19445, MPFR_RNDN); mpfr_set_d(r19447, U, MPFR_RNDN); mpfr_mul(r19448, r19446, r19447, MPFR_RNDN); mpfr_set_d(r19449, t, MPFR_RNDN); mpfr_set_d(r19450, l, MPFR_RNDN); mpfr_sqr(r19451, r19450, MPFR_RNDN); mpfr_set_d(r19452, Om, MPFR_RNDN); mpfr_div(r19453, r19451, r19452, MPFR_RNDN); mpfr_mul(r19454, r19444, r19453, MPFR_RNDN); mpfr_sub(r19455, r19449, r19454, MPFR_RNDN); mpfr_div(r19456, r19450, r19452, MPFR_RNDN); mpfr_sqr(r19457, r19456, MPFR_RNDN); mpfr_mul(r19458, r19445, r19457, MPFR_RNDN); mpfr_set_d(r19459, U_, MPFR_RNDN); mpfr_sub(r19460, r19447, r19459, MPFR_RNDN); mpfr_mul(r19461, r19458, r19460, MPFR_RNDN); mpfr_sub(r19462, r19455, r19461, MPFR_RNDN); mpfr_mul(r19463, r19448, r19462, MPFR_RNDN); mpfr_sqrt(r19464, r19463, MPFR_RNDN); return mpfr_get_d(r19464, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r19465); mpfr_init_set_str(r19466, "4.5602004537906414e-307", 10, MPFR_RNDN); mpfr_init(r19467); mpfr_init_set_str(r19468, "2", 10, MPFR_RNDN); mpfr_init(r19469); mpfr_init(r19470); mpfr_init(r19471); mpfr_init(r19472); mpfr_init(r19473); mpfr_init(r19474); mpfr_init(r19475); mpfr_init(r19476); mpfr_init(r19477); mpfr_init(r19478); mpfr_init(r19479); mpfr_init(r19480); mpfr_init(r19481); mpfr_init(r19482); mpfr_init(r19483); mpfr_init(r19484); mpfr_init(r19485); mpfr_init(r19486); mpfr_init(r19487); mpfr_init(r19488); mpfr_init(r19489); mpfr_init(r19490); mpfr_init(r19491); mpfr_init(r19492); mpfr_init(r19493); } double f_fm(double n, double U, double t, double l, double Om, double U_) { mpfr_set_d(r19465, n, MPFR_RNDN); ; mpfr_set_si(r19467, mpfr_cmp(r19465, r19466) <= 0, MPFR_RNDN); ; mpfr_mul(r19469, r19468, r19465, MPFR_RNDN); mpfr_set_d(r19470, U, MPFR_RNDN); mpfr_set_d(r19471, t, MPFR_RNDN); mpfr_set_d(r19472, l, MPFR_RNDN); mpfr_set_d(r19473, Om, MPFR_RNDN); mpfr_div(r19474, r19473, r19472, MPFR_RNDN); mpfr_div(r19475, r19472, r19474, MPFR_RNDN); mpfr_mul(r19476, r19468, r19475, MPFR_RNDN); mpfr_sub(r19477, r19471, r19476, MPFR_RNDN); mpfr_div(r19478, r19472, r19473, MPFR_RNDN); mpfr_sqr(r19479, r19478, MPFR_RNDN); mpfr_mul(r19480, r19465, r19479, MPFR_RNDN); mpfr_set_d(r19481, U_, MPFR_RNDN); mpfr_sub(r19482, r19470, r19481, MPFR_RNDN); mpfr_mul(r19483, r19480, r19482, MPFR_RNDN); mpfr_sub(r19484, r19477, r19483, MPFR_RNDN); mpfr_mul(r19485, r19470, r19484, MPFR_RNDN); mpfr_mul(r19486, r19469, r19485, MPFR_RNDN); mpfr_sqrt(r19487, r19486, MPFR_RNDN); mpfr_cbrt(r19488, r19487, MPFR_RNDN); mpfr_mul(r19489, r19488, r19488, MPFR_RNDN); mpfr_mul(r19489, r19489, r19488, MPFR_RNDN); mpfr_sqrt(r19490, r19469, MPFR_RNDN); mpfr_sqrt(r19491, r19485, MPFR_RNDN); mpfr_mul(r19492, r19490, r19491, MPFR_RNDN); if (mpfr_get_si(r19467, MPFR_RNDN)) { mpfr_set(r19493, r19489, MPFR_RNDN); } else { mpfr_set(r19493, r19492, MPFR_RNDN); }; return mpfr_get_d(r19493, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r19494); mpfr_init_set_str(r19495, "4.5602004537906414e-307", 10, MPFR_RNDN); mpfr_init(r19496); mpfr_init_set_str(r19497, "2", 10, MPFR_RNDN); mpfr_init(r19498); mpfr_init(r19499); mpfr_init(r19500); mpfr_init(r19501); mpfr_init(r19502); mpfr_init(r19503); mpfr_init(r19504); mpfr_init(r19505); mpfr_init(r19506); mpfr_init(r19507); mpfr_init(r19508); mpfr_init(r19509); mpfr_init(r19510); mpfr_init(r19511); mpfr_init(r19512); mpfr_init(r19513); mpfr_init(r19514); mpfr_init(r19515); mpfr_init(r19516); mpfr_init(r19517); mpfr_init(r19518); mpfr_init(r19519); mpfr_init(r19520); mpfr_init(r19521); mpfr_init(r19522); } double f_dm(double n, double U, double t, double l, double Om, double U_) { mpfr_set_d(r19494, n, MPFR_RNDN); ; mpfr_set_si(r19496, mpfr_cmp(r19494, r19495) <= 0, MPFR_RNDN); ; mpfr_mul(r19498, r19497, r19494, MPFR_RNDN); mpfr_set_d(r19499, U, MPFR_RNDN); mpfr_set_d(r19500, t, MPFR_RNDN); mpfr_set_d(r19501, l, MPFR_RNDN); mpfr_set_d(r19502, Om, MPFR_RNDN); mpfr_div(r19503, r19502, r19501, MPFR_RNDN); mpfr_div(r19504, r19501, r19503, MPFR_RNDN); mpfr_mul(r19505, r19497, r19504, MPFR_RNDN); mpfr_sub(r19506, r19500, r19505, MPFR_RNDN); mpfr_div(r19507, r19501, r19502, MPFR_RNDN); mpfr_sqr(r19508, r19507, MPFR_RNDN); mpfr_mul(r19509, r19494, r19508, MPFR_RNDN); mpfr_set_d(r19510, U_, MPFR_RNDN); mpfr_sub(r19511, r19499, r19510, MPFR_RNDN); mpfr_mul(r19512, r19509, r19511, MPFR_RNDN); mpfr_sub(r19513, r19506, r19512, MPFR_RNDN); mpfr_mul(r19514, r19499, r19513, MPFR_RNDN); mpfr_mul(r19515, r19498, r19514, MPFR_RNDN); mpfr_sqrt(r19516, r19515, MPFR_RNDN); mpfr_cbrt(r19517, r19516, MPFR_RNDN); mpfr_mul(r19518, r19517, r19517, MPFR_RNDN); mpfr_mul(r19518, r19518, r19517, MPFR_RNDN); mpfr_sqrt(r19519, r19498, MPFR_RNDN); mpfr_sqrt(r19520, r19514, MPFR_RNDN); mpfr_mul(r19521, r19519, r19520, MPFR_RNDN); if (mpfr_get_si(r19496, MPFR_RNDN)) { mpfr_set(r19522, r19518, MPFR_RNDN); } else { mpfr_set(r19522, r19521, MPFR_RNDN); }; return mpfr_get_d(r19522, MPFR_RNDN); }