#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Toniolo and Linder, Equation (2)";

double f_if(float t, float l, float Om, float Omc) {
        float r22151 = 1;
        float r22152 = Om;
        float r22153 = Omc;
        float r22154 = r22152 / r22153;
        float r22155 = 2;
        float r22156 = pow(r22154, r22155);
        float r22157 = r22151 - r22156;
        float r22158 = t;
        float r22159 = l;
        float r22160 = r22158 / r22159;
        float r22161 = pow(r22160, r22155);
        float r22162 = r22155 * r22161;
        float r22163 = r22151 + r22162;
        float r22164 = r22157 / r22163;
        float r22165 = sqrt(r22164);
        float r22166 = asin(r22165);
        return r22166;
}

double f_id(double t, double l, double Om, double Omc) {
        double r22167 = 1;
        double r22168 = Om;
        double r22169 = Omc;
        double r22170 = r22168 / r22169;
        double r22171 = 2;
        double r22172 = pow(r22170, r22171);
        double r22173 = r22167 - r22172;
        double r22174 = t;
        double r22175 = l;
        double r22176 = r22174 / r22175;
        double r22177 = pow(r22176, r22171);
        double r22178 = r22171 * r22177;
        double r22179 = r22167 + r22178;
        double r22180 = r22173 / r22179;
        double r22181 = sqrt(r22180);
        double r22182 = asin(r22181);
        return r22182;
}


double f_of(float t, float l, float Om, float Omc) {
        float r22183 = 1;
        float r22184 = Om;
        float r22185 = Omc;
        float r22186 = r22184 / r22185;
        float r22187 = 2;
        float r22188 = pow(r22186, r22187);
        float r22189 = r22183 - r22188;
        float r22190 = t;
        float r22191 = l;
        float r22192 = r22190 / r22191;
        float r22193 = pow(r22192, r22187);
        float r22194 = r22187 * r22193;
        float r22195 = r22183 + r22194;
        float r22196 = r22189 / r22195;
        float r22197 = sqrt(r22196);
        float r22198 = asin(r22197);
        float r22199 = log1p(r22198);
        float r22200 = expm1(r22199);
        return r22200;
}

double f_od(double t, double l, double Om, double Omc) {
        double r22201 = 1;
        double r22202 = Om;
        double r22203 = Omc;
        double r22204 = r22202 / r22203;
        double r22205 = 2;
        double r22206 = pow(r22204, r22205);
        double r22207 = r22201 - r22206;
        double r22208 = t;
        double r22209 = l;
        double r22210 = r22208 / r22209;
        double r22211 = pow(r22210, r22205);
        double r22212 = r22205 * r22211;
        double r22213 = r22201 + r22212;
        double r22214 = r22207 / r22213;
        double r22215 = sqrt(r22214);
        double r22216 = asin(r22215);
        double r22217 = log1p(r22216);
        double r22218 = expm1(r22217);
        return r22218;
}

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 r22219, r22220, r22221, r22222, r22223, r22224, r22225, r22226, r22227, r22228, r22229, r22230, r22231, r22232, r22233, r22234;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r22219, "1", 10, MPFR_RNDN);
        mpfr_init(r22220);
        mpfr_init(r22221);
        mpfr_init(r22222);
        mpfr_init_set_str(r22223, "2", 10, MPFR_RNDN);
        mpfr_init(r22224);
        mpfr_init(r22225);
        mpfr_init(r22226);
        mpfr_init(r22227);
        mpfr_init(r22228);
        mpfr_init(r22229);
        mpfr_init(r22230);
        mpfr_init(r22231);
        mpfr_init(r22232);
        mpfr_init(r22233);
        mpfr_init(r22234);
}

double f_im(double t, double l, double Om, double Omc) {
        ;
        mpfr_set_d(r22220, Om, MPFR_RNDN);
        mpfr_set_d(r22221, Omc, MPFR_RNDN);
        mpfr_div(r22222, r22220, r22221, MPFR_RNDN);
        ;
        mpfr_pow(r22224, r22222, r22223, MPFR_RNDN);
        mpfr_sub(r22225, r22219, r22224, MPFR_RNDN);
        mpfr_set_d(r22226, t, MPFR_RNDN);
        mpfr_set_d(r22227, l, MPFR_RNDN);
        mpfr_div(r22228, r22226, r22227, MPFR_RNDN);
        mpfr_pow(r22229, r22228, r22223, MPFR_RNDN);
        mpfr_mul(r22230, r22223, r22229, MPFR_RNDN);
        mpfr_add(r22231, r22219, r22230, MPFR_RNDN);
        mpfr_div(r22232, r22225, r22231, MPFR_RNDN);
        mpfr_sqrt(r22233, r22232, MPFR_RNDN);
        mpfr_asin(r22234, r22233, MPFR_RNDN);
        return mpfr_get_d(r22234, MPFR_RNDN);
}

static mpfr_t r22235, r22236, r22237, r22238, r22239, r22240, r22241, r22242, r22243, r22244, r22245, r22246, r22247, r22248, r22249, r22250, r22251, r22252;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r22235, "1", 10, MPFR_RNDN);
        mpfr_init(r22236);
        mpfr_init(r22237);
        mpfr_init(r22238);
        mpfr_init_set_str(r22239, "2", 10, MPFR_RNDN);
        mpfr_init(r22240);
        mpfr_init(r22241);
        mpfr_init(r22242);
        mpfr_init(r22243);
        mpfr_init(r22244);
        mpfr_init(r22245);
        mpfr_init(r22246);
        mpfr_init(r22247);
        mpfr_init(r22248);
        mpfr_init(r22249);
        mpfr_init(r22250);
        mpfr_init(r22251);
        mpfr_init(r22252);
}

double f_fm(double t, double l, double Om, double Omc) {
        ;
        mpfr_set_d(r22236, Om, MPFR_RNDN);
        mpfr_set_d(r22237, Omc, MPFR_RNDN);
        mpfr_div(r22238, r22236, r22237, MPFR_RNDN);
        ;
        mpfr_pow(r22240, r22238, r22239, MPFR_RNDN);
        mpfr_sub(r22241, r22235, r22240, MPFR_RNDN);
        mpfr_set_d(r22242, t, MPFR_RNDN);
        mpfr_set_d(r22243, l, MPFR_RNDN);
        mpfr_div(r22244, r22242, r22243, MPFR_RNDN);
        mpfr_pow(r22245, r22244, r22239, MPFR_RNDN);
        mpfr_mul(r22246, r22239, r22245, MPFR_RNDN);
        mpfr_add(r22247, r22235, r22246, MPFR_RNDN);
        mpfr_div(r22248, r22241, r22247, MPFR_RNDN);
        mpfr_sqrt(r22249, r22248, MPFR_RNDN);
        mpfr_asin(r22250, r22249, MPFR_RNDN);
        mpfr_log1p(r22251, r22250, MPFR_RNDN);
        mpfr_expm1(r22252, r22251, MPFR_RNDN);
        return mpfr_get_d(r22252, MPFR_RNDN);
}

static mpfr_t r22253, r22254, r22255, r22256, r22257, r22258, r22259, r22260, r22261, r22262, r22263, r22264, r22265, r22266, r22267, r22268, r22269, r22270;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r22253, "1", 10, MPFR_RNDN);
        mpfr_init(r22254);
        mpfr_init(r22255);
        mpfr_init(r22256);
        mpfr_init_set_str(r22257, "2", 10, MPFR_RNDN);
        mpfr_init(r22258);
        mpfr_init(r22259);
        mpfr_init(r22260);
        mpfr_init(r22261);
        mpfr_init(r22262);
        mpfr_init(r22263);
        mpfr_init(r22264);
        mpfr_init(r22265);
        mpfr_init(r22266);
        mpfr_init(r22267);
        mpfr_init(r22268);
        mpfr_init(r22269);
        mpfr_init(r22270);
}

double f_dm(double t, double l, double Om, double Omc) {
        ;
        mpfr_set_d(r22254, Om, MPFR_RNDN);
        mpfr_set_d(r22255, Omc, MPFR_RNDN);
        mpfr_div(r22256, r22254, r22255, MPFR_RNDN);
        ;
        mpfr_pow(r22258, r22256, r22257, MPFR_RNDN);
        mpfr_sub(r22259, r22253, r22258, MPFR_RNDN);
        mpfr_set_d(r22260, t, MPFR_RNDN);
        mpfr_set_d(r22261, l, MPFR_RNDN);
        mpfr_div(r22262, r22260, r22261, MPFR_RNDN);
        mpfr_pow(r22263, r22262, r22257, MPFR_RNDN);
        mpfr_mul(r22264, r22257, r22263, MPFR_RNDN);
        mpfr_add(r22265, r22253, r22264, MPFR_RNDN);
        mpfr_div(r22266, r22259, r22265, MPFR_RNDN);
        mpfr_sqrt(r22267, r22266, MPFR_RNDN);
        mpfr_asin(r22268, r22267, MPFR_RNDN);
        mpfr_log1p(r22269, r22268, MPFR_RNDN);
        mpfr_expm1(r22270, r22269, MPFR_RNDN);
        return mpfr_get_d(r22270, MPFR_RNDN);
}