#include #include #include #include #include char *name = "FRP.Yampa.Vector3:vector3Rho from Yampa-0.10.2"; double f_if(float x, float y, float z) { float r51824 = x; float r51825 = r51824 * r51824; float r51826 = y; float r51827 = r51826 * r51826; float r51828 = r51825 + r51827; float r51829 = z; float r51830 = r51829 * r51829; float r51831 = r51828 + r51830; float r51832 = sqrt(r51831); return r51832; } double f_id(double x, double y, double z) { double r51833 = x; double r51834 = r51833 * r51833; double r51835 = y; double r51836 = r51835 * r51835; double r51837 = r51834 + r51836; double r51838 = z; double r51839 = r51838 * r51838; double r51840 = r51837 + r51839; double r51841 = sqrt(r51840); return r51841; } double f_of(float x, float y, float z) { float r51842 = x; float r51843 = -2.013881564242057e+151; bool r51844 = r51842 <= r51843; float r51845 = -r51842; float r51846 = 1.6698180397332825e+161; bool r51847 = r51842 <= r51846; float r51848 = r51842 * r51842; float r51849 = y; float r51850 = r51849 * r51849; float r51851 = r51848 + r51850; float r51852 = z; float r51853 = r51852 * r51852; float r51854 = r51851 + r51853; float r51855 = sqrt(r51854); float r51856 = r51847 ? r51855 : r51842; float r51857 = r51844 ? r51845 : r51856; return r51857; } double f_od(double x, double y, double z) { double r51858 = x; double r51859 = -2.013881564242057e+151; bool r51860 = r51858 <= r51859; double r51861 = -r51858; double r51862 = 1.6698180397332825e+161; bool r51863 = r51858 <= r51862; double r51864 = r51858 * r51858; double r51865 = y; double r51866 = r51865 * r51865; double r51867 = r51864 + r51866; double r51868 = z; double r51869 = r51868 * r51868; double r51870 = r51867 + r51869; double r51871 = sqrt(r51870); double r51872 = r51863 ? r51871 : r51858; double r51873 = r51860 ? r51861 : r51872; return r51873; } 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 r51874, r51875, r51876, r51877, r51878, r51879, r51880, r51881, r51882; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r51874); mpfr_init(r51875); mpfr_init(r51876); mpfr_init(r51877); mpfr_init(r51878); mpfr_init(r51879); mpfr_init(r51880); mpfr_init(r51881); mpfr_init(r51882); } double f_im(double x, double y, double z) { mpfr_set_d(r51874, x, MPFR_RNDN); mpfr_mul(r51875, r51874, r51874, MPFR_RNDN); mpfr_set_d(r51876, y, MPFR_RNDN); mpfr_mul(r51877, r51876, r51876, MPFR_RNDN); mpfr_add(r51878, r51875, r51877, MPFR_RNDN); mpfr_set_d(r51879, z, MPFR_RNDN); mpfr_mul(r51880, r51879, r51879, MPFR_RNDN); mpfr_add(r51881, r51878, r51880, MPFR_RNDN); mpfr_sqrt(r51882, r51881, MPFR_RNDN); return mpfr_get_d(r51882, MPFR_RNDN); } static mpfr_t r51883, r51884, r51885, r51886, r51887, r51888, r51889, r51890, r51891, r51892, r51893, r51894, r51895, r51896, r51897, r51898; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r51883); mpfr_init_set_str(r51884, "-2.013881564242057e+151", 10, MPFR_RNDN); mpfr_init(r51885); mpfr_init(r51886); mpfr_init_set_str(r51887, "1.6698180397332825e+161", 10, MPFR_RNDN); mpfr_init(r51888); mpfr_init(r51889); mpfr_init(r51890); mpfr_init(r51891); mpfr_init(r51892); mpfr_init(r51893); mpfr_init(r51894); mpfr_init(r51895); mpfr_init(r51896); mpfr_init(r51897); mpfr_init(r51898); } double f_fm(double x, double y, double z) { mpfr_set_d(r51883, x, MPFR_RNDN); ; mpfr_set_si(r51885, mpfr_cmp(r51883, r51884) <= 0, MPFR_RNDN); mpfr_neg(r51886, r51883, MPFR_RNDN); ; mpfr_set_si(r51888, mpfr_cmp(r51883, r51887) <= 0, MPFR_RNDN); mpfr_mul(r51889, r51883, r51883, MPFR_RNDN); mpfr_set_d(r51890, y, MPFR_RNDN); mpfr_mul(r51891, r51890, r51890, MPFR_RNDN); mpfr_add(r51892, r51889, r51891, MPFR_RNDN); mpfr_set_d(r51893, z, MPFR_RNDN); mpfr_mul(r51894, r51893, r51893, MPFR_RNDN); mpfr_add(r51895, r51892, r51894, MPFR_RNDN); mpfr_sqrt(r51896, r51895, MPFR_RNDN); if (mpfr_get_si(r51888, MPFR_RNDN)) { mpfr_set(r51897, r51896, MPFR_RNDN); } else { mpfr_set(r51897, r51883, MPFR_RNDN); }; if (mpfr_get_si(r51885, MPFR_RNDN)) { mpfr_set(r51898, r51886, MPFR_RNDN); } else { mpfr_set(r51898, r51897, MPFR_RNDN); }; return mpfr_get_d(r51898, MPFR_RNDN); } static mpfr_t r51899, r51900, r51901, r51902, r51903, r51904, r51905, r51906, r51907, r51908, r51909, r51910, r51911, r51912, r51913, r51914; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r51899); mpfr_init_set_str(r51900, "-2.013881564242057e+151", 10, MPFR_RNDN); mpfr_init(r51901); mpfr_init(r51902); mpfr_init_set_str(r51903, "1.6698180397332825e+161", 10, MPFR_RNDN); mpfr_init(r51904); mpfr_init(r51905); mpfr_init(r51906); mpfr_init(r51907); mpfr_init(r51908); mpfr_init(r51909); mpfr_init(r51910); mpfr_init(r51911); mpfr_init(r51912); mpfr_init(r51913); mpfr_init(r51914); } double f_dm(double x, double y, double z) { mpfr_set_d(r51899, x, MPFR_RNDN); ; mpfr_set_si(r51901, mpfr_cmp(r51899, r51900) <= 0, MPFR_RNDN); mpfr_neg(r51902, r51899, MPFR_RNDN); ; mpfr_set_si(r51904, mpfr_cmp(r51899, r51903) <= 0, MPFR_RNDN); mpfr_mul(r51905, r51899, r51899, MPFR_RNDN); mpfr_set_d(r51906, y, MPFR_RNDN); mpfr_mul(r51907, r51906, r51906, MPFR_RNDN); mpfr_add(r51908, r51905, r51907, MPFR_RNDN); mpfr_set_d(r51909, z, MPFR_RNDN); mpfr_mul(r51910, r51909, r51909, MPFR_RNDN); mpfr_add(r51911, r51908, r51910, MPFR_RNDN); mpfr_sqrt(r51912, r51911, MPFR_RNDN); if (mpfr_get_si(r51904, MPFR_RNDN)) { mpfr_set(r51913, r51912, MPFR_RNDN); } else { mpfr_set(r51913, r51899, MPFR_RNDN); }; if (mpfr_get_si(r51901, MPFR_RNDN)) { mpfr_set(r51914, r51902, MPFR_RNDN); } else { mpfr_set(r51914, r51913, MPFR_RNDN); }; return mpfr_get_d(r51914, MPFR_RNDN); }