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

char *name = "Data.Colour.Matrix:inverse from colour-2.3.3, B";

double f_if(float x, float y, float z, float t, float a) {
        float r56909 = x;
        float r56910 = y;
        float r56911 = r56909 * r56910;
        float r56912 = z;
        float r56913 = t;
        float r56914 = r56912 * r56913;
        float r56915 = r56911 - r56914;
        float r56916 = a;
        float r56917 = r56915 / r56916;
        return r56917;
}

double f_id(double x, double y, double z, double t, double a) {
        double r56918 = x;
        double r56919 = y;
        double r56920 = r56918 * r56919;
        double r56921 = z;
        double r56922 = t;
        double r56923 = r56921 * r56922;
        double r56924 = r56920 - r56923;
        double r56925 = a;
        double r56926 = r56924 / r56925;
        return r56926;
}


double f_of(float x, float y, float z, float t, float a) {
        float r56927 = a;
        float r56928 = 6.625574720452708e+86;
        bool r56929 = r56927 <= r56928;
        float r56930 = x;
        float r56931 = y;
        float r56932 = r56930 * r56931;
        float r56933 = z;
        float r56934 = t;
        float r56935 = r56933 * r56934;
        float r56936 = r56932 - r56935;
        float r56937 = 1;
        float r56938 = r56937 / r56927;
        float r56939 = r56936 * r56938;
        float r56940 = 2.7445378437622284e+195;
        bool r56941 = r56927 <= r56940;
        float r56942 = -1;
        float r56943 = cbrt(r56942);
        float r56944 = r56933 / r56927;
        float r56945 = r56934 * r56944;
        float r56946 = r56931 / r56927;
        float r56947 = r56930 * r56946;
        float r56948 = r56945 - r56947;
        float r56949 = r56943 * r56948;
        float r56950 = r56941 ? r56949 : r56939;
        float r56951 = r56929 ? r56939 : r56950;
        return r56951;
}

double f_od(double x, double y, double z, double t, double a) {
        double r56952 = a;
        double r56953 = 6.625574720452708e+86;
        bool r56954 = r56952 <= r56953;
        double r56955 = x;
        double r56956 = y;
        double r56957 = r56955 * r56956;
        double r56958 = z;
        double r56959 = t;
        double r56960 = r56958 * r56959;
        double r56961 = r56957 - r56960;
        double r56962 = 1;
        double r56963 = r56962 / r56952;
        double r56964 = r56961 * r56963;
        double r56965 = 2.7445378437622284e+195;
        bool r56966 = r56952 <= r56965;
        double r56967 = -1;
        double r56968 = cbrt(r56967);
        double r56969 = r56958 / r56952;
        double r56970 = r56959 * r56969;
        double r56971 = r56956 / r56952;
        double r56972 = r56955 * r56971;
        double r56973 = r56970 - r56972;
        double r56974 = r56968 * r56973;
        double r56975 = r56966 ? r56974 : r56964;
        double r56976 = r56954 ? r56964 : r56975;
        return r56976;
}

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 r56977, r56978, r56979, r56980, r56981, r56982, r56983, r56984, r56985;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r56977);
        mpfr_init(r56978);
        mpfr_init(r56979);
        mpfr_init(r56980);
        mpfr_init(r56981);
        mpfr_init(r56982);
        mpfr_init(r56983);
        mpfr_init(r56984);
        mpfr_init(r56985);
}

double f_im(double x, double y, double z, double t, double a) {
        mpfr_set_d(r56977, x, MPFR_RNDN);
        mpfr_set_d(r56978, y, MPFR_RNDN);
        mpfr_mul(r56979, r56977, r56978, MPFR_RNDN);
        mpfr_set_d(r56980, z, MPFR_RNDN);
        mpfr_set_d(r56981, t, MPFR_RNDN);
        mpfr_mul(r56982, r56980, r56981, MPFR_RNDN);
        mpfr_sub(r56983, r56979, r56982, MPFR_RNDN);
        mpfr_set_d(r56984, a, MPFR_RNDN);
        mpfr_div(r56985, r56983, r56984, MPFR_RNDN);
        return mpfr_get_d(r56985, MPFR_RNDN);
}

static mpfr_t r56986, r56987, r56988, r56989, r56990, r56991, r56992, r56993, r56994, r56995, r56996, r56997, r56998, r56999, r57000, r57001, r57002, r57003, r57004, r57005, r57006, r57007, r57008, r57009, r57010;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r56986);
        mpfr_init_set_str(r56987, "6.625574720452708e+86", 10, MPFR_RNDN);
        mpfr_init(r56988);
        mpfr_init(r56989);
        mpfr_init(r56990);
        mpfr_init(r56991);
        mpfr_init(r56992);
        mpfr_init(r56993);
        mpfr_init(r56994);
        mpfr_init(r56995);
        mpfr_init_set_str(r56996, "1", 10, MPFR_RNDN);
        mpfr_init(r56997);
        mpfr_init(r56998);
        mpfr_init_set_str(r56999, "2.7445378437622284e+195", 10, MPFR_RNDN);
        mpfr_init(r57000);
        mpfr_init_set_str(r57001, "-1", 10, MPFR_RNDN);
        mpfr_init(r57002);
        mpfr_init(r57003);
        mpfr_init(r57004);
        mpfr_init(r57005);
        mpfr_init(r57006);
        mpfr_init(r57007);
        mpfr_init(r57008);
        mpfr_init(r57009);
        mpfr_init(r57010);
}

double f_fm(double x, double y, double z, double t, double a) {
        mpfr_set_d(r56986, a, MPFR_RNDN);
        ;
        mpfr_set_si(r56988, mpfr_cmp(r56986, r56987) <= 0, MPFR_RNDN);
        mpfr_set_d(r56989, x, MPFR_RNDN);
        mpfr_set_d(r56990, y, MPFR_RNDN);
        mpfr_mul(r56991, r56989, r56990, MPFR_RNDN);
        mpfr_set_d(r56992, z, MPFR_RNDN);
        mpfr_set_d(r56993, t, MPFR_RNDN);
        mpfr_mul(r56994, r56992, r56993, MPFR_RNDN);
        mpfr_sub(r56995, r56991, r56994, MPFR_RNDN);
        ;
        mpfr_div(r56997, r56996, r56986, MPFR_RNDN);
        mpfr_mul(r56998, r56995, r56997, MPFR_RNDN);
        ;
        mpfr_set_si(r57000, mpfr_cmp(r56986, r56999) <= 0, MPFR_RNDN);
        ;
        mpfr_cbrt(r57002, r57001, MPFR_RNDN);
        mpfr_div(r57003, r56992, r56986, MPFR_RNDN);
        mpfr_mul(r57004, r56993, r57003, MPFR_RNDN);
        mpfr_div(r57005, r56990, r56986, MPFR_RNDN);
        mpfr_mul(r57006, r56989, r57005, MPFR_RNDN);
        mpfr_sub(r57007, r57004, r57006, MPFR_RNDN);
        mpfr_mul(r57008, r57002, r57007, MPFR_RNDN);
        if (mpfr_get_si(r57000, MPFR_RNDN)) { mpfr_set(r57009, r57008, MPFR_RNDN); } else { mpfr_set(r57009, r56998, MPFR_RNDN); };
        if (mpfr_get_si(r56988, MPFR_RNDN)) { mpfr_set(r57010, r56998, MPFR_RNDN); } else { mpfr_set(r57010, r57009, MPFR_RNDN); };
        return mpfr_get_d(r57010, MPFR_RNDN);
}

static mpfr_t r57011, r57012, r57013, r57014, r57015, r57016, r57017, r57018, r57019, r57020, r57021, r57022, r57023, r57024, r57025, r57026, r57027, r57028, r57029, r57030, r57031, r57032, r57033, r57034, r57035;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r57011);
        mpfr_init_set_str(r57012, "6.625574720452708e+86", 10, MPFR_RNDN);
        mpfr_init(r57013);
        mpfr_init(r57014);
        mpfr_init(r57015);
        mpfr_init(r57016);
        mpfr_init(r57017);
        mpfr_init(r57018);
        mpfr_init(r57019);
        mpfr_init(r57020);
        mpfr_init_set_str(r57021, "1", 10, MPFR_RNDN);
        mpfr_init(r57022);
        mpfr_init(r57023);
        mpfr_init_set_str(r57024, "2.7445378437622284e+195", 10, MPFR_RNDN);
        mpfr_init(r57025);
        mpfr_init_set_str(r57026, "-1", 10, MPFR_RNDN);
        mpfr_init(r57027);
        mpfr_init(r57028);
        mpfr_init(r57029);
        mpfr_init(r57030);
        mpfr_init(r57031);
        mpfr_init(r57032);
        mpfr_init(r57033);
        mpfr_init(r57034);
        mpfr_init(r57035);
}

double f_dm(double x, double y, double z, double t, double a) {
        mpfr_set_d(r57011, a, MPFR_RNDN);
        ;
        mpfr_set_si(r57013, mpfr_cmp(r57011, r57012) <= 0, MPFR_RNDN);
        mpfr_set_d(r57014, x, MPFR_RNDN);
        mpfr_set_d(r57015, y, MPFR_RNDN);
        mpfr_mul(r57016, r57014, r57015, MPFR_RNDN);
        mpfr_set_d(r57017, z, MPFR_RNDN);
        mpfr_set_d(r57018, t, MPFR_RNDN);
        mpfr_mul(r57019, r57017, r57018, MPFR_RNDN);
        mpfr_sub(r57020, r57016, r57019, MPFR_RNDN);
        ;
        mpfr_div(r57022, r57021, r57011, MPFR_RNDN);
        mpfr_mul(r57023, r57020, r57022, MPFR_RNDN);
        ;
        mpfr_set_si(r57025, mpfr_cmp(r57011, r57024) <= 0, MPFR_RNDN);
        ;
        mpfr_cbrt(r57027, r57026, MPFR_RNDN);
        mpfr_div(r57028, r57017, r57011, MPFR_RNDN);
        mpfr_mul(r57029, r57018, r57028, MPFR_RNDN);
        mpfr_div(r57030, r57015, r57011, MPFR_RNDN);
        mpfr_mul(r57031, r57014, r57030, MPFR_RNDN);
        mpfr_sub(r57032, r57029, r57031, MPFR_RNDN);
        mpfr_mul(r57033, r57027, r57032, MPFR_RNDN);
        if (mpfr_get_si(r57025, MPFR_RNDN)) { mpfr_set(r57034, r57033, MPFR_RNDN); } else { mpfr_set(r57034, r57023, MPFR_RNDN); };
        if (mpfr_get_si(r57013, MPFR_RNDN)) { mpfr_set(r57035, r57023, MPFR_RNDN); } else { mpfr_set(r57035, r57034, MPFR_RNDN); };
        return mpfr_get_d(r57035, MPFR_RNDN);
}