#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 r21959 = x;
        float r21960 = y;
        float r21961 = r21959 + r21960;
        float r21962 = log(r21961);
        float r21963 = z;
        float r21964 = log(r21963);
        float r21965 = r21962 + r21964;
        float r21966 = t;
        float r21967 = r21965 - r21966;
        float r21968 = a;
        float r21969 = 0.5;
        float r21970 = r21968 - r21969;
        float r21971 = log(r21966);
        float r21972 = r21970 * r21971;
        float r21973 = r21967 + r21972;
        return r21973;
}

double f_id(double x, double y, double z, double t, double a) {
        double r21974 = x;
        double r21975 = y;
        double r21976 = r21974 + r21975;
        double r21977 = log(r21976);
        double r21978 = z;
        double r21979 = log(r21978);
        double r21980 = r21977 + r21979;
        double r21981 = t;
        double r21982 = r21980 - r21981;
        double r21983 = a;
        double r21984 = 0.5;
        double r21985 = r21983 - r21984;
        double r21986 = log(r21981);
        double r21987 = r21985 * r21986;
        double r21988 = r21982 + r21987;
        return r21988;
}


double f_of(float x, float y, float z, float t, float a) {
        float r21989 = z;
        float r21990 = log(r21989);
        float r21991 = y;
        float r21992 = x;
        float r21993 = r21991 + r21992;
        float r21994 = log(r21993);
        float r21995 = r21990 + r21994;
        float r21996 = 3;
        float r21997 = pow(r21995, r21996);
        float r21998 = cbrt(r21997);
        float r21999 = t;
        float r22000 = r21998 - r21999;
        float r22001 = a;
        float r22002 = 0.5;
        float r22003 = r22001 - r22002;
        float r22004 = log(r21999);
        float r22005 = r22003 * r22004;
        float r22006 = r22000 + r22005;
        return r22006;
}

double f_od(double x, double y, double z, double t, double a) {
        double r22007 = z;
        double r22008 = log(r22007);
        double r22009 = y;
        double r22010 = x;
        double r22011 = r22009 + r22010;
        double r22012 = log(r22011);
        double r22013 = r22008 + r22012;
        double r22014 = 3;
        double r22015 = pow(r22013, r22014);
        double r22016 = cbrt(r22015);
        double r22017 = t;
        double r22018 = r22016 - r22017;
        double r22019 = a;
        double r22020 = 0.5;
        double r22021 = r22019 - r22020;
        double r22022 = log(r22017);
        double r22023 = r22021 * r22022;
        double r22024 = r22018 + r22023;
        return r22024;
}

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 r22025, r22026, r22027, r22028, r22029, r22030, r22031, r22032, r22033, r22034, r22035, r22036, r22037, r22038, r22039;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r22025);
        mpfr_init(r22026);
        mpfr_init(r22027);
        mpfr_init(r22028);
        mpfr_init(r22029);
        mpfr_init(r22030);
        mpfr_init(r22031);
        mpfr_init(r22032);
        mpfr_init(r22033);
        mpfr_init(r22034);
        mpfr_init_set_str(r22035, "0.5", 10, MPFR_RNDN);
        mpfr_init(r22036);
        mpfr_init(r22037);
        mpfr_init(r22038);
        mpfr_init(r22039);
}

double f_im(double x, double y, double z, double t, double a) {
        mpfr_set_d(r22025, x, MPFR_RNDN);
        mpfr_set_d(r22026, y, MPFR_RNDN);
        mpfr_add(r22027, r22025, r22026, MPFR_RNDN);
        mpfr_log(r22028, r22027, MPFR_RNDN);
        mpfr_set_d(r22029, z, MPFR_RNDN);
        mpfr_log(r22030, r22029, MPFR_RNDN);
        mpfr_add(r22031, r22028, r22030, MPFR_RNDN);
        mpfr_set_d(r22032, t, MPFR_RNDN);
        mpfr_sub(r22033, r22031, r22032, MPFR_RNDN);
        mpfr_set_d(r22034, a, MPFR_RNDN);
        ;
        mpfr_sub(r22036, r22034, r22035, MPFR_RNDN);
        mpfr_log(r22037, r22032, MPFR_RNDN);
        mpfr_mul(r22038, r22036, r22037, MPFR_RNDN);
        mpfr_add(r22039, r22033, r22038, MPFR_RNDN);
        return mpfr_get_d(r22039, MPFR_RNDN);
}

static mpfr_t r22040, r22041, r22042, r22043, r22044, r22045, r22046, r22047, r22048, r22049, r22050, r22051, r22052, r22053, r22054, r22055, r22056, r22057;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r22040);
        mpfr_init(r22041);
        mpfr_init(r22042);
        mpfr_init(r22043);
        mpfr_init(r22044);
        mpfr_init(r22045);
        mpfr_init(r22046);
        mpfr_init_set_str(r22047, "3", 10, MPFR_RNDN);
        mpfr_init(r22048);
        mpfr_init(r22049);
        mpfr_init(r22050);
        mpfr_init(r22051);
        mpfr_init(r22052);
        mpfr_init_set_str(r22053, "0.5", 10, MPFR_RNDN);
        mpfr_init(r22054);
        mpfr_init(r22055);
        mpfr_init(r22056);
        mpfr_init(r22057);
}

double f_fm(double x, double y, double z, double t, double a) {
        mpfr_set_d(r22040, z, MPFR_RNDN);
        mpfr_log(r22041, r22040, MPFR_RNDN);
        mpfr_set_d(r22042, y, MPFR_RNDN);
        mpfr_set_d(r22043, x, MPFR_RNDN);
        mpfr_add(r22044, r22042, r22043, MPFR_RNDN);
        mpfr_log(r22045, r22044, MPFR_RNDN);
        mpfr_add(r22046, r22041, r22045, MPFR_RNDN);
        ;
        mpfr_pow(r22048, r22046, r22047, MPFR_RNDN);
        mpfr_cbrt(r22049, r22048, MPFR_RNDN);
        mpfr_set_d(r22050, t, MPFR_RNDN);
        mpfr_sub(r22051, r22049, r22050, MPFR_RNDN);
        mpfr_set_d(r22052, a, MPFR_RNDN);
        ;
        mpfr_sub(r22054, r22052, r22053, MPFR_RNDN);
        mpfr_log(r22055, r22050, MPFR_RNDN);
        mpfr_mul(r22056, r22054, r22055, MPFR_RNDN);
        mpfr_add(r22057, r22051, r22056, MPFR_RNDN);
        return mpfr_get_d(r22057, MPFR_RNDN);
}

static mpfr_t r22058, r22059, r22060, r22061, r22062, r22063, r22064, r22065, r22066, r22067, r22068, r22069, r22070, r22071, r22072, r22073, r22074, r22075;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r22058);
        mpfr_init(r22059);
        mpfr_init(r22060);
        mpfr_init(r22061);
        mpfr_init(r22062);
        mpfr_init(r22063);
        mpfr_init(r22064);
        mpfr_init_set_str(r22065, "3", 10, MPFR_RNDN);
        mpfr_init(r22066);
        mpfr_init(r22067);
        mpfr_init(r22068);
        mpfr_init(r22069);
        mpfr_init(r22070);
        mpfr_init_set_str(r22071, "0.5", 10, MPFR_RNDN);
        mpfr_init(r22072);
        mpfr_init(r22073);
        mpfr_init(r22074);
        mpfr_init(r22075);
}

double f_dm(double x, double y, double z, double t, double a) {
        mpfr_set_d(r22058, z, MPFR_RNDN);
        mpfr_log(r22059, r22058, MPFR_RNDN);
        mpfr_set_d(r22060, y, MPFR_RNDN);
        mpfr_set_d(r22061, x, MPFR_RNDN);
        mpfr_add(r22062, r22060, r22061, MPFR_RNDN);
        mpfr_log(r22063, r22062, MPFR_RNDN);
        mpfr_add(r22064, r22059, r22063, MPFR_RNDN);
        ;
        mpfr_pow(r22066, r22064, r22065, MPFR_RNDN);
        mpfr_cbrt(r22067, r22066, MPFR_RNDN);
        mpfr_set_d(r22068, t, MPFR_RNDN);
        mpfr_sub(r22069, r22067, r22068, MPFR_RNDN);
        mpfr_set_d(r22070, a, MPFR_RNDN);
        ;
        mpfr_sub(r22072, r22070, r22071, MPFR_RNDN);
        mpfr_log(r22073, r22068, MPFR_RNDN);
        mpfr_mul(r22074, r22072, r22073, MPFR_RNDN);
        mpfr_add(r22075, r22069, r22074, MPFR_RNDN);
        return mpfr_get_d(r22075, MPFR_RNDN);
}