#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r25999 = 1; float r26000 = 0.1049934947; float r26001 = x; float r26002 = r26001 * r26001; float r26003 = r26000 * r26002; float r26004 = r25999 + r26003; float r26005 = 0.0424060604; float r26006 = r26002 * r26002; float r26007 = r26005 * r26006; float r26008 = r26004 + r26007; float r26009 = 0.0072644182; float r26010 = r26006 * r26002; float r26011 = r26009 * r26010; float r26012 = r26008 + r26011; float r26013 = 0.0005064034; float r26014 = r26010 * r26002; float r26015 = r26013 * r26014; float r26016 = r26012 + r26015; float r26017 = 0.0001789971; float r26018 = r26014 * r26002; float r26019 = r26017 * r26018; float r26020 = r26016 + r26019; float r26021 = 0.7715471019; float r26022 = r26021 * r26002; float r26023 = r25999 + r26022; float r26024 = 0.2909738639; float r26025 = r26024 * r26006; float r26026 = r26023 + r26025; float r26027 = 0.0694555761; float r26028 = r26027 * r26010; float r26029 = r26026 + r26028; float r26030 = 0.0140005442; float r26031 = r26030 * r26014; float r26032 = r26029 + r26031; float r26033 = 0.0008327945; float r26034 = r26033 * r26018; float r26035 = r26032 + r26034; float r26036 = 2; float r26037 = r26036 * r26017; float r26038 = r26018 * r26002; float r26039 = r26037 * r26038; float r26040 = r26035 + r26039; float r26041 = r26020 / r26040; float r26042 = r26041 * r26001; return r26042; } double f_id(double x) { double r26043 = 1; double r26044 = 0.1049934947; double r26045 = x; double r26046 = r26045 * r26045; double r26047 = r26044 * r26046; double r26048 = r26043 + r26047; double r26049 = 0.0424060604; double r26050 = r26046 * r26046; double r26051 = r26049 * r26050; double r26052 = r26048 + r26051; double r26053 = 0.0072644182; double r26054 = r26050 * r26046; double r26055 = r26053 * r26054; double r26056 = r26052 + r26055; double r26057 = 0.0005064034; double r26058 = r26054 * r26046; double r26059 = r26057 * r26058; double r26060 = r26056 + r26059; double r26061 = 0.0001789971; double r26062 = r26058 * r26046; double r26063 = r26061 * r26062; double r26064 = r26060 + r26063; double r26065 = 0.7715471019; double r26066 = r26065 * r26046; double r26067 = r26043 + r26066; double r26068 = 0.2909738639; double r26069 = r26068 * r26050; double r26070 = r26067 + r26069; double r26071 = 0.0694555761; double r26072 = r26071 * r26054; double r26073 = r26070 + r26072; double r26074 = 0.0140005442; double r26075 = r26074 * r26058; double r26076 = r26073 + r26075; double r26077 = 0.0008327945; double r26078 = r26077 * r26062; double r26079 = r26076 + r26078; double r26080 = 2; double r26081 = r26080 * r26061; double r26082 = r26062 * r26046; double r26083 = r26081 * r26082; double r26084 = r26079 + r26083; double r26085 = r26064 / r26084; double r26086 = r26085 * r26045; return r26086; } double f_of(float x) { float r26087 = x; float r26088 = -1.4994409406823423e+25; bool r26089 = r26087 <= r26088; float r26090 = 6155815.439411731; bool r26091 = r26087 <= r26090; float r26092 = !r26091; bool r26093 = r26089 || r26092; float r26094 = 0.2514179000665375; float r26095 = r26094 * r26087; float r26096 = 4; float r26097 = pow(r26087, r26096); float r26098 = r26095 / r26097; float r26099 = 0.15298196345929327; float r26100 = r26087 * r26099; float r26101 = 6; float r26102 = pow(r26087, r26101); float r26103 = r26100 / r26102; float r26104 = r26098 + r26103; float r26105 = 0.5; float r26106 = r26105 / r26087; float r26107 = r26104 + r26106; float r26108 = r26087 * r26087; float r26109 = r26108 * r26108; float r26110 = r26108 * r26109; float r26111 = r26108 * r26110; float r26112 = r26108 * r26111; float r26113 = 0.0001789971; float r26114 = r26112 * r26113; float r26115 = 0.0005064034; float r26116 = r26111 * r26115; float r26117 = 0.0424060604; float r26118 = r26109 * r26117; float r26119 = 1; float r26120 = 0.1049934947; float r26121 = r26120 * r26108; float r26122 = r26119 + r26121; float r26123 = r26118 + r26122; float r26124 = 0.0072644182; float r26125 = r26110 * r26124; float r26126 = r26123 + r26125; float r26127 = r26116 + r26126; float r26128 = r26114 + r26127; float r26129 = 0.0008327945; float r26130 = r26112 * r26129; float r26131 = 0.0140005442; float r26132 = r26131 * r26111; float r26133 = 0.0694555761; float r26134 = r26133 * r26110; float r26135 = 0.7715471019; float r26136 = r26108 * r26135; float r26137 = r26136 + r26119; float r26138 = 0.2909738639; float r26139 = r26109 * r26138; float r26140 = r26137 + r26139; float r26141 = r26134 + r26140; float r26142 = r26132 + r26141; float r26143 = r26130 + r26142; float r26144 = 2; float r26145 = r26113 * r26144; float r26146 = r26112 * r26108; float r26147 = r26145 * r26146; float r26148 = r26143 + r26147; float r26149 = r26128 / r26148; float r26150 = r26087 * r26149; float r26151 = r26093 ? r26107 : r26150; return r26151; } double f_od(double x) { double r26152 = x; double r26153 = -1.4994409406823423e+25; bool r26154 = r26152 <= r26153; double r26155 = 6155815.439411731; bool r26156 = r26152 <= r26155; double r26157 = !r26156; bool r26158 = r26154 || r26157; double r26159 = 0.2514179000665375; double r26160 = r26159 * r26152; double r26161 = 4; double r26162 = pow(r26152, r26161); double r26163 = r26160 / r26162; double r26164 = 0.15298196345929327; double r26165 = r26152 * r26164; double r26166 = 6; double r26167 = pow(r26152, r26166); double r26168 = r26165 / r26167; double r26169 = r26163 + r26168; double r26170 = 0.5; double r26171 = r26170 / r26152; double r26172 = r26169 + r26171; double r26173 = r26152 * r26152; double r26174 = r26173 * r26173; double r26175 = r26173 * r26174; double r26176 = r26173 * r26175; double r26177 = r26173 * r26176; double r26178 = 0.0001789971; double r26179 = r26177 * r26178; double r26180 = 0.0005064034; double r26181 = r26176 * r26180; double r26182 = 0.0424060604; double r26183 = r26174 * r26182; double r26184 = 1; double r26185 = 0.1049934947; double r26186 = r26185 * r26173; double r26187 = r26184 + r26186; double r26188 = r26183 + r26187; double r26189 = 0.0072644182; double r26190 = r26175 * r26189; double r26191 = r26188 + r26190; double r26192 = r26181 + r26191; double r26193 = r26179 + r26192; double r26194 = 0.0008327945; double r26195 = r26177 * r26194; double r26196 = 0.0140005442; double r26197 = r26196 * r26176; double r26198 = 0.0694555761; double r26199 = r26198 * r26175; double r26200 = 0.7715471019; double r26201 = r26173 * r26200; double r26202 = r26201 + r26184; double r26203 = 0.2909738639; double r26204 = r26174 * r26203; double r26205 = r26202 + r26204; double r26206 = r26199 + r26205; double r26207 = r26197 + r26206; double r26208 = r26195 + r26207; double r26209 = 2; double r26210 = r26178 * r26209; double r26211 = r26177 * r26173; double r26212 = r26210 * r26211; double r26213 = r26208 + r26212; double r26214 = r26193 / r26213; double r26215 = r26152 * r26214; double r26216 = r26158 ? r26172 : r26215; return r26216; } 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 r26217, r26218, r26219, r26220, r26221, r26222, r26223, r26224, r26225, r26226, r26227, r26228, r26229, r26230, r26231, r26232, r26233, r26234, r26235, r26236, r26237, r26238, r26239, r26240, r26241, r26242, r26243, r26244, r26245, r26246, r26247, r26248, r26249, r26250, r26251, r26252, r26253, r26254, r26255, r26256, r26257, r26258, r26259, r26260; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r26217, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26218, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r26219); mpfr_init(r26220); mpfr_init(r26221); mpfr_init(r26222); mpfr_init_set_str(r26223, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r26224); mpfr_init(r26225); mpfr_init(r26226); mpfr_init_set_str(r26227, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r26228); mpfr_init(r26229); mpfr_init(r26230); mpfr_init_set_str(r26231, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r26232); mpfr_init(r26233); mpfr_init(r26234); mpfr_init_set_str(r26235, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r26236); mpfr_init(r26237); mpfr_init(r26238); mpfr_init_set_str(r26239, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r26240); mpfr_init(r26241); mpfr_init_set_str(r26242, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r26243); mpfr_init(r26244); mpfr_init_set_str(r26245, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r26246); mpfr_init(r26247); mpfr_init_set_str(r26248, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r26249); mpfr_init(r26250); mpfr_init_set_str(r26251, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r26252); mpfr_init(r26253); mpfr_init_set_str(r26254, "2", 10, MPFR_RNDN); mpfr_init(r26255); mpfr_init(r26256); mpfr_init(r26257); mpfr_init(r26258); mpfr_init(r26259); mpfr_init(r26260); } double f_im(double x) { ; ; mpfr_set_d(r26219, x, MPFR_RNDN); mpfr_mul(r26220, r26219, r26219, MPFR_RNDN); mpfr_mul(r26221, r26218, r26220, MPFR_RNDN); mpfr_add(r26222, r26217, r26221, MPFR_RNDN); ; mpfr_mul(r26224, r26220, r26220, MPFR_RNDN); mpfr_mul(r26225, r26223, r26224, MPFR_RNDN); mpfr_add(r26226, r26222, r26225, MPFR_RNDN); ; mpfr_mul(r26228, r26224, r26220, MPFR_RNDN); mpfr_mul(r26229, r26227, r26228, MPFR_RNDN); mpfr_add(r26230, r26226, r26229, MPFR_RNDN); ; mpfr_mul(r26232, r26228, r26220, MPFR_RNDN); mpfr_mul(r26233, r26231, r26232, MPFR_RNDN); mpfr_add(r26234, r26230, r26233, MPFR_RNDN); ; mpfr_mul(r26236, r26232, r26220, MPFR_RNDN); mpfr_mul(r26237, r26235, r26236, MPFR_RNDN); mpfr_add(r26238, r26234, r26237, MPFR_RNDN); ; mpfr_mul(r26240, r26239, r26220, MPFR_RNDN); mpfr_add(r26241, r26217, r26240, MPFR_RNDN); ; mpfr_mul(r26243, r26242, r26224, MPFR_RNDN); mpfr_add(r26244, r26241, r26243, MPFR_RNDN); ; mpfr_mul(r26246, r26245, r26228, MPFR_RNDN); mpfr_add(r26247, r26244, r26246, MPFR_RNDN); ; mpfr_mul(r26249, r26248, r26232, MPFR_RNDN); mpfr_add(r26250, r26247, r26249, MPFR_RNDN); ; mpfr_mul(r26252, r26251, r26236, MPFR_RNDN); mpfr_add(r26253, r26250, r26252, MPFR_RNDN); ; mpfr_mul(r26255, r26254, r26235, MPFR_RNDN); mpfr_mul(r26256, r26236, r26220, MPFR_RNDN); mpfr_mul(r26257, r26255, r26256, MPFR_RNDN); mpfr_add(r26258, r26253, r26257, MPFR_RNDN); mpfr_div(r26259, r26238, r26258, MPFR_RNDN); mpfr_mul(r26260, r26259, r26219, MPFR_RNDN); return mpfr_get_d(r26260, MPFR_RNDN); } static mpfr_t r26261, r26262, r26263, r26264, r26265, r26266, r26267, r26268, r26269, r26270, r26271, r26272, r26273, r26274, r26275, r26276, r26277, r26278, r26279, r26280, r26281, r26282, r26283, r26284, r26285, r26286, r26287, r26288, r26289, r26290, r26291, r26292, r26293, r26294, r26295, r26296, r26297, r26298, r26299, r26300, r26301, r26302, r26303, r26304, r26305, r26306, r26307, r26308, r26309, r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318, r26319, r26320, r26321, r26322, r26323, r26324, r26325; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26261); mpfr_init_set_str(r26262, "-1.4994409406823423e+25", 10, MPFR_RNDN); mpfr_init(r26263); mpfr_init_set_str(r26264, "6155815.439411731", 10, MPFR_RNDN); mpfr_init(r26265); mpfr_init(r26266); mpfr_init(r26267); mpfr_init_set_str(r26268, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r26269); mpfr_init_set_str(r26270, "4", 10, MPFR_RNDN); mpfr_init(r26271); mpfr_init(r26272); mpfr_init_set_str(r26273, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r26274); mpfr_init_set_str(r26275, "6", 10, MPFR_RNDN); mpfr_init(r26276); mpfr_init(r26277); mpfr_init(r26278); mpfr_init_set_str(r26279, "0.5", 10, MPFR_RNDN); mpfr_init(r26280); mpfr_init(r26281); mpfr_init(r26282); mpfr_init(r26283); mpfr_init(r26284); mpfr_init(r26285); mpfr_init(r26286); mpfr_init_set_str(r26287, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r26288); mpfr_init_set_str(r26289, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r26290); mpfr_init_set_str(r26291, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r26292); mpfr_init_set_str(r26293, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26294, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r26295); mpfr_init(r26296); mpfr_init(r26297); mpfr_init_set_str(r26298, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r26299); mpfr_init(r26300); mpfr_init(r26301); mpfr_init(r26302); mpfr_init_set_str(r26303, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r26304); mpfr_init_set_str(r26305, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r26306); mpfr_init_set_str(r26307, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r26308); mpfr_init_set_str(r26309, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r26310); mpfr_init(r26311); mpfr_init_set_str(r26312, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r26313); mpfr_init(r26314); mpfr_init(r26315); mpfr_init(r26316); mpfr_init(r26317); mpfr_init_set_str(r26318, "2", 10, MPFR_RNDN); mpfr_init(r26319); mpfr_init(r26320); mpfr_init(r26321); mpfr_init(r26322); mpfr_init(r26323); mpfr_init(r26324); mpfr_init(r26325); } double f_fm(double x) { mpfr_set_d(r26261, x, MPFR_RNDN); ; mpfr_set_si(r26263, mpfr_cmp(r26261, r26262) <= 0, MPFR_RNDN); ; mpfr_set_si(r26265, mpfr_cmp(r26261, r26264) <= 0, MPFR_RNDN); mpfr_set_si(r26266, !mpfr_get_si(r26265, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r26267, mpfr_get_si(r26263, MPFR_RNDN) || mpfr_get_si(r26266, MPFR_RNDN), MPFR_RNDN); ; mpfr_mul(r26269, r26268, r26261, MPFR_RNDN); ; mpfr_pow(r26271, r26261, r26270, MPFR_RNDN); mpfr_div(r26272, r26269, r26271, MPFR_RNDN); ; mpfr_mul(r26274, r26261, r26273, MPFR_RNDN); ; mpfr_pow(r26276, r26261, r26275, MPFR_RNDN); mpfr_div(r26277, r26274, r26276, MPFR_RNDN); mpfr_add(r26278, r26272, r26277, MPFR_RNDN); ; mpfr_div(r26280, r26279, r26261, MPFR_RNDN); mpfr_add(r26281, r26278, r26280, MPFR_RNDN); mpfr_mul(r26282, r26261, r26261, MPFR_RNDN); mpfr_mul(r26283, r26282, r26282, MPFR_RNDN); mpfr_mul(r26284, r26282, r26283, MPFR_RNDN); mpfr_mul(r26285, r26282, r26284, MPFR_RNDN); mpfr_mul(r26286, r26282, r26285, MPFR_RNDN); ; mpfr_mul(r26288, r26286, r26287, MPFR_RNDN); ; mpfr_mul(r26290, r26285, r26289, MPFR_RNDN); ; mpfr_mul(r26292, r26283, r26291, MPFR_RNDN); ; ; mpfr_mul(r26295, r26294, r26282, MPFR_RNDN); mpfr_add(r26296, r26293, r26295, MPFR_RNDN); mpfr_add(r26297, r26292, r26296, MPFR_RNDN); ; mpfr_mul(r26299, r26284, r26298, MPFR_RNDN); mpfr_add(r26300, r26297, r26299, MPFR_RNDN); mpfr_add(r26301, r26290, r26300, MPFR_RNDN); mpfr_add(r26302, r26288, r26301, MPFR_RNDN); ; mpfr_mul(r26304, r26286, r26303, MPFR_RNDN); ; mpfr_mul(r26306, r26305, r26285, MPFR_RNDN); ; mpfr_mul(r26308, r26307, r26284, MPFR_RNDN); ; mpfr_mul(r26310, r26282, r26309, MPFR_RNDN); mpfr_add(r26311, r26310, r26293, MPFR_RNDN); ; mpfr_mul(r26313, r26283, r26312, MPFR_RNDN); mpfr_add(r26314, r26311, r26313, MPFR_RNDN); mpfr_add(r26315, r26308, r26314, MPFR_RNDN); mpfr_add(r26316, r26306, r26315, MPFR_RNDN); mpfr_add(r26317, r26304, r26316, MPFR_RNDN); ; mpfr_mul(r26319, r26287, r26318, MPFR_RNDN); mpfr_mul(r26320, r26286, r26282, MPFR_RNDN); mpfr_mul(r26321, r26319, r26320, MPFR_RNDN); mpfr_add(r26322, r26317, r26321, MPFR_RNDN); mpfr_div(r26323, r26302, r26322, MPFR_RNDN); mpfr_mul(r26324, r26261, r26323, MPFR_RNDN); if (mpfr_get_si(r26267, MPFR_RNDN)) { mpfr_set(r26325, r26281, MPFR_RNDN); } else { mpfr_set(r26325, r26324, MPFR_RNDN); }; return mpfr_get_d(r26325, MPFR_RNDN); } static mpfr_t r26326, r26327, r26328, r26329, r26330, r26331, r26332, r26333, r26334, r26335, r26336, r26337, r26338, r26339, r26340, r26341, r26342, r26343, r26344, r26345, r26346, r26347, r26348, r26349, r26350, r26351, r26352, r26353, r26354, r26355, r26356, r26357, r26358, r26359, r26360, r26361, r26362, r26363, r26364, r26365, r26366, r26367, r26368, r26369, r26370, r26371, r26372, r26373, r26374, r26375, r26376, r26377, r26378, r26379, r26380, r26381, r26382, r26383, r26384, r26385, r26386, r26387, r26388, r26389, r26390; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26326); mpfr_init_set_str(r26327, "-1.4994409406823423e+25", 10, MPFR_RNDN); mpfr_init(r26328); mpfr_init_set_str(r26329, "6155815.439411731", 10, MPFR_RNDN); mpfr_init(r26330); mpfr_init(r26331); mpfr_init(r26332); mpfr_init_set_str(r26333, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r26334); mpfr_init_set_str(r26335, "4", 10, MPFR_RNDN); mpfr_init(r26336); mpfr_init(r26337); mpfr_init_set_str(r26338, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r26339); mpfr_init_set_str(r26340, "6", 10, MPFR_RNDN); mpfr_init(r26341); mpfr_init(r26342); mpfr_init(r26343); mpfr_init_set_str(r26344, "0.5", 10, MPFR_RNDN); mpfr_init(r26345); mpfr_init(r26346); mpfr_init(r26347); mpfr_init(r26348); mpfr_init(r26349); mpfr_init(r26350); mpfr_init(r26351); mpfr_init_set_str(r26352, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r26353); mpfr_init_set_str(r26354, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r26355); mpfr_init_set_str(r26356, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r26357); mpfr_init_set_str(r26358, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26359, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r26360); mpfr_init(r26361); mpfr_init(r26362); mpfr_init_set_str(r26363, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r26364); mpfr_init(r26365); mpfr_init(r26366); mpfr_init(r26367); mpfr_init_set_str(r26368, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r26369); mpfr_init_set_str(r26370, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r26371); mpfr_init_set_str(r26372, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r26373); mpfr_init_set_str(r26374, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r26375); mpfr_init(r26376); mpfr_init_set_str(r26377, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r26378); mpfr_init(r26379); mpfr_init(r26380); mpfr_init(r26381); mpfr_init(r26382); mpfr_init_set_str(r26383, "2", 10, MPFR_RNDN); mpfr_init(r26384); mpfr_init(r26385); mpfr_init(r26386); mpfr_init(r26387); mpfr_init(r26388); mpfr_init(r26389); mpfr_init(r26390); } double f_dm(double x) { mpfr_set_d(r26326, x, MPFR_RNDN); ; mpfr_set_si(r26328, mpfr_cmp(r26326, r26327) <= 0, MPFR_RNDN); ; mpfr_set_si(r26330, mpfr_cmp(r26326, r26329) <= 0, MPFR_RNDN); mpfr_set_si(r26331, !mpfr_get_si(r26330, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r26332, mpfr_get_si(r26328, MPFR_RNDN) || mpfr_get_si(r26331, MPFR_RNDN), MPFR_RNDN); ; mpfr_mul(r26334, r26333, r26326, MPFR_RNDN); ; mpfr_pow(r26336, r26326, r26335, MPFR_RNDN); mpfr_div(r26337, r26334, r26336, MPFR_RNDN); ; mpfr_mul(r26339, r26326, r26338, MPFR_RNDN); ; mpfr_pow(r26341, r26326, r26340, MPFR_RNDN); mpfr_div(r26342, r26339, r26341, MPFR_RNDN); mpfr_add(r26343, r26337, r26342, MPFR_RNDN); ; mpfr_div(r26345, r26344, r26326, MPFR_RNDN); mpfr_add(r26346, r26343, r26345, MPFR_RNDN); mpfr_mul(r26347, r26326, r26326, MPFR_RNDN); mpfr_mul(r26348, r26347, r26347, MPFR_RNDN); mpfr_mul(r26349, r26347, r26348, MPFR_RNDN); mpfr_mul(r26350, r26347, r26349, MPFR_RNDN); mpfr_mul(r26351, r26347, r26350, MPFR_RNDN); ; mpfr_mul(r26353, r26351, r26352, MPFR_RNDN); ; mpfr_mul(r26355, r26350, r26354, MPFR_RNDN); ; mpfr_mul(r26357, r26348, r26356, MPFR_RNDN); ; ; mpfr_mul(r26360, r26359, r26347, MPFR_RNDN); mpfr_add(r26361, r26358, r26360, MPFR_RNDN); mpfr_add(r26362, r26357, r26361, MPFR_RNDN); ; mpfr_mul(r26364, r26349, r26363, MPFR_RNDN); mpfr_add(r26365, r26362, r26364, MPFR_RNDN); mpfr_add(r26366, r26355, r26365, MPFR_RNDN); mpfr_add(r26367, r26353, r26366, MPFR_RNDN); ; mpfr_mul(r26369, r26351, r26368, MPFR_RNDN); ; mpfr_mul(r26371, r26370, r26350, MPFR_RNDN); ; mpfr_mul(r26373, r26372, r26349, MPFR_RNDN); ; mpfr_mul(r26375, r26347, r26374, MPFR_RNDN); mpfr_add(r26376, r26375, r26358, MPFR_RNDN); ; mpfr_mul(r26378, r26348, r26377, MPFR_RNDN); mpfr_add(r26379, r26376, r26378, MPFR_RNDN); mpfr_add(r26380, r26373, r26379, MPFR_RNDN); mpfr_add(r26381, r26371, r26380, MPFR_RNDN); mpfr_add(r26382, r26369, r26381, MPFR_RNDN); ; mpfr_mul(r26384, r26352, r26383, MPFR_RNDN); mpfr_mul(r26385, r26351, r26347, MPFR_RNDN); mpfr_mul(r26386, r26384, r26385, MPFR_RNDN); mpfr_add(r26387, r26382, r26386, MPFR_RNDN); mpfr_div(r26388, r26367, r26387, MPFR_RNDN); mpfr_mul(r26389, r26326, r26388, MPFR_RNDN); if (mpfr_get_si(r26332, MPFR_RNDN)) { mpfr_set(r26390, r26346, MPFR_RNDN); } else { mpfr_set(r26390, r26389, MPFR_RNDN); }; return mpfr_get_d(r26390, MPFR_RNDN); }