#include #include #include #include #include char *name = "Spherical law of cosines"; double f_if(float R, float lambda1, float lambda2, float phi1, float phi2) { float r20693 = phi1; float r20694 = sin(r20693); float r20695 = phi2; float r20696 = sin(r20695); float r20697 = r20694 * r20696; float r20698 = cos(r20693); float r20699 = cos(r20695); float r20700 = r20698 * r20699; float r20701 = lambda1; float r20702 = lambda2; float r20703 = r20701 - r20702; float r20704 = cos(r20703); float r20705 = r20700 * r20704; float r20706 = r20697 + r20705; float r20707 = acos(r20706); float r20708 = R; float r20709 = r20707 * r20708; return r20709; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r20710 = phi1; double r20711 = sin(r20710); double r20712 = phi2; double r20713 = sin(r20712); double r20714 = r20711 * r20713; double r20715 = cos(r20710); double r20716 = cos(r20712); double r20717 = r20715 * r20716; double r20718 = lambda1; double r20719 = lambda2; double r20720 = r20718 - r20719; double r20721 = cos(r20720); double r20722 = r20717 * r20721; double r20723 = r20714 + r20722; double r20724 = acos(r20723); double r20725 = R; double r20726 = r20724 * r20725; return r20726; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r20727 = phi1; float r20728 = cos(r20727); float r20729 = phi2; float r20730 = cos(r20729); float r20731 = r20728 * r20730; float r20732 = lambda1; float r20733 = cos(r20732); float r20734 = lambda2; float r20735 = cos(r20734); float r20736 = sin(r20734); float r20737 = sin(r20732); float r20738 = r20736 * r20737; float r20739 = fma(r20733, r20735, r20738); float r20740 = sin(r20727); float r20741 = sin(r20729); float r20742 = r20740 * r20741; float r20743 = fma(r20731, r20739, r20742); float r20744 = acos(r20743); float r20745 = log(r20744); float r20746 = exp(r20745); float r20747 = exp(r20746); float r20748 = log(r20747); float r20749 = R; float r20750 = r20748 * r20749; return r20750; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r20751 = phi1; double r20752 = cos(r20751); double r20753 = phi2; double r20754 = cos(r20753); double r20755 = r20752 * r20754; double r20756 = lambda1; double r20757 = cos(r20756); double r20758 = lambda2; double r20759 = cos(r20758); double r20760 = sin(r20758); double r20761 = sin(r20756); double r20762 = r20760 * r20761; double r20763 = fma(r20757, r20759, r20762); double r20764 = sin(r20751); double r20765 = sin(r20753); double r20766 = r20764 * r20765; double r20767 = fma(r20755, r20763, r20766); double r20768 = acos(r20767); double r20769 = log(r20768); double r20770 = exp(r20769); double r20771 = exp(r20770); double r20772 = log(r20771); double r20773 = R; double r20774 = r20772 * r20773; return r20774; } 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 r20775, r20776, r20777, r20778, r20779, r20780, r20781, r20782, r20783, r20784, r20785, r20786, r20787, r20788, r20789, r20790, r20791; void setup_mpfr_f_im() { mpfr_set_default_prec(1936); mpfr_init(r20775); mpfr_init(r20776); mpfr_init(r20777); mpfr_init(r20778); mpfr_init(r20779); mpfr_init(r20780); mpfr_init(r20781); mpfr_init(r20782); mpfr_init(r20783); mpfr_init(r20784); mpfr_init(r20785); mpfr_init(r20786); mpfr_init(r20787); mpfr_init(r20788); mpfr_init(r20789); mpfr_init(r20790); mpfr_init(r20791); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r20775, phi1, MPFR_RNDN); mpfr_sin(r20776, r20775, MPFR_RNDN); mpfr_set_d(r20777, phi2, MPFR_RNDN); mpfr_sin(r20778, r20777, MPFR_RNDN); mpfr_mul(r20779, r20776, r20778, MPFR_RNDN); mpfr_cos(r20780, r20775, MPFR_RNDN); mpfr_cos(r20781, r20777, MPFR_RNDN); mpfr_mul(r20782, r20780, r20781, MPFR_RNDN); mpfr_set_d(r20783, lambda1, MPFR_RNDN); mpfr_set_d(r20784, lambda2, MPFR_RNDN); mpfr_sub(r20785, r20783, r20784, MPFR_RNDN); mpfr_cos(r20786, r20785, MPFR_RNDN); mpfr_mul(r20787, r20782, r20786, MPFR_RNDN); mpfr_add(r20788, r20779, r20787, MPFR_RNDN); mpfr_acos(r20789, r20788, MPFR_RNDN); mpfr_set_d(r20790, R, MPFR_RNDN); mpfr_mul(r20791, r20789, r20790, MPFR_RNDN); return mpfr_get_d(r20791, MPFR_RNDN); } static mpfr_t r20792, r20793, r20794, r20795, r20796, r20797, r20798, r20799, r20800, r20801, r20802, r20803, r20804, r20805, r20806, r20807, r20808, r20809, r20810, r20811, r20812, r20813, r20814, r20815; void setup_mpfr_f_fm() { mpfr_set_default_prec(1936); mpfr_init(r20792); mpfr_init(r20793); mpfr_init(r20794); mpfr_init(r20795); mpfr_init(r20796); mpfr_init(r20797); mpfr_init(r20798); mpfr_init(r20799); mpfr_init(r20800); mpfr_init(r20801); mpfr_init(r20802); mpfr_init(r20803); mpfr_init(r20804); mpfr_init(r20805); mpfr_init(r20806); mpfr_init(r20807); mpfr_init(r20808); mpfr_init(r20809); mpfr_init(r20810); mpfr_init(r20811); mpfr_init(r20812); mpfr_init(r20813); mpfr_init(r20814); mpfr_init(r20815); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r20792, phi1, MPFR_RNDN); mpfr_cos(r20793, r20792, MPFR_RNDN); mpfr_set_d(r20794, phi2, MPFR_RNDN); mpfr_cos(r20795, r20794, MPFR_RNDN); mpfr_mul(r20796, r20793, r20795, MPFR_RNDN); mpfr_set_d(r20797, lambda1, MPFR_RNDN); mpfr_cos(r20798, r20797, MPFR_RNDN); mpfr_set_d(r20799, lambda2, MPFR_RNDN); mpfr_cos(r20800, r20799, MPFR_RNDN); mpfr_sin(r20801, r20799, MPFR_RNDN); mpfr_sin(r20802, r20797, MPFR_RNDN); mpfr_mul(r20803, r20801, r20802, MPFR_RNDN); mpfr_fma(r20804, r20798, r20800, r20803, MPFR_RNDN); mpfr_sin(r20805, r20792, MPFR_RNDN); mpfr_sin(r20806, r20794, MPFR_RNDN); mpfr_mul(r20807, r20805, r20806, MPFR_RNDN); mpfr_fma(r20808, r20796, r20804, r20807, MPFR_RNDN); mpfr_acos(r20809, r20808, MPFR_RNDN); mpfr_log(r20810, r20809, MPFR_RNDN); mpfr_exp(r20811, r20810, MPFR_RNDN); mpfr_exp(r20812, r20811, MPFR_RNDN); mpfr_log(r20813, r20812, MPFR_RNDN); mpfr_set_d(r20814, R, MPFR_RNDN); mpfr_mul(r20815, r20813, r20814, MPFR_RNDN); return mpfr_get_d(r20815, MPFR_RNDN); } static mpfr_t r20816, r20817, r20818, r20819, r20820, r20821, r20822, r20823, r20824, r20825, r20826, r20827, r20828, r20829, r20830, r20831, r20832, r20833, r20834, r20835, r20836, r20837, r20838, r20839; void setup_mpfr_f_dm() { mpfr_set_default_prec(1936); mpfr_init(r20816); mpfr_init(r20817); mpfr_init(r20818); mpfr_init(r20819); mpfr_init(r20820); mpfr_init(r20821); mpfr_init(r20822); mpfr_init(r20823); mpfr_init(r20824); mpfr_init(r20825); mpfr_init(r20826); mpfr_init(r20827); mpfr_init(r20828); mpfr_init(r20829); mpfr_init(r20830); mpfr_init(r20831); mpfr_init(r20832); mpfr_init(r20833); mpfr_init(r20834); mpfr_init(r20835); mpfr_init(r20836); mpfr_init(r20837); mpfr_init(r20838); mpfr_init(r20839); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r20816, phi1, MPFR_RNDN); mpfr_cos(r20817, r20816, MPFR_RNDN); mpfr_set_d(r20818, phi2, MPFR_RNDN); mpfr_cos(r20819, r20818, MPFR_RNDN); mpfr_mul(r20820, r20817, r20819, MPFR_RNDN); mpfr_set_d(r20821, lambda1, MPFR_RNDN); mpfr_cos(r20822, r20821, MPFR_RNDN); mpfr_set_d(r20823, lambda2, MPFR_RNDN); mpfr_cos(r20824, r20823, MPFR_RNDN); mpfr_sin(r20825, r20823, MPFR_RNDN); mpfr_sin(r20826, r20821, MPFR_RNDN); mpfr_mul(r20827, r20825, r20826, MPFR_RNDN); mpfr_fma(r20828, r20822, r20824, r20827, MPFR_RNDN); mpfr_sin(r20829, r20816, MPFR_RNDN); mpfr_sin(r20830, r20818, MPFR_RNDN); mpfr_mul(r20831, r20829, r20830, MPFR_RNDN); mpfr_fma(r20832, r20820, r20828, r20831, MPFR_RNDN); mpfr_acos(r20833, r20832, MPFR_RNDN); mpfr_log(r20834, r20833, MPFR_RNDN); mpfr_exp(r20835, r20834, MPFR_RNDN); mpfr_exp(r20836, r20835, MPFR_RNDN); mpfr_log(r20837, r20836, MPFR_RNDN); mpfr_set_d(r20838, R, MPFR_RNDN); mpfr_mul(r20839, r20837, r20838, MPFR_RNDN); return mpfr_get_d(r20839, MPFR_RNDN); }