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

char *name = "Linear.Matrix:det33 from linear-1.19.1.3";

double f_if(float x, float y, float z, float t, float a, float b, float c, float i, float j) {
        float r10873 = x;
        float r10874 = y;
        float r10875 = z;
        float r10876 = r10874 * r10875;
        float r10877 = t;
        float r10878 = a;
        float r10879 = r10877 * r10878;
        float r10880 = r10876 - r10879;
        float r10881 = r10873 * r10880;
        float r10882 = b;
        float r10883 = c;
        float r10884 = r10883 * r10875;
        float r10885 = i;
        float r10886 = r10885 * r10878;
        float r10887 = r10884 - r10886;
        float r10888 = r10882 * r10887;
        float r10889 = r10881 - r10888;
        float r10890 = j;
        float r10891 = r10883 * r10877;
        float r10892 = r10885 * r10874;
        float r10893 = r10891 - r10892;
        float r10894 = r10890 * r10893;
        float r10895 = r10889 + r10894;
        return r10895;
}

double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
        double r10896 = x;
        double r10897 = y;
        double r10898 = z;
        double r10899 = r10897 * r10898;
        double r10900 = t;
        double r10901 = a;
        double r10902 = r10900 * r10901;
        double r10903 = r10899 - r10902;
        double r10904 = r10896 * r10903;
        double r10905 = b;
        double r10906 = c;
        double r10907 = r10906 * r10898;
        double r10908 = i;
        double r10909 = r10908 * r10901;
        double r10910 = r10907 - r10909;
        double r10911 = r10905 * r10910;
        double r10912 = r10904 - r10911;
        double r10913 = j;
        double r10914 = r10906 * r10900;
        double r10915 = r10908 * r10897;
        double r10916 = r10914 - r10915;
        double r10917 = r10913 * r10916;
        double r10918 = r10912 + r10917;
        return r10918;
}


double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j) {
        float r10919 = t;
        float r10920 = c;
        float r10921 = r10919 * r10920;
        float r10922 = y;
        float r10923 = i;
        float r10924 = r10922 * r10923;
        float r10925 = r10921 - r10924;
        float r10926 = j;
        float r10927 = r10925 * r10926;
        float r10928 = z;
        float r10929 = r10928 * r10922;
        float r10930 = a;
        float r10931 = r10919 * r10930;
        float r10932 = r10929 - r10931;
        float r10933 = x;
        float r10934 = r10932 * r10933;
        float r10935 = r10928 * r10920;
        float r10936 = r10923 * r10930;
        float r10937 = r10935 - r10936;
        float r10938 = b;
        float r10939 = r10937 * r10938;
        float r10940 = r10934 - r10939;
        float r10941 = r10927 + r10940;
        float r10942 = -7.233071798759434e+304;
        bool r10943 = r10941 <= r10942;
        float r10944 = 1.7836295536525312e+308;
        bool r10945 = r10941 <= r10944;
        float r10946 = !r10945;
        bool r10947 = r10943 || r10946;
        float r10948 = r10933 * r10922;
        float r10949 = r10938 * r10920;
        float r10950 = r10948 - r10949;
        float r10951 = r10928 * r10950;
        float r10952 = r10938 * r10923;
        float r10953 = r10919 * r10933;
        float r10954 = r10952 - r10953;
        float r10955 = r10954 * r10930;
        float r10956 = r10951 + r10955;
        float r10957 = cbrt(r10934);
        float r10958 = r10957 * r10957;
        float r10959 = r10957 * r10958;
        float r10960 = r10959 - r10939;
        float r10961 = r10960 + r10927;
        float r10962 = r10947 ? r10956 : r10961;
        return r10962;
}

double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
        double r10963 = t;
        double r10964 = c;
        double r10965 = r10963 * r10964;
        double r10966 = y;
        double r10967 = i;
        double r10968 = r10966 * r10967;
        double r10969 = r10965 - r10968;
        double r10970 = j;
        double r10971 = r10969 * r10970;
        double r10972 = z;
        double r10973 = r10972 * r10966;
        double r10974 = a;
        double r10975 = r10963 * r10974;
        double r10976 = r10973 - r10975;
        double r10977 = x;
        double r10978 = r10976 * r10977;
        double r10979 = r10972 * r10964;
        double r10980 = r10967 * r10974;
        double r10981 = r10979 - r10980;
        double r10982 = b;
        double r10983 = r10981 * r10982;
        double r10984 = r10978 - r10983;
        double r10985 = r10971 + r10984;
        double r10986 = -7.233071798759434e+304;
        bool r10987 = r10985 <= r10986;
        double r10988 = 1.7836295536525312e+308;
        bool r10989 = r10985 <= r10988;
        double r10990 = !r10989;
        bool r10991 = r10987 || r10990;
        double r10992 = r10977 * r10966;
        double r10993 = r10982 * r10964;
        double r10994 = r10992 - r10993;
        double r10995 = r10972 * r10994;
        double r10996 = r10982 * r10967;
        double r10997 = r10963 * r10977;
        double r10998 = r10996 - r10997;
        double r10999 = r10998 * r10974;
        double r11000 = r10995 + r10999;
        double r11001 = cbrt(r10978);
        double r11002 = r11001 * r11001;
        double r11003 = r11001 * r11002;
        double r11004 = r11003 - r10983;
        double r11005 = r11004 + r10971;
        double r11006 = r10991 ? r11000 : r11005;
        return r11006;
}

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 r11007, r11008, r11009, r11010, r11011, r11012, r11013, r11014, r11015, r11016, r11017, r11018, r11019, r11020, r11021, r11022, r11023, r11024, r11025, r11026, r11027, r11028, r11029;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r11007);
        mpfr_init(r11008);
        mpfr_init(r11009);
        mpfr_init(r11010);
        mpfr_init(r11011);
        mpfr_init(r11012);
        mpfr_init(r11013);
        mpfr_init(r11014);
        mpfr_init(r11015);
        mpfr_init(r11016);
        mpfr_init(r11017);
        mpfr_init(r11018);
        mpfr_init(r11019);
        mpfr_init(r11020);
        mpfr_init(r11021);
        mpfr_init(r11022);
        mpfr_init(r11023);
        mpfr_init(r11024);
        mpfr_init(r11025);
        mpfr_init(r11026);
        mpfr_init(r11027);
        mpfr_init(r11028);
        mpfr_init(r11029);
}

double f_im(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
        mpfr_set_d(r11007, x, MPFR_RNDN);
        mpfr_set_d(r11008, y, MPFR_RNDN);
        mpfr_set_d(r11009, z, MPFR_RNDN);
        mpfr_mul(r11010, r11008, r11009, MPFR_RNDN);
        mpfr_set_d(r11011, t, MPFR_RNDN);
        mpfr_set_d(r11012, a, MPFR_RNDN);
        mpfr_mul(r11013, r11011, r11012, MPFR_RNDN);
        mpfr_sub(r11014, r11010, r11013, MPFR_RNDN);
        mpfr_mul(r11015, r11007, r11014, MPFR_RNDN);
        mpfr_set_d(r11016, b, MPFR_RNDN);
        mpfr_set_d(r11017, c, MPFR_RNDN);
        mpfr_mul(r11018, r11017, r11009, MPFR_RNDN);
        mpfr_set_d(r11019, i, MPFR_RNDN);
        mpfr_mul(r11020, r11019, r11012, MPFR_RNDN);
        mpfr_sub(r11021, r11018, r11020, MPFR_RNDN);
        mpfr_mul(r11022, r11016, r11021, MPFR_RNDN);
        mpfr_sub(r11023, r11015, r11022, MPFR_RNDN);
        mpfr_set_d(r11024, j, MPFR_RNDN);
        mpfr_mul(r11025, r11017, r11011, MPFR_RNDN);
        mpfr_mul(r11026, r11019, r11008, MPFR_RNDN);
        mpfr_sub(r11027, r11025, r11026, MPFR_RNDN);
        mpfr_mul(r11028, r11024, r11027, MPFR_RNDN);
        mpfr_add(r11029, r11023, r11028, MPFR_RNDN);
        return mpfr_get_d(r11029, MPFR_RNDN);
}

static mpfr_t r11030, r11031, r11032, r11033, r11034, r11035, r11036, r11037, r11038, r11039, r11040, r11041, r11042, r11043, r11044, r11045, r11046, r11047, r11048, r11049, r11050, r11051, r11052, r11053, r11054, r11055, r11056, r11057, r11058, r11059, r11060, r11061, r11062, r11063, r11064, r11065, r11066, r11067, r11068, r11069, r11070, r11071, r11072, r11073;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r11030);
        mpfr_init(r11031);
        mpfr_init(r11032);
        mpfr_init(r11033);
        mpfr_init(r11034);
        mpfr_init(r11035);
        mpfr_init(r11036);
        mpfr_init(r11037);
        mpfr_init(r11038);
        mpfr_init(r11039);
        mpfr_init(r11040);
        mpfr_init(r11041);
        mpfr_init(r11042);
        mpfr_init(r11043);
        mpfr_init(r11044);
        mpfr_init(r11045);
        mpfr_init(r11046);
        mpfr_init(r11047);
        mpfr_init(r11048);
        mpfr_init(r11049);
        mpfr_init(r11050);
        mpfr_init(r11051);
        mpfr_init(r11052);
        mpfr_init_set_str(r11053, "-7.233071798759434e+304", 10, MPFR_RNDN);
        mpfr_init(r11054);
        mpfr_init_set_str(r11055, "1.7836295536525312e+308", 10, MPFR_RNDN);
        mpfr_init(r11056);
        mpfr_init(r11057);
        mpfr_init(r11058);
        mpfr_init(r11059);
        mpfr_init(r11060);
        mpfr_init(r11061);
        mpfr_init(r11062);
        mpfr_init(r11063);
        mpfr_init(r11064);
        mpfr_init(r11065);
        mpfr_init(r11066);
        mpfr_init(r11067);
        mpfr_init(r11068);
        mpfr_init(r11069);
        mpfr_init(r11070);
        mpfr_init(r11071);
        mpfr_init(r11072);
        mpfr_init(r11073);
}

double f_fm(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
        mpfr_set_d(r11030, t, MPFR_RNDN);
        mpfr_set_d(r11031, c, MPFR_RNDN);
        mpfr_mul(r11032, r11030, r11031, MPFR_RNDN);
        mpfr_set_d(r11033, y, MPFR_RNDN);
        mpfr_set_d(r11034, i, MPFR_RNDN);
        mpfr_mul(r11035, r11033, r11034, MPFR_RNDN);
        mpfr_sub(r11036, r11032, r11035, MPFR_RNDN);
        mpfr_set_d(r11037, j, MPFR_RNDN);
        mpfr_mul(r11038, r11036, r11037, MPFR_RNDN);
        mpfr_set_d(r11039, z, MPFR_RNDN);
        mpfr_mul(r11040, r11039, r11033, MPFR_RNDN);
        mpfr_set_d(r11041, a, MPFR_RNDN);
        mpfr_mul(r11042, r11030, r11041, MPFR_RNDN);
        mpfr_sub(r11043, r11040, r11042, MPFR_RNDN);
        mpfr_set_d(r11044, x, MPFR_RNDN);
        mpfr_mul(r11045, r11043, r11044, MPFR_RNDN);
        mpfr_mul(r11046, r11039, r11031, MPFR_RNDN);
        mpfr_mul(r11047, r11034, r11041, MPFR_RNDN);
        mpfr_sub(r11048, r11046, r11047, MPFR_RNDN);
        mpfr_set_d(r11049, b, MPFR_RNDN);
        mpfr_mul(r11050, r11048, r11049, MPFR_RNDN);
        mpfr_sub(r11051, r11045, r11050, MPFR_RNDN);
        mpfr_add(r11052, r11038, r11051, MPFR_RNDN);
        ;
        mpfr_set_si(r11054, mpfr_cmp(r11052, r11053) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r11056, mpfr_cmp(r11052, r11055) <= 0, MPFR_RNDN);
        mpfr_set_si(r11057, !mpfr_get_si(r11056, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r11058, mpfr_get_si(r11054, MPFR_RNDN) || mpfr_get_si(r11057, MPFR_RNDN), MPFR_RNDN);
        mpfr_mul(r11059, r11044, r11033, MPFR_RNDN);
        mpfr_mul(r11060, r11049, r11031, MPFR_RNDN);
        mpfr_sub(r11061, r11059, r11060, MPFR_RNDN);
        mpfr_mul(r11062, r11039, r11061, MPFR_RNDN);
        mpfr_mul(r11063, r11049, r11034, MPFR_RNDN);
        mpfr_mul(r11064, r11030, r11044, MPFR_RNDN);
        mpfr_sub(r11065, r11063, r11064, MPFR_RNDN);
        mpfr_mul(r11066, r11065, r11041, MPFR_RNDN);
        mpfr_add(r11067, r11062, r11066, MPFR_RNDN);
        mpfr_cbrt(r11068, r11045, MPFR_RNDN);
        mpfr_mul(r11069, r11068, r11068, MPFR_RNDN);
        mpfr_mul(r11070, r11068, r11069, MPFR_RNDN);
        mpfr_sub(r11071, r11070, r11050, MPFR_RNDN);
        mpfr_add(r11072, r11071, r11038, MPFR_RNDN);
        if (mpfr_get_si(r11058, MPFR_RNDN)) { mpfr_set(r11073, r11067, MPFR_RNDN); } else { mpfr_set(r11073, r11072, MPFR_RNDN); };
        return mpfr_get_d(r11073, MPFR_RNDN);
}

static mpfr_t r11074, r11075, r11076, r11077, r11078, r11079, r11080, r11081, r11082, r11083, r11084, r11085, r11086, r11087, r11088, r11089, r11090, r11091, r11092, r11093, r11094, r11095, r11096, r11097, r11098, r11099, r11100, r11101, r11102, r11103, r11104, r11105, r11106, r11107, r11108, r11109, r11110, r11111, r11112, r11113, r11114, r11115, r11116, r11117;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r11074);
        mpfr_init(r11075);
        mpfr_init(r11076);
        mpfr_init(r11077);
        mpfr_init(r11078);
        mpfr_init(r11079);
        mpfr_init(r11080);
        mpfr_init(r11081);
        mpfr_init(r11082);
        mpfr_init(r11083);
        mpfr_init(r11084);
        mpfr_init(r11085);
        mpfr_init(r11086);
        mpfr_init(r11087);
        mpfr_init(r11088);
        mpfr_init(r11089);
        mpfr_init(r11090);
        mpfr_init(r11091);
        mpfr_init(r11092);
        mpfr_init(r11093);
        mpfr_init(r11094);
        mpfr_init(r11095);
        mpfr_init(r11096);
        mpfr_init_set_str(r11097, "-7.233071798759434e+304", 10, MPFR_RNDN);
        mpfr_init(r11098);
        mpfr_init_set_str(r11099, "1.7836295536525312e+308", 10, MPFR_RNDN);
        mpfr_init(r11100);
        mpfr_init(r11101);
        mpfr_init(r11102);
        mpfr_init(r11103);
        mpfr_init(r11104);
        mpfr_init(r11105);
        mpfr_init(r11106);
        mpfr_init(r11107);
        mpfr_init(r11108);
        mpfr_init(r11109);
        mpfr_init(r11110);
        mpfr_init(r11111);
        mpfr_init(r11112);
        mpfr_init(r11113);
        mpfr_init(r11114);
        mpfr_init(r11115);
        mpfr_init(r11116);
        mpfr_init(r11117);
}

double f_dm(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
        mpfr_set_d(r11074, t, MPFR_RNDN);
        mpfr_set_d(r11075, c, MPFR_RNDN);
        mpfr_mul(r11076, r11074, r11075, MPFR_RNDN);
        mpfr_set_d(r11077, y, MPFR_RNDN);
        mpfr_set_d(r11078, i, MPFR_RNDN);
        mpfr_mul(r11079, r11077, r11078, MPFR_RNDN);
        mpfr_sub(r11080, r11076, r11079, MPFR_RNDN);
        mpfr_set_d(r11081, j, MPFR_RNDN);
        mpfr_mul(r11082, r11080, r11081, MPFR_RNDN);
        mpfr_set_d(r11083, z, MPFR_RNDN);
        mpfr_mul(r11084, r11083, r11077, MPFR_RNDN);
        mpfr_set_d(r11085, a, MPFR_RNDN);
        mpfr_mul(r11086, r11074, r11085, MPFR_RNDN);
        mpfr_sub(r11087, r11084, r11086, MPFR_RNDN);
        mpfr_set_d(r11088, x, MPFR_RNDN);
        mpfr_mul(r11089, r11087, r11088, MPFR_RNDN);
        mpfr_mul(r11090, r11083, r11075, MPFR_RNDN);
        mpfr_mul(r11091, r11078, r11085, MPFR_RNDN);
        mpfr_sub(r11092, r11090, r11091, MPFR_RNDN);
        mpfr_set_d(r11093, b, MPFR_RNDN);
        mpfr_mul(r11094, r11092, r11093, MPFR_RNDN);
        mpfr_sub(r11095, r11089, r11094, MPFR_RNDN);
        mpfr_add(r11096, r11082, r11095, MPFR_RNDN);
        ;
        mpfr_set_si(r11098, mpfr_cmp(r11096, r11097) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r11100, mpfr_cmp(r11096, r11099) <= 0, MPFR_RNDN);
        mpfr_set_si(r11101, !mpfr_get_si(r11100, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r11102, mpfr_get_si(r11098, MPFR_RNDN) || mpfr_get_si(r11101, MPFR_RNDN), MPFR_RNDN);
        mpfr_mul(r11103, r11088, r11077, MPFR_RNDN);
        mpfr_mul(r11104, r11093, r11075, MPFR_RNDN);
        mpfr_sub(r11105, r11103, r11104, MPFR_RNDN);
        mpfr_mul(r11106, r11083, r11105, MPFR_RNDN);
        mpfr_mul(r11107, r11093, r11078, MPFR_RNDN);
        mpfr_mul(r11108, r11074, r11088, MPFR_RNDN);
        mpfr_sub(r11109, r11107, r11108, MPFR_RNDN);
        mpfr_mul(r11110, r11109, r11085, MPFR_RNDN);
        mpfr_add(r11111, r11106, r11110, MPFR_RNDN);
        mpfr_cbrt(r11112, r11089, MPFR_RNDN);
        mpfr_mul(r11113, r11112, r11112, MPFR_RNDN);
        mpfr_mul(r11114, r11112, r11113, MPFR_RNDN);
        mpfr_sub(r11115, r11114, r11094, MPFR_RNDN);
        mpfr_add(r11116, r11115, r11082, MPFR_RNDN);
        if (mpfr_get_si(r11102, MPFR_RNDN)) { mpfr_set(r11117, r11111, MPFR_RNDN); } else { mpfr_set(r11117, r11116, MPFR_RNDN); };
        return mpfr_get_d(r11117, MPFR_RNDN);
}