#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 r10821 = x; float r10822 = y; float r10823 = z; float r10824 = r10822 * r10823; float r10825 = t; float r10826 = a; float r10827 = r10825 * r10826; float r10828 = r10824 - r10827; float r10829 = r10821 * r10828; float r10830 = b; float r10831 = c; float r10832 = r10831 * r10823; float r10833 = i; float r10834 = r10833 * r10826; float r10835 = r10832 - r10834; float r10836 = r10830 * r10835; float r10837 = r10829 - r10836; float r10838 = j; float r10839 = r10831 * r10825; float r10840 = r10833 * r10822; float r10841 = r10839 - r10840; float r10842 = r10838 * r10841; float r10843 = r10837 + r10842; return r10843; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j) { double r10844 = x; double r10845 = y; double r10846 = z; double r10847 = r10845 * r10846; double r10848 = t; double r10849 = a; double r10850 = r10848 * r10849; double r10851 = r10847 - r10850; double r10852 = r10844 * r10851; double r10853 = b; double r10854 = c; double r10855 = r10854 * r10846; double r10856 = i; double r10857 = r10856 * r10849; double r10858 = r10855 - r10857; double r10859 = r10853 * r10858; double r10860 = r10852 - r10859; double r10861 = j; double r10862 = r10854 * r10848; double r10863 = r10856 * r10845; double r10864 = r10862 - r10863; double r10865 = r10861 * r10864; double r10866 = r10860 + r10865; return r10866; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j) { float r10867 = x; float r10868 = z; float r10869 = y; float r10870 = r10868 * r10869; float r10871 = t; float r10872 = a; float r10873 = r10871 * r10872; float r10874 = r10870 - r10873; float r10875 = c; float r10876 = j; float r10877 = r10876 * r10871; float r10878 = r10875 * r10877; float r10879 = fma(r10867, r10874, r10878); float r10880 = r10875 * r10868; float r10881 = i; float r10882 = r10872 * r10881; float r10883 = r10880 - r10882; float r10884 = b; float r10885 = r10876 * r10869; float r10886 = r10885 * r10881; float r10887 = fma(r10883, r10884, r10886); float r10888 = r10879 - r10887; float r10889 = -1.9612421800185972e+298; bool r10890 = r10888 <= r10889; float r10891 = r10871 * r10875; float r10892 = r10869 * r10881; float r10893 = r10891 - r10892; float r10894 = r10893 * r10876; float r10895 = r10884 * r10875; float r10896 = r10868 * r10895; float r10897 = r10881 * r10884; float r10898 = r10872 * r10897; float r10899 = r10896 - r10898; float r10900 = r10894 - r10899; float r10901 = fma(r10874, r10867, r10900); float r10902 = 7.055696161566341e+305; bool r10903 = r10888 <= r10902; float r10904 = r10868 * r10867; float r10905 = r10869 * r10904; float r10906 = r10875 * r10871; float r10907 = r10881 * r10869; float r10908 = r10906 - r10907; float r10909 = r10876 * r10908; float r10910 = r10881 * r10872; float r10911 = r10880 - r10910; float r10912 = r10872 * r10867; float r10913 = r10912 * r10871; float r10914 = fma(r10884, r10911, r10913); float r10915 = r10909 - r10914; float r10916 = r10905 + r10915; float r10917 = r10903 ? r10888 : r10916; float r10918 = r10890 ? r10901 : r10917; return r10918; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j) { double r10919 = x; double r10920 = z; double r10921 = y; double r10922 = r10920 * r10921; double r10923 = t; double r10924 = a; double r10925 = r10923 * r10924; double r10926 = r10922 - r10925; double r10927 = c; double r10928 = j; double r10929 = r10928 * r10923; double r10930 = r10927 * r10929; double r10931 = fma(r10919, r10926, r10930); double r10932 = r10927 * r10920; double r10933 = i; double r10934 = r10924 * r10933; double r10935 = r10932 - r10934; double r10936 = b; double r10937 = r10928 * r10921; double r10938 = r10937 * r10933; double r10939 = fma(r10935, r10936, r10938); double r10940 = r10931 - r10939; double r10941 = -1.9612421800185972e+298; bool r10942 = r10940 <= r10941; double r10943 = r10923 * r10927; double r10944 = r10921 * r10933; double r10945 = r10943 - r10944; double r10946 = r10945 * r10928; double r10947 = r10936 * r10927; double r10948 = r10920 * r10947; double r10949 = r10933 * r10936; double r10950 = r10924 * r10949; double r10951 = r10948 - r10950; double r10952 = r10946 - r10951; double r10953 = fma(r10926, r10919, r10952); double r10954 = 7.055696161566341e+305; bool r10955 = r10940 <= r10954; double r10956 = r10920 * r10919; double r10957 = r10921 * r10956; double r10958 = r10927 * r10923; double r10959 = r10933 * r10921; double r10960 = r10958 - r10959; double r10961 = r10928 * r10960; double r10962 = r10933 * r10924; double r10963 = r10932 - r10962; double r10964 = r10924 * r10919; double r10965 = r10964 * r10923; double r10966 = fma(r10936, r10963, r10965); double r10967 = r10961 - r10966; double r10968 = r10957 + r10967; double r10969 = r10955 ? r10940 : r10968; double r10970 = r10942 ? r10953 : r10969; return r10970; } 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 r10971, r10972, r10973, r10974, r10975, r10976, r10977, r10978, r10979, r10980, r10981, r10982, r10983, r10984, r10985, r10986, r10987, r10988, r10989, r10990, r10991, r10992, r10993; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r10971); mpfr_init(r10972); mpfr_init(r10973); mpfr_init(r10974); mpfr_init(r10975); mpfr_init(r10976); mpfr_init(r10977); mpfr_init(r10978); mpfr_init(r10979); mpfr_init(r10980); mpfr_init(r10981); mpfr_init(r10982); mpfr_init(r10983); mpfr_init(r10984); mpfr_init(r10985); mpfr_init(r10986); mpfr_init(r10987); mpfr_init(r10988); mpfr_init(r10989); mpfr_init(r10990); mpfr_init(r10991); mpfr_init(r10992); mpfr_init(r10993); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r10971, x, MPFR_RNDN); mpfr_set_d(r10972, y, MPFR_RNDN); mpfr_set_d(r10973, z, MPFR_RNDN); mpfr_mul(r10974, r10972, r10973, MPFR_RNDN); mpfr_set_d(r10975, t, MPFR_RNDN); mpfr_set_d(r10976, a, MPFR_RNDN); mpfr_mul(r10977, r10975, r10976, MPFR_RNDN); mpfr_sub(r10978, r10974, r10977, MPFR_RNDN); mpfr_mul(r10979, r10971, r10978, MPFR_RNDN); mpfr_set_d(r10980, b, MPFR_RNDN); mpfr_set_d(r10981, c, MPFR_RNDN); mpfr_mul(r10982, r10981, r10973, MPFR_RNDN); mpfr_set_d(r10983, i, MPFR_RNDN); mpfr_mul(r10984, r10983, r10976, MPFR_RNDN); mpfr_sub(r10985, r10982, r10984, MPFR_RNDN); mpfr_mul(r10986, r10980, r10985, MPFR_RNDN); mpfr_sub(r10987, r10979, r10986, MPFR_RNDN); mpfr_set_d(r10988, j, MPFR_RNDN); mpfr_mul(r10989, r10981, r10975, MPFR_RNDN); mpfr_mul(r10990, r10983, r10972, MPFR_RNDN); mpfr_sub(r10991, r10989, r10990, MPFR_RNDN); mpfr_mul(r10992, r10988, r10991, MPFR_RNDN); mpfr_add(r10993, r10987, r10992, MPFR_RNDN); return mpfr_get_d(r10993, MPFR_RNDN); } static mpfr_t r10994, r10995, r10996, r10997, r10998, r10999, r11000, r11001, r11002, r11003, r11004, r11005, r11006, r11007, r11008, r11009, r11010, r11011, r11012, r11013, r11014, r11015, r11016, r11017, r11018, r11019, r11020, r11021, r11022, r11023, r11024, r11025, r11026, r11027, r11028, r11029, r11030, r11031, r11032, r11033, r11034, r11035, r11036, r11037, r11038, r11039, r11040, r11041, r11042, r11043, r11044, r11045; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r10994); mpfr_init(r10995); mpfr_init(r10996); mpfr_init(r10997); mpfr_init(r10998); mpfr_init(r10999); mpfr_init(r11000); mpfr_init(r11001); mpfr_init(r11002); mpfr_init(r11003); mpfr_init(r11004); mpfr_init(r11005); mpfr_init(r11006); 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_set_str(r11016, "-1.9612421800185972e+298", 10, MPFR_RNDN); 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_set_str(r11029, "7.055696161566341e+305", 10, MPFR_RNDN); 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); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r10994, x, MPFR_RNDN); mpfr_set_d(r10995, z, MPFR_RNDN); mpfr_set_d(r10996, y, MPFR_RNDN); mpfr_mul(r10997, r10995, r10996, MPFR_RNDN); mpfr_set_d(r10998, t, MPFR_RNDN); mpfr_set_d(r10999, a, MPFR_RNDN); mpfr_mul(r11000, r10998, r10999, MPFR_RNDN); mpfr_sub(r11001, r10997, r11000, MPFR_RNDN); mpfr_set_d(r11002, c, MPFR_RNDN); mpfr_set_d(r11003, j, MPFR_RNDN); mpfr_mul(r11004, r11003, r10998, MPFR_RNDN); mpfr_mul(r11005, r11002, r11004, MPFR_RNDN); mpfr_fma(r11006, r10994, r11001, r11005, MPFR_RNDN); mpfr_mul(r11007, r11002, r10995, MPFR_RNDN); mpfr_set_d(r11008, i, MPFR_RNDN); mpfr_mul(r11009, r10999, r11008, MPFR_RNDN); mpfr_sub(r11010, r11007, r11009, MPFR_RNDN); mpfr_set_d(r11011, b, MPFR_RNDN); mpfr_mul(r11012, r11003, r10996, MPFR_RNDN); mpfr_mul(r11013, r11012, r11008, MPFR_RNDN); mpfr_fma(r11014, r11010, r11011, r11013, MPFR_RNDN); mpfr_sub(r11015, r11006, r11014, MPFR_RNDN); ; mpfr_set_si(r11017, mpfr_cmp(r11015, r11016) <= 0, MPFR_RNDN); mpfr_mul(r11018, r10998, r11002, MPFR_RNDN); mpfr_mul(r11019, r10996, r11008, MPFR_RNDN); mpfr_sub(r11020, r11018, r11019, MPFR_RNDN); mpfr_mul(r11021, r11020, r11003, MPFR_RNDN); mpfr_mul(r11022, r11011, r11002, MPFR_RNDN); mpfr_mul(r11023, r10995, r11022, MPFR_RNDN); mpfr_mul(r11024, r11008, r11011, MPFR_RNDN); mpfr_mul(r11025, r10999, r11024, MPFR_RNDN); mpfr_sub(r11026, r11023, r11025, MPFR_RNDN); mpfr_sub(r11027, r11021, r11026, MPFR_RNDN); mpfr_fma(r11028, r11001, r10994, r11027, MPFR_RNDN); ; mpfr_set_si(r11030, mpfr_cmp(r11015, r11029) <= 0, MPFR_RNDN); mpfr_mul(r11031, r10995, r10994, MPFR_RNDN); mpfr_mul(r11032, r10996, r11031, MPFR_RNDN); mpfr_mul(r11033, r11002, r10998, MPFR_RNDN); mpfr_mul(r11034, r11008, r10996, MPFR_RNDN); mpfr_sub(r11035, r11033, r11034, MPFR_RNDN); mpfr_mul(r11036, r11003, r11035, MPFR_RNDN); mpfr_mul(r11037, r11008, r10999, MPFR_RNDN); mpfr_sub(r11038, r11007, r11037, MPFR_RNDN); mpfr_mul(r11039, r10999, r10994, MPFR_RNDN); mpfr_mul(r11040, r11039, r10998, MPFR_RNDN); mpfr_fma(r11041, r11011, r11038, r11040, MPFR_RNDN); mpfr_sub(r11042, r11036, r11041, MPFR_RNDN); mpfr_add(r11043, r11032, r11042, MPFR_RNDN); if (mpfr_get_si(r11030, MPFR_RNDN)) { mpfr_set(r11044, r11015, MPFR_RNDN); } else { mpfr_set(r11044, r11043, MPFR_RNDN); }; if (mpfr_get_si(r11017, MPFR_RNDN)) { mpfr_set(r11045, r11028, MPFR_RNDN); } else { mpfr_set(r11045, r11044, MPFR_RNDN); }; return mpfr_get_d(r11045, MPFR_RNDN); } static mpfr_t 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, r11074, r11075, r11076, r11077, r11078, r11079, r11080, r11081, r11082, r11083, r11084, r11085, r11086, r11087, r11088, r11089, r11090, r11091, r11092, r11093, r11094, r11095, r11096, r11097; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r11046); mpfr_init(r11047); mpfr_init(r11048); mpfr_init(r11049); mpfr_init(r11050); mpfr_init(r11051); mpfr_init(r11052); mpfr_init(r11053); mpfr_init(r11054); mpfr_init(r11055); 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_set_str(r11068, "-1.9612421800185972e+298", 10, MPFR_RNDN); mpfr_init(r11069); mpfr_init(r11070); mpfr_init(r11071); mpfr_init(r11072); mpfr_init(r11073); mpfr_init(r11074); mpfr_init(r11075); mpfr_init(r11076); mpfr_init(r11077); mpfr_init(r11078); mpfr_init(r11079); mpfr_init(r11080); mpfr_init_set_str(r11081, "7.055696161566341e+305", 10, MPFR_RNDN); 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(r11097); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r11046, x, MPFR_RNDN); mpfr_set_d(r11047, z, MPFR_RNDN); mpfr_set_d(r11048, y, MPFR_RNDN); mpfr_mul(r11049, r11047, r11048, MPFR_RNDN); mpfr_set_d(r11050, t, MPFR_RNDN); mpfr_set_d(r11051, a, MPFR_RNDN); mpfr_mul(r11052, r11050, r11051, MPFR_RNDN); mpfr_sub(r11053, r11049, r11052, MPFR_RNDN); mpfr_set_d(r11054, c, MPFR_RNDN); mpfr_set_d(r11055, j, MPFR_RNDN); mpfr_mul(r11056, r11055, r11050, MPFR_RNDN); mpfr_mul(r11057, r11054, r11056, MPFR_RNDN); mpfr_fma(r11058, r11046, r11053, r11057, MPFR_RNDN); mpfr_mul(r11059, r11054, r11047, MPFR_RNDN); mpfr_set_d(r11060, i, MPFR_RNDN); mpfr_mul(r11061, r11051, r11060, MPFR_RNDN); mpfr_sub(r11062, r11059, r11061, MPFR_RNDN); mpfr_set_d(r11063, b, MPFR_RNDN); mpfr_mul(r11064, r11055, r11048, MPFR_RNDN); mpfr_mul(r11065, r11064, r11060, MPFR_RNDN); mpfr_fma(r11066, r11062, r11063, r11065, MPFR_RNDN); mpfr_sub(r11067, r11058, r11066, MPFR_RNDN); ; mpfr_set_si(r11069, mpfr_cmp(r11067, r11068) <= 0, MPFR_RNDN); mpfr_mul(r11070, r11050, r11054, MPFR_RNDN); mpfr_mul(r11071, r11048, r11060, MPFR_RNDN); mpfr_sub(r11072, r11070, r11071, MPFR_RNDN); mpfr_mul(r11073, r11072, r11055, MPFR_RNDN); mpfr_mul(r11074, r11063, r11054, MPFR_RNDN); mpfr_mul(r11075, r11047, r11074, MPFR_RNDN); mpfr_mul(r11076, r11060, r11063, MPFR_RNDN); mpfr_mul(r11077, r11051, r11076, MPFR_RNDN); mpfr_sub(r11078, r11075, r11077, MPFR_RNDN); mpfr_sub(r11079, r11073, r11078, MPFR_RNDN); mpfr_fma(r11080, r11053, r11046, r11079, MPFR_RNDN); ; mpfr_set_si(r11082, mpfr_cmp(r11067, r11081) <= 0, MPFR_RNDN); mpfr_mul(r11083, r11047, r11046, MPFR_RNDN); mpfr_mul(r11084, r11048, r11083, MPFR_RNDN); mpfr_mul(r11085, r11054, r11050, MPFR_RNDN); mpfr_mul(r11086, r11060, r11048, MPFR_RNDN); mpfr_sub(r11087, r11085, r11086, MPFR_RNDN); mpfr_mul(r11088, r11055, r11087, MPFR_RNDN); mpfr_mul(r11089, r11060, r11051, MPFR_RNDN); mpfr_sub(r11090, r11059, r11089, MPFR_RNDN); mpfr_mul(r11091, r11051, r11046, MPFR_RNDN); mpfr_mul(r11092, r11091, r11050, MPFR_RNDN); mpfr_fma(r11093, r11063, r11090, r11092, MPFR_RNDN); mpfr_sub(r11094, r11088, r11093, MPFR_RNDN); mpfr_add(r11095, r11084, r11094, MPFR_RNDN); if (mpfr_get_si(r11082, MPFR_RNDN)) { mpfr_set(r11096, r11067, MPFR_RNDN); } else { mpfr_set(r11096, r11095, MPFR_RNDN); }; if (mpfr_get_si(r11069, MPFR_RNDN)) { mpfr_set(r11097, r11080, MPFR_RNDN); } else { mpfr_set(r11097, r11096, MPFR_RNDN); }; return mpfr_get_d(r11097, MPFR_RNDN); }