#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 r26643 = x; float r26644 = y; float r26645 = z; float r26646 = r26644 * r26645; float r26647 = t; float r26648 = a; float r26649 = r26647 * r26648; float r26650 = r26646 - r26649; float r26651 = r26643 * r26650; float r26652 = b; float r26653 = c; float r26654 = r26653 * r26645; float r26655 = i; float r26656 = r26655 * r26648; float r26657 = r26654 - r26656; float r26658 = r26652 * r26657; float r26659 = r26651 - r26658; float r26660 = j; float r26661 = r26653 * r26647; float r26662 = r26655 * r26644; float r26663 = r26661 - r26662; float r26664 = r26660 * r26663; float r26665 = r26659 + r26664; return r26665; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j) { double r26666 = x; double r26667 = y; double r26668 = z; double r26669 = r26667 * r26668; double r26670 = t; double r26671 = a; double r26672 = r26670 * r26671; double r26673 = r26669 - r26672; double r26674 = r26666 * r26673; double r26675 = b; double r26676 = c; double r26677 = r26676 * r26668; double r26678 = i; double r26679 = r26678 * r26671; double r26680 = r26677 - r26679; double r26681 = r26675 * r26680; double r26682 = r26674 - r26681; double r26683 = j; double r26684 = r26676 * r26670; double r26685 = r26678 * r26667; double r26686 = r26684 - r26685; double r26687 = r26683 * r26686; double r26688 = r26682 + r26687; return r26688; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j) { float r26689 = x; float r26690 = y; float r26691 = z; float r26692 = r26690 * r26691; float r26693 = t; float r26694 = a; float r26695 = r26693 * r26694; float r26696 = r26692 - r26695; float r26697 = r26689 * r26696; float r26698 = b; float r26699 = c; float r26700 = r26699 * r26691; float r26701 = i; float r26702 = r26701 * r26694; float r26703 = r26700 - r26702; float r26704 = r26698 * r26703; float r26705 = r26697 - r26704; float r26706 = j; float r26707 = r26699 * r26693; float r26708 = r26701 * r26690; float r26709 = r26707 - r26708; float r26710 = r26706 * r26709; float r26711 = cbrt(r26710); float r26712 = r26711 * r26711; float r26713 = cbrt(r26706); float r26714 = cbrt(r26709); float r26715 = r26713 * r26714; float r26716 = r26712 * r26715; float r26717 = r26705 + r26716; return r26717; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j) { double r26718 = x; double r26719 = y; double r26720 = z; double r26721 = r26719 * r26720; double r26722 = t; double r26723 = a; double r26724 = r26722 * r26723; double r26725 = r26721 - r26724; double r26726 = r26718 * r26725; double r26727 = b; double r26728 = c; double r26729 = r26728 * r26720; double r26730 = i; double r26731 = r26730 * r26723; double r26732 = r26729 - r26731; double r26733 = r26727 * r26732; double r26734 = r26726 - r26733; double r26735 = j; double r26736 = r26728 * r26722; double r26737 = r26730 * r26719; double r26738 = r26736 - r26737; double r26739 = r26735 * r26738; double r26740 = cbrt(r26739); double r26741 = r26740 * r26740; double r26742 = cbrt(r26735); double r26743 = cbrt(r26738); double r26744 = r26742 * r26743; double r26745 = r26741 * r26744; double r26746 = r26734 + r26745; return r26746; } 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 r26747, r26748, r26749, r26750, r26751, r26752, r26753, r26754, r26755, r26756, r26757, r26758, r26759, r26760, r26761, r26762, r26763, r26764, r26765, r26766, r26767, r26768, r26769; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26747); mpfr_init(r26748); mpfr_init(r26749); mpfr_init(r26750); mpfr_init(r26751); mpfr_init(r26752); mpfr_init(r26753); mpfr_init(r26754); mpfr_init(r26755); mpfr_init(r26756); mpfr_init(r26757); mpfr_init(r26758); mpfr_init(r26759); mpfr_init(r26760); mpfr_init(r26761); mpfr_init(r26762); mpfr_init(r26763); mpfr_init(r26764); mpfr_init(r26765); mpfr_init(r26766); mpfr_init(r26767); mpfr_init(r26768); mpfr_init(r26769); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r26747, x, MPFR_RNDN); mpfr_set_d(r26748, y, MPFR_RNDN); mpfr_set_d(r26749, z, MPFR_RNDN); mpfr_mul(r26750, r26748, r26749, MPFR_RNDN); mpfr_set_d(r26751, t, MPFR_RNDN); mpfr_set_d(r26752, a, MPFR_RNDN); mpfr_mul(r26753, r26751, r26752, MPFR_RNDN); mpfr_sub(r26754, r26750, r26753, MPFR_RNDN); mpfr_mul(r26755, r26747, r26754, MPFR_RNDN); mpfr_set_d(r26756, b, MPFR_RNDN); mpfr_set_d(r26757, c, MPFR_RNDN); mpfr_mul(r26758, r26757, r26749, MPFR_RNDN); mpfr_set_d(r26759, i, MPFR_RNDN); mpfr_mul(r26760, r26759, r26752, MPFR_RNDN); mpfr_sub(r26761, r26758, r26760, MPFR_RNDN); mpfr_mul(r26762, r26756, r26761, MPFR_RNDN); mpfr_sub(r26763, r26755, r26762, MPFR_RNDN); mpfr_set_d(r26764, j, MPFR_RNDN); mpfr_mul(r26765, r26757, r26751, MPFR_RNDN); mpfr_mul(r26766, r26759, r26748, MPFR_RNDN); mpfr_sub(r26767, r26765, r26766, MPFR_RNDN); mpfr_mul(r26768, r26764, r26767, MPFR_RNDN); mpfr_add(r26769, r26763, r26768, MPFR_RNDN); return mpfr_get_d(r26769, MPFR_RNDN); } static mpfr_t r26770, r26771, r26772, r26773, r26774, r26775, r26776, r26777, r26778, r26779, r26780, r26781, r26782, r26783, r26784, r26785, r26786, r26787, r26788, r26789, r26790, r26791, r26792, r26793, r26794, r26795, r26796, r26797, r26798; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26770); mpfr_init(r26771); mpfr_init(r26772); mpfr_init(r26773); mpfr_init(r26774); mpfr_init(r26775); mpfr_init(r26776); mpfr_init(r26777); mpfr_init(r26778); mpfr_init(r26779); mpfr_init(r26780); mpfr_init(r26781); mpfr_init(r26782); mpfr_init(r26783); mpfr_init(r26784); mpfr_init(r26785); mpfr_init(r26786); mpfr_init(r26787); mpfr_init(r26788); mpfr_init(r26789); mpfr_init(r26790); mpfr_init(r26791); mpfr_init(r26792); mpfr_init(r26793); mpfr_init(r26794); mpfr_init(r26795); mpfr_init(r26796); mpfr_init(r26797); mpfr_init(r26798); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r26770, x, MPFR_RNDN); mpfr_set_d(r26771, y, MPFR_RNDN); mpfr_set_d(r26772, z, MPFR_RNDN); mpfr_mul(r26773, r26771, r26772, MPFR_RNDN); mpfr_set_d(r26774, t, MPFR_RNDN); mpfr_set_d(r26775, a, MPFR_RNDN); mpfr_mul(r26776, r26774, r26775, MPFR_RNDN); mpfr_sub(r26777, r26773, r26776, MPFR_RNDN); mpfr_mul(r26778, r26770, r26777, MPFR_RNDN); mpfr_set_d(r26779, b, MPFR_RNDN); mpfr_set_d(r26780, c, MPFR_RNDN); mpfr_mul(r26781, r26780, r26772, MPFR_RNDN); mpfr_set_d(r26782, i, MPFR_RNDN); mpfr_mul(r26783, r26782, r26775, MPFR_RNDN); mpfr_sub(r26784, r26781, r26783, MPFR_RNDN); mpfr_mul(r26785, r26779, r26784, MPFR_RNDN); mpfr_sub(r26786, r26778, r26785, MPFR_RNDN); mpfr_set_d(r26787, j, MPFR_RNDN); mpfr_mul(r26788, r26780, r26774, MPFR_RNDN); mpfr_mul(r26789, r26782, r26771, MPFR_RNDN); mpfr_sub(r26790, r26788, r26789, MPFR_RNDN); mpfr_mul(r26791, r26787, r26790, MPFR_RNDN); mpfr_cbrt(r26792, r26791, MPFR_RNDN); mpfr_mul(r26793, r26792, r26792, MPFR_RNDN); mpfr_cbrt(r26794, r26787, MPFR_RNDN); mpfr_cbrt(r26795, r26790, MPFR_RNDN); mpfr_mul(r26796, r26794, r26795, MPFR_RNDN); mpfr_mul(r26797, r26793, r26796, MPFR_RNDN); mpfr_add(r26798, r26786, r26797, MPFR_RNDN); return mpfr_get_d(r26798, MPFR_RNDN); } static mpfr_t r26799, r26800, r26801, r26802, r26803, r26804, r26805, r26806, r26807, r26808, r26809, r26810, r26811, r26812, r26813, r26814, r26815, r26816, r26817, r26818, r26819, r26820, r26821, r26822, r26823, r26824, r26825, r26826, r26827; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26799); mpfr_init(r26800); mpfr_init(r26801); mpfr_init(r26802); mpfr_init(r26803); mpfr_init(r26804); mpfr_init(r26805); mpfr_init(r26806); mpfr_init(r26807); mpfr_init(r26808); mpfr_init(r26809); mpfr_init(r26810); mpfr_init(r26811); mpfr_init(r26812); mpfr_init(r26813); mpfr_init(r26814); mpfr_init(r26815); mpfr_init(r26816); mpfr_init(r26817); mpfr_init(r26818); mpfr_init(r26819); mpfr_init(r26820); mpfr_init(r26821); mpfr_init(r26822); mpfr_init(r26823); mpfr_init(r26824); mpfr_init(r26825); mpfr_init(r26826); mpfr_init(r26827); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r26799, x, MPFR_RNDN); mpfr_set_d(r26800, y, MPFR_RNDN); mpfr_set_d(r26801, z, MPFR_RNDN); mpfr_mul(r26802, r26800, r26801, MPFR_RNDN); mpfr_set_d(r26803, t, MPFR_RNDN); mpfr_set_d(r26804, a, MPFR_RNDN); mpfr_mul(r26805, r26803, r26804, MPFR_RNDN); mpfr_sub(r26806, r26802, r26805, MPFR_RNDN); mpfr_mul(r26807, r26799, r26806, MPFR_RNDN); mpfr_set_d(r26808, b, MPFR_RNDN); mpfr_set_d(r26809, c, MPFR_RNDN); mpfr_mul(r26810, r26809, r26801, MPFR_RNDN); mpfr_set_d(r26811, i, MPFR_RNDN); mpfr_mul(r26812, r26811, r26804, MPFR_RNDN); mpfr_sub(r26813, r26810, r26812, MPFR_RNDN); mpfr_mul(r26814, r26808, r26813, MPFR_RNDN); mpfr_sub(r26815, r26807, r26814, MPFR_RNDN); mpfr_set_d(r26816, j, MPFR_RNDN); mpfr_mul(r26817, r26809, r26803, MPFR_RNDN); mpfr_mul(r26818, r26811, r26800, MPFR_RNDN); mpfr_sub(r26819, r26817, r26818, MPFR_RNDN); mpfr_mul(r26820, r26816, r26819, MPFR_RNDN); mpfr_cbrt(r26821, r26820, MPFR_RNDN); mpfr_mul(r26822, r26821, r26821, MPFR_RNDN); mpfr_cbrt(r26823, r26816, MPFR_RNDN); mpfr_cbrt(r26824, r26819, MPFR_RNDN); mpfr_mul(r26825, r26823, r26824, MPFR_RNDN); mpfr_mul(r26826, r26822, r26825, MPFR_RNDN); mpfr_add(r26827, r26815, r26826, MPFR_RNDN); return mpfr_get_d(r26827, MPFR_RNDN); }