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

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

double f_if(float kx, float ky, float th) {
        float r20194 = ky;
        float r20195 = sin(r20194);
        float r20196 = kx;
        float r20197 = sin(r20196);
        float r20198 = 2;
        float r20199 = pow(r20197, r20198);
        float r20200 = pow(r20195, r20198);
        float r20201 = r20199 + r20200;
        float r20202 = sqrt(r20201);
        float r20203 = r20195 / r20202;
        float r20204 = th;
        float r20205 = sin(r20204);
        float r20206 = r20203 * r20205;
        return r20206;
}

double f_id(double kx, double ky, double th) {
        double r20207 = ky;
        double r20208 = sin(r20207);
        double r20209 = kx;
        double r20210 = sin(r20209);
        double r20211 = 2;
        double r20212 = pow(r20210, r20211);
        double r20213 = pow(r20208, r20211);
        double r20214 = r20212 + r20213;
        double r20215 = sqrt(r20214);
        double r20216 = r20208 / r20215;
        double r20217 = th;
        double r20218 = sin(r20217);
        double r20219 = r20216 * r20218;
        return r20219;
}


double f_of(float kx, float ky, float th) {
        float r20220 = ky;
        float r20221 = sin(r20220);
        float r20222 = 1;
        float r20223 = kx;
        float r20224 = sin(r20223);
        float r20225 = 2;
        float r20226 = pow(r20224, r20225);
        float r20227 = pow(r20221, r20225);
        float r20228 = r20226 + r20227;
        float r20229 = sqrt(r20228);
        float r20230 = r20222 / r20229;
        float r20231 = th;
        float r20232 = sin(r20231);
        float r20233 = r20230 * r20232;
        float r20234 = r20221 * r20233;
        return r20234;
}

double f_od(double kx, double ky, double th) {
        double r20235 = ky;
        double r20236 = sin(r20235);
        double r20237 = 1;
        double r20238 = kx;
        double r20239 = sin(r20238);
        double r20240 = 2;
        double r20241 = pow(r20239, r20240);
        double r20242 = pow(r20236, r20240);
        double r20243 = r20241 + r20242;
        double r20244 = sqrt(r20243);
        double r20245 = r20237 / r20244;
        double r20246 = th;
        double r20247 = sin(r20246);
        double r20248 = r20245 * r20247;
        double r20249 = r20236 * r20248;
        return r20249;
}

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 r20250, r20251, r20252, r20253, r20254, r20255, r20256, r20257, r20258, r20259, r20260, r20261, r20262;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r20250);
        mpfr_init(r20251);
        mpfr_init(r20252);
        mpfr_init(r20253);
        mpfr_init_set_str(r20254, "2", 10, MPFR_RNDN);
        mpfr_init(r20255);
        mpfr_init(r20256);
        mpfr_init(r20257);
        mpfr_init(r20258);
        mpfr_init(r20259);
        mpfr_init(r20260);
        mpfr_init(r20261);
        mpfr_init(r20262);
}

double f_im(double kx, double ky, double th) {
        mpfr_set_d(r20250, ky, MPFR_RNDN);
        mpfr_sin(r20251, r20250, MPFR_RNDN);
        mpfr_set_d(r20252, kx, MPFR_RNDN);
        mpfr_sin(r20253, r20252, MPFR_RNDN);
        ;
        mpfr_pow(r20255, r20253, r20254, MPFR_RNDN);
        mpfr_pow(r20256, r20251, r20254, MPFR_RNDN);
        mpfr_add(r20257, r20255, r20256, MPFR_RNDN);
        mpfr_sqrt(r20258, r20257, MPFR_RNDN);
        mpfr_div(r20259, r20251, r20258, MPFR_RNDN);
        mpfr_set_d(r20260, th, MPFR_RNDN);
        mpfr_sin(r20261, r20260, MPFR_RNDN);
        mpfr_mul(r20262, r20259, r20261, MPFR_RNDN);
        return mpfr_get_d(r20262, MPFR_RNDN);
}

static mpfr_t r20263, r20264, r20265, r20266, r20267, r20268, r20269, r20270, r20271, r20272, r20273, r20274, r20275, r20276, r20277;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20263);
        mpfr_init(r20264);
        mpfr_init_set_str(r20265, "1", 10, MPFR_RNDN);
        mpfr_init(r20266);
        mpfr_init(r20267);
        mpfr_init_set_str(r20268, "2", 10, MPFR_RNDN);
        mpfr_init(r20269);
        mpfr_init(r20270);
        mpfr_init(r20271);
        mpfr_init(r20272);
        mpfr_init(r20273);
        mpfr_init(r20274);
        mpfr_init(r20275);
        mpfr_init(r20276);
        mpfr_init(r20277);
}

double f_fm(double kx, double ky, double th) {
        mpfr_set_d(r20263, ky, MPFR_RNDN);
        mpfr_sin(r20264, r20263, MPFR_RNDN);
        ;
        mpfr_set_d(r20266, kx, MPFR_RNDN);
        mpfr_sin(r20267, r20266, MPFR_RNDN);
        ;
        mpfr_pow(r20269, r20267, r20268, MPFR_RNDN);
        mpfr_pow(r20270, r20264, r20268, MPFR_RNDN);
        mpfr_add(r20271, r20269, r20270, MPFR_RNDN);
        mpfr_sqrt(r20272, r20271, MPFR_RNDN);
        mpfr_div(r20273, r20265, r20272, MPFR_RNDN);
        mpfr_set_d(r20274, th, MPFR_RNDN);
        mpfr_sin(r20275, r20274, MPFR_RNDN);
        mpfr_mul(r20276, r20273, r20275, MPFR_RNDN);
        mpfr_mul(r20277, r20264, r20276, MPFR_RNDN);
        return mpfr_get_d(r20277, MPFR_RNDN);
}

static mpfr_t r20278, r20279, r20280, r20281, r20282, r20283, r20284, r20285, r20286, r20287, r20288, r20289, r20290, r20291, r20292;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20278);
        mpfr_init(r20279);
        mpfr_init_set_str(r20280, "1", 10, MPFR_RNDN);
        mpfr_init(r20281);
        mpfr_init(r20282);
        mpfr_init_set_str(r20283, "2", 10, MPFR_RNDN);
        mpfr_init(r20284);
        mpfr_init(r20285);
        mpfr_init(r20286);
        mpfr_init(r20287);
        mpfr_init(r20288);
        mpfr_init(r20289);
        mpfr_init(r20290);
        mpfr_init(r20291);
        mpfr_init(r20292);
}

double f_dm(double kx, double ky, double th) {
        mpfr_set_d(r20278, ky, MPFR_RNDN);
        mpfr_sin(r20279, r20278, MPFR_RNDN);
        ;
        mpfr_set_d(r20281, kx, MPFR_RNDN);
        mpfr_sin(r20282, r20281, MPFR_RNDN);
        ;
        mpfr_pow(r20284, r20282, r20283, MPFR_RNDN);
        mpfr_pow(r20285, r20279, r20283, MPFR_RNDN);
        mpfr_add(r20286, r20284, r20285, MPFR_RNDN);
        mpfr_sqrt(r20287, r20286, MPFR_RNDN);
        mpfr_div(r20288, r20280, r20287, MPFR_RNDN);
        mpfr_set_d(r20289, th, MPFR_RNDN);
        mpfr_sin(r20290, r20289, MPFR_RNDN);
        mpfr_mul(r20291, r20288, r20290, MPFR_RNDN);
        mpfr_mul(r20292, r20279, r20291, MPFR_RNDN);
        return mpfr_get_d(r20292, MPFR_RNDN);
}