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

char *name = "Numeric.SpecFunctions:logGammaL from math-functions-0.1.5.2";

double f_if(float x, float y, float z, float t, float a) {
        float r9199 = x;
        float r9200 = y;
        float r9201 = r9199 + r9200;
        float r9202 = log(r9201);
        float r9203 = z;
        float r9204 = log(r9203);
        float r9205 = r9202 + r9204;
        float r9206 = t;
        float r9207 = r9205 - r9206;
        float r9208 = a;
        float r9209 = 0.5;
        float r9210 = r9208 - r9209;
        float r9211 = log(r9206);
        float r9212 = r9210 * r9211;
        float r9213 = r9207 + r9212;
        return r9213;
}

double f_id(double x, double y, double z, double t, double a) {
        double r9214 = x;
        double r9215 = y;
        double r9216 = r9214 + r9215;
        double r9217 = log(r9216);
        double r9218 = z;
        double r9219 = log(r9218);
        double r9220 = r9217 + r9219;
        double r9221 = t;
        double r9222 = r9220 - r9221;
        double r9223 = a;
        double r9224 = 0.5;
        double r9225 = r9223 - r9224;
        double r9226 = log(r9221);
        double r9227 = r9225 * r9226;
        double r9228 = r9222 + r9227;
        return r9228;
}


double f_of(float x, float y, float z, float t, float a) {
        float r9229 = x;
        float r9230 = y;
        float r9231 = r9229 + r9230;
        float r9232 = log(r9231);
        float r9233 = z;
        float r9234 = log(r9233);
        float r9235 = r9232 + r9234;
        float r9236 = t;
        float r9237 = r9235 - r9236;
        float r9238 = a;
        float r9239 = 0.5;
        float r9240 = r9238 - r9239;
        float r9241 = log(r9236);
        float r9242 = r9240 * r9241;
        float r9243 = r9237 + r9242;
        return r9243;
}

double f_od(double x, double y, double z, double t, double a) {
        double r9244 = x;
        double r9245 = y;
        double r9246 = r9244 + r9245;
        double r9247 = log(r9246);
        double r9248 = z;
        double r9249 = log(r9248);
        double r9250 = r9247 + r9249;
        double r9251 = t;
        double r9252 = r9250 - r9251;
        double r9253 = a;
        double r9254 = 0.5;
        double r9255 = r9253 - r9254;
        double r9256 = log(r9251);
        double r9257 = r9255 * r9256;
        double r9258 = r9252 + r9257;
        return r9258;
}

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 r9259, r9260, r9261, r9262, r9263, r9264, r9265, r9266, r9267, r9268, r9269, r9270, r9271, r9272, r9273;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r9259);
        mpfr_init(r9260);
        mpfr_init(r9261);
        mpfr_init(r9262);
        mpfr_init(r9263);
        mpfr_init(r9264);
        mpfr_init(r9265);
        mpfr_init(r9266);
        mpfr_init(r9267);
        mpfr_init(r9268);
        mpfr_init_set_str(r9269, "0.5", 10, MPFR_RNDN);
        mpfr_init(r9270);
        mpfr_init(r9271);
        mpfr_init(r9272);
        mpfr_init(r9273);
}

double f_im(double x, double y, double z, double t, double a) {
        mpfr_set_d(r9259, x, MPFR_RNDN);
        mpfr_set_d(r9260, y, MPFR_RNDN);
        mpfr_add(r9261, r9259, r9260, MPFR_RNDN);
        mpfr_log(r9262, r9261, MPFR_RNDN);
        mpfr_set_d(r9263, z, MPFR_RNDN);
        mpfr_log(r9264, r9263, MPFR_RNDN);
        mpfr_add(r9265, r9262, r9264, MPFR_RNDN);
        mpfr_set_d(r9266, t, MPFR_RNDN);
        mpfr_sub(r9267, r9265, r9266, MPFR_RNDN);
        mpfr_set_d(r9268, a, MPFR_RNDN);
        ;
        mpfr_sub(r9270, r9268, r9269, MPFR_RNDN);
        mpfr_log(r9271, r9266, MPFR_RNDN);
        mpfr_mul(r9272, r9270, r9271, MPFR_RNDN);
        mpfr_add(r9273, r9267, r9272, MPFR_RNDN);
        return mpfr_get_d(r9273, MPFR_RNDN);
}

static mpfr_t r9274, r9275, r9276, r9277, r9278, r9279, r9280, r9281, r9282, r9283, r9284, r9285, r9286, r9287, r9288;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r9274);
        mpfr_init(r9275);
        mpfr_init(r9276);
        mpfr_init(r9277);
        mpfr_init(r9278);
        mpfr_init(r9279);
        mpfr_init(r9280);
        mpfr_init(r9281);
        mpfr_init(r9282);
        mpfr_init(r9283);
        mpfr_init_set_str(r9284, "0.5", 10, MPFR_RNDN);
        mpfr_init(r9285);
        mpfr_init(r9286);
        mpfr_init(r9287);
        mpfr_init(r9288);
}

double f_fm(double x, double y, double z, double t, double a) {
        mpfr_set_d(r9274, x, MPFR_RNDN);
        mpfr_set_d(r9275, y, MPFR_RNDN);
        mpfr_add(r9276, r9274, r9275, MPFR_RNDN);
        mpfr_log(r9277, r9276, MPFR_RNDN);
        mpfr_set_d(r9278, z, MPFR_RNDN);
        mpfr_log(r9279, r9278, MPFR_RNDN);
        mpfr_add(r9280, r9277, r9279, MPFR_RNDN);
        mpfr_set_d(r9281, t, MPFR_RNDN);
        mpfr_sub(r9282, r9280, r9281, MPFR_RNDN);
        mpfr_set_d(r9283, a, MPFR_RNDN);
        ;
        mpfr_sub(r9285, r9283, r9284, MPFR_RNDN);
        mpfr_log(r9286, r9281, MPFR_RNDN);
        mpfr_mul(r9287, r9285, r9286, MPFR_RNDN);
        mpfr_add(r9288, r9282, r9287, MPFR_RNDN);
        return mpfr_get_d(r9288, MPFR_RNDN);
}

static mpfr_t r9289, r9290, r9291, r9292, r9293, r9294, r9295, r9296, r9297, r9298, r9299, r9300, r9301, r9302, r9303;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r9289);
        mpfr_init(r9290);
        mpfr_init(r9291);
        mpfr_init(r9292);
        mpfr_init(r9293);
        mpfr_init(r9294);
        mpfr_init(r9295);
        mpfr_init(r9296);
        mpfr_init(r9297);
        mpfr_init(r9298);
        mpfr_init_set_str(r9299, "0.5", 10, MPFR_RNDN);
        mpfr_init(r9300);
        mpfr_init(r9301);
        mpfr_init(r9302);
        mpfr_init(r9303);
}

double f_dm(double x, double y, double z, double t, double a) {
        mpfr_set_d(r9289, x, MPFR_RNDN);
        mpfr_set_d(r9290, y, MPFR_RNDN);
        mpfr_add(r9291, r9289, r9290, MPFR_RNDN);
        mpfr_log(r9292, r9291, MPFR_RNDN);
        mpfr_set_d(r9293, z, MPFR_RNDN);
        mpfr_log(r9294, r9293, MPFR_RNDN);
        mpfr_add(r9295, r9292, r9294, MPFR_RNDN);
        mpfr_set_d(r9296, t, MPFR_RNDN);
        mpfr_sub(r9297, r9295, r9296, MPFR_RNDN);
        mpfr_set_d(r9298, a, MPFR_RNDN);
        ;
        mpfr_sub(r9300, r9298, r9299, MPFR_RNDN);
        mpfr_log(r9301, r9296, MPFR_RNDN);
        mpfr_mul(r9302, r9300, r9301, MPFR_RNDN);
        mpfr_add(r9303, r9297, r9302, MPFR_RNDN);
        return mpfr_get_d(r9303, MPFR_RNDN);
}