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

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

double f_if(float l, float Om, float kx, float ky) {
        float r22141 = 1;
        float r22142 = 2;
        float r22143 = r22141 / r22142;
        float r22144 = l;
        float r22145 = r22142 * r22144;
        float r22146 = Om;
        float r22147 = r22145 / r22146;
        float r22148 = pow(r22147, r22142);
        float r22149 = kx;
        float r22150 = sin(r22149);
        float r22151 = pow(r22150, r22142);
        float r22152 = ky;
        float r22153 = sin(r22152);
        float r22154 = pow(r22153, r22142);
        float r22155 = r22151 + r22154;
        float r22156 = r22148 * r22155;
        float r22157 = r22141 + r22156;
        float r22158 = sqrt(r22157);
        float r22159 = r22141 / r22158;
        float r22160 = r22141 + r22159;
        float r22161 = r22143 * r22160;
        float r22162 = sqrt(r22161);
        return r22162;
}

double f_id(double l, double Om, double kx, double ky) {
        double r22163 = 1;
        double r22164 = 2;
        double r22165 = r22163 / r22164;
        double r22166 = l;
        double r22167 = r22164 * r22166;
        double r22168 = Om;
        double r22169 = r22167 / r22168;
        double r22170 = pow(r22169, r22164);
        double r22171 = kx;
        double r22172 = sin(r22171);
        double r22173 = pow(r22172, r22164);
        double r22174 = ky;
        double r22175 = sin(r22174);
        double r22176 = pow(r22175, r22164);
        double r22177 = r22173 + r22176;
        double r22178 = r22170 * r22177;
        double r22179 = r22163 + r22178;
        double r22180 = sqrt(r22179);
        double r22181 = r22163 / r22180;
        double r22182 = r22163 + r22181;
        double r22183 = r22165 * r22182;
        double r22184 = sqrt(r22183);
        return r22184;
}


double f_of(float l, float Om, float kx, float ky) {
        float r22185 = 1;
        float r22186 = 2;
        float r22187 = r22185 / r22186;
        float r22188 = l;
        float r22189 = r22186 * r22188;
        float r22190 = Om;
        float r22191 = r22189 / r22190;
        float r22192 = cbrt(r22191);
        float r22193 = r22192 * r22192;
        float r22194 = pow(r22193, r22186);
        float r22195 = pow(r22192, r22186);
        float r22196 = kx;
        float r22197 = sin(r22196);
        float r22198 = pow(r22197, r22186);
        float r22199 = ky;
        float r22200 = sin(r22199);
        float r22201 = pow(r22200, r22186);
        float r22202 = r22198 + r22201;
        float r22203 = r22195 * r22202;
        float r22204 = r22194 * r22203;
        float r22205 = r22185 + r22204;
        float r22206 = sqrt(r22205);
        float r22207 = r22185 / r22206;
        float r22208 = r22185 + r22207;
        float r22209 = r22187 * r22208;
        float r22210 = sqrt(r22209);
        return r22210;
}

double f_od(double l, double Om, double kx, double ky) {
        double r22211 = 1;
        double r22212 = 2;
        double r22213 = r22211 / r22212;
        double r22214 = l;
        double r22215 = r22212 * r22214;
        double r22216 = Om;
        double r22217 = r22215 / r22216;
        double r22218 = cbrt(r22217);
        double r22219 = r22218 * r22218;
        double r22220 = pow(r22219, r22212);
        double r22221 = pow(r22218, r22212);
        double r22222 = kx;
        double r22223 = sin(r22222);
        double r22224 = pow(r22223, r22212);
        double r22225 = ky;
        double r22226 = sin(r22225);
        double r22227 = pow(r22226, r22212);
        double r22228 = r22224 + r22227;
        double r22229 = r22221 * r22228;
        double r22230 = r22220 * r22229;
        double r22231 = r22211 + r22230;
        double r22232 = sqrt(r22231);
        double r22233 = r22211 / r22232;
        double r22234 = r22211 + r22233;
        double r22235 = r22213 * r22234;
        double r22236 = sqrt(r22235);
        return r22236;
}

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 r22237, r22238, r22239, r22240, r22241, r22242, r22243, r22244, r22245, r22246, r22247, r22248, r22249, r22250, r22251, r22252, r22253, r22254, r22255, r22256, r22257, r22258;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r22237, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22238, "2", 10, MPFR_RNDN);
        mpfr_init(r22239);
        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);
        mpfr_init(r22253);
        mpfr_init(r22254);
        mpfr_init(r22255);
        mpfr_init(r22256);
        mpfr_init(r22257);
        mpfr_init(r22258);
}

double f_im(double l, double Om, double kx, double ky) {
        ;
        ;
        mpfr_div(r22239, r22237, r22238, MPFR_RNDN);
        mpfr_set_d(r22240, l, MPFR_RNDN);
        mpfr_mul(r22241, r22238, r22240, MPFR_RNDN);
        mpfr_set_d(r22242, Om, MPFR_RNDN);
        mpfr_div(r22243, r22241, r22242, MPFR_RNDN);
        mpfr_pow(r22244, r22243, r22238, MPFR_RNDN);
        mpfr_set_d(r22245, kx, MPFR_RNDN);
        mpfr_sin(r22246, r22245, MPFR_RNDN);
        mpfr_pow(r22247, r22246, r22238, MPFR_RNDN);
        mpfr_set_d(r22248, ky, MPFR_RNDN);
        mpfr_sin(r22249, r22248, MPFR_RNDN);
        mpfr_pow(r22250, r22249, r22238, MPFR_RNDN);
        mpfr_add(r22251, r22247, r22250, MPFR_RNDN);
        mpfr_mul(r22252, r22244, r22251, MPFR_RNDN);
        mpfr_add(r22253, r22237, r22252, MPFR_RNDN);
        mpfr_sqrt(r22254, r22253, MPFR_RNDN);
        mpfr_div(r22255, r22237, r22254, MPFR_RNDN);
        mpfr_add(r22256, r22237, r22255, MPFR_RNDN);
        mpfr_mul(r22257, r22239, r22256, MPFR_RNDN);
        mpfr_sqrt(r22258, r22257, MPFR_RNDN);
        return mpfr_get_d(r22258, MPFR_RNDN);
}

static mpfr_t r22259, r22260, r22261, r22262, r22263, r22264, r22265, r22266, r22267, r22268, r22269, r22270, r22271, r22272, r22273, r22274, r22275, r22276, r22277, r22278, r22279, r22280, r22281, r22282, r22283, r22284;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r22259, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22260, "2", 10, MPFR_RNDN);
        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);
        mpfr_init(r22271);
        mpfr_init(r22272);
        mpfr_init(r22273);
        mpfr_init(r22274);
        mpfr_init(r22275);
        mpfr_init(r22276);
        mpfr_init(r22277);
        mpfr_init(r22278);
        mpfr_init(r22279);
        mpfr_init(r22280);
        mpfr_init(r22281);
        mpfr_init(r22282);
        mpfr_init(r22283);
        mpfr_init(r22284);
}

double f_fm(double l, double Om, double kx, double ky) {
        ;
        ;
        mpfr_div(r22261, r22259, r22260, MPFR_RNDN);
        mpfr_set_d(r22262, l, MPFR_RNDN);
        mpfr_mul(r22263, r22260, r22262, MPFR_RNDN);
        mpfr_set_d(r22264, Om, MPFR_RNDN);
        mpfr_div(r22265, r22263, r22264, MPFR_RNDN);
        mpfr_cbrt(r22266, r22265, MPFR_RNDN);
        mpfr_mul(r22267, r22266, r22266, MPFR_RNDN);
        mpfr_pow(r22268, r22267, r22260, MPFR_RNDN);
        mpfr_pow(r22269, r22266, r22260, MPFR_RNDN);
        mpfr_set_d(r22270, kx, MPFR_RNDN);
        mpfr_sin(r22271, r22270, MPFR_RNDN);
        mpfr_pow(r22272, r22271, r22260, MPFR_RNDN);
        mpfr_set_d(r22273, ky, MPFR_RNDN);
        mpfr_sin(r22274, r22273, MPFR_RNDN);
        mpfr_pow(r22275, r22274, r22260, MPFR_RNDN);
        mpfr_add(r22276, r22272, r22275, MPFR_RNDN);
        mpfr_mul(r22277, r22269, r22276, MPFR_RNDN);
        mpfr_mul(r22278, r22268, r22277, MPFR_RNDN);
        mpfr_add(r22279, r22259, r22278, MPFR_RNDN);
        mpfr_sqrt(r22280, r22279, MPFR_RNDN);
        mpfr_div(r22281, r22259, r22280, MPFR_RNDN);
        mpfr_add(r22282, r22259, r22281, MPFR_RNDN);
        mpfr_mul(r22283, r22261, r22282, MPFR_RNDN);
        mpfr_sqrt(r22284, r22283, MPFR_RNDN);
        return mpfr_get_d(r22284, MPFR_RNDN);
}

static mpfr_t r22285, r22286, r22287, r22288, r22289, r22290, r22291, r22292, r22293, r22294, r22295, r22296, r22297, r22298, r22299, r22300, r22301, r22302, r22303, r22304, r22305, r22306, r22307, r22308, r22309, r22310;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r22285, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22286, "2", 10, MPFR_RNDN);
        mpfr_init(r22287);
        mpfr_init(r22288);
        mpfr_init(r22289);
        mpfr_init(r22290);
        mpfr_init(r22291);
        mpfr_init(r22292);
        mpfr_init(r22293);
        mpfr_init(r22294);
        mpfr_init(r22295);
        mpfr_init(r22296);
        mpfr_init(r22297);
        mpfr_init(r22298);
        mpfr_init(r22299);
        mpfr_init(r22300);
        mpfr_init(r22301);
        mpfr_init(r22302);
        mpfr_init(r22303);
        mpfr_init(r22304);
        mpfr_init(r22305);
        mpfr_init(r22306);
        mpfr_init(r22307);
        mpfr_init(r22308);
        mpfr_init(r22309);
        mpfr_init(r22310);
}

double f_dm(double l, double Om, double kx, double ky) {
        ;
        ;
        mpfr_div(r22287, r22285, r22286, MPFR_RNDN);
        mpfr_set_d(r22288, l, MPFR_RNDN);
        mpfr_mul(r22289, r22286, r22288, MPFR_RNDN);
        mpfr_set_d(r22290, Om, MPFR_RNDN);
        mpfr_div(r22291, r22289, r22290, MPFR_RNDN);
        mpfr_cbrt(r22292, r22291, MPFR_RNDN);
        mpfr_mul(r22293, r22292, r22292, MPFR_RNDN);
        mpfr_pow(r22294, r22293, r22286, MPFR_RNDN);
        mpfr_pow(r22295, r22292, r22286, MPFR_RNDN);
        mpfr_set_d(r22296, kx, MPFR_RNDN);
        mpfr_sin(r22297, r22296, MPFR_RNDN);
        mpfr_pow(r22298, r22297, r22286, MPFR_RNDN);
        mpfr_set_d(r22299, ky, MPFR_RNDN);
        mpfr_sin(r22300, r22299, MPFR_RNDN);
        mpfr_pow(r22301, r22300, r22286, MPFR_RNDN);
        mpfr_add(r22302, r22298, r22301, MPFR_RNDN);
        mpfr_mul(r22303, r22295, r22302, MPFR_RNDN);
        mpfr_mul(r22304, r22294, r22303, MPFR_RNDN);
        mpfr_add(r22305, r22285, r22304, MPFR_RNDN);
        mpfr_sqrt(r22306, r22305, MPFR_RNDN);
        mpfr_div(r22307, r22285, r22306, MPFR_RNDN);
        mpfr_add(r22308, r22285, r22307, MPFR_RNDN);
        mpfr_mul(r22309, r22287, r22308, MPFR_RNDN);
        mpfr_sqrt(r22310, r22309, MPFR_RNDN);
        return mpfr_get_d(r22310, MPFR_RNDN);
}