#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 r24123 = phi1; float r24124 = sin(r24123); float r24125 = phi2; float r24126 = sin(r24125); float r24127 = r24124 * r24126; float r24128 = cos(r24123); float r24129 = cos(r24125); float r24130 = r24128 * r24129; float r24131 = lambda1; float r24132 = lambda2; float r24133 = r24131 - r24132; float r24134 = cos(r24133); float r24135 = r24130 * r24134; float r24136 = r24127 + r24135; float r24137 = acos(r24136); float r24138 = R; float r24139 = r24137 * r24138; return r24139; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r24140 = phi1; double r24141 = sin(r24140); double r24142 = phi2; double r24143 = sin(r24142); double r24144 = r24141 * r24143; double r24145 = cos(r24140); double r24146 = cos(r24142); double r24147 = r24145 * r24146; double r24148 = lambda1; double r24149 = lambda2; double r24150 = r24148 - r24149; double r24151 = cos(r24150); double r24152 = r24147 * r24151; double r24153 = r24144 + r24152; double r24154 = acos(r24153); double r24155 = R; double r24156 = r24154 * r24155; return r24156; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r24157 = phi1; float r24158 = sin(r24157); float r24159 = phi2; float r24160 = sin(r24159); float r24161 = r24158 * r24160; float r24162 = cos(r24159); float r24163 = cos(r24157); float r24164 = r24162 * r24163; float r24165 = exp(r24164); float r24166 = lambda2; float r24167 = cos(r24166); float r24168 = lambda1; float r24169 = cos(r24168); float r24170 = r24167 * r24169; float r24171 = sin(r24168); float r24172 = sin(r24166); float r24173 = r24171 * r24172; float r24174 = r24170 + r24173; float r24175 = pow(r24165, r24174); float r24176 = log(r24175); float r24177 = r24161 + r24176; float r24178 = acos(r24177); float r24179 = R; float r24180 = r24178 * r24179; return r24180; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r24181 = phi1; double r24182 = sin(r24181); double r24183 = phi2; double r24184 = sin(r24183); double r24185 = r24182 * r24184; double r24186 = cos(r24183); double r24187 = cos(r24181); double r24188 = r24186 * r24187; double r24189 = exp(r24188); double r24190 = lambda2; double r24191 = cos(r24190); double r24192 = lambda1; double r24193 = cos(r24192); double r24194 = r24191 * r24193; double r24195 = sin(r24192); double r24196 = sin(r24190); double r24197 = r24195 * r24196; double r24198 = r24194 + r24197; double r24199 = pow(r24189, r24198); double r24200 = log(r24199); double r24201 = r24185 + r24200; double r24202 = acos(r24201); double r24203 = R; double r24204 = r24202 * r24203; return r24204; } 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 r24205, r24206, r24207, r24208, r24209, r24210, r24211, r24212, r24213, r24214, r24215, r24216, r24217, r24218, r24219, r24220, r24221; void setup_mpfr_f_im() { mpfr_set_default_prec(2128); mpfr_init(r24205); mpfr_init(r24206); mpfr_init(r24207); mpfr_init(r24208); mpfr_init(r24209); mpfr_init(r24210); mpfr_init(r24211); mpfr_init(r24212); mpfr_init(r24213); mpfr_init(r24214); mpfr_init(r24215); mpfr_init(r24216); mpfr_init(r24217); mpfr_init(r24218); mpfr_init(r24219); mpfr_init(r24220); mpfr_init(r24221); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24205, phi1, MPFR_RNDN); mpfr_sin(r24206, r24205, MPFR_RNDN); mpfr_set_d(r24207, phi2, MPFR_RNDN); mpfr_sin(r24208, r24207, MPFR_RNDN); mpfr_mul(r24209, r24206, r24208, MPFR_RNDN); mpfr_cos(r24210, r24205, MPFR_RNDN); mpfr_cos(r24211, r24207, MPFR_RNDN); mpfr_mul(r24212, r24210, r24211, MPFR_RNDN); mpfr_set_d(r24213, lambda1, MPFR_RNDN); mpfr_set_d(r24214, lambda2, MPFR_RNDN); mpfr_sub(r24215, r24213, r24214, MPFR_RNDN); mpfr_cos(r24216, r24215, MPFR_RNDN); mpfr_mul(r24217, r24212, r24216, MPFR_RNDN); mpfr_add(r24218, r24209, r24217, MPFR_RNDN); mpfr_acos(r24219, r24218, MPFR_RNDN); mpfr_set_d(r24220, R, MPFR_RNDN); mpfr_mul(r24221, r24219, r24220, MPFR_RNDN); return mpfr_get_d(r24221, MPFR_RNDN); } static mpfr_t r24222, r24223, r24224, r24225, r24226, r24227, r24228, r24229, r24230, r24231, r24232, r24233, r24234, r24235, r24236, r24237, r24238, r24239, r24240, r24241, r24242, r24243, r24244, r24245; void setup_mpfr_f_fm() { mpfr_set_default_prec(2128); mpfr_init(r24222); mpfr_init(r24223); mpfr_init(r24224); mpfr_init(r24225); mpfr_init(r24226); mpfr_init(r24227); mpfr_init(r24228); mpfr_init(r24229); mpfr_init(r24230); mpfr_init(r24231); mpfr_init(r24232); mpfr_init(r24233); mpfr_init(r24234); mpfr_init(r24235); mpfr_init(r24236); mpfr_init(r24237); mpfr_init(r24238); mpfr_init(r24239); mpfr_init(r24240); mpfr_init(r24241); mpfr_init(r24242); mpfr_init(r24243); mpfr_init(r24244); mpfr_init(r24245); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24222, phi1, MPFR_RNDN); mpfr_sin(r24223, r24222, MPFR_RNDN); mpfr_set_d(r24224, phi2, MPFR_RNDN); mpfr_sin(r24225, r24224, MPFR_RNDN); mpfr_mul(r24226, r24223, r24225, MPFR_RNDN); mpfr_cos(r24227, r24224, MPFR_RNDN); mpfr_cos(r24228, r24222, MPFR_RNDN); mpfr_mul(r24229, r24227, r24228, MPFR_RNDN); mpfr_exp(r24230, r24229, MPFR_RNDN); mpfr_set_d(r24231, lambda2, MPFR_RNDN); mpfr_cos(r24232, r24231, MPFR_RNDN); mpfr_set_d(r24233, lambda1, MPFR_RNDN); mpfr_cos(r24234, r24233, MPFR_RNDN); mpfr_mul(r24235, r24232, r24234, MPFR_RNDN); mpfr_sin(r24236, r24233, MPFR_RNDN); mpfr_sin(r24237, r24231, MPFR_RNDN); mpfr_mul(r24238, r24236, r24237, MPFR_RNDN); mpfr_add(r24239, r24235, r24238, MPFR_RNDN); mpfr_pow(r24240, r24230, r24239, MPFR_RNDN); mpfr_log(r24241, r24240, MPFR_RNDN); mpfr_add(r24242, r24226, r24241, MPFR_RNDN); mpfr_acos(r24243, r24242, MPFR_RNDN); mpfr_set_d(r24244, R, MPFR_RNDN); mpfr_mul(r24245, r24243, r24244, MPFR_RNDN); return mpfr_get_d(r24245, MPFR_RNDN); } static mpfr_t r24246, r24247, r24248, r24249, r24250, r24251, r24252, r24253, r24254, r24255, r24256, r24257, r24258, r24259, r24260, r24261, r24262, r24263, r24264, r24265, r24266, r24267, r24268, r24269; void setup_mpfr_f_dm() { mpfr_set_default_prec(2128); mpfr_init(r24246); mpfr_init(r24247); mpfr_init(r24248); mpfr_init(r24249); mpfr_init(r24250); mpfr_init(r24251); mpfr_init(r24252); mpfr_init(r24253); mpfr_init(r24254); mpfr_init(r24255); mpfr_init(r24256); mpfr_init(r24257); mpfr_init(r24258); mpfr_init(r24259); mpfr_init(r24260); mpfr_init(r24261); mpfr_init(r24262); mpfr_init(r24263); mpfr_init(r24264); mpfr_init(r24265); mpfr_init(r24266); mpfr_init(r24267); mpfr_init(r24268); mpfr_init(r24269); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24246, phi1, MPFR_RNDN); mpfr_sin(r24247, r24246, MPFR_RNDN); mpfr_set_d(r24248, phi2, MPFR_RNDN); mpfr_sin(r24249, r24248, MPFR_RNDN); mpfr_mul(r24250, r24247, r24249, MPFR_RNDN); mpfr_cos(r24251, r24248, MPFR_RNDN); mpfr_cos(r24252, r24246, MPFR_RNDN); mpfr_mul(r24253, r24251, r24252, MPFR_RNDN); mpfr_exp(r24254, r24253, MPFR_RNDN); mpfr_set_d(r24255, lambda2, MPFR_RNDN); mpfr_cos(r24256, r24255, MPFR_RNDN); mpfr_set_d(r24257, lambda1, MPFR_RNDN); mpfr_cos(r24258, r24257, MPFR_RNDN); mpfr_mul(r24259, r24256, r24258, MPFR_RNDN); mpfr_sin(r24260, r24257, MPFR_RNDN); mpfr_sin(r24261, r24255, MPFR_RNDN); mpfr_mul(r24262, r24260, r24261, MPFR_RNDN); mpfr_add(r24263, r24259, r24262, MPFR_RNDN); mpfr_pow(r24264, r24254, r24263, MPFR_RNDN); mpfr_log(r24265, r24264, MPFR_RNDN); mpfr_add(r24266, r24250, r24265, MPFR_RNDN); mpfr_acos(r24267, r24266, MPFR_RNDN); mpfr_set_d(r24268, R, MPFR_RNDN); mpfr_mul(r24269, r24267, r24268, MPFR_RNDN); return mpfr_get_d(r24269, MPFR_RNDN); }