#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 r8468 = ky;
        float r8469 = sin(r8468);
        float r8470 = kx;
        float r8471 = sin(r8470);
        float r8472 = 2;
        float r8473 = pow(r8471, r8472);
        float r8474 = pow(r8469, r8472);
        float r8475 = r8473 + r8474;
        float r8476 = sqrt(r8475);
        float r8477 = r8469 / r8476;
        float r8478 = th;
        float r8479 = sin(r8478);
        float r8480 = r8477 * r8479;
        return r8480;
}

double f_id(double kx, double ky, double th) {
        double r8481 = ky;
        double r8482 = sin(r8481);
        double r8483 = kx;
        double r8484 = sin(r8483);
        double r8485 = 2;
        double r8486 = pow(r8484, r8485);
        double r8487 = pow(r8482, r8485);
        double r8488 = r8486 + r8487;
        double r8489 = sqrt(r8488);
        double r8490 = r8482 / r8489;
        double r8491 = th;
        double r8492 = sin(r8491);
        double r8493 = r8490 * r8492;
        return r8493;
}


double f_of(float kx, float ky, float th) {
        float r8494 = 1;
        float r8495 = kx;
        float r8496 = sin(r8495);
        float r8497 = 2;
        float r8498 = pow(r8496, r8497);
        float r8499 = ky;
        float r8500 = sin(r8499);
        float r8501 = pow(r8500, r8497);
        float r8502 = r8498 + r8501;
        float r8503 = sqrt(r8502);
        float r8504 = r8503 / r8500;
        float r8505 = r8494 / r8504;
        float r8506 = th;
        float r8507 = sin(r8506);
        float r8508 = r8505 * r8507;
        return r8508;
}

double f_od(double kx, double ky, double th) {
        double r8509 = 1;
        double r8510 = kx;
        double r8511 = sin(r8510);
        double r8512 = 2;
        double r8513 = pow(r8511, r8512);
        double r8514 = ky;
        double r8515 = sin(r8514);
        double r8516 = pow(r8515, r8512);
        double r8517 = r8513 + r8516;
        double r8518 = sqrt(r8517);
        double r8519 = r8518 / r8515;
        double r8520 = r8509 / r8519;
        double r8521 = th;
        double r8522 = sin(r8521);
        double r8523 = r8520 * r8522;
        return r8523;
}

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 r8524, r8525, r8526, r8527, r8528, r8529, r8530, r8531, r8532, r8533, r8534, r8535, r8536;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r8524);
        mpfr_init(r8525);
        mpfr_init(r8526);
        mpfr_init(r8527);
        mpfr_init_set_str(r8528, "2", 10, MPFR_RNDN);
        mpfr_init(r8529);
        mpfr_init(r8530);
        mpfr_init(r8531);
        mpfr_init(r8532);
        mpfr_init(r8533);
        mpfr_init(r8534);
        mpfr_init(r8535);
        mpfr_init(r8536);
}

double f_im(double kx, double ky, double th) {
        mpfr_set_d(r8524, ky, MPFR_RNDN);
        mpfr_sin(r8525, r8524, MPFR_RNDN);
        mpfr_set_d(r8526, kx, MPFR_RNDN);
        mpfr_sin(r8527, r8526, MPFR_RNDN);
        ;
        mpfr_pow(r8529, r8527, r8528, MPFR_RNDN);
        mpfr_pow(r8530, r8525, r8528, MPFR_RNDN);
        mpfr_add(r8531, r8529, r8530, MPFR_RNDN);
        mpfr_sqrt(r8532, r8531, MPFR_RNDN);
        mpfr_div(r8533, r8525, r8532, MPFR_RNDN);
        mpfr_set_d(r8534, th, MPFR_RNDN);
        mpfr_sin(r8535, r8534, MPFR_RNDN);
        mpfr_mul(r8536, r8533, r8535, MPFR_RNDN);
        return mpfr_get_d(r8536, MPFR_RNDN);
}

static mpfr_t r8537, r8538, r8539, r8540, r8541, r8542, r8543, r8544, r8545, r8546, r8547, r8548, r8549, r8550, r8551;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r8537, "1", 10, MPFR_RNDN);
        mpfr_init(r8538);
        mpfr_init(r8539);
        mpfr_init_set_str(r8540, "2", 10, MPFR_RNDN);
        mpfr_init(r8541);
        mpfr_init(r8542);
        mpfr_init(r8543);
        mpfr_init(r8544);
        mpfr_init(r8545);
        mpfr_init(r8546);
        mpfr_init(r8547);
        mpfr_init(r8548);
        mpfr_init(r8549);
        mpfr_init(r8550);
        mpfr_init(r8551);
}

double f_fm(double kx, double ky, double th) {
        ;
        mpfr_set_d(r8538, kx, MPFR_RNDN);
        mpfr_sin(r8539, r8538, MPFR_RNDN);
        ;
        mpfr_pow(r8541, r8539, r8540, MPFR_RNDN);
        mpfr_set_d(r8542, ky, MPFR_RNDN);
        mpfr_sin(r8543, r8542, MPFR_RNDN);
        mpfr_pow(r8544, r8543, r8540, MPFR_RNDN);
        mpfr_add(r8545, r8541, r8544, MPFR_RNDN);
        mpfr_sqrt(r8546, r8545, MPFR_RNDN);
        mpfr_div(r8547, r8546, r8543, MPFR_RNDN);
        mpfr_div(r8548, r8537, r8547, MPFR_RNDN);
        mpfr_set_d(r8549, th, MPFR_RNDN);
        mpfr_sin(r8550, r8549, MPFR_RNDN);
        mpfr_mul(r8551, r8548, r8550, MPFR_RNDN);
        return mpfr_get_d(r8551, MPFR_RNDN);
}

static mpfr_t r8552, r8553, r8554, r8555, r8556, r8557, r8558, r8559, r8560, r8561, r8562, r8563, r8564, r8565, r8566;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r8552, "1", 10, MPFR_RNDN);
        mpfr_init(r8553);
        mpfr_init(r8554);
        mpfr_init_set_str(r8555, "2", 10, MPFR_RNDN);
        mpfr_init(r8556);
        mpfr_init(r8557);
        mpfr_init(r8558);
        mpfr_init(r8559);
        mpfr_init(r8560);
        mpfr_init(r8561);
        mpfr_init(r8562);
        mpfr_init(r8563);
        mpfr_init(r8564);
        mpfr_init(r8565);
        mpfr_init(r8566);
}

double f_dm(double kx, double ky, double th) {
        ;
        mpfr_set_d(r8553, kx, MPFR_RNDN);
        mpfr_sin(r8554, r8553, MPFR_RNDN);
        ;
        mpfr_pow(r8556, r8554, r8555, MPFR_RNDN);
        mpfr_set_d(r8557, ky, MPFR_RNDN);
        mpfr_sin(r8558, r8557, MPFR_RNDN);
        mpfr_pow(r8559, r8558, r8555, MPFR_RNDN);
        mpfr_add(r8560, r8556, r8559, MPFR_RNDN);
        mpfr_sqrt(r8561, r8560, MPFR_RNDN);
        mpfr_div(r8562, r8561, r8558, MPFR_RNDN);
        mpfr_div(r8563, r8552, r8562, MPFR_RNDN);
        mpfr_set_d(r8564, th, MPFR_RNDN);
        mpfr_sin(r8565, r8564, MPFR_RNDN);
        mpfr_mul(r8566, r8563, r8565, MPFR_RNDN);
        return mpfr_get_d(r8566, MPFR_RNDN);
}