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