#include #include #include #include #include char *name = "Toniolo and Linder, Equation (2)"; double f_if(float t, float l, float Om, float Omc) { float r8357 = 1.0f; float r8358 = Om; float r8359 = Omc; float r8360 = r8358 / r8359; float r8361 = 2.0f; float r8362 = pow(r8360, r8361); float r8363 = r8357 - r8362; float r8364 = t; float r8365 = l; float r8366 = r8364 / r8365; float r8367 = pow(r8366, r8361); float r8368 = r8361 * r8367; float r8369 = r8357 + r8368; float r8370 = r8363 / r8369; float r8371 = sqrt(r8370); float r8372 = asin(r8371); return r8372; } double f_id(double t, double l, double Om, double Omc) { double r8373 = 1.0; double r8374 = Om; double r8375 = Omc; double r8376 = r8374 / r8375; double r8377 = 2.0; double r8378 = pow(r8376, r8377); double r8379 = r8373 - r8378; double r8380 = t; double r8381 = l; double r8382 = r8380 / r8381; double r8383 = pow(r8382, r8377); double r8384 = r8377 * r8383; double r8385 = r8373 + r8384; double r8386 = r8379 / r8385; double r8387 = sqrt(r8386); double r8388 = asin(r8387); return r8388; } double f_of(float t, float l, float Om, float Omc) { float r8389 = 1.0f; float r8390 = Om; float r8391 = Omc; float r8392 = r8390 / r8391; float r8393 = 2.0f; float r8394 = pow(r8392, r8393); float r8395 = r8389 - r8394; float r8396 = t; float r8397 = l; float r8398 = r8396 / r8397; float r8399 = pow(r8398, r8393); float r8400 = r8393 * r8399; float r8401 = r8389 + r8400; float r8402 = r8395 / r8401; float r8403 = log1p(r8402); float r8404 = expm1(r8403); float r8405 = sqrt(r8404); float r8406 = asin(r8405); return r8406; } double f_od(double t, double l, double Om, double Omc) { double r8407 = 1.0; double r8408 = Om; double r8409 = Omc; double r8410 = r8408 / r8409; double r8411 = 2.0; double r8412 = pow(r8410, r8411); double r8413 = r8407 - r8412; double r8414 = t; double r8415 = l; double r8416 = r8414 / r8415; double r8417 = pow(r8416, r8411); double r8418 = r8411 * r8417; double r8419 = r8407 + r8418; double r8420 = r8413 / r8419; double r8421 = log1p(r8420); double r8422 = expm1(r8421); double r8423 = sqrt(r8422); double r8424 = asin(r8423); return r8424; } 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 r8425, r8426, r8427, r8428, r8429, r8430, r8431, r8432, r8433, r8434, r8435, r8436, r8437, r8438, r8439, r8440; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r8425, "1", 10, MPFR_RNDN); mpfr_init(r8426); mpfr_init(r8427); mpfr_init(r8428); mpfr_init_set_str(r8429, "2", 10, MPFR_RNDN); mpfr_init(r8430); mpfr_init(r8431); mpfr_init(r8432); mpfr_init(r8433); mpfr_init(r8434); mpfr_init(r8435); mpfr_init(r8436); mpfr_init(r8437); mpfr_init(r8438); mpfr_init(r8439); mpfr_init(r8440); } double f_im(double t, double l, double Om, double Omc) { ; mpfr_set_d(r8426, Om, MPFR_RNDN); mpfr_set_d(r8427, Omc, MPFR_RNDN); mpfr_div(r8428, r8426, r8427, MPFR_RNDN); ; mpfr_pow(r8430, r8428, r8429, MPFR_RNDN); mpfr_sub(r8431, r8425, r8430, MPFR_RNDN); mpfr_set_d(r8432, t, MPFR_RNDN); mpfr_set_d(r8433, l, MPFR_RNDN); mpfr_div(r8434, r8432, r8433, MPFR_RNDN); mpfr_pow(r8435, r8434, r8429, MPFR_RNDN); mpfr_mul(r8436, r8429, r8435, MPFR_RNDN); mpfr_add(r8437, r8425, r8436, MPFR_RNDN); mpfr_div(r8438, r8431, r8437, MPFR_RNDN); mpfr_sqrt(r8439, r8438, MPFR_RNDN); mpfr_asin(r8440, r8439, MPFR_RNDN); return mpfr_get_d(r8440, MPFR_RNDN); } static mpfr_t r8441, r8442, r8443, r8444, r8445, r8446, r8447, r8448, r8449, r8450, r8451, r8452, r8453, r8454, r8455, r8456, r8457, r8458; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r8441, "1", 10, MPFR_RNDN); mpfr_init(r8442); mpfr_init(r8443); mpfr_init(r8444); mpfr_init_set_str(r8445, "2", 10, MPFR_RNDN); mpfr_init(r8446); mpfr_init(r8447); mpfr_init(r8448); mpfr_init(r8449); mpfr_init(r8450); mpfr_init(r8451); mpfr_init(r8452); mpfr_init(r8453); mpfr_init(r8454); mpfr_init(r8455); mpfr_init(r8456); mpfr_init(r8457); mpfr_init(r8458); } double f_fm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r8442, Om, MPFR_RNDN); mpfr_set_d(r8443, Omc, MPFR_RNDN); mpfr_div(r8444, r8442, r8443, MPFR_RNDN); ; mpfr_pow(r8446, r8444, r8445, MPFR_RNDN); mpfr_sub(r8447, r8441, r8446, MPFR_RNDN); mpfr_set_d(r8448, t, MPFR_RNDN); mpfr_set_d(r8449, l, MPFR_RNDN); mpfr_div(r8450, r8448, r8449, MPFR_RNDN); mpfr_pow(r8451, r8450, r8445, MPFR_RNDN); mpfr_mul(r8452, r8445, r8451, MPFR_RNDN); mpfr_add(r8453, r8441, r8452, MPFR_RNDN); mpfr_div(r8454, r8447, r8453, MPFR_RNDN); mpfr_log1p(r8455, r8454, MPFR_RNDN); mpfr_expm1(r8456, r8455, MPFR_RNDN); mpfr_sqrt(r8457, r8456, MPFR_RNDN); mpfr_asin(r8458, r8457, MPFR_RNDN); return mpfr_get_d(r8458, MPFR_RNDN); } static mpfr_t r8459, r8460, r8461, r8462, r8463, r8464, r8465, r8466, r8467, r8468, r8469, r8470, r8471, r8472, r8473, r8474, r8475, r8476; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r8459, "1", 10, MPFR_RNDN); mpfr_init(r8460); mpfr_init(r8461); mpfr_init(r8462); mpfr_init_set_str(r8463, "2", 10, MPFR_RNDN); mpfr_init(r8464); mpfr_init(r8465); mpfr_init(r8466); mpfr_init(r8467); mpfr_init(r8468); mpfr_init(r8469); mpfr_init(r8470); mpfr_init(r8471); mpfr_init(r8472); mpfr_init(r8473); mpfr_init(r8474); mpfr_init(r8475); mpfr_init(r8476); } double f_dm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r8460, Om, MPFR_RNDN); mpfr_set_d(r8461, Omc, MPFR_RNDN); mpfr_div(r8462, r8460, r8461, MPFR_RNDN); ; mpfr_pow(r8464, r8462, r8463, MPFR_RNDN); mpfr_sub(r8465, r8459, r8464, MPFR_RNDN); mpfr_set_d(r8466, t, MPFR_RNDN); mpfr_set_d(r8467, l, MPFR_RNDN); mpfr_div(r8468, r8466, r8467, MPFR_RNDN); mpfr_pow(r8469, r8468, r8463, MPFR_RNDN); mpfr_mul(r8470, r8463, r8469, MPFR_RNDN); mpfr_add(r8471, r8459, r8470, MPFR_RNDN); mpfr_div(r8472, r8465, r8471, MPFR_RNDN); mpfr_log1p(r8473, r8472, MPFR_RNDN); mpfr_expm1(r8474, r8473, MPFR_RNDN); mpfr_sqrt(r8475, r8474, MPFR_RNDN); mpfr_asin(r8476, r8475, MPFR_RNDN); return mpfr_get_d(r8476, MPFR_RNDN); }