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

char *name = "cos(2*x)/(cos^2(x)*sin^2(x))";

double f_if(float x, float cos, float sin) {
        float r24767 = 2;
        float r24768 = x;
        float r24769 = r24767 * r24768;
        float r24770 = cos(r24769);
        float r24771 = cos;
        float r24772 = pow(r24771, r24767);
        float r24773 = sin;
        float r24774 = pow(r24773, r24767);
        float r24775 = r24768 * r24774;
        float r24776 = r24775 * r24768;
        float r24777 = r24772 * r24776;
        float r24778 = r24770 / r24777;
        return r24778;
}

double f_id(double x, double cos, double sin) {
        double r24779 = 2;
        double r24780 = x;
        double r24781 = r24779 * r24780;
        double r24782 = cos(r24781);
        double r24783 = cos;
        double r24784 = pow(r24783, r24779);
        double r24785 = sin;
        double r24786 = pow(r24785, r24779);
        double r24787 = r24780 * r24786;
        double r24788 = r24787 * r24780;
        double r24789 = r24784 * r24788;
        double r24790 = r24782 / r24789;
        return r24790;
}


double f_of(float x, float cos, float sin) {
        float r24791 = x;
        float r24792 = 2;
        float r24793 = r24791 * r24792;
        float r24794 = cos(r24793);
        float r24795 = cos;
        float r24796 = sin;
        float r24797 = r24796 * r24791;
        float r24798 = r24795 * r24797;
        float r24799 = r24794 / r24798;
        float r24800 = r24799 / r24798;
        float r24801 = 1.3649867799669e-317;
        bool r24802 = r24800 <= r24801;
        float r24803 = cos(r24791);
        float r24804 = r24796 * r24795;
        float r24805 = r24804 * r24791;
        float r24806 = fabs(r24805);
        float r24807 = r24803 / r24806;
        float r24808 = sin(r24791);
        float r24809 = -r24808;
        float r24810 = r24808 * r24809;
        float r24811 = r24805 * r24805;
        float r24812 = r24810 / r24811;
        float r24813 = fma(r24807, r24807, r24812);
        float r24814 = 1.2738224146149126e+242;
        bool r24815 = r24800 <= r24814;
        float r24816 = r24815 ? r24800 : r24813;
        float r24817 = r24802 ? r24813 : r24816;
        return r24817;
}

double f_od(double x, double cos, double sin) {
        double r24818 = x;
        double r24819 = 2;
        double r24820 = r24818 * r24819;
        double r24821 = cos(r24820);
        double r24822 = cos;
        double r24823 = sin;
        double r24824 = r24823 * r24818;
        double r24825 = r24822 * r24824;
        double r24826 = r24821 / r24825;
        double r24827 = r24826 / r24825;
        double r24828 = 1.3649867799669e-317;
        bool r24829 = r24827 <= r24828;
        double r24830 = cos(r24818);
        double r24831 = r24823 * r24822;
        double r24832 = r24831 * r24818;
        double r24833 = fabs(r24832);
        double r24834 = r24830 / r24833;
        double r24835 = sin(r24818);
        double r24836 = -r24835;
        double r24837 = r24835 * r24836;
        double r24838 = r24832 * r24832;
        double r24839 = r24837 / r24838;
        double r24840 = fma(r24834, r24834, r24839);
        double r24841 = 1.2738224146149126e+242;
        bool r24842 = r24827 <= r24841;
        double r24843 = r24842 ? r24827 : r24840;
        double r24844 = r24829 ? r24840 : r24843;
        return r24844;
}

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 r24845, r24846, r24847, r24848, r24849, r24850, r24851, r24852, r24853, r24854, r24855, r24856;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r24845, "2", 10, MPFR_RNDN);
        mpfr_init(r24846);
        mpfr_init(r24847);
        mpfr_init(r24848);
        mpfr_init(r24849);
        mpfr_init(r24850);
        mpfr_init(r24851);
        mpfr_init(r24852);
        mpfr_init(r24853);
        mpfr_init(r24854);
        mpfr_init(r24855);
        mpfr_init(r24856);
}

double f_im(double x, double cos, double sin) {
        ;
        mpfr_set_d(r24846, x, MPFR_RNDN);
        mpfr_mul(r24847, r24845, r24846, MPFR_RNDN);
        mpfr_cos(r24848, r24847, MPFR_RNDN);
        mpfr_set_d(r24849, cos, MPFR_RNDN);
        mpfr_pow(r24850, r24849, r24845, MPFR_RNDN);
        mpfr_set_d(r24851, sin, MPFR_RNDN);
        mpfr_pow(r24852, r24851, r24845, MPFR_RNDN);
        mpfr_mul(r24853, r24846, r24852, MPFR_RNDN);
        mpfr_mul(r24854, r24853, r24846, MPFR_RNDN);
        mpfr_mul(r24855, r24850, r24854, MPFR_RNDN);
        mpfr_div(r24856, r24848, r24855, MPFR_RNDN);
        return mpfr_get_d(r24856, MPFR_RNDN);
}

static mpfr_t r24857, r24858, r24859, r24860, r24861, r24862, r24863, r24864, r24865, r24866, r24867, r24868, r24869, r24870, r24871, r24872, r24873, r24874, r24875, r24876, r24877, r24878, r24879, r24880, r24881, r24882, r24883;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r24857);
        mpfr_init_set_str(r24858, "2", 10, MPFR_RNDN);
        mpfr_init(r24859);
        mpfr_init(r24860);
        mpfr_init(r24861);
        mpfr_init(r24862);
        mpfr_init(r24863);
        mpfr_init(r24864);
        mpfr_init(r24865);
        mpfr_init(r24866);
        mpfr_init_set_str(r24867, "1.3649867799669e-317", 10, MPFR_RNDN);
        mpfr_init(r24868);
        mpfr_init(r24869);
        mpfr_init(r24870);
        mpfr_init(r24871);
        mpfr_init(r24872);
        mpfr_init(r24873);
        mpfr_init(r24874);
        mpfr_init(r24875);
        mpfr_init(r24876);
        mpfr_init(r24877);
        mpfr_init(r24878);
        mpfr_init(r24879);
        mpfr_init_set_str(r24880, "1.2738224146149126e+242", 10, MPFR_RNDN);
        mpfr_init(r24881);
        mpfr_init(r24882);
        mpfr_init(r24883);
}

double f_fm(double x, double cos, double sin) {
        mpfr_set_d(r24857, x, MPFR_RNDN);
        ;
        mpfr_mul(r24859, r24857, r24858, MPFR_RNDN);
        mpfr_cos(r24860, r24859, MPFR_RNDN);
        mpfr_set_d(r24861, cos, MPFR_RNDN);
        mpfr_set_d(r24862, sin, MPFR_RNDN);
        mpfr_mul(r24863, r24862, r24857, MPFR_RNDN);
        mpfr_mul(r24864, r24861, r24863, MPFR_RNDN);
        mpfr_div(r24865, r24860, r24864, MPFR_RNDN);
        mpfr_div(r24866, r24865, r24864, MPFR_RNDN);
        ;
        mpfr_set_si(r24868, mpfr_cmp(r24866, r24867) <= 0, MPFR_RNDN);
        mpfr_cos(r24869, r24857, MPFR_RNDN);
        mpfr_mul(r24870, r24862, r24861, MPFR_RNDN);
        mpfr_mul(r24871, r24870, r24857, MPFR_RNDN);
        mpfr_abs(r24872, r24871, MPFR_RNDN);
        mpfr_div(r24873, r24869, r24872, MPFR_RNDN);
        mpfr_sin(r24874, r24857, MPFR_RNDN);
        mpfr_neg(r24875, r24874, MPFR_RNDN);
        mpfr_mul(r24876, r24874, r24875, MPFR_RNDN);
        mpfr_mul(r24877, r24871, r24871, MPFR_RNDN);
        mpfr_div(r24878, r24876, r24877, MPFR_RNDN);
        mpfr_fma(r24879, r24873, r24873, r24878, MPFR_RNDN);
        ;
        mpfr_set_si(r24881, mpfr_cmp(r24866, r24880) <= 0, MPFR_RNDN);
        if (mpfr_get_si(r24881, MPFR_RNDN)) { mpfr_set(r24882, r24866, MPFR_RNDN); } else { mpfr_set(r24882, r24879, MPFR_RNDN); };
        if (mpfr_get_si(r24868, MPFR_RNDN)) { mpfr_set(r24883, r24879, MPFR_RNDN); } else { mpfr_set(r24883, r24882, MPFR_RNDN); };
        return mpfr_get_d(r24883, MPFR_RNDN);
}

static mpfr_t r24884, r24885, r24886, r24887, r24888, r24889, r24890, r24891, r24892, r24893, r24894, r24895, r24896, r24897, r24898, r24899, r24900, r24901, r24902, r24903, r24904, r24905, r24906, r24907, r24908, r24909, r24910;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r24884);
        mpfr_init_set_str(r24885, "2", 10, MPFR_RNDN);
        mpfr_init(r24886);
        mpfr_init(r24887);
        mpfr_init(r24888);
        mpfr_init(r24889);
        mpfr_init(r24890);
        mpfr_init(r24891);
        mpfr_init(r24892);
        mpfr_init(r24893);
        mpfr_init_set_str(r24894, "1.3649867799669e-317", 10, MPFR_RNDN);
        mpfr_init(r24895);
        mpfr_init(r24896);
        mpfr_init(r24897);
        mpfr_init(r24898);
        mpfr_init(r24899);
        mpfr_init(r24900);
        mpfr_init(r24901);
        mpfr_init(r24902);
        mpfr_init(r24903);
        mpfr_init(r24904);
        mpfr_init(r24905);
        mpfr_init(r24906);
        mpfr_init_set_str(r24907, "1.2738224146149126e+242", 10, MPFR_RNDN);
        mpfr_init(r24908);
        mpfr_init(r24909);
        mpfr_init(r24910);
}

double f_dm(double x, double cos, double sin) {
        mpfr_set_d(r24884, x, MPFR_RNDN);
        ;
        mpfr_mul(r24886, r24884, r24885, MPFR_RNDN);
        mpfr_cos(r24887, r24886, MPFR_RNDN);
        mpfr_set_d(r24888, cos, MPFR_RNDN);
        mpfr_set_d(r24889, sin, MPFR_RNDN);
        mpfr_mul(r24890, r24889, r24884, MPFR_RNDN);
        mpfr_mul(r24891, r24888, r24890, MPFR_RNDN);
        mpfr_div(r24892, r24887, r24891, MPFR_RNDN);
        mpfr_div(r24893, r24892, r24891, MPFR_RNDN);
        ;
        mpfr_set_si(r24895, mpfr_cmp(r24893, r24894) <= 0, MPFR_RNDN);
        mpfr_cos(r24896, r24884, MPFR_RNDN);
        mpfr_mul(r24897, r24889, r24888, MPFR_RNDN);
        mpfr_mul(r24898, r24897, r24884, MPFR_RNDN);
        mpfr_abs(r24899, r24898, MPFR_RNDN);
        mpfr_div(r24900, r24896, r24899, MPFR_RNDN);
        mpfr_sin(r24901, r24884, MPFR_RNDN);
        mpfr_neg(r24902, r24901, MPFR_RNDN);
        mpfr_mul(r24903, r24901, r24902, MPFR_RNDN);
        mpfr_mul(r24904, r24898, r24898, MPFR_RNDN);
        mpfr_div(r24905, r24903, r24904, MPFR_RNDN);
        mpfr_fma(r24906, r24900, r24900, r24905, MPFR_RNDN);
        ;
        mpfr_set_si(r24908, mpfr_cmp(r24893, r24907) <= 0, MPFR_RNDN);
        if (mpfr_get_si(r24908, MPFR_RNDN)) { mpfr_set(r24909, r24893, MPFR_RNDN); } else { mpfr_set(r24909, r24906, MPFR_RNDN); };
        if (mpfr_get_si(r24895, MPFR_RNDN)) { mpfr_set(r24910, r24906, MPFR_RNDN); } else { mpfr_set(r24910, r24909, MPFR_RNDN); };
        return mpfr_get_d(r24910, MPFR_RNDN);
}