#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 r24184 = phi1; float r24185 = sin(r24184); float r24186 = phi2; float r24187 = sin(r24186); float r24188 = r24185 * r24187; float r24189 = cos(r24184); float r24190 = cos(r24186); float r24191 = r24189 * r24190; float r24192 = lambda1; float r24193 = lambda2; float r24194 = r24192 - r24193; float r24195 = cos(r24194); float r24196 = r24191 * r24195; float r24197 = r24188 + r24196; float r24198 = acos(r24197); float r24199 = R; float r24200 = r24198 * r24199; return r24200; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r24201 = phi1; double r24202 = sin(r24201); double r24203 = phi2; double r24204 = sin(r24203); double r24205 = r24202 * r24204; double r24206 = cos(r24201); double r24207 = cos(r24203); double r24208 = r24206 * r24207; double r24209 = lambda1; double r24210 = lambda2; double r24211 = r24209 - r24210; double r24212 = cos(r24211); double r24213 = r24208 * r24212; double r24214 = r24205 + r24213; double r24215 = acos(r24214); double r24216 = R; double r24217 = r24215 * r24216; return r24217; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r24218 = lambda2; float r24219 = sin(r24218); float r24220 = lambda1; float r24221 = sin(r24220); float r24222 = cos(r24218); float r24223 = cos(r24220); float r24224 = r24222 * r24223; float r24225 = fma(r24219, r24221, r24224); float r24226 = phi2; float r24227 = cos(r24226); float r24228 = phi1; float r24229 = cos(r24228); float r24230 = r24227 * r24229; float r24231 = sin(r24226); float r24232 = sin(r24228); float r24233 = r24231 * r24232; float r24234 = fma(r24225, r24230, r24233); float r24235 = acos(r24234); float r24236 = expm1(r24235); float r24237 = exp(r24236); float r24238 = log(r24237); float r24239 = log1p(r24238); float r24240 = R; float r24241 = r24239 * r24240; return r24241; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r24242 = lambda2; double r24243 = sin(r24242); double r24244 = lambda1; double r24245 = sin(r24244); double r24246 = cos(r24242); double r24247 = cos(r24244); double r24248 = r24246 * r24247; double r24249 = fma(r24243, r24245, r24248); double r24250 = phi2; double r24251 = cos(r24250); double r24252 = phi1; double r24253 = cos(r24252); double r24254 = r24251 * r24253; double r24255 = sin(r24250); double r24256 = sin(r24252); double r24257 = r24255 * r24256; double r24258 = fma(r24249, r24254, r24257); double r24259 = acos(r24258); double r24260 = expm1(r24259); double r24261 = exp(r24260); double r24262 = log(r24261); double r24263 = log1p(r24262); double r24264 = R; double r24265 = r24263 * r24264; return r24265; } 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 r24266, r24267, r24268, r24269, r24270, r24271, r24272, r24273, r24274, r24275, r24276, r24277, r24278, r24279, r24280, r24281, r24282; void setup_mpfr_f_im() { mpfr_set_default_prec(1872); mpfr_init(r24266); mpfr_init(r24267); mpfr_init(r24268); mpfr_init(r24269); mpfr_init(r24270); mpfr_init(r24271); mpfr_init(r24272); mpfr_init(r24273); mpfr_init(r24274); mpfr_init(r24275); mpfr_init(r24276); mpfr_init(r24277); mpfr_init(r24278); mpfr_init(r24279); mpfr_init(r24280); mpfr_init(r24281); mpfr_init(r24282); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24266, phi1, MPFR_RNDN); mpfr_sin(r24267, r24266, MPFR_RNDN); mpfr_set_d(r24268, phi2, MPFR_RNDN); mpfr_sin(r24269, r24268, MPFR_RNDN); mpfr_mul(r24270, r24267, r24269, MPFR_RNDN); mpfr_cos(r24271, r24266, MPFR_RNDN); mpfr_cos(r24272, r24268, MPFR_RNDN); mpfr_mul(r24273, r24271, r24272, MPFR_RNDN); mpfr_set_d(r24274, lambda1, MPFR_RNDN); mpfr_set_d(r24275, lambda2, MPFR_RNDN); mpfr_sub(r24276, r24274, r24275, MPFR_RNDN); mpfr_cos(r24277, r24276, MPFR_RNDN); mpfr_mul(r24278, r24273, r24277, MPFR_RNDN); mpfr_add(r24279, r24270, r24278, MPFR_RNDN); mpfr_acos(r24280, r24279, MPFR_RNDN); mpfr_set_d(r24281, R, MPFR_RNDN); mpfr_mul(r24282, r24280, r24281, MPFR_RNDN); return mpfr_get_d(r24282, MPFR_RNDN); } static mpfr_t r24283, r24284, r24285, r24286, r24287, r24288, r24289, r24290, r24291, r24292, r24293, r24294, r24295, r24296, r24297, r24298, r24299, r24300, r24301, r24302, r24303, r24304, r24305, r24306; void setup_mpfr_f_fm() { mpfr_set_default_prec(1872); mpfr_init(r24283); mpfr_init(r24284); mpfr_init(r24285); mpfr_init(r24286); mpfr_init(r24287); mpfr_init(r24288); mpfr_init(r24289); mpfr_init(r24290); mpfr_init(r24291); mpfr_init(r24292); mpfr_init(r24293); mpfr_init(r24294); mpfr_init(r24295); mpfr_init(r24296); mpfr_init(r24297); mpfr_init(r24298); mpfr_init(r24299); mpfr_init(r24300); mpfr_init(r24301); mpfr_init(r24302); mpfr_init(r24303); mpfr_init(r24304); mpfr_init(r24305); mpfr_init(r24306); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24283, lambda2, MPFR_RNDN); mpfr_sin(r24284, r24283, MPFR_RNDN); mpfr_set_d(r24285, lambda1, MPFR_RNDN); mpfr_sin(r24286, r24285, MPFR_RNDN); mpfr_cos(r24287, r24283, MPFR_RNDN); mpfr_cos(r24288, r24285, MPFR_RNDN); mpfr_mul(r24289, r24287, r24288, MPFR_RNDN); mpfr_fma(r24290, r24284, r24286, r24289, MPFR_RNDN); mpfr_set_d(r24291, phi2, MPFR_RNDN); mpfr_cos(r24292, r24291, MPFR_RNDN); mpfr_set_d(r24293, phi1, MPFR_RNDN); mpfr_cos(r24294, r24293, MPFR_RNDN); mpfr_mul(r24295, r24292, r24294, MPFR_RNDN); mpfr_sin(r24296, r24291, MPFR_RNDN); mpfr_sin(r24297, r24293, MPFR_RNDN); mpfr_mul(r24298, r24296, r24297, MPFR_RNDN); mpfr_fma(r24299, r24290, r24295, r24298, MPFR_RNDN); mpfr_acos(r24300, r24299, MPFR_RNDN); mpfr_expm1(r24301, r24300, MPFR_RNDN); mpfr_exp(r24302, r24301, MPFR_RNDN); mpfr_log(r24303, r24302, MPFR_RNDN); mpfr_log1p(r24304, r24303, MPFR_RNDN); mpfr_set_d(r24305, R, MPFR_RNDN); mpfr_mul(r24306, r24304, r24305, MPFR_RNDN); return mpfr_get_d(r24306, MPFR_RNDN); } static mpfr_t r24307, r24308, r24309, r24310, r24311, r24312, r24313, r24314, r24315, r24316, r24317, r24318, r24319, r24320, r24321, r24322, r24323, r24324, r24325, r24326, r24327, r24328, r24329, r24330; void setup_mpfr_f_dm() { mpfr_set_default_prec(1872); mpfr_init(r24307); mpfr_init(r24308); mpfr_init(r24309); mpfr_init(r24310); mpfr_init(r24311); mpfr_init(r24312); mpfr_init(r24313); mpfr_init(r24314); mpfr_init(r24315); mpfr_init(r24316); mpfr_init(r24317); mpfr_init(r24318); mpfr_init(r24319); mpfr_init(r24320); mpfr_init(r24321); mpfr_init(r24322); mpfr_init(r24323); mpfr_init(r24324); mpfr_init(r24325); mpfr_init(r24326); mpfr_init(r24327); mpfr_init(r24328); mpfr_init(r24329); mpfr_init(r24330); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24307, lambda2, MPFR_RNDN); mpfr_sin(r24308, r24307, MPFR_RNDN); mpfr_set_d(r24309, lambda1, MPFR_RNDN); mpfr_sin(r24310, r24309, MPFR_RNDN); mpfr_cos(r24311, r24307, MPFR_RNDN); mpfr_cos(r24312, r24309, MPFR_RNDN); mpfr_mul(r24313, r24311, r24312, MPFR_RNDN); mpfr_fma(r24314, r24308, r24310, r24313, MPFR_RNDN); mpfr_set_d(r24315, phi2, MPFR_RNDN); mpfr_cos(r24316, r24315, MPFR_RNDN); mpfr_set_d(r24317, phi1, MPFR_RNDN); mpfr_cos(r24318, r24317, MPFR_RNDN); mpfr_mul(r24319, r24316, r24318, MPFR_RNDN); mpfr_sin(r24320, r24315, MPFR_RNDN); mpfr_sin(r24321, r24317, MPFR_RNDN); mpfr_mul(r24322, r24320, r24321, MPFR_RNDN); mpfr_fma(r24323, r24314, r24319, r24322, MPFR_RNDN); mpfr_acos(r24324, r24323, MPFR_RNDN); mpfr_expm1(r24325, r24324, MPFR_RNDN); mpfr_exp(r24326, r24325, MPFR_RNDN); mpfr_log(r24327, r24326, MPFR_RNDN); mpfr_log1p(r24328, r24327, MPFR_RNDN); mpfr_set_d(r24329, R, MPFR_RNDN); mpfr_mul(r24330, r24328, r24329, MPFR_RNDN); return mpfr_get_d(r24330, MPFR_RNDN); }