#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 r20669 = phi1; float r20670 = sin(r20669); float r20671 = phi2; float r20672 = sin(r20671); float r20673 = r20670 * r20672; float r20674 = cos(r20669); float r20675 = cos(r20671); float r20676 = r20674 * r20675; float r20677 = lambda1; float r20678 = lambda2; float r20679 = r20677 - r20678; float r20680 = cos(r20679); float r20681 = r20676 * r20680; float r20682 = r20673 + r20681; float r20683 = acos(r20682); float r20684 = R; float r20685 = r20683 * r20684; return r20685; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r20686 = phi1; double r20687 = sin(r20686); double r20688 = phi2; double r20689 = sin(r20688); double r20690 = r20687 * r20689; double r20691 = cos(r20686); double r20692 = cos(r20688); double r20693 = r20691 * r20692; double r20694 = lambda1; double r20695 = lambda2; double r20696 = r20694 - r20695; double r20697 = cos(r20696); double r20698 = r20693 * r20697; double r20699 = r20690 + r20698; double r20700 = acos(r20699); double r20701 = R; double r20702 = r20700 * r20701; return r20702; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r20703 = phi2; float r20704 = cos(r20703); float r20705 = phi1; float r20706 = cos(r20705); float r20707 = r20704 * r20706; float r20708 = lambda1; float r20709 = cos(r20708); float r20710 = lambda2; float r20711 = cos(r20710); float r20712 = sin(r20710); float r20713 = sin(r20708); float r20714 = r20712 * r20713; float r20715 = fma(r20709, r20711, r20714); float r20716 = sin(r20703); float r20717 = sin(r20705); float r20718 = r20716 * r20717; float r20719 = fma(r20707, r20715, r20718); float r20720 = log1p(r20719); float r20721 = expm1(r20720); float r20722 = acos(r20721); float r20723 = exp(r20722); float r20724 = log(r20723); float r20725 = R; float r20726 = r20724 * r20725; return r20726; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r20727 = phi2; double r20728 = cos(r20727); double r20729 = phi1; double r20730 = cos(r20729); double r20731 = r20728 * r20730; double r20732 = lambda1; double r20733 = cos(r20732); double r20734 = lambda2; double r20735 = cos(r20734); double r20736 = sin(r20734); double r20737 = sin(r20732); double r20738 = r20736 * r20737; double r20739 = fma(r20733, r20735, r20738); double r20740 = sin(r20727); double r20741 = sin(r20729); double r20742 = r20740 * r20741; double r20743 = fma(r20731, r20739, r20742); double r20744 = log1p(r20743); double r20745 = expm1(r20744); double r20746 = acos(r20745); double r20747 = exp(r20746); double r20748 = log(r20747); double r20749 = R; double r20750 = r20748 * r20749; return r20750; } 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 r20751, r20752, r20753, r20754, r20755, r20756, r20757, r20758, r20759, r20760, r20761, r20762, r20763, r20764, r20765, r20766, r20767; void setup_mpfr_f_im() { mpfr_set_default_prec(2192); mpfr_init(r20751); mpfr_init(r20752); mpfr_init(r20753); mpfr_init(r20754); mpfr_init(r20755); mpfr_init(r20756); mpfr_init(r20757); mpfr_init(r20758); mpfr_init(r20759); mpfr_init(r20760); mpfr_init(r20761); mpfr_init(r20762); mpfr_init(r20763); mpfr_init(r20764); mpfr_init(r20765); mpfr_init(r20766); mpfr_init(r20767); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r20751, phi1, MPFR_RNDN); mpfr_sin(r20752, r20751, MPFR_RNDN); mpfr_set_d(r20753, phi2, MPFR_RNDN); mpfr_sin(r20754, r20753, MPFR_RNDN); mpfr_mul(r20755, r20752, r20754, MPFR_RNDN); mpfr_cos(r20756, r20751, MPFR_RNDN); mpfr_cos(r20757, r20753, MPFR_RNDN); mpfr_mul(r20758, r20756, r20757, MPFR_RNDN); mpfr_set_d(r20759, lambda1, MPFR_RNDN); mpfr_set_d(r20760, lambda2, MPFR_RNDN); mpfr_sub(r20761, r20759, r20760, MPFR_RNDN); mpfr_cos(r20762, r20761, MPFR_RNDN); mpfr_mul(r20763, r20758, r20762, MPFR_RNDN); mpfr_add(r20764, r20755, r20763, MPFR_RNDN); mpfr_acos(r20765, r20764, MPFR_RNDN); mpfr_set_d(r20766, R, MPFR_RNDN); mpfr_mul(r20767, r20765, r20766, MPFR_RNDN); return mpfr_get_d(r20767, MPFR_RNDN); } static mpfr_t r20768, r20769, r20770, r20771, r20772, r20773, r20774, r20775, r20776, r20777, r20778, r20779, r20780, r20781, r20782, r20783, r20784, r20785, r20786, r20787, r20788, r20789, r20790, r20791; void setup_mpfr_f_fm() { mpfr_set_default_prec(2192); mpfr_init(r20768); mpfr_init(r20769); mpfr_init(r20770); mpfr_init(r20771); mpfr_init(r20772); mpfr_init(r20773); mpfr_init(r20774); 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_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r20768, phi2, MPFR_RNDN); mpfr_cos(r20769, r20768, MPFR_RNDN); mpfr_set_d(r20770, phi1, MPFR_RNDN); mpfr_cos(r20771, r20770, MPFR_RNDN); mpfr_mul(r20772, r20769, r20771, MPFR_RNDN); mpfr_set_d(r20773, lambda1, MPFR_RNDN); mpfr_cos(r20774, r20773, MPFR_RNDN); mpfr_set_d(r20775, lambda2, MPFR_RNDN); mpfr_cos(r20776, r20775, MPFR_RNDN); mpfr_sin(r20777, r20775, MPFR_RNDN); mpfr_sin(r20778, r20773, MPFR_RNDN); mpfr_mul(r20779, r20777, r20778, MPFR_RNDN); mpfr_fma(r20780, r20774, r20776, r20779, MPFR_RNDN); mpfr_sin(r20781, r20768, MPFR_RNDN); mpfr_sin(r20782, r20770, MPFR_RNDN); mpfr_mul(r20783, r20781, r20782, MPFR_RNDN); mpfr_fma(r20784, r20772, r20780, r20783, MPFR_RNDN); mpfr_log1p(r20785, r20784, MPFR_RNDN); mpfr_expm1(r20786, r20785, MPFR_RNDN); mpfr_acos(r20787, r20786, MPFR_RNDN); mpfr_exp(r20788, r20787, MPFR_RNDN); mpfr_log(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_dm() { mpfr_set_default_prec(2192); 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_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r20792, phi2, MPFR_RNDN); mpfr_cos(r20793, r20792, MPFR_RNDN); mpfr_set_d(r20794, phi1, 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_log1p(r20809, r20808, MPFR_RNDN); mpfr_expm1(r20810, r20809, MPFR_RNDN); mpfr_acos(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); }