#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 r24095 = phi1; float r24096 = sin(r24095); float r24097 = phi2; float r24098 = sin(r24097); float r24099 = r24096 * r24098; float r24100 = cos(r24095); float r24101 = cos(r24097); float r24102 = r24100 * r24101; float r24103 = lambda1; float r24104 = lambda2; float r24105 = r24103 - r24104; float r24106 = cos(r24105); float r24107 = r24102 * r24106; float r24108 = r24099 + r24107; float r24109 = acos(r24108); float r24110 = R; float r24111 = r24109 * r24110; return r24111; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r24112 = phi1; double r24113 = sin(r24112); double r24114 = phi2; double r24115 = sin(r24114); double r24116 = r24113 * r24115; double r24117 = cos(r24112); double r24118 = cos(r24114); double r24119 = r24117 * r24118; double r24120 = lambda1; double r24121 = lambda2; double r24122 = r24120 - r24121; double r24123 = cos(r24122); double r24124 = r24119 * r24123; double r24125 = r24116 + r24124; double r24126 = acos(r24125); double r24127 = R; double r24128 = r24126 * r24127; return r24128; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r24129 = lambda1; float r24130 = cos(r24129); float r24131 = lambda2; float r24132 = cos(r24131); float r24133 = sin(r24129); float r24134 = sin(r24131); float r24135 = r24133 * r24134; float r24136 = fma(r24130, r24132, r24135); float r24137 = phi1; float r24138 = cos(r24137); float r24139 = phi2; float r24140 = cos(r24139); float r24141 = r24138 * r24140; float r24142 = sin(r24137); float r24143 = sin(r24139); float r24144 = r24142 * r24143; float r24145 = fma(r24136, r24141, r24144); float r24146 = acos(r24145); float r24147 = log(r24146); float r24148 = expm1(r24147); float r24149 = log1p(r24148); float r24150 = exp(r24149); float r24151 = R; float r24152 = r24150 * r24151; return r24152; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r24153 = lambda1; double r24154 = cos(r24153); double r24155 = lambda2; double r24156 = cos(r24155); double r24157 = sin(r24153); double r24158 = sin(r24155); double r24159 = r24157 * r24158; double r24160 = fma(r24154, r24156, r24159); double r24161 = phi1; double r24162 = cos(r24161); double r24163 = phi2; double r24164 = cos(r24163); double r24165 = r24162 * r24164; double r24166 = sin(r24161); double r24167 = sin(r24163); double r24168 = r24166 * r24167; double r24169 = fma(r24160, r24165, r24168); double r24170 = acos(r24169); double r24171 = log(r24170); double r24172 = expm1(r24171); double r24173 = log1p(r24172); double r24174 = exp(r24173); double r24175 = R; double r24176 = r24174 * r24175; return r24176; } 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 r24177, r24178, r24179, r24180, r24181, r24182, r24183, r24184, r24185, r24186, r24187, r24188, r24189, r24190, r24191, r24192, r24193; void setup_mpfr_f_im() { mpfr_set_default_prec(1872); mpfr_init(r24177); mpfr_init(r24178); mpfr_init(r24179); mpfr_init(r24180); mpfr_init(r24181); mpfr_init(r24182); mpfr_init(r24183); mpfr_init(r24184); mpfr_init(r24185); mpfr_init(r24186); mpfr_init(r24187); mpfr_init(r24188); mpfr_init(r24189); mpfr_init(r24190); mpfr_init(r24191); mpfr_init(r24192); mpfr_init(r24193); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24177, phi1, MPFR_RNDN); mpfr_sin(r24178, r24177, MPFR_RNDN); mpfr_set_d(r24179, phi2, MPFR_RNDN); mpfr_sin(r24180, r24179, MPFR_RNDN); mpfr_mul(r24181, r24178, r24180, MPFR_RNDN); mpfr_cos(r24182, r24177, MPFR_RNDN); mpfr_cos(r24183, r24179, MPFR_RNDN); mpfr_mul(r24184, r24182, r24183, MPFR_RNDN); mpfr_set_d(r24185, lambda1, MPFR_RNDN); mpfr_set_d(r24186, lambda2, MPFR_RNDN); mpfr_sub(r24187, r24185, r24186, MPFR_RNDN); mpfr_cos(r24188, r24187, MPFR_RNDN); mpfr_mul(r24189, r24184, r24188, MPFR_RNDN); mpfr_add(r24190, r24181, r24189, MPFR_RNDN); mpfr_acos(r24191, r24190, MPFR_RNDN); mpfr_set_d(r24192, R, MPFR_RNDN); mpfr_mul(r24193, r24191, r24192, MPFR_RNDN); return mpfr_get_d(r24193, MPFR_RNDN); } static mpfr_t r24194, r24195, r24196, r24197, r24198, r24199, r24200, r24201, r24202, r24203, r24204, r24205, r24206, r24207, r24208, r24209, r24210, r24211, r24212, r24213, r24214, r24215, r24216, r24217; void setup_mpfr_f_fm() { mpfr_set_default_prec(1872); mpfr_init(r24194); mpfr_init(r24195); mpfr_init(r24196); mpfr_init(r24197); mpfr_init(r24198); mpfr_init(r24199); mpfr_init(r24200); mpfr_init(r24201); mpfr_init(r24202); mpfr_init(r24203); mpfr_init(r24204); 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); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24194, lambda1, MPFR_RNDN); mpfr_cos(r24195, r24194, MPFR_RNDN); mpfr_set_d(r24196, lambda2, MPFR_RNDN); mpfr_cos(r24197, r24196, MPFR_RNDN); mpfr_sin(r24198, r24194, MPFR_RNDN); mpfr_sin(r24199, r24196, MPFR_RNDN); mpfr_mul(r24200, r24198, r24199, MPFR_RNDN); mpfr_fma(r24201, r24195, r24197, r24200, MPFR_RNDN); mpfr_set_d(r24202, phi1, MPFR_RNDN); mpfr_cos(r24203, r24202, MPFR_RNDN); mpfr_set_d(r24204, phi2, MPFR_RNDN); mpfr_cos(r24205, r24204, MPFR_RNDN); mpfr_mul(r24206, r24203, r24205, MPFR_RNDN); mpfr_sin(r24207, r24202, MPFR_RNDN); mpfr_sin(r24208, r24204, MPFR_RNDN); mpfr_mul(r24209, r24207, r24208, MPFR_RNDN); mpfr_fma(r24210, r24201, r24206, r24209, MPFR_RNDN); mpfr_acos(r24211, r24210, MPFR_RNDN); mpfr_log(r24212, r24211, MPFR_RNDN); mpfr_expm1(r24213, r24212, MPFR_RNDN); mpfr_log1p(r24214, r24213, MPFR_RNDN); mpfr_exp(r24215, r24214, MPFR_RNDN); mpfr_set_d(r24216, R, MPFR_RNDN); mpfr_mul(r24217, r24215, r24216, MPFR_RNDN); return mpfr_get_d(r24217, MPFR_RNDN); } static mpfr_t r24218, r24219, r24220, r24221, r24222, r24223, r24224, r24225, r24226, r24227, r24228, r24229, r24230, r24231, r24232, r24233, r24234, r24235, r24236, r24237, r24238, r24239, r24240, r24241; void setup_mpfr_f_dm() { mpfr_set_default_prec(1872); mpfr_init(r24218); mpfr_init(r24219); mpfr_init(r24220); mpfr_init(r24221); 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); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24218, lambda1, MPFR_RNDN); mpfr_cos(r24219, r24218, MPFR_RNDN); mpfr_set_d(r24220, lambda2, MPFR_RNDN); mpfr_cos(r24221, r24220, MPFR_RNDN); mpfr_sin(r24222, r24218, MPFR_RNDN); mpfr_sin(r24223, r24220, MPFR_RNDN); mpfr_mul(r24224, r24222, r24223, MPFR_RNDN); mpfr_fma(r24225, r24219, r24221, r24224, MPFR_RNDN); mpfr_set_d(r24226, phi1, MPFR_RNDN); mpfr_cos(r24227, r24226, MPFR_RNDN); mpfr_set_d(r24228, phi2, MPFR_RNDN); mpfr_cos(r24229, r24228, MPFR_RNDN); mpfr_mul(r24230, r24227, r24229, MPFR_RNDN); mpfr_sin(r24231, r24226, MPFR_RNDN); mpfr_sin(r24232, r24228, MPFR_RNDN); mpfr_mul(r24233, r24231, r24232, MPFR_RNDN); mpfr_fma(r24234, r24225, r24230, r24233, MPFR_RNDN); mpfr_acos(r24235, r24234, MPFR_RNDN); mpfr_log(r24236, r24235, MPFR_RNDN); mpfr_expm1(r24237, r24236, MPFR_RNDN); mpfr_log1p(r24238, r24237, MPFR_RNDN); mpfr_exp(r24239, r24238, MPFR_RNDN); mpfr_set_d(r24240, R, MPFR_RNDN); mpfr_mul(r24241, r24239, r24240, MPFR_RNDN); return mpfr_get_d(r24241, MPFR_RNDN); }