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

char *name = "Octave 3.8, jcobi/1";

double f_if(float alpha, float beta) {
        float r17793 = beta;
        float r17794 = alpha;
        float r17795 = r17793 - r17794;
        float r17796 = r17794 + r17793;
        float r17797 = 2.0f;
        float r17798 = r17796 + r17797;
        float r17799 = r17795 / r17798;
        float r17800 = 1.0f;
        float r17801 = r17799 + r17800;
        float r17802 = r17801 / r17797;
        return r17802;
}

double f_id(double alpha, double beta) {
        double r17803 = beta;
        double r17804 = alpha;
        double r17805 = r17803 - r17804;
        double r17806 = r17804 + r17803;
        double r17807 = 2.0;
        double r17808 = r17806 + r17807;
        double r17809 = r17805 / r17808;
        double r17810 = 1.0;
        double r17811 = r17809 + r17810;
        double r17812 = r17811 / r17807;
        return r17812;
}


double f_of(float alpha, float beta) {
        float r17813 = beta;
        float r17814 = alpha;
        float r17815 = r17813 - r17814;
        float r17816 = r17814 + r17813;
        float r17817 = 2.0f;
        float r17818 = r17816 + r17817;
        float r17819 = r17815 / r17818;
        float r17820 = -0.9997063657335743f;
        bool r17821 = r17819 <= r17820;
        float r17822 = 8.0f;
        float r17823 = r17822 / r17814;
        float r17824 = r17823 / r17814;
        float r17825 = r17817 + r17824;
        float r17826 = r17817 * r17814;
        float r17827 = r17825 / r17826;
        float r17828 = r17813 / r17817;
        float r17829 = r17817 + r17813;
        float r17830 = r17814 + r17829;
        float r17831 = r17828 / r17830;
        float r17832 = 4.0f;
        float r17833 = r17814 * r17814;
        float r17834 = r17832 / r17833;
        float r17835 = r17834 / r17817;
        float r17836 = r17831 - r17835;
        float r17837 = r17827 + r17836;
        float r17838 = r17813 / r17818;
        float r17839 = 1.0f;
        float r17840 = r17839 / r17818;
        float r17841 = r17814 * r17840;
        float r17842 = 1.0f;
        float r17843 = r17841 - r17842;
        float r17844 = r17838 - r17843;
        float r17845 = r17844 / r17817;
        float r17846 = r17821 ? r17837 : r17845;
        return r17846;
}

double f_od(double alpha, double beta) {
        double r17847 = beta;
        double r17848 = alpha;
        double r17849 = r17847 - r17848;
        double r17850 = r17848 + r17847;
        double r17851 = 2.0;
        double r17852 = r17850 + r17851;
        double r17853 = r17849 / r17852;
        double r17854 = -0.9997063657335743;
        bool r17855 = r17853 <= r17854;
        double r17856 = 8.0;
        double r17857 = r17856 / r17848;
        double r17858 = r17857 / r17848;
        double r17859 = r17851 + r17858;
        double r17860 = r17851 * r17848;
        double r17861 = r17859 / r17860;
        double r17862 = r17847 / r17851;
        double r17863 = r17851 + r17847;
        double r17864 = r17848 + r17863;
        double r17865 = r17862 / r17864;
        double r17866 = 4.0;
        double r17867 = r17848 * r17848;
        double r17868 = r17866 / r17867;
        double r17869 = r17868 / r17851;
        double r17870 = r17865 - r17869;
        double r17871 = r17861 + r17870;
        double r17872 = r17847 / r17852;
        double r17873 = 1.0;
        double r17874 = r17873 / r17852;
        double r17875 = r17848 * r17874;
        double r17876 = 1.0;
        double r17877 = r17875 - r17876;
        double r17878 = r17872 - r17877;
        double r17879 = r17878 / r17851;
        double r17880 = r17855 ? r17871 : r17879;
        return r17880;
}

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 r17881, r17882, r17883, r17884, r17885, r17886, r17887, r17888, r17889, r17890;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17881);
        mpfr_init(r17882);
        mpfr_init(r17883);
        mpfr_init(r17884);
        mpfr_init_set_str(r17885, "2.0", 10, MPFR_RNDN);
        mpfr_init(r17886);
        mpfr_init(r17887);
        mpfr_init_set_str(r17888, "1.0", 10, MPFR_RNDN);
        mpfr_init(r17889);
        mpfr_init(r17890);
}

double f_im(double alpha, double beta) {
        mpfr_set_d(r17881, beta, MPFR_RNDN);
        mpfr_set_d(r17882, alpha, MPFR_RNDN);
        mpfr_sub(r17883, r17881, r17882, MPFR_RNDN);
        mpfr_add(r17884, r17882, r17881, MPFR_RNDN);
        ;
        mpfr_add(r17886, r17884, r17885, MPFR_RNDN);
        mpfr_div(r17887, r17883, r17886, MPFR_RNDN);
        ;
        mpfr_add(r17889, r17887, r17888, MPFR_RNDN);
        mpfr_div(r17890, r17889, r17885, MPFR_RNDN);
        return mpfr_get_d(r17890, MPFR_RNDN);
}

static mpfr_t r17891, r17892, r17893, r17894, r17895, r17896, r17897, r17898, r17899, r17900, r17901, r17902, r17903, r17904, r17905, r17906, r17907, r17908, r17909, r17910, r17911, r17912, r17913, r17914, r17915, r17916, r17917, r17918, r17919, r17920, r17921, r17922, r17923, r17924;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17891);
        mpfr_init(r17892);
        mpfr_init(r17893);
        mpfr_init(r17894);
        mpfr_init_set_str(r17895, "2.0", 10, MPFR_RNDN);
        mpfr_init(r17896);
        mpfr_init(r17897);
        mpfr_init_set_str(r17898, "-0.9997063657335743", 10, MPFR_RNDN);
        mpfr_init(r17899);
        mpfr_init_set_str(r17900, "8.0", 10, MPFR_RNDN);
        mpfr_init(r17901);
        mpfr_init(r17902);
        mpfr_init(r17903);
        mpfr_init(r17904);
        mpfr_init(r17905);
        mpfr_init(r17906);
        mpfr_init(r17907);
        mpfr_init(r17908);
        mpfr_init(r17909);
        mpfr_init_set_str(r17910, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17911);
        mpfr_init(r17912);
        mpfr_init(r17913);
        mpfr_init(r17914);
        mpfr_init(r17915);
        mpfr_init(r17916);
        mpfr_init_set_str(r17917, "1", 10, MPFR_RNDN);
        mpfr_init(r17918);
        mpfr_init(r17919);
        mpfr_init_set_str(r17920, "1.0", 10, MPFR_RNDN);
        mpfr_init(r17921);
        mpfr_init(r17922);
        mpfr_init(r17923);
        mpfr_init(r17924);
}

double f_fm(double alpha, double beta) {
        mpfr_set_d(r17891, beta, MPFR_RNDN);
        mpfr_set_d(r17892, alpha, MPFR_RNDN);
        mpfr_sub(r17893, r17891, r17892, MPFR_RNDN);
        mpfr_add(r17894, r17892, r17891, MPFR_RNDN);
        ;
        mpfr_add(r17896, r17894, r17895, MPFR_RNDN);
        mpfr_div(r17897, r17893, r17896, MPFR_RNDN);
        ;
        mpfr_set_si(r17899, mpfr_cmp(r17897, r17898) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r17901, r17900, r17892, MPFR_RNDN);
        mpfr_div(r17902, r17901, r17892, MPFR_RNDN);
        mpfr_add(r17903, r17895, r17902, MPFR_RNDN);
        mpfr_mul(r17904, r17895, r17892, MPFR_RNDN);
        mpfr_div(r17905, r17903, r17904, MPFR_RNDN);
        mpfr_div(r17906, r17891, r17895, MPFR_RNDN);
        mpfr_add(r17907, r17895, r17891, MPFR_RNDN);
        mpfr_add(r17908, r17892, r17907, MPFR_RNDN);
        mpfr_div(r17909, r17906, r17908, MPFR_RNDN);
        ;
        mpfr_mul(r17911, r17892, r17892, MPFR_RNDN);
        mpfr_div(r17912, r17910, r17911, MPFR_RNDN);
        mpfr_div(r17913, r17912, r17895, MPFR_RNDN);
        mpfr_sub(r17914, r17909, r17913, MPFR_RNDN);
        mpfr_add(r17915, r17905, r17914, MPFR_RNDN);
        mpfr_div(r17916, r17891, r17896, MPFR_RNDN);
        ;
        mpfr_div(r17918, r17917, r17896, MPFR_RNDN);
        mpfr_mul(r17919, r17892, r17918, MPFR_RNDN);
        ;
        mpfr_sub(r17921, r17919, r17920, MPFR_RNDN);
        mpfr_sub(r17922, r17916, r17921, MPFR_RNDN);
        mpfr_div(r17923, r17922, r17895, MPFR_RNDN);
        if (mpfr_get_si(r17899, MPFR_RNDN)) { mpfr_set(r17924, r17915, MPFR_RNDN); } else { mpfr_set(r17924, r17923, MPFR_RNDN); };
        return mpfr_get_d(r17924, MPFR_RNDN);
}

static mpfr_t r17925, r17926, r17927, r17928, r17929, r17930, r17931, r17932, r17933, r17934, r17935, r17936, r17937, r17938, r17939, r17940, r17941, r17942, r17943, r17944, r17945, r17946, r17947, r17948, r17949, r17950, r17951, r17952, r17953, r17954, r17955, r17956, r17957, r17958;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17925);
        mpfr_init(r17926);
        mpfr_init(r17927);
        mpfr_init(r17928);
        mpfr_init_set_str(r17929, "2.0", 10, MPFR_RNDN);
        mpfr_init(r17930);
        mpfr_init(r17931);
        mpfr_init_set_str(r17932, "-0.9997063657335743", 10, MPFR_RNDN);
        mpfr_init(r17933);
        mpfr_init_set_str(r17934, "8.0", 10, MPFR_RNDN);
        mpfr_init(r17935);
        mpfr_init(r17936);
        mpfr_init(r17937);
        mpfr_init(r17938);
        mpfr_init(r17939);
        mpfr_init(r17940);
        mpfr_init(r17941);
        mpfr_init(r17942);
        mpfr_init(r17943);
        mpfr_init_set_str(r17944, "4.0", 10, MPFR_RNDN);
        mpfr_init(r17945);
        mpfr_init(r17946);
        mpfr_init(r17947);
        mpfr_init(r17948);
        mpfr_init(r17949);
        mpfr_init(r17950);
        mpfr_init_set_str(r17951, "1", 10, MPFR_RNDN);
        mpfr_init(r17952);
        mpfr_init(r17953);
        mpfr_init_set_str(r17954, "1.0", 10, MPFR_RNDN);
        mpfr_init(r17955);
        mpfr_init(r17956);
        mpfr_init(r17957);
        mpfr_init(r17958);
}

double f_dm(double alpha, double beta) {
        mpfr_set_d(r17925, beta, MPFR_RNDN);
        mpfr_set_d(r17926, alpha, MPFR_RNDN);
        mpfr_sub(r17927, r17925, r17926, MPFR_RNDN);
        mpfr_add(r17928, r17926, r17925, MPFR_RNDN);
        ;
        mpfr_add(r17930, r17928, r17929, MPFR_RNDN);
        mpfr_div(r17931, r17927, r17930, MPFR_RNDN);
        ;
        mpfr_set_si(r17933, mpfr_cmp(r17931, r17932) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r17935, r17934, r17926, MPFR_RNDN);
        mpfr_div(r17936, r17935, r17926, MPFR_RNDN);
        mpfr_add(r17937, r17929, r17936, MPFR_RNDN);
        mpfr_mul(r17938, r17929, r17926, MPFR_RNDN);
        mpfr_div(r17939, r17937, r17938, MPFR_RNDN);
        mpfr_div(r17940, r17925, r17929, MPFR_RNDN);
        mpfr_add(r17941, r17929, r17925, MPFR_RNDN);
        mpfr_add(r17942, r17926, r17941, MPFR_RNDN);
        mpfr_div(r17943, r17940, r17942, MPFR_RNDN);
        ;
        mpfr_mul(r17945, r17926, r17926, MPFR_RNDN);
        mpfr_div(r17946, r17944, r17945, MPFR_RNDN);
        mpfr_div(r17947, r17946, r17929, MPFR_RNDN);
        mpfr_sub(r17948, r17943, r17947, MPFR_RNDN);
        mpfr_add(r17949, r17939, r17948, MPFR_RNDN);
        mpfr_div(r17950, r17925, r17930, MPFR_RNDN);
        ;
        mpfr_div(r17952, r17951, r17930, MPFR_RNDN);
        mpfr_mul(r17953, r17926, r17952, MPFR_RNDN);
        ;
        mpfr_sub(r17955, r17953, r17954, MPFR_RNDN);
        mpfr_sub(r17956, r17950, r17955, MPFR_RNDN);
        mpfr_div(r17957, r17956, r17929, MPFR_RNDN);
        if (mpfr_get_si(r17933, MPFR_RNDN)) { mpfr_set(r17958, r17949, MPFR_RNDN); } else { mpfr_set(r17958, r17957, MPFR_RNDN); };
        return mpfr_get_d(r17958, MPFR_RNDN);
}