#include #include #include #include #include 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 r17033 = x; float r17034 = y; float r17035 = z; float r17036 = r17034 * r17035; float r17037 = t; float r17038 = a; float r17039 = r17037 * r17038; float r17040 = r17036 - r17039; float r17041 = r17033 * r17040; float r17042 = b; float r17043 = c; float r17044 = r17043 * r17035; float r17045 = i; float r17046 = r17045 * r17038; float r17047 = r17044 - r17046; float r17048 = r17042 * r17047; float r17049 = r17041 - r17048; float r17050 = j; float r17051 = r17043 * r17037; float r17052 = r17045 * r17034; float r17053 = r17051 - r17052; float r17054 = r17050 * r17053; float r17055 = r17049 + r17054; return r17055; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j) { double r17056 = x; double r17057 = y; double r17058 = z; double r17059 = r17057 * r17058; double r17060 = t; double r17061 = a; double r17062 = r17060 * r17061; double r17063 = r17059 - r17062; double r17064 = r17056 * r17063; double r17065 = b; double r17066 = c; double r17067 = r17066 * r17058; double r17068 = i; double r17069 = r17068 * r17061; double r17070 = r17067 - r17069; double r17071 = r17065 * r17070; double r17072 = r17064 - r17071; double r17073 = j; double r17074 = r17066 * r17060; double r17075 = r17068 * r17057; double r17076 = r17074 - r17075; double r17077 = r17073 * r17076; double r17078 = r17072 + r17077; return r17078; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j) { float r17079 = x; float r17080 = y; float r17081 = z; float r17082 = r17080 * r17081; float r17083 = t; float r17084 = a; float r17085 = r17083 * r17084; float r17086 = r17082 - r17085; float r17087 = r17079 * r17086; float r17088 = b; float r17089 = c; float r17090 = r17089 * r17081; float r17091 = i; float r17092 = r17091 * r17084; float r17093 = r17090 - r17092; float r17094 = r17088 * r17093; float r17095 = r17087 - r17094; float r17096 = j; float r17097 = r17089 * r17083; float r17098 = r17091 * r17080; float r17099 = r17097 - r17098; float r17100 = r17096 * r17099; float r17101 = r17095 + r17100; float r17102 = cbrt(r17101); float r17103 = r17102 * (r17102 * r17102); return r17103; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j) { double r17104 = x; double r17105 = y; double r17106 = z; double r17107 = r17105 * r17106; double r17108 = t; double r17109 = a; double r17110 = r17108 * r17109; double r17111 = r17107 - r17110; double r17112 = r17104 * r17111; double r17113 = b; double r17114 = c; double r17115 = r17114 * r17106; double r17116 = i; double r17117 = r17116 * r17109; double r17118 = r17115 - r17117; double r17119 = r17113 * r17118; double r17120 = r17112 - r17119; double r17121 = j; double r17122 = r17114 * r17108; double r17123 = r17116 * r17105; double r17124 = r17122 - r17123; double r17125 = r17121 * r17124; double r17126 = r17120 + r17125; double r17127 = cbrt(r17126); double r17128 = r17127 * (r17127 * r17127); return r17128; } 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 r17129, r17130, r17131, r17132, r17133, r17134, r17135, r17136, r17137, r17138, r17139, r17140, r17141, r17142, r17143, r17144, r17145, r17146, r17147, r17148, r17149, r17150, r17151; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17129); mpfr_init(r17130); mpfr_init(r17131); mpfr_init(r17132); mpfr_init(r17133); mpfr_init(r17134); mpfr_init(r17135); mpfr_init(r17136); mpfr_init(r17137); mpfr_init(r17138); mpfr_init(r17139); mpfr_init(r17140); mpfr_init(r17141); mpfr_init(r17142); mpfr_init(r17143); mpfr_init(r17144); mpfr_init(r17145); mpfr_init(r17146); mpfr_init(r17147); mpfr_init(r17148); mpfr_init(r17149); mpfr_init(r17150); mpfr_init(r17151); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r17129, x, MPFR_RNDN); mpfr_set_d(r17130, y, MPFR_RNDN); mpfr_set_d(r17131, z, MPFR_RNDN); mpfr_mul(r17132, r17130, r17131, MPFR_RNDN); mpfr_set_d(r17133, t, MPFR_RNDN); mpfr_set_d(r17134, a, MPFR_RNDN); mpfr_mul(r17135, r17133, r17134, MPFR_RNDN); mpfr_sub(r17136, r17132, r17135, MPFR_RNDN); mpfr_mul(r17137, r17129, r17136, MPFR_RNDN); mpfr_set_d(r17138, b, MPFR_RNDN); mpfr_set_d(r17139, c, MPFR_RNDN); mpfr_mul(r17140, r17139, r17131, MPFR_RNDN); mpfr_set_d(r17141, i, MPFR_RNDN); mpfr_mul(r17142, r17141, r17134, MPFR_RNDN); mpfr_sub(r17143, r17140, r17142, MPFR_RNDN); mpfr_mul(r17144, r17138, r17143, MPFR_RNDN); mpfr_sub(r17145, r17137, r17144, MPFR_RNDN); mpfr_set_d(r17146, j, MPFR_RNDN); mpfr_mul(r17147, r17139, r17133, MPFR_RNDN); mpfr_mul(r17148, r17141, r17130, MPFR_RNDN); mpfr_sub(r17149, r17147, r17148, MPFR_RNDN); mpfr_mul(r17150, r17146, r17149, MPFR_RNDN); mpfr_add(r17151, r17145, r17150, MPFR_RNDN); return mpfr_get_d(r17151, MPFR_RNDN); } static mpfr_t r17152, r17153, r17154, r17155, r17156, r17157, r17158, r17159, r17160, r17161, r17162, r17163, r17164, r17165, r17166, r17167, r17168, r17169, r17170, r17171, r17172, r17173, r17174, r17175, r17176; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17152); mpfr_init(r17153); mpfr_init(r17154); mpfr_init(r17155); mpfr_init(r17156); mpfr_init(r17157); mpfr_init(r17158); mpfr_init(r17159); mpfr_init(r17160); mpfr_init(r17161); mpfr_init(r17162); mpfr_init(r17163); mpfr_init(r17164); mpfr_init(r17165); mpfr_init(r17166); mpfr_init(r17167); mpfr_init(r17168); mpfr_init(r17169); mpfr_init(r17170); mpfr_init(r17171); mpfr_init(r17172); mpfr_init(r17173); mpfr_init(r17174); mpfr_init(r17175); mpfr_init(r17176); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r17152, x, MPFR_RNDN); mpfr_set_d(r17153, y, MPFR_RNDN); mpfr_set_d(r17154, z, MPFR_RNDN); mpfr_mul(r17155, r17153, r17154, MPFR_RNDN); mpfr_set_d(r17156, t, MPFR_RNDN); mpfr_set_d(r17157, a, MPFR_RNDN); mpfr_mul(r17158, r17156, r17157, MPFR_RNDN); mpfr_sub(r17159, r17155, r17158, MPFR_RNDN); mpfr_mul(r17160, r17152, r17159, MPFR_RNDN); mpfr_set_d(r17161, b, MPFR_RNDN); mpfr_set_d(r17162, c, MPFR_RNDN); mpfr_mul(r17163, r17162, r17154, MPFR_RNDN); mpfr_set_d(r17164, i, MPFR_RNDN); mpfr_mul(r17165, r17164, r17157, MPFR_RNDN); mpfr_sub(r17166, r17163, r17165, MPFR_RNDN); mpfr_mul(r17167, r17161, r17166, MPFR_RNDN); mpfr_sub(r17168, r17160, r17167, MPFR_RNDN); mpfr_set_d(r17169, j, MPFR_RNDN); mpfr_mul(r17170, r17162, r17156, MPFR_RNDN); mpfr_mul(r17171, r17164, r17153, MPFR_RNDN); mpfr_sub(r17172, r17170, r17171, MPFR_RNDN); mpfr_mul(r17173, r17169, r17172, MPFR_RNDN); mpfr_add(r17174, r17168, r17173, MPFR_RNDN); mpfr_cbrt(r17175, r17174, MPFR_RNDN); mpfr_mul(r17176, r17175, r17175, MPFR_RNDN); mpfr_mul(r17176, r17176, r17175, MPFR_RNDN); return mpfr_get_d(r17176, MPFR_RNDN); } static mpfr_t r17177, r17178, r17179, r17180, r17181, r17182, r17183, r17184, r17185, r17186, r17187, r17188, r17189, r17190, r17191, r17192, r17193, r17194, r17195, r17196, r17197, r17198, r17199, r17200, r17201; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17177); mpfr_init(r17178); mpfr_init(r17179); mpfr_init(r17180); mpfr_init(r17181); mpfr_init(r17182); mpfr_init(r17183); mpfr_init(r17184); mpfr_init(r17185); mpfr_init(r17186); mpfr_init(r17187); mpfr_init(r17188); mpfr_init(r17189); mpfr_init(r17190); mpfr_init(r17191); mpfr_init(r17192); mpfr_init(r17193); mpfr_init(r17194); mpfr_init(r17195); mpfr_init(r17196); mpfr_init(r17197); mpfr_init(r17198); mpfr_init(r17199); mpfr_init(r17200); mpfr_init(r17201); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r17177, x, MPFR_RNDN); mpfr_set_d(r17178, y, MPFR_RNDN); mpfr_set_d(r17179, z, MPFR_RNDN); mpfr_mul(r17180, r17178, r17179, MPFR_RNDN); mpfr_set_d(r17181, t, MPFR_RNDN); mpfr_set_d(r17182, a, MPFR_RNDN); mpfr_mul(r17183, r17181, r17182, MPFR_RNDN); mpfr_sub(r17184, r17180, r17183, MPFR_RNDN); mpfr_mul(r17185, r17177, r17184, MPFR_RNDN); mpfr_set_d(r17186, b, MPFR_RNDN); mpfr_set_d(r17187, c, MPFR_RNDN); mpfr_mul(r17188, r17187, r17179, MPFR_RNDN); mpfr_set_d(r17189, i, MPFR_RNDN); mpfr_mul(r17190, r17189, r17182, MPFR_RNDN); mpfr_sub(r17191, r17188, r17190, MPFR_RNDN); mpfr_mul(r17192, r17186, r17191, MPFR_RNDN); mpfr_sub(r17193, r17185, r17192, MPFR_RNDN); mpfr_set_d(r17194, j, MPFR_RNDN); mpfr_mul(r17195, r17187, r17181, MPFR_RNDN); mpfr_mul(r17196, r17189, r17178, MPFR_RNDN); mpfr_sub(r17197, r17195, r17196, MPFR_RNDN); mpfr_mul(r17198, r17194, r17197, MPFR_RNDN); mpfr_add(r17199, r17193, r17198, MPFR_RNDN); mpfr_cbrt(r17200, r17199, MPFR_RNDN); mpfr_mul(r17201, r17200, r17200, MPFR_RNDN); mpfr_mul(r17201, r17201, r17200, MPFR_RNDN); return mpfr_get_d(r17201, MPFR_RNDN); }