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

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

double f_if(float alpha, float beta) {
        float r26188 = alpha;
        float r26189 = beta;
        float r26190 = r26188 + r26189;
        float r26191 = r26189 * r26188;
        float r26192 = r26190 + r26191;
        float r26193 = 1.0;
        float r26194 = r26192 + r26193;
        float r26195 = 2;
        float r26196 = 1;
        float r26197 = r26195 * r26196;
        float r26198 = r26190 + r26197;
        float r26199 = r26194 / r26198;
        float r26200 = r26199 / r26198;
        float r26201 = r26198 + r26193;
        float r26202 = r26200 / r26201;
        return r26202;
}

double f_id(double alpha, double beta) {
        double r26203 = alpha;
        double r26204 = beta;
        double r26205 = r26203 + r26204;
        double r26206 = r26204 * r26203;
        double r26207 = r26205 + r26206;
        double r26208 = 1.0;
        double r26209 = r26207 + r26208;
        double r26210 = 2;
        double r26211 = 1;
        double r26212 = r26210 * r26211;
        double r26213 = r26205 + r26212;
        double r26214 = r26209 / r26213;
        double r26215 = r26214 / r26213;
        double r26216 = r26213 + r26208;
        double r26217 = r26215 / r26216;
        return r26217;
}


double f_of(float alpha, float beta) {
        float r26218 = beta;
        float r26219 = alpha;
        float r26220 = r26218 * r26219;
        float r26221 = r26218 + r26219;
        float r26222 = r26220 + r26221;
        float r26223 = 1.0;
        float r26224 = r26222 + r26223;
        float r26225 = 2;
        float r26226 = r26225 + r26221;
        float r26227 = r26224 / r26226;
        float r26228 = r26227 / r26226;
        float r26229 = r26223 + r26226;
        float r26230 = r26228 / r26229;
        float r26231 = 1.9688390490121056e+227;
        bool r26232 = r26230 <= r26231;
        float r26233 = 1;
        float r26234 = r26233 / r26226;
        float r26235 = r26224 * r26234;
        float r26236 = r26235 / r26226;
        float r26237 = r26236 / r26229;
        float r26238 = 0.25;
        float r26239 = r26238 * r26221;
        float r26240 = 0.5;
        float r26241 = r26239 + r26240;
        float r26242 = r26223 + r26225;
        float r26243 = r26242 + r26221;
        float r26244 = r26226 * r26243;
        float r26245 = r26241 / r26244;
        float r26246 = r26232 ? r26237 : r26245;
        return r26246;
}

double f_od(double alpha, double beta) {
        double r26247 = beta;
        double r26248 = alpha;
        double r26249 = r26247 * r26248;
        double r26250 = r26247 + r26248;
        double r26251 = r26249 + r26250;
        double r26252 = 1.0;
        double r26253 = r26251 + r26252;
        double r26254 = 2;
        double r26255 = r26254 + r26250;
        double r26256 = r26253 / r26255;
        double r26257 = r26256 / r26255;
        double r26258 = r26252 + r26255;
        double r26259 = r26257 / r26258;
        double r26260 = 1.9688390490121056e+227;
        bool r26261 = r26259 <= r26260;
        double r26262 = 1;
        double r26263 = r26262 / r26255;
        double r26264 = r26253 * r26263;
        double r26265 = r26264 / r26255;
        double r26266 = r26265 / r26258;
        double r26267 = 0.25;
        double r26268 = r26267 * r26250;
        double r26269 = 0.5;
        double r26270 = r26268 + r26269;
        double r26271 = r26252 + r26254;
        double r26272 = r26271 + r26250;
        double r26273 = r26255 * r26272;
        double r26274 = r26270 / r26273;
        double r26275 = r26261 ? r26266 : r26274;
        return r26275;
}

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 r26276, r26277, r26278, r26279, r26280, r26281, r26282, r26283, r26284, r26285, r26286, r26287, r26288, r26289, r26290;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26276);
        mpfr_init(r26277);
        mpfr_init(r26278);
        mpfr_init(r26279);
        mpfr_init(r26280);
        mpfr_init_set_str(r26281, "1.0", 10, MPFR_RNDN);
        mpfr_init(r26282);
        mpfr_init_set_str(r26283, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r26284, "1", 10, MPFR_RNDN);
        mpfr_init(r26285);
        mpfr_init(r26286);
        mpfr_init(r26287);
        mpfr_init(r26288);
        mpfr_init(r26289);
        mpfr_init(r26290);
}

double f_im(double alpha, double beta) {
        mpfr_set_d(r26276, alpha, MPFR_RNDN);
        mpfr_set_d(r26277, beta, MPFR_RNDN);
        mpfr_add(r26278, r26276, r26277, MPFR_RNDN);
        mpfr_mul(r26279, r26277, r26276, MPFR_RNDN);
        mpfr_add(r26280, r26278, r26279, MPFR_RNDN);
        ;
        mpfr_add(r26282, r26280, r26281, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r26285, r26283, r26284, MPFR_RNDN);
        mpfr_add(r26286, r26278, r26285, MPFR_RNDN);
        mpfr_div(r26287, r26282, r26286, MPFR_RNDN);
        mpfr_div(r26288, r26287, r26286, MPFR_RNDN);
        mpfr_add(r26289, r26286, r26281, MPFR_RNDN);
        mpfr_div(r26290, r26288, r26289, MPFR_RNDN);
        return mpfr_get_d(r26290, MPFR_RNDN);
}

static mpfr_t r26291, r26292, r26293, r26294, r26295, r26296, r26297, r26298, r26299, r26300, r26301, r26302, r26303, r26304, r26305, r26306, r26307, r26308, r26309, r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318, r26319;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26291);
        mpfr_init(r26292);
        mpfr_init(r26293);
        mpfr_init(r26294);
        mpfr_init(r26295);
        mpfr_init_set_str(r26296, "1.0", 10, MPFR_RNDN);
        mpfr_init(r26297);
        mpfr_init_set_str(r26298, "2", 10, MPFR_RNDN);
        mpfr_init(r26299);
        mpfr_init(r26300);
        mpfr_init(r26301);
        mpfr_init(r26302);
        mpfr_init(r26303);
        mpfr_init_set_str(r26304, "1.9688390490121056e+227", 10, MPFR_RNDN);
        mpfr_init(r26305);
        mpfr_init_set_str(r26306, "1", 10, MPFR_RNDN);
        mpfr_init(r26307);
        mpfr_init(r26308);
        mpfr_init(r26309);
        mpfr_init(r26310);
        mpfr_init_set_str(r26311, "0.25", 10, MPFR_RNDN);
        mpfr_init(r26312);
        mpfr_init_set_str(r26313, "0.5", 10, MPFR_RNDN);
        mpfr_init(r26314);
        mpfr_init(r26315);
        mpfr_init(r26316);
        mpfr_init(r26317);
        mpfr_init(r26318);
        mpfr_init(r26319);
}

double f_fm(double alpha, double beta) {
        mpfr_set_d(r26291, beta, MPFR_RNDN);
        mpfr_set_d(r26292, alpha, MPFR_RNDN);
        mpfr_mul(r26293, r26291, r26292, MPFR_RNDN);
        mpfr_add(r26294, r26291, r26292, MPFR_RNDN);
        mpfr_add(r26295, r26293, r26294, MPFR_RNDN);
        ;
        mpfr_add(r26297, r26295, r26296, MPFR_RNDN);
        ;
        mpfr_add(r26299, r26298, r26294, MPFR_RNDN);
        mpfr_div(r26300, r26297, r26299, MPFR_RNDN);
        mpfr_div(r26301, r26300, r26299, MPFR_RNDN);
        mpfr_add(r26302, r26296, r26299, MPFR_RNDN);
        mpfr_div(r26303, r26301, r26302, MPFR_RNDN);
        ;
        mpfr_set_si(r26305, mpfr_cmp(r26303, r26304) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r26307, r26306, r26299, MPFR_RNDN);
        mpfr_mul(r26308, r26297, r26307, MPFR_RNDN);
        mpfr_div(r26309, r26308, r26299, MPFR_RNDN);
        mpfr_div(r26310, r26309, r26302, MPFR_RNDN);
        ;
        mpfr_mul(r26312, r26311, r26294, MPFR_RNDN);
        ;
        mpfr_add(r26314, r26312, r26313, MPFR_RNDN);
        mpfr_add(r26315, r26296, r26298, MPFR_RNDN);
        mpfr_add(r26316, r26315, r26294, MPFR_RNDN);
        mpfr_mul(r26317, r26299, r26316, MPFR_RNDN);
        mpfr_div(r26318, r26314, r26317, MPFR_RNDN);
        if (mpfr_get_si(r26305, MPFR_RNDN)) { mpfr_set(r26319, r26310, MPFR_RNDN); } else { mpfr_set(r26319, r26318, MPFR_RNDN); };
        return mpfr_get_d(r26319, MPFR_RNDN);
}

static mpfr_t r26320, r26321, r26322, r26323, r26324, r26325, r26326, r26327, r26328, r26329, r26330, r26331, r26332, r26333, r26334, r26335, r26336, r26337, r26338, r26339, r26340, r26341, r26342, r26343, r26344, r26345, r26346, r26347, r26348;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26320);
        mpfr_init(r26321);
        mpfr_init(r26322);
        mpfr_init(r26323);
        mpfr_init(r26324);
        mpfr_init_set_str(r26325, "1.0", 10, MPFR_RNDN);
        mpfr_init(r26326);
        mpfr_init_set_str(r26327, "2", 10, MPFR_RNDN);
        mpfr_init(r26328);
        mpfr_init(r26329);
        mpfr_init(r26330);
        mpfr_init(r26331);
        mpfr_init(r26332);
        mpfr_init_set_str(r26333, "1.9688390490121056e+227", 10, MPFR_RNDN);
        mpfr_init(r26334);
        mpfr_init_set_str(r26335, "1", 10, MPFR_RNDN);
        mpfr_init(r26336);
        mpfr_init(r26337);
        mpfr_init(r26338);
        mpfr_init(r26339);
        mpfr_init_set_str(r26340, "0.25", 10, MPFR_RNDN);
        mpfr_init(r26341);
        mpfr_init_set_str(r26342, "0.5", 10, MPFR_RNDN);
        mpfr_init(r26343);
        mpfr_init(r26344);
        mpfr_init(r26345);
        mpfr_init(r26346);
        mpfr_init(r26347);
        mpfr_init(r26348);
}

double f_dm(double alpha, double beta) {
        mpfr_set_d(r26320, beta, MPFR_RNDN);
        mpfr_set_d(r26321, alpha, MPFR_RNDN);
        mpfr_mul(r26322, r26320, r26321, MPFR_RNDN);
        mpfr_add(r26323, r26320, r26321, MPFR_RNDN);
        mpfr_add(r26324, r26322, r26323, MPFR_RNDN);
        ;
        mpfr_add(r26326, r26324, r26325, MPFR_RNDN);
        ;
        mpfr_add(r26328, r26327, r26323, MPFR_RNDN);
        mpfr_div(r26329, r26326, r26328, MPFR_RNDN);
        mpfr_div(r26330, r26329, r26328, MPFR_RNDN);
        mpfr_add(r26331, r26325, r26328, MPFR_RNDN);
        mpfr_div(r26332, r26330, r26331, MPFR_RNDN);
        ;
        mpfr_set_si(r26334, mpfr_cmp(r26332, r26333) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r26336, r26335, r26328, MPFR_RNDN);
        mpfr_mul(r26337, r26326, r26336, MPFR_RNDN);
        mpfr_div(r26338, r26337, r26328, MPFR_RNDN);
        mpfr_div(r26339, r26338, r26331, MPFR_RNDN);
        ;
        mpfr_mul(r26341, r26340, r26323, MPFR_RNDN);
        ;
        mpfr_add(r26343, r26341, r26342, MPFR_RNDN);
        mpfr_add(r26344, r26325, r26327, MPFR_RNDN);
        mpfr_add(r26345, r26344, r26323, MPFR_RNDN);
        mpfr_mul(r26346, r26328, r26345, MPFR_RNDN);
        mpfr_div(r26347, r26343, r26346, MPFR_RNDN);
        if (mpfr_get_si(r26334, MPFR_RNDN)) { mpfr_set(r26348, r26339, MPFR_RNDN); } else { mpfr_set(r26348, r26347, MPFR_RNDN); };
        return mpfr_get_d(r26348, MPFR_RNDN);
}