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