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

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

double f_if(float t, float l, float k) {
        float r22622 = 2;
        float r22623 = t;
        float r22624 = 3;
        float r22625 = pow(r22623, r22624);
        float r22626 = l;
        float r22627 = r22626 * r22626;
        float r22628 = r22625 / r22627;
        float r22629 = k;
        float r22630 = sin(r22629);
        float r22631 = r22628 * r22630;
        float r22632 = tan(r22629);
        float r22633 = r22631 * r22632;
        float r22634 = 1;
        float r22635 = r22629 / r22623;
        float r22636 = pow(r22635, r22622);
        float r22637 = r22634 + r22636;
        float r22638 = r22637 - r22634;
        float r22639 = r22633 * r22638;
        float r22640 = r22622 / r22639;
        return r22640;
}

double f_id(double t, double l, double k) {
        double r22641 = 2;
        double r22642 = t;
        double r22643 = 3;
        double r22644 = pow(r22642, r22643);
        double r22645 = l;
        double r22646 = r22645 * r22645;
        double r22647 = r22644 / r22646;
        double r22648 = k;
        double r22649 = sin(r22648);
        double r22650 = r22647 * r22649;
        double r22651 = tan(r22648);
        double r22652 = r22650 * r22651;
        double r22653 = 1;
        double r22654 = r22648 / r22642;
        double r22655 = pow(r22654, r22641);
        double r22656 = r22653 + r22655;
        double r22657 = r22656 - r22653;
        double r22658 = r22652 * r22657;
        double r22659 = r22641 / r22658;
        return r22659;
}


double f_of(float t, float l, float k) {
        float r22660 = k;
        float r22661 = -9.809847864742864e-106;
        bool r22662 = r22660 <= r22661;
        float r22663 = 2;
        float r22664 = l;
        float r22665 = r22660 / r22664;
        float r22666 = r22663 / r22665;
        float r22667 = t;
        float r22668 = cos(r22660);
        float r22669 = r22667 / r22668;
        float r22670 = r22666 / r22669;
        float r22671 = sin(r22660);
        float r22672 = r22664 / r22671;
        float r22673 = r22671 * r22660;
        float r22674 = r22672 / r22673;
        float r22675 = r22670 * r22674;
        float r22676 = 5.037389137920242e-150;
        bool r22677 = r22660 <= r22676;
        float r22678 = 1;
        float r22679 = r22678 / r22660;
        float r22680 = r22667 / r22663;
        float r22681 = r22680 * r22665;
        float r22682 = r22671 / r22664;
        float r22683 = r22681 * r22682;
        float r22684 = r22671 / r22668;
        float r22685 = r22683 * r22684;
        float r22686 = r22679 / r22685;
        float r22687 = r22677 ? r22686 : r22675;
        float r22688 = r22662 ? r22675 : r22687;
        return r22688;
}

double f_od(double t, double l, double k) {
        double r22689 = k;
        double r22690 = -9.809847864742864e-106;
        bool r22691 = r22689 <= r22690;
        double r22692 = 2;
        double r22693 = l;
        double r22694 = r22689 / r22693;
        double r22695 = r22692 / r22694;
        double r22696 = t;
        double r22697 = cos(r22689);
        double r22698 = r22696 / r22697;
        double r22699 = r22695 / r22698;
        double r22700 = sin(r22689);
        double r22701 = r22693 / r22700;
        double r22702 = r22700 * r22689;
        double r22703 = r22701 / r22702;
        double r22704 = r22699 * r22703;
        double r22705 = 5.037389137920242e-150;
        bool r22706 = r22689 <= r22705;
        double r22707 = 1;
        double r22708 = r22707 / r22689;
        double r22709 = r22696 / r22692;
        double r22710 = r22709 * r22694;
        double r22711 = r22700 / r22693;
        double r22712 = r22710 * r22711;
        double r22713 = r22700 / r22697;
        double r22714 = r22712 * r22713;
        double r22715 = r22708 / r22714;
        double r22716 = r22706 ? r22715 : r22704;
        double r22717 = r22691 ? r22704 : r22716;
        return r22717;
}

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 r22718, r22719, r22720, r22721, r22722, r22723, r22724, r22725, r22726, r22727, r22728, r22729, r22730, r22731, r22732, r22733, r22734, r22735, r22736;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(4240);
        mpfr_init_set_str(r22718, "2", 10, MPFR_RNDN);
        mpfr_init(r22719);
        mpfr_init_set_str(r22720, "3", 10, MPFR_RNDN);
        mpfr_init(r22721);
        mpfr_init(r22722);
        mpfr_init(r22723);
        mpfr_init(r22724);
        mpfr_init(r22725);
        mpfr_init(r22726);
        mpfr_init(r22727);
        mpfr_init(r22728);
        mpfr_init(r22729);
        mpfr_init_set_str(r22730, "1", 10, MPFR_RNDN);
        mpfr_init(r22731);
        mpfr_init(r22732);
        mpfr_init(r22733);
        mpfr_init(r22734);
        mpfr_init(r22735);
        mpfr_init(r22736);
}

double f_im(double t, double l, double k) {
        ;
        mpfr_set_d(r22719, t, MPFR_RNDN);
        ;
        mpfr_pow(r22721, r22719, r22720, MPFR_RNDN);
        mpfr_set_d(r22722, l, MPFR_RNDN);
        mpfr_mul(r22723, r22722, r22722, MPFR_RNDN);
        mpfr_div(r22724, r22721, r22723, MPFR_RNDN);
        mpfr_set_d(r22725, k, MPFR_RNDN);
        mpfr_sin(r22726, r22725, MPFR_RNDN);
        mpfr_mul(r22727, r22724, r22726, MPFR_RNDN);
        mpfr_tan(r22728, r22725, MPFR_RNDN);
        mpfr_mul(r22729, r22727, r22728, MPFR_RNDN);
        ;
        mpfr_div(r22731, r22725, r22719, MPFR_RNDN);
        mpfr_pow(r22732, r22731, r22718, MPFR_RNDN);
        mpfr_add(r22733, r22730, r22732, MPFR_RNDN);
        mpfr_sub(r22734, r22733, r22730, MPFR_RNDN);
        mpfr_mul(r22735, r22729, r22734, MPFR_RNDN);
        mpfr_div(r22736, r22718, r22735, MPFR_RNDN);
        return mpfr_get_d(r22736, MPFR_RNDN);
}

static mpfr_t r22737, r22738, r22739, r22740, r22741, r22742, r22743, r22744, r22745, r22746, r22747, r22748, r22749, r22750, r22751, r22752, r22753, r22754, r22755, r22756, r22757, r22758, r22759, r22760, r22761, r22762, r22763, r22764, r22765;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(4240);
        mpfr_init(r22737);
        mpfr_init_set_str(r22738, "-9.809847864742864e-106", 10, MPFR_RNDN);
        mpfr_init(r22739);
        mpfr_init_set_str(r22740, "2", 10, MPFR_RNDN);
        mpfr_init(r22741);
        mpfr_init(r22742);
        mpfr_init(r22743);
        mpfr_init(r22744);
        mpfr_init(r22745);
        mpfr_init(r22746);
        mpfr_init(r22747);
        mpfr_init(r22748);
        mpfr_init(r22749);
        mpfr_init(r22750);
        mpfr_init(r22751);
        mpfr_init(r22752);
        mpfr_init_set_str(r22753, "5.037389137920242e-150", 10, MPFR_RNDN);
        mpfr_init(r22754);
        mpfr_init_set_str(r22755, "1", 10, MPFR_RNDN);
        mpfr_init(r22756);
        mpfr_init(r22757);
        mpfr_init(r22758);
        mpfr_init(r22759);
        mpfr_init(r22760);
        mpfr_init(r22761);
        mpfr_init(r22762);
        mpfr_init(r22763);
        mpfr_init(r22764);
        mpfr_init(r22765);
}

double f_fm(double t, double l, double k) {
        mpfr_set_d(r22737, k, MPFR_RNDN);
        ;
        mpfr_set_si(r22739, mpfr_cmp(r22737, r22738) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r22741, l, MPFR_RNDN);
        mpfr_div(r22742, r22737, r22741, MPFR_RNDN);
        mpfr_div(r22743, r22740, r22742, MPFR_RNDN);
        mpfr_set_d(r22744, t, MPFR_RNDN);
        mpfr_cos(r22745, r22737, MPFR_RNDN);
        mpfr_div(r22746, r22744, r22745, MPFR_RNDN);
        mpfr_div(r22747, r22743, r22746, MPFR_RNDN);
        mpfr_sin(r22748, r22737, MPFR_RNDN);
        mpfr_div(r22749, r22741, r22748, MPFR_RNDN);
        mpfr_mul(r22750, r22748, r22737, MPFR_RNDN);
        mpfr_div(r22751, r22749, r22750, MPFR_RNDN);
        mpfr_mul(r22752, r22747, r22751, MPFR_RNDN);
        ;
        mpfr_set_si(r22754, mpfr_cmp(r22737, r22753) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r22756, r22755, r22737, MPFR_RNDN);
        mpfr_div(r22757, r22744, r22740, MPFR_RNDN);
        mpfr_mul(r22758, r22757, r22742, MPFR_RNDN);
        mpfr_div(r22759, r22748, r22741, MPFR_RNDN);
        mpfr_mul(r22760, r22758, r22759, MPFR_RNDN);
        mpfr_div(r22761, r22748, r22745, MPFR_RNDN);
        mpfr_mul(r22762, r22760, r22761, MPFR_RNDN);
        mpfr_div(r22763, r22756, r22762, MPFR_RNDN);
        if (mpfr_get_si(r22754, MPFR_RNDN)) { mpfr_set(r22764, r22763, MPFR_RNDN); } else { mpfr_set(r22764, r22752, MPFR_RNDN); };
        if (mpfr_get_si(r22739, MPFR_RNDN)) { mpfr_set(r22765, r22752, MPFR_RNDN); } else { mpfr_set(r22765, r22764, MPFR_RNDN); };
        return mpfr_get_d(r22765, MPFR_RNDN);
}

static mpfr_t r22766, r22767, r22768, r22769, r22770, r22771, r22772, r22773, r22774, r22775, r22776, r22777, r22778, r22779, r22780, r22781, r22782, r22783, r22784, r22785, r22786, r22787, r22788, r22789, r22790, r22791, r22792, r22793, r22794;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(4240);
        mpfr_init(r22766);
        mpfr_init_set_str(r22767, "-9.809847864742864e-106", 10, MPFR_RNDN);
        mpfr_init(r22768);
        mpfr_init_set_str(r22769, "2", 10, MPFR_RNDN);
        mpfr_init(r22770);
        mpfr_init(r22771);
        mpfr_init(r22772);
        mpfr_init(r22773);
        mpfr_init(r22774);
        mpfr_init(r22775);
        mpfr_init(r22776);
        mpfr_init(r22777);
        mpfr_init(r22778);
        mpfr_init(r22779);
        mpfr_init(r22780);
        mpfr_init(r22781);
        mpfr_init_set_str(r22782, "5.037389137920242e-150", 10, MPFR_RNDN);
        mpfr_init(r22783);
        mpfr_init_set_str(r22784, "1", 10, MPFR_RNDN);
        mpfr_init(r22785);
        mpfr_init(r22786);
        mpfr_init(r22787);
        mpfr_init(r22788);
        mpfr_init(r22789);
        mpfr_init(r22790);
        mpfr_init(r22791);
        mpfr_init(r22792);
        mpfr_init(r22793);
        mpfr_init(r22794);
}

double f_dm(double t, double l, double k) {
        mpfr_set_d(r22766, k, MPFR_RNDN);
        ;
        mpfr_set_si(r22768, mpfr_cmp(r22766, r22767) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r22770, l, MPFR_RNDN);
        mpfr_div(r22771, r22766, r22770, MPFR_RNDN);
        mpfr_div(r22772, r22769, r22771, MPFR_RNDN);
        mpfr_set_d(r22773, t, MPFR_RNDN);
        mpfr_cos(r22774, r22766, MPFR_RNDN);
        mpfr_div(r22775, r22773, r22774, MPFR_RNDN);
        mpfr_div(r22776, r22772, r22775, MPFR_RNDN);
        mpfr_sin(r22777, r22766, MPFR_RNDN);
        mpfr_div(r22778, r22770, r22777, MPFR_RNDN);
        mpfr_mul(r22779, r22777, r22766, MPFR_RNDN);
        mpfr_div(r22780, r22778, r22779, MPFR_RNDN);
        mpfr_mul(r22781, r22776, r22780, MPFR_RNDN);
        ;
        mpfr_set_si(r22783, mpfr_cmp(r22766, r22782) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r22785, r22784, r22766, MPFR_RNDN);
        mpfr_div(r22786, r22773, r22769, MPFR_RNDN);
        mpfr_mul(r22787, r22786, r22771, MPFR_RNDN);
        mpfr_div(r22788, r22777, r22770, MPFR_RNDN);
        mpfr_mul(r22789, r22787, r22788, MPFR_RNDN);
        mpfr_div(r22790, r22777, r22774, MPFR_RNDN);
        mpfr_mul(r22791, r22789, r22790, MPFR_RNDN);
        mpfr_div(r22792, r22785, r22791, MPFR_RNDN);
        if (mpfr_get_si(r22783, MPFR_RNDN)) { mpfr_set(r22793, r22792, MPFR_RNDN); } else { mpfr_set(r22793, r22781, MPFR_RNDN); };
        if (mpfr_get_si(r22768, MPFR_RNDN)) { mpfr_set(r22794, r22781, MPFR_RNDN); } else { mpfr_set(r22794, r22793, MPFR_RNDN); };
        return mpfr_get_d(r22794, MPFR_RNDN);
}