#include #include #include #include #include char *name = "Henrywood and Agarwal, Equation (9a)"; double f_if(float w0, float M, float D, float h, float l, float d) { float r10133 = w0; float r10134 = 1; float r10135 = M; float r10136 = D; float r10137 = r10135 * r10136; float r10138 = 2; float r10139 = d; float r10140 = r10138 * r10139; float r10141 = r10137 / r10140; float r10142 = pow(r10141, r10138); float r10143 = h; float r10144 = l; float r10145 = r10143 / r10144; float r10146 = r10142 * r10145; float r10147 = r10134 - r10146; float r10148 = sqrt(r10147); float r10149 = r10133 * r10148; return r10149; } double f_id(double w0, double M, double D, double h, double l, double d) { double r10150 = w0; double r10151 = 1; double r10152 = M; double r10153 = D; double r10154 = r10152 * r10153; double r10155 = 2; double r10156 = d; double r10157 = r10155 * r10156; double r10158 = r10154 / r10157; double r10159 = pow(r10158, r10155); double r10160 = h; double r10161 = l; double r10162 = r10160 / r10161; double r10163 = r10159 * r10162; double r10164 = r10151 - r10163; double r10165 = sqrt(r10164); double r10166 = r10150 * r10165; return r10166; } double f_of(float w0, float M, float D, float h, float l, float d) { float r10167 = M; float r10168 = D; float r10169 = r10167 * r10168; float r10170 = 2; float r10171 = d; float r10172 = r10170 * r10171; float r10173 = r10169 / r10172; float r10174 = h; float r10175 = r10174 * r10173; float r10176 = r10173 * r10175; float r10177 = -3.4790670620477553e+245; bool r10178 = r10176 <= r10177; float r10179 = 1.629810084747972e+111; bool r10180 = r10176 <= r10179; float r10181 = !r10180; bool r10182 = r10178 || r10181; float r10183 = w0; float r10184 = 1; float r10185 = r10167 / r10171; float r10186 = 1/2; float r10187 = r10186 * r10168; float r10188 = r10185 * r10187; float r10189 = l; float r10190 = r10188 / r10189; float r10191 = r10168 / r10170; float r10192 = r10191 * r10174; float r10193 = r10185 * r10192; float r10194 = r10190 * r10193; float r10195 = r10184 - r10194; float r10196 = sqrt(r10195); float r10197 = r10183 * r10196; float r10198 = r10184 / r10189; float r10199 = r10198 * r10176; float r10200 = r10184 - r10199; float r10201 = sqrt(r10200); float r10202 = r10201 * r10183; float r10203 = r10182 ? r10197 : r10202; return r10203; } double f_od(double w0, double M, double D, double h, double l, double d) { double r10204 = M; double r10205 = D; double r10206 = r10204 * r10205; double r10207 = 2; double r10208 = d; double r10209 = r10207 * r10208; double r10210 = r10206 / r10209; double r10211 = h; double r10212 = r10211 * r10210; double r10213 = r10210 * r10212; double r10214 = -3.4790670620477553e+245; bool r10215 = r10213 <= r10214; double r10216 = 1.629810084747972e+111; bool r10217 = r10213 <= r10216; double r10218 = !r10217; bool r10219 = r10215 || r10218; double r10220 = w0; double r10221 = 1; double r10222 = r10204 / r10208; double r10223 = 1/2; double r10224 = r10223 * r10205; double r10225 = r10222 * r10224; double r10226 = l; double r10227 = r10225 / r10226; double r10228 = r10205 / r10207; double r10229 = r10228 * r10211; double r10230 = r10222 * r10229; double r10231 = r10227 * r10230; double r10232 = r10221 - r10231; double r10233 = sqrt(r10232); double r10234 = r10220 * r10233; double r10235 = r10221 / r10226; double r10236 = r10235 * r10213; double r10237 = r10221 - r10236; double r10238 = sqrt(r10237); double r10239 = r10238 * r10220; double r10240 = r10219 ? r10234 : r10239; return r10240; } 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 r10241, r10242, r10243, r10244, r10245, r10246, r10247, r10248, r10249, r10250, r10251, r10252, r10253, r10254, r10255, r10256, r10257; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r10241); mpfr_init_set_str(r10242, "1", 10, MPFR_RNDN); mpfr_init(r10243); mpfr_init(r10244); mpfr_init(r10245); mpfr_init_set_str(r10246, "2", 10, MPFR_RNDN); mpfr_init(r10247); mpfr_init(r10248); mpfr_init(r10249); mpfr_init(r10250); mpfr_init(r10251); mpfr_init(r10252); mpfr_init(r10253); mpfr_init(r10254); mpfr_init(r10255); mpfr_init(r10256); mpfr_init(r10257); } double f_im(double w0, double M, double D, double h, double l, double d) { mpfr_set_d(r10241, w0, MPFR_RNDN); ; mpfr_set_d(r10243, M, MPFR_RNDN); mpfr_set_d(r10244, D, MPFR_RNDN); mpfr_mul(r10245, r10243, r10244, MPFR_RNDN); ; mpfr_set_d(r10247, d, MPFR_RNDN); mpfr_mul(r10248, r10246, r10247, MPFR_RNDN); mpfr_div(r10249, r10245, r10248, MPFR_RNDN); mpfr_pow(r10250, r10249, r10246, MPFR_RNDN); mpfr_set_d(r10251, h, MPFR_RNDN); mpfr_set_d(r10252, l, MPFR_RNDN); mpfr_div(r10253, r10251, r10252, MPFR_RNDN); mpfr_mul(r10254, r10250, r10253, MPFR_RNDN); mpfr_sub(r10255, r10242, r10254, MPFR_RNDN); mpfr_sqrt(r10256, r10255, MPFR_RNDN); mpfr_mul(r10257, r10241, r10256, MPFR_RNDN); return mpfr_get_d(r10257, MPFR_RNDN); } static mpfr_t r10258, r10259, r10260, r10261, r10262, r10263, r10264, r10265, r10266, r10267, r10268, r10269, r10270, r10271, r10272, r10273, r10274, r10275, r10276, r10277, r10278, r10279, r10280, r10281, r10282, r10283, r10284, r10285, r10286, r10287, r10288, r10289, r10290, r10291, r10292, r10293, r10294; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r10258); mpfr_init(r10259); mpfr_init(r10260); mpfr_init_set_str(r10261, "2", 10, MPFR_RNDN); mpfr_init(r10262); mpfr_init(r10263); mpfr_init(r10264); mpfr_init(r10265); mpfr_init(r10266); mpfr_init(r10267); mpfr_init_set_str(r10268, "-3.4790670620477553e+245", 10, MPFR_RNDN); mpfr_init(r10269); mpfr_init_set_str(r10270, "1.629810084747972e+111", 10, MPFR_RNDN); mpfr_init(r10271); mpfr_init(r10272); mpfr_init(r10273); mpfr_init(r10274); mpfr_init_set_str(r10275, "1", 10, MPFR_RNDN); mpfr_init(r10276); mpfr_init_set_str(r10277, "1/2", 10, MPFR_RNDN); mpfr_init(r10278); mpfr_init(r10279); mpfr_init(r10280); mpfr_init(r10281); mpfr_init(r10282); mpfr_init(r10283); mpfr_init(r10284); mpfr_init(r10285); mpfr_init(r10286); mpfr_init(r10287); mpfr_init(r10288); mpfr_init(r10289); mpfr_init(r10290); mpfr_init(r10291); mpfr_init(r10292); mpfr_init(r10293); mpfr_init(r10294); } double f_fm(double w0, double M, double D, double h, double l, double d) { mpfr_set_d(r10258, M, MPFR_RNDN); mpfr_set_d(r10259, D, MPFR_RNDN); mpfr_mul(r10260, r10258, r10259, MPFR_RNDN); ; mpfr_set_d(r10262, d, MPFR_RNDN); mpfr_mul(r10263, r10261, r10262, MPFR_RNDN); mpfr_div(r10264, r10260, r10263, MPFR_RNDN); mpfr_set_d(r10265, h, MPFR_RNDN); mpfr_mul(r10266, r10265, r10264, MPFR_RNDN); mpfr_mul(r10267, r10264, r10266, MPFR_RNDN); ; mpfr_set_si(r10269, mpfr_cmp(r10267, r10268) <= 0, MPFR_RNDN); ; mpfr_set_si(r10271, mpfr_cmp(r10267, r10270) <= 0, MPFR_RNDN); mpfr_set_si(r10272, !mpfr_get_si(r10271, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r10273, mpfr_get_si(r10269, MPFR_RNDN) || mpfr_get_si(r10272, MPFR_RNDN), MPFR_RNDN); mpfr_set_d(r10274, w0, MPFR_RNDN); ; mpfr_div(r10276, r10258, r10262, MPFR_RNDN); ; mpfr_mul(r10278, r10277, r10259, MPFR_RNDN); mpfr_mul(r10279, r10276, r10278, MPFR_RNDN); mpfr_set_d(r10280, l, MPFR_RNDN); mpfr_div(r10281, r10279, r10280, MPFR_RNDN); mpfr_div(r10282, r10259, r10261, MPFR_RNDN); mpfr_mul(r10283, r10282, r10265, MPFR_RNDN); mpfr_mul(r10284, r10276, r10283, MPFR_RNDN); mpfr_mul(r10285, r10281, r10284, MPFR_RNDN); mpfr_sub(r10286, r10275, r10285, MPFR_RNDN); mpfr_sqrt(r10287, r10286, MPFR_RNDN); mpfr_mul(r10288, r10274, r10287, MPFR_RNDN); mpfr_div(r10289, r10275, r10280, MPFR_RNDN); mpfr_mul(r10290, r10289, r10267, MPFR_RNDN); mpfr_sub(r10291, r10275, r10290, MPFR_RNDN); mpfr_sqrt(r10292, r10291, MPFR_RNDN); mpfr_mul(r10293, r10292, r10274, MPFR_RNDN); if (mpfr_get_si(r10273, MPFR_RNDN)) { mpfr_set(r10294, r10288, MPFR_RNDN); } else { mpfr_set(r10294, r10293, MPFR_RNDN); }; return mpfr_get_d(r10294, MPFR_RNDN); } static mpfr_t r10295, r10296, r10297, r10298, r10299, r10300, r10301, r10302, r10303, r10304, r10305, r10306, r10307, r10308, r10309, r10310, r10311, r10312, r10313, r10314, r10315, r10316, r10317, r10318, r10319, r10320, r10321, r10322, r10323, r10324, r10325, r10326, r10327, r10328, r10329, r10330, r10331; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r10295); mpfr_init(r10296); mpfr_init(r10297); mpfr_init_set_str(r10298, "2", 10, MPFR_RNDN); mpfr_init(r10299); mpfr_init(r10300); mpfr_init(r10301); mpfr_init(r10302); mpfr_init(r10303); mpfr_init(r10304); mpfr_init_set_str(r10305, "-3.4790670620477553e+245", 10, MPFR_RNDN); mpfr_init(r10306); mpfr_init_set_str(r10307, "1.629810084747972e+111", 10, MPFR_RNDN); mpfr_init(r10308); mpfr_init(r10309); mpfr_init(r10310); mpfr_init(r10311); mpfr_init_set_str(r10312, "1", 10, MPFR_RNDN); mpfr_init(r10313); mpfr_init_set_str(r10314, "1/2", 10, MPFR_RNDN); mpfr_init(r10315); mpfr_init(r10316); mpfr_init(r10317); mpfr_init(r10318); mpfr_init(r10319); mpfr_init(r10320); mpfr_init(r10321); mpfr_init(r10322); mpfr_init(r10323); mpfr_init(r10324); mpfr_init(r10325); mpfr_init(r10326); mpfr_init(r10327); mpfr_init(r10328); mpfr_init(r10329); mpfr_init(r10330); mpfr_init(r10331); } double f_dm(double w0, double M, double D, double h, double l, double d) { mpfr_set_d(r10295, M, MPFR_RNDN); mpfr_set_d(r10296, D, MPFR_RNDN); mpfr_mul(r10297, r10295, r10296, MPFR_RNDN); ; mpfr_set_d(r10299, d, MPFR_RNDN); mpfr_mul(r10300, r10298, r10299, MPFR_RNDN); mpfr_div(r10301, r10297, r10300, MPFR_RNDN); mpfr_set_d(r10302, h, MPFR_RNDN); mpfr_mul(r10303, r10302, r10301, MPFR_RNDN); mpfr_mul(r10304, r10301, r10303, MPFR_RNDN); ; mpfr_set_si(r10306, mpfr_cmp(r10304, r10305) <= 0, MPFR_RNDN); ; mpfr_set_si(r10308, mpfr_cmp(r10304, r10307) <= 0, MPFR_RNDN); mpfr_set_si(r10309, !mpfr_get_si(r10308, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r10310, mpfr_get_si(r10306, MPFR_RNDN) || mpfr_get_si(r10309, MPFR_RNDN), MPFR_RNDN); mpfr_set_d(r10311, w0, MPFR_RNDN); ; mpfr_div(r10313, r10295, r10299, MPFR_RNDN); ; mpfr_mul(r10315, r10314, r10296, MPFR_RNDN); mpfr_mul(r10316, r10313, r10315, MPFR_RNDN); mpfr_set_d(r10317, l, MPFR_RNDN); mpfr_div(r10318, r10316, r10317, MPFR_RNDN); mpfr_div(r10319, r10296, r10298, MPFR_RNDN); mpfr_mul(r10320, r10319, r10302, MPFR_RNDN); mpfr_mul(r10321, r10313, r10320, MPFR_RNDN); mpfr_mul(r10322, r10318, r10321, MPFR_RNDN); mpfr_sub(r10323, r10312, r10322, MPFR_RNDN); mpfr_sqrt(r10324, r10323, MPFR_RNDN); mpfr_mul(r10325, r10311, r10324, MPFR_RNDN); mpfr_div(r10326, r10312, r10317, MPFR_RNDN); mpfr_mul(r10327, r10326, r10304, MPFR_RNDN); mpfr_sub(r10328, r10312, r10327, MPFR_RNDN); mpfr_sqrt(r10329, r10328, MPFR_RNDN); mpfr_mul(r10330, r10329, r10311, MPFR_RNDN); if (mpfr_get_si(r10310, MPFR_RNDN)) { mpfr_set(r10331, r10325, MPFR_RNDN); } else { mpfr_set(r10331, r10330, MPFR_RNDN); }; return mpfr_get_d(r10331, MPFR_RNDN); }