#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 r10130 = x;
        float r10131 = y;
        float r10132 = z;
        float r10133 = r10131 * r10132;
        float r10134 = t;
        float r10135 = a;
        float r10136 = r10134 * r10135;
        float r10137 = r10133 - r10136;
        float r10138 = r10130 * r10137;
        float r10139 = b;
        float r10140 = c;
        float r10141 = r10140 * r10132;
        float r10142 = i;
        float r10143 = r10142 * r10135;
        float r10144 = r10141 - r10143;
        float r10145 = r10139 * r10144;
        float r10146 = r10138 - r10145;
        float r10147 = j;
        float r10148 = r10140 * r10134;
        float r10149 = r10142 * r10131;
        float r10150 = r10148 - r10149;
        float r10151 = r10147 * r10150;
        float r10152 = r10146 + r10151;
        return r10152;
}

double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
        double r10153 = x;
        double r10154 = y;
        double r10155 = z;
        double r10156 = r10154 * r10155;
        double r10157 = t;
        double r10158 = a;
        double r10159 = r10157 * r10158;
        double r10160 = r10156 - r10159;
        double r10161 = r10153 * r10160;
        double r10162 = b;
        double r10163 = c;
        double r10164 = r10163 * r10155;
        double r10165 = i;
        double r10166 = r10165 * r10158;
        double r10167 = r10164 - r10166;
        double r10168 = r10162 * r10167;
        double r10169 = r10161 - r10168;
        double r10170 = j;
        double r10171 = r10163 * r10157;
        double r10172 = r10165 * r10154;
        double r10173 = r10171 - r10172;
        double r10174 = r10170 * r10173;
        double r10175 = r10169 + r10174;
        return r10175;
}


double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j) {
        float r10176 = x;
        float r10177 = -1.0412696191654271e-33f;
        bool r10178 = r10176 <= r10177;
        float r10179 = 5.914453084821963e-159f;
        bool r10180 = r10176 <= r10179;
        float r10181 = !r10180;
        bool r10182 = r10178 || r10181;
        float r10183 = y;
        float r10184 = z;
        float r10185 = r10183 * r10184;
        float r10186 = a;
        float r10187 = t;
        float r10188 = r10186 * r10187;
        float r10189 = r10185 - r10188;
        float r10190 = r10189 * r10176;
        float r10191 = b;
        float r10192 = c;
        float r10193 = r10184 * r10192;
        float r10194 = i;
        float r10195 = r10194 * r10186;
        float r10196 = r10193 - r10195;
        float r10197 = r10191 * r10196;
        float r10198 = r10190 - r10197;
        float r10199 = r10192 * r10187;
        float r10200 = r10194 * r10183;
        float r10201 = r10199 - r10200;
        float r10202 = cbrt(r10201);
        float r10203 = r10202 * r10202;
        float r10204 = j;
        float r10205 = r10203 * r10204;
        float r10206 = r10205 * r10202;
        float r10207 = r10198 + r10206;
        float r10208 = r10195 - r10193;
        float r10209 = r10191 * r10208;
        float r10210 = r10176 * r10187;
        float r10211 = r10210 * r10186;
        float r10212 = r10209 - r10211;
        float r10213 = r10204 * r10201;
        float r10214 = r10212 + r10213;
        float r10215 = r10182 ? r10207 : r10214;
        return r10215;
}

double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
        double r10216 = x;
        double r10217 = -1.0412696191654271e-33;
        bool r10218 = r10216 <= r10217;
        double r10219 = 5.914453084821963e-159;
        bool r10220 = r10216 <= r10219;
        double r10221 = !r10220;
        bool r10222 = r10218 || r10221;
        double r10223 = y;
        double r10224 = z;
        double r10225 = r10223 * r10224;
        double r10226 = a;
        double r10227 = t;
        double r10228 = r10226 * r10227;
        double r10229 = r10225 - r10228;
        double r10230 = r10229 * r10216;
        double r10231 = b;
        double r10232 = c;
        double r10233 = r10224 * r10232;
        double r10234 = i;
        double r10235 = r10234 * r10226;
        double r10236 = r10233 - r10235;
        double r10237 = r10231 * r10236;
        double r10238 = r10230 - r10237;
        double r10239 = r10232 * r10227;
        double r10240 = r10234 * r10223;
        double r10241 = r10239 - r10240;
        double r10242 = cbrt(r10241);
        double r10243 = r10242 * r10242;
        double r10244 = j;
        double r10245 = r10243 * r10244;
        double r10246 = r10245 * r10242;
        double r10247 = r10238 + r10246;
        double r10248 = r10235 - r10233;
        double r10249 = r10231 * r10248;
        double r10250 = r10216 * r10227;
        double r10251 = r10250 * r10226;
        double r10252 = r10249 - r10251;
        double r10253 = r10244 * r10241;
        double r10254 = r10252 + r10253;
        double r10255 = r10222 ? r10247 : r10254;
        return r10255;
}

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 r10256, r10257, r10258, r10259, r10260, r10261, r10262, r10263, r10264, r10265, r10266, r10267, r10268, r10269, r10270, r10271, r10272, r10273, r10274, r10275, r10276, r10277, r10278;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r10256);
        mpfr_init(r10257);
        mpfr_init(r10258);
        mpfr_init(r10259);
        mpfr_init(r10260);
        mpfr_init(r10261);
        mpfr_init(r10262);
        mpfr_init(r10263);
        mpfr_init(r10264);
        mpfr_init(r10265);
        mpfr_init(r10266);
        mpfr_init(r10267);
        mpfr_init(r10268);
        mpfr_init(r10269);
        mpfr_init(r10270);
        mpfr_init(r10271);
        mpfr_init(r10272);
        mpfr_init(r10273);
        mpfr_init(r10274);
        mpfr_init(r10275);
        mpfr_init(r10276);
        mpfr_init(r10277);
        mpfr_init(r10278);
}

double f_im(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
        mpfr_set_d(r10256, x, MPFR_RNDN);
        mpfr_set_d(r10257, y, MPFR_RNDN);
        mpfr_set_d(r10258, z, MPFR_RNDN);
        mpfr_mul(r10259, r10257, r10258, MPFR_RNDN);
        mpfr_set_d(r10260, t, MPFR_RNDN);
        mpfr_set_d(r10261, a, MPFR_RNDN);
        mpfr_mul(r10262, r10260, r10261, MPFR_RNDN);
        mpfr_sub(r10263, r10259, r10262, MPFR_RNDN);
        mpfr_mul(r10264, r10256, r10263, MPFR_RNDN);
        mpfr_set_d(r10265, b, MPFR_RNDN);
        mpfr_set_d(r10266, c, MPFR_RNDN);
        mpfr_mul(r10267, r10266, r10258, MPFR_RNDN);
        mpfr_set_d(r10268, i, MPFR_RNDN);
        mpfr_mul(r10269, r10268, r10261, MPFR_RNDN);
        mpfr_sub(r10270, r10267, r10269, MPFR_RNDN);
        mpfr_mul(r10271, r10265, r10270, MPFR_RNDN);
        mpfr_sub(r10272, r10264, r10271, MPFR_RNDN);
        mpfr_set_d(r10273, j, MPFR_RNDN);
        mpfr_mul(r10274, r10266, r10260, MPFR_RNDN);
        mpfr_mul(r10275, r10268, r10257, MPFR_RNDN);
        mpfr_sub(r10276, r10274, r10275, MPFR_RNDN);
        mpfr_mul(r10277, r10273, r10276, MPFR_RNDN);
        mpfr_add(r10278, r10272, r10277, MPFR_RNDN);
        return mpfr_get_d(r10278, MPFR_RNDN);
}

static mpfr_t r10279, r10280, r10281, r10282, r10283, r10284, r10285, r10286, r10287, r10288, r10289, r10290, r10291, r10292, r10293, r10294, r10295, r10296, r10297, r10298, r10299, r10300, r10301, r10302, r10303, r10304, r10305, r10306, r10307, r10308, r10309, r10310, r10311, r10312, r10313, r10314, r10315, r10316, r10317, r10318;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r10279);
        mpfr_init_set_str(r10280, "-1.0412696191654271e-33", 10, MPFR_RNDN);
        mpfr_init(r10281);
        mpfr_init_set_str(r10282, "5.914453084821963e-159", 10, MPFR_RNDN);
        mpfr_init(r10283);
        mpfr_init(r10284);
        mpfr_init(r10285);
        mpfr_init(r10286);
        mpfr_init(r10287);
        mpfr_init(r10288);
        mpfr_init(r10289);
        mpfr_init(r10290);
        mpfr_init(r10291);
        mpfr_init(r10292);
        mpfr_init(r10293);
        mpfr_init(r10294);
        mpfr_init(r10295);
        mpfr_init(r10296);
        mpfr_init(r10297);
        mpfr_init(r10298);
        mpfr_init(r10299);
        mpfr_init(r10300);
        mpfr_init(r10301);
        mpfr_init(r10302);
        mpfr_init(r10303);
        mpfr_init(r10304);
        mpfr_init(r10305);
        mpfr_init(r10306);
        mpfr_init(r10307);
        mpfr_init(r10308);
        mpfr_init(r10309);
        mpfr_init(r10310);
        mpfr_init(r10311);
        mpfr_init(r10312);
        mpfr_init(r10313);
        mpfr_init(r10314);
        mpfr_init(r10315);
        mpfr_init(r10316);
        mpfr_init(r10317);
        mpfr_init(r10318);
}

double f_fm(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
        mpfr_set_d(r10279, x, MPFR_RNDN);
        ;
        mpfr_set_si(r10281, mpfr_cmp(r10279, r10280) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r10283, mpfr_cmp(r10279, r10282) <= 0, MPFR_RNDN);
        mpfr_set_si(r10284, !mpfr_get_si(r10283, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r10285, mpfr_get_si(r10281, MPFR_RNDN) || mpfr_get_si(r10284, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_d(r10286, y, MPFR_RNDN);
        mpfr_set_d(r10287, z, MPFR_RNDN);
        mpfr_mul(r10288, r10286, r10287, MPFR_RNDN);
        mpfr_set_d(r10289, a, MPFR_RNDN);
        mpfr_set_d(r10290, t, MPFR_RNDN);
        mpfr_mul(r10291, r10289, r10290, MPFR_RNDN);
        mpfr_sub(r10292, r10288, r10291, MPFR_RNDN);
        mpfr_mul(r10293, r10292, r10279, MPFR_RNDN);
        mpfr_set_d(r10294, b, MPFR_RNDN);
        mpfr_set_d(r10295, c, MPFR_RNDN);
        mpfr_mul(r10296, r10287, r10295, MPFR_RNDN);
        mpfr_set_d(r10297, i, MPFR_RNDN);
        mpfr_mul(r10298, r10297, r10289, MPFR_RNDN);
        mpfr_sub(r10299, r10296, r10298, MPFR_RNDN);
        mpfr_mul(r10300, r10294, r10299, MPFR_RNDN);
        mpfr_sub(r10301, r10293, r10300, MPFR_RNDN);
        mpfr_mul(r10302, r10295, r10290, MPFR_RNDN);
        mpfr_mul(r10303, r10297, r10286, MPFR_RNDN);
        mpfr_sub(r10304, r10302, r10303, MPFR_RNDN);
        mpfr_cbrt(r10305, r10304, MPFR_RNDN);
        mpfr_mul(r10306, r10305, r10305, MPFR_RNDN);
        mpfr_set_d(r10307, j, MPFR_RNDN);
        mpfr_mul(r10308, r10306, r10307, MPFR_RNDN);
        mpfr_mul(r10309, r10308, r10305, MPFR_RNDN);
        mpfr_add(r10310, r10301, r10309, MPFR_RNDN);
        mpfr_sub(r10311, r10298, r10296, MPFR_RNDN);
        mpfr_mul(r10312, r10294, r10311, MPFR_RNDN);
        mpfr_mul(r10313, r10279, r10290, MPFR_RNDN);
        mpfr_mul(r10314, r10313, r10289, MPFR_RNDN);
        mpfr_sub(r10315, r10312, r10314, MPFR_RNDN);
        mpfr_mul(r10316, r10307, r10304, MPFR_RNDN);
        mpfr_add(r10317, r10315, r10316, MPFR_RNDN);
        if (mpfr_get_si(r10285, MPFR_RNDN)) { mpfr_set(r10318, r10310, MPFR_RNDN); } else { mpfr_set(r10318, r10317, MPFR_RNDN); };
        return mpfr_get_d(r10318, MPFR_RNDN);
}

static mpfr_t r10319, r10320, r10321, r10322, r10323, r10324, r10325, r10326, r10327, r10328, r10329, r10330, r10331, r10332, r10333, r10334, r10335, r10336, r10337, r10338, r10339, r10340, r10341, r10342, r10343, r10344, r10345, r10346, r10347, r10348, r10349, r10350, r10351, r10352, r10353, r10354, r10355, r10356, r10357, r10358;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r10319);
        mpfr_init_set_str(r10320, "-1.0412696191654271e-33", 10, MPFR_RNDN);
        mpfr_init(r10321);
        mpfr_init_set_str(r10322, "5.914453084821963e-159", 10, MPFR_RNDN);
        mpfr_init(r10323);
        mpfr_init(r10324);
        mpfr_init(r10325);
        mpfr_init(r10326);
        mpfr_init(r10327);
        mpfr_init(r10328);
        mpfr_init(r10329);
        mpfr_init(r10330);
        mpfr_init(r10331);
        mpfr_init(r10332);
        mpfr_init(r10333);
        mpfr_init(r10334);
        mpfr_init(r10335);
        mpfr_init(r10336);
        mpfr_init(r10337);
        mpfr_init(r10338);
        mpfr_init(r10339);
        mpfr_init(r10340);
        mpfr_init(r10341);
        mpfr_init(r10342);
        mpfr_init(r10343);
        mpfr_init(r10344);
        mpfr_init(r10345);
        mpfr_init(r10346);
        mpfr_init(r10347);
        mpfr_init(r10348);
        mpfr_init(r10349);
        mpfr_init(r10350);
        mpfr_init(r10351);
        mpfr_init(r10352);
        mpfr_init(r10353);
        mpfr_init(r10354);
        mpfr_init(r10355);
        mpfr_init(r10356);
        mpfr_init(r10357);
        mpfr_init(r10358);
}

double f_dm(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
        mpfr_set_d(r10319, x, MPFR_RNDN);
        ;
        mpfr_set_si(r10321, mpfr_cmp(r10319, r10320) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r10323, mpfr_cmp(r10319, r10322) <= 0, MPFR_RNDN);
        mpfr_set_si(r10324, !mpfr_get_si(r10323, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r10325, mpfr_get_si(r10321, MPFR_RNDN) || mpfr_get_si(r10324, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_d(r10326, y, MPFR_RNDN);
        mpfr_set_d(r10327, z, MPFR_RNDN);
        mpfr_mul(r10328, r10326, r10327, MPFR_RNDN);
        mpfr_set_d(r10329, a, MPFR_RNDN);
        mpfr_set_d(r10330, t, MPFR_RNDN);
        mpfr_mul(r10331, r10329, r10330, MPFR_RNDN);
        mpfr_sub(r10332, r10328, r10331, MPFR_RNDN);
        mpfr_mul(r10333, r10332, r10319, MPFR_RNDN);
        mpfr_set_d(r10334, b, MPFR_RNDN);
        mpfr_set_d(r10335, c, MPFR_RNDN);
        mpfr_mul(r10336, r10327, r10335, MPFR_RNDN);
        mpfr_set_d(r10337, i, MPFR_RNDN);
        mpfr_mul(r10338, r10337, r10329, MPFR_RNDN);
        mpfr_sub(r10339, r10336, r10338, MPFR_RNDN);
        mpfr_mul(r10340, r10334, r10339, MPFR_RNDN);
        mpfr_sub(r10341, r10333, r10340, MPFR_RNDN);
        mpfr_mul(r10342, r10335, r10330, MPFR_RNDN);
        mpfr_mul(r10343, r10337, r10326, MPFR_RNDN);
        mpfr_sub(r10344, r10342, r10343, MPFR_RNDN);
        mpfr_cbrt(r10345, r10344, MPFR_RNDN);
        mpfr_mul(r10346, r10345, r10345, MPFR_RNDN);
        mpfr_set_d(r10347, j, MPFR_RNDN);
        mpfr_mul(r10348, r10346, r10347, MPFR_RNDN);
        mpfr_mul(r10349, r10348, r10345, MPFR_RNDN);
        mpfr_add(r10350, r10341, r10349, MPFR_RNDN);
        mpfr_sub(r10351, r10338, r10336, MPFR_RNDN);
        mpfr_mul(r10352, r10334, r10351, MPFR_RNDN);
        mpfr_mul(r10353, r10319, r10330, MPFR_RNDN);
        mpfr_mul(r10354, r10353, r10329, MPFR_RNDN);
        mpfr_sub(r10355, r10352, r10354, MPFR_RNDN);
        mpfr_mul(r10356, r10347, r10344, MPFR_RNDN);
        mpfr_add(r10357, r10355, r10356, MPFR_RNDN);
        if (mpfr_get_si(r10325, MPFR_RNDN)) { mpfr_set(r10358, r10350, MPFR_RNDN); } else { mpfr_set(r10358, r10357, MPFR_RNDN); };
        return mpfr_get_d(r10358, MPFR_RNDN);
}