diff --git a/networks/metal_chem/actual_rhs.H b/networks/metal_chem/actual_rhs.H index 7d4fdb8c3..228dd22c1 100644 --- a/networks/metal_chem/actual_rhs.H +++ b/networks/metal_chem/actual_rhs.H @@ -54,13 +54,13 @@ void rhs_specie(const burn_t& state, Real x10 = 1.0/x4; - Real x11 = 5.7884371785482823e-10*x10*std::exp((-1.7524)*std::log(std::abs(0.00060040841663220993*x4 + 1.0)))*std::exp((-0.24759999999999999)*std::log(std::abs(0.32668576019240059*x4 + 1.0))) + 1.0e-99; + Real x11 = 5.7884371785482823e-10*x10*std::pow(0.00060040841663220993*x4 + 1.0, -1.7524)*std::pow(0.32668576019240059*x4 + 1.0, -0.24759999999999999) + 1.0e-99; Real x12 = X(13)*X(2); Real x13 = -x11*x12; - Real x14 = std::exp((-0.69999999999999996)*std::log(std::abs(x1))); + Real x14 = std::pow(x1, -0.69999999999999996); Real x15 = X(2)*X(31); @@ -72,7 +72,7 @@ void rhs_specie(const burn_t& state, Real x19 = -x17*x18; - Real x20 = 9.7355109978586414e-8*std::exp((-0.47999999999999998)*std::log(std::abs(x1))) + 1.0e-99; + Real x20 = 9.7355109978586414e-8*std::pow(x1, -0.47999999999999998) + 1.0e-99; Real x21 = X(2)*X(24); @@ -96,7 +96,7 @@ void rhs_specie(const burn_t& state, Real x31 = std::exp((8)*std::log(std::abs(x24))); - Real x32 = std::exp((23.915965629999999)*std::log(std::abs(x1)))*std::exp(-941.91483008144996*x25 + 2506.9866529060901*x26 - 4361.9927099007555*x27 + 4879.7345146260486*x28 - 3366.4639698826941*x29 + 1300.3028484326148*x30 - 214.82451513312137*x31); + Real x32 = std::pow(x1, 23.915965629999999)*std::exp(-941.91483008144996*x25 + 2506.9866529060901*x26 - 4361.9927099007555*x27 + 4879.7345146260486*x28 - 3366.4639698826941*x29 + 1300.3028484326148*x30 - 214.82451513312137*x31); Real x33 = 4.3524079114767552e-117*x32 + 1.0e-99; @@ -104,7 +104,7 @@ void rhs_specie(const burn_t& state, Real x35 = -x33*x34; - Real x36 = std::exp((-0.75)*std::log(std::abs(x1))); + Real x36 = std::pow(x1, -0.75); Real x37 = 2.5950363272655348e-10*x36 + 1.0e-99; @@ -112,13 +112,13 @@ void rhs_specie(const burn_t& state, Real x39 = -x37*x38; - Real x40 = 6.3350403285625223e-6*std::exp((-0.55000000000000004)*std::log(std::abs(x1))) + 1.0e-99; + Real x40 = 6.3350403285625223e-6*std::pow(x1, -0.55000000000000004) + 1.0e-99; Real x41 = X(2)*X(29); Real x42 = -x40*x41; - Real x43 = std::exp((2.360852208681)*std::log(std::abs(x1)))*std::exp(-24.766609674457612*x25 + 13.307984239358756*x26 - 258.18559308467116*x27 + 846.15238706523724*x28 - 1113.0879095147111*x29 + 671.95094388835207*x30 - 154.90262957142161*x31); + Real x43 = std::pow(x1, 2.360852208681)*std::exp(-24.766609674457612*x25 + 13.307984239358756*x26 - 258.18559308467116*x27 + 846.15238706523724*x28 - 1113.0879095147111*x29 + 671.95094388835207*x30 - 154.90262957142161*x31); Real x44 = 3.7903999274394518e-18*x43 + 1.0e-99; @@ -126,13 +126,13 @@ void rhs_specie(const burn_t& state, Real x46 = -x44*x45; - Real x47 = 1.3300135414628029e-18*std::exp((0.94999999999999996)*std::log(std::abs(x1)))*std::exp(-0.00010729613733905579*x1) + 1.0e-99; + Real x47 = 1.3300135414628029e-18*std::pow(x1, 0.94999999999999996)*std::exp(-0.00010729613733905579*x1) + 1.0e-99; Real x48 = X(2)*X(7); Real x49 = -x47*x48; - Real x50 = std::exp((0.25)*std::log(std::abs(x1))); + Real x50 = std::pow(x1, 0.25); Real x51 = 1.0/x50; @@ -148,7 +148,7 @@ void rhs_specie(const burn_t& state, Real x57 = -x55*x56; - Real x58 = std::exp((-0.34000000000000002)*std::log(std::abs(x1))); + Real x58 = std::pow(x1, -0.34000000000000002); Real x59 = 157821.26811013557*x52; @@ -160,22 +160,22 @@ void rhs_specie(const burn_t& state, Real x63 = -x61*x62; - Real x64 = std::exp((0.34999999999999998)*std::log(std::abs(x1)))*std::exp(-102000.0*x52); + Real x64 = std::pow(x1, 0.34999999999999998)*std::exp(-102000.0*x52); Real x65 = 4.3799999999999999e-10*x64 + 1.0e-99; Real x66 = X(10)*X(2); - Real x67 = std::exp((43.933476326349997)*std::log(std::abs(x1)))*std::exp(-1618.789587733125*x25 + 3854.4033653120221*x26 - 5902.1601240760483*x27 + 5825.9326359379538*x28 - 3578.1439181805955*x29 + 1242.7294446825149*x30 - 186.35635455381879*x31); + Real x67 = std::pow(x1, 43.933476326349997)*std::exp(-1618.789587733125*x25 + 3854.4033653120221*x26 - 5902.1601240760483*x27 + 5825.9326359379538*x28 - 3578.1439181805955*x29 + 1242.7294446825149*x30 - 186.35635455381879*x31); Real x68 = X(14)*X(2); - Real x69 = std::exp((13.536555999999999)*std::log(std::abs(x1)))*std::exp(-502.72883252679094*x25 + 1281.477767828706*x26 - 2207.4643501257692*x27 + 2500.8077583366976*x28 - 1768.8867461266502*x29 + 704.19926629500367*x30 - 120.0438480494693*x31); + Real x69 = std::pow(x1, 13.536555999999999)*std::exp(-502.72883252679094*x25 + 1281.477767828706*x26 - 2207.4643501257692*x27 + 2500.8077583366976*x28 - 1768.8867461266502*x29 + 704.19926629500367*x30 - 120.0438480494693*x31); Real x70 = X(2)*X(4); Real x71 = ((x1 <= 1160.0) ? ( - 1.4643482606109061e-16*std::exp((1.78186)*std::log(std::abs(x1))) + 1.4643482606109061e-16*std::pow(x1, 1.78186) ) : ( 0 @@ -188,7 +188,7 @@ void rhs_specie(const burn_t& state, Real x74 = ((x24)*(x24)*(x24)*(x24)*(x24)*(x24)*(x24)*(x24)*(x24)); Real x75 = ((x1 > 1160.0) ? ( - 3.3178155742407601e-14*std::exp((1.139449335841631)*std::log(std::abs(x1)))*std::exp(-12.447178055372776*x25 + 6.9391784778399117*x26 - 10.993097527150175*x27 + 14.449862906216714*x28 + 58.228375789703179*x29 - 162.59852239006702*x30 + 144.55426734953477*x31 - 44.454280878123605*x74) + 3.3178155742407601e-14*std::pow(x1, 1.139449335841631)*std::exp(-12.447178055372776*x25 + 6.9391784778399117*x26 - 10.993097527150175*x27 + 14.449862906216714*x28 + 58.228375789703179*x29 - 162.59852239006702*x30 + 144.55426734953477*x31 - 44.454280878123605*x74) ) : ( 0 @@ -227,14 +227,14 @@ void rhs_specie(const burn_t& state, Real x89 = X(2)*X(8); - Real x90 = 1.0e-8*std::exp((-0.40000000000000002)*std::log(std::abs(x1))) + 1.0e-99; + Real x90 = 1.0e-8*std::pow(x1, -0.40000000000000002) + 1.0e-99; Real x91 = X(3)*X(5); Real x92 = x89*(-x88 - 1.0e-99) + x90*x91; Real x93 = ((x1 <= 400.0) ? ( - 1.2999999999999999e-10*std::exp((-0.64000000000000001)*std::log(std::abs(x1))) + 1.2999999999999999e-10*std::pow(x1, -0.64000000000000001) ) : ( 0 @@ -242,10 +242,10 @@ void rhs_specie(const burn_t& state, Real x94 = X(2)*X(21); - Real x95 = std::exp((-1.5)*std::log(std::abs(x1))); + Real x95 = std::pow(x1, -1.5); Real x96 = ((x1 > 400.0) ? ( - 1.41e-10*std::exp((-0.66000000000000003)*std::log(std::abs(x1))) + 0.00073999999999999999*x95*(1.0 + 0.062*std::exp(-145000.0*x52))*std::exp(-175000.0*x52) + 1.41e-10*std::pow(x1, -0.66000000000000003) + 0.00073999999999999999*x95*(1.0 + 0.062*std::exp(-145000.0*x52))*std::exp(-175000.0*x52) ) : ( 0 @@ -253,7 +253,7 @@ void rhs_specie(const burn_t& state, Real x97 = -x93*x94 - x94*x96; - Real x98 = 1.4981088130721367e-10*std::exp((-0.63529999999999998)*std::log(std::abs(x1))); + Real x98 = 1.4981088130721367e-10*std::pow(x1, -0.63529999999999998); Real x99 = x1 > 2.73; @@ -277,7 +277,7 @@ void rhs_specie(const burn_t& state, Real x104 = -x100*x68 - x102*x68 - x103*x68; - Real x105 = 2.6534040307116387e-9*std::exp((-0.10000000000000001)*std::log(std::abs(x1))) + 1.0e-99; + Real x105 = 2.6534040307116387e-9*std::pow(x1, -0.10000000000000001) + 1.0e-99; Real x106 = X(4)*X(9); @@ -286,7 +286,7 @@ void rhs_specie(const burn_t& state, Real x108 = x105*x106 + x105*x107; Real x109 = ((x1 <= 7950.0) ? ( - 1.4308352583277889e-10*std::exp((-0.59999999999999998)*std::log(std::abs(x1))) + 1.4308352583277889e-10*std::pow(x1, -0.59999999999999998) ) : ( 0 @@ -295,14 +295,14 @@ void rhs_specie(const burn_t& state, Real x110 = X(16)*X(2); Real x111 = ((x1 > 21140.0) ? ( - 0.00023813936486848612*std::exp((-1.3700000000000001)*std::log(std::abs(x1)))*std::exp(-115786.2*x52) + 0.00023813936486848612*std::pow(x1, -1.3700000000000001)*std::exp(-115786.2*x52) ) : ( 0 )) + 1.0e-99; Real x112 = ((x1 <= 21140.0 && x1 > 7950.0) ? ( - 8.3535906003219973e-24*std::exp((2.4900000000000002)*std::log(std::abs(x1)))*std::exp(21845.599999999999*x52) + 8.3535906003219973e-24*std::pow(x1, 2.4900000000000002)*std::exp(21845.599999999999*x52) ) : ( 0 @@ -320,7 +320,7 @@ void rhs_specie(const burn_t& state, Real x115 = X(2)*X(3); Real x116 = ((x101 && x1 > 5500.0) ? ( - 3.2867337024382733e-10*std::exp((-0.7241125657826851)*std::log(std::abs(x1)))*std::exp(-1.7746868094247411*x25 - 1.9518356165136789*x26 - 2.4649195146505537*x27 - 1.0207737270119371*x28 + 3.3530579587656565*x29 + 3.6203127646377788*x30 - 1.0930705283186734*x31 - 1.6921001126637108*x74) + 3.2867337024382733e-10*std::pow(x1, -0.7241125657826851)*std::exp(-1.7746868094247411*x25 - 1.9518356165136789*x26 - 2.4649195146505537*x27 - 1.0207737270119371*x28 + 3.3530579587656565*x29 + 3.6203127646377788*x30 - 1.0930705283186734*x31 - 1.6921001126637108*x74) ) : ( 0 @@ -328,15 +328,15 @@ void rhs_specie(const burn_t& state, Real x117 = -x114*x115 - x115*x116; - Real x118 = 35.5*std::exp((-2.2799999999999998)*std::log(std::abs(x1)))*std::exp(-46707.0*x52) + 1.0e-99; + Real x118 = 35.5*std::pow(x1, -2.2799999999999998)*std::exp(-46707.0*x52) + 1.0e-99; - Real x119 = (1.3500000000000001e-9*std::exp((0.098492999999999997)*std::log(std::abs(x1))) + 4.4350199999999998e-10*std::exp((0.55610000000000004)*std::log(std::abs(x1))) + 3.7408500000000004e-16*std::exp((2.1825999999999999)*std::log(std::abs(x1))))/(0.0061910000000000003*std::exp((1.0461)*std::log(std::abs(x1))) + 8.9711999999999997e-11*std::exp((3.0424000000000002)*std::log(std::abs(x1))) + 3.2575999999999999e-14*std::exp((3.7740999999999998)*std::log(std::abs(x1))) + 1.0) + 1.0e-99; + Real x119 = (1.3500000000000001e-9*std::pow(x1, 0.098492999999999997) + 4.4350199999999998e-10*std::pow(x1, 0.55610000000000004) + 3.7408500000000004e-16*std::pow(x1, 2.1825999999999999))/(0.0061910000000000003*std::pow(x1, 1.0461) + 8.9711999999999997e-11*std::pow(x1, 3.0424000000000002) + 3.2575999999999999e-14*std::pow(x1, 3.7740999999999998) + 1.0) + 1.0e-99; Real x120 = -x118*x66 + x119*x73 - x65*x66; Real x121 = 5.9082438637265071e-70*x69 + 1.0e-99; - Real x122 = 1.4000000000000001e-18*std::exp((0.92800000000000005)*std::log(std::abs(x1)))*std::exp(-6.1728395061728397e-5*x1) + 1.0e-99; + Real x122 = 1.4000000000000001e-18*std::pow(x1, 0.92800000000000005)*std::exp(-6.1728395061728397e-5*x1) + 1.0e-99; Real x123 = -x121*x70 - x122*x70; @@ -357,7 +357,7 @@ void rhs_specie(const burn_t& state, Real x131 = X(14)*X(4); Real x132 = ((x1 > 10000.0) ? ( - 4.0000000000000003e-37*std::exp((4.7400000000000002)*std::log(std::abs(x1))) + 4.0000000000000003e-37*std::pow(x1, 4.7400000000000002) ) : ( 0 @@ -429,7 +429,7 @@ void rhs_specie(const burn_t& state, Real x155 = -x90*x91 + x91*(-x154 - 1.0e-99); - Real x156 = std::exp((-0.12690000000000001*std::exp((-3.0)*std::log(std::abs(x137)))*((x23)*(x23)*(x23)) + 1.1180000000000001*std::exp((-2.0)*std::log(std::abs(x137)))*((x23)*(x23)) - 1.5229999999999999*x139 - 19.379999999999999)*std::log(std::abs(10.0))) + 1.0e-99; + Real x156 = std::pow(10.0, -0.12690000000000001*amrex::Math::powi<-3>(x137)*((x23)*(x23)*(x23)) + 1.1180000000000001*amrex::Math::powi<-2>(x137)*((x23)*(x23)) - 1.5229999999999999*x139 - 19.379999999999999) + 1.0e-99; Real x157 = -x156; @@ -438,14 +438,14 @@ void rhs_specie(const burn_t& state, Real x159 = x1 > 50.0; Real x160 = ((x159) ? ( - 2.0000000000000001e-10*std::exp((0.40200000000000002)*std::log(std::abs(x1)))*std::exp(-37.100000000000001*x52) - 3.3099999999999998e-17*std::exp((1.48)*std::log(std::abs(x1))) + 2.0000000000000001e-10*std::pow(x1, 0.40200000000000002)*std::exp(-37.100000000000001*x52) - 3.3099999999999998e-17*std::pow(x1, 1.48) ) : ( 0 )) + 1.0e-99; Real x161 = ((x159) ? ( - 2.0299999999999998e-9*std::exp((-0.33200000000000002)*std::log(std::abs(x1))) + 2.0600000000000001e-10*std::exp((0.39600000000000002)*std::log(std::abs(x1)))*std::exp(-33.0*x52) + 2.0299999999999998e-9*std::pow(x1, -0.33200000000000002) + 2.0600000000000001e-10*std::pow(x1, 0.39600000000000002)*std::exp(-33.0*x52) ) : ( 0 @@ -455,11 +455,11 @@ void rhs_specie(const burn_t& state, Real x163 = x157*x158 - x158*x160 + x161*x162; - Real x164 = 8.7882738061352424e-22*std::exp((1.96)*std::log(std::abs(x1)))*std::exp(-170000.0*x52) + 1.0e-99; + Real x164 = 8.7882738061352424e-22*std::pow(x1, 1.96)*std::exp(-170000.0*x52) + 1.0e-99; Real x165 = X(16)*X(4); - Real x166 = 3.8999999999999998e-16*std::exp((0.21299999999999999)*std::log(std::abs(x1))) + 1.0e-99; + Real x166 = 3.8999999999999998e-16*std::pow(x1, 0.21299999999999999) + 1.0e-99; Real x167 = X(17)*X(3); @@ -483,11 +483,11 @@ void rhs_specie(const burn_t& state, Real x172 = -x169*x170 - x170*x171; - Real x173 = 7.5399999999999998e-10*std::exp((-0.45800000000000002)*std::log(std::abs(x1))) + 4.9899999999999997e-11*std::exp((0.40500000000000003)*std::log(std::abs(x1))) + 1.0e-99; + Real x173 = 7.5399999999999998e-10*std::pow(x1, -0.45800000000000002) + 4.9899999999999997e-11*std::pow(x1, 0.40500000000000003) + 1.0e-99; Real x174 = X(21)*X(4); - Real x175 = (4.0000000000000001e-10*std::exp((0.0066899999999999998)*std::log(std::abs(x1))) + 1.0799999999999999e-11*std::exp((0.51700000000000002)*std::log(std::abs(x1))))*std::exp(-227.0*x52) + 1.0e-99; + Real x175 = (4.0000000000000001e-10*std::pow(x1, 0.0066899999999999998) + 1.0799999999999999e-11*std::pow(x1, 0.51700000000000002))*std::exp(-227.0*x52) + 1.0e-99; Real x176 = X(22)*X(3); @@ -515,7 +515,7 @@ void rhs_specie(const burn_t& state, Real x185 = x179/((x137)*(x137)*(x137)); Real x186 = ((x1 > 30.0) ? ( - std::exp((-3.194*x139 + 1.786*x140*x141 - 0.2072*x185 - 18.199999999999999)*std::log(std::abs(10))) + std::pow(10, -3.194*x139 + 1.786*x140*x141 - 0.2072*x185 - 18.199999999999999) ) : ( 0 @@ -523,7 +523,7 @@ void rhs_specie(const burn_t& state, Real x187 = X(3)*X(4); - Real x188 = std::exp((-0.14999999999999999)*std::log(std::abs(x1))); + Real x188 = std::pow(x1, -0.14999999999999999); Real x189 = ((x1 <= 30.0) ? ( 3.4977396723747635e-20*x188 @@ -550,9 +550,9 @@ void rhs_specie(const burn_t& state, Real x198 = 2.0*X(10) + X(11) + X(12) + X(18) + 2.0*X(19) + 3.0*X(20) + 4.0*X(23) + X(24) + X(25) + 2.0*X(26) + 2.0*X(27) + 3.0*X(28) + X(3) + X(4) + X(5) + 2.0*X(8); - Real x199 = 1.0/(std::exp((1.3*x138*(x197 - 9.2103403719761836) - 137.42519902360013*x141*((0.10857362047581294*x197 - 1)*(0.10857362047581294*x197 - 1)) - 4.8449999999999998)*std::log(std::abs(10.0)))*x198 + 1.0); + Real x199 = 1.0/(std::pow(10.0, 1.3*x138*(x197 - 9.2103403719761836) - 137.42519902360013*x141*((0.10857362047581294*x197 - 1)*(0.10857362047581294*x197 - 1)) - 4.8449999999999998)*x198 + 1.0); - Real x200 = std::exp((x199)*std::log(std::abs(1.1800000000000001e-10*std::exp(-69500.0*x52))))*std::exp((1.0 - x199)*std::log(std::abs(8.1250000000000003e-8*std::exp((-0.5)*std::log(std::abs(T)))*(1.0 - std::exp(-6000.0*x52))*std::exp(-52000.0*x52)))); + Real x200 = std::pow(1.1800000000000001e-10*std::exp(-69500.0*x52), x199)*std::pow(8.1250000000000003e-8*std::pow(T, -0.5)*(1.0 - std::exp(-6000.0*x52))*std::exp(-52000.0*x52), 1.0 - x199); Real x201 = 2*x200 + 2.0e-99; @@ -570,7 +570,7 @@ void rhs_specie(const burn_t& state, Real x208 = 2.9375070000000001*x206 + 0.23588480000000001*x207 + 0.75022860000000002; - Real x209 = std::exp((21.360939999999999*x138*x204 + 42.707410000000003*x138*x23 - 2.0273650000000001*x142 - x203 - 21467.790000000001*x52 - 1657.4099999999999*x52/(std::exp((x208)*std::log(std::abs(std::exp((-x205 - 8.1313220000000008)*std::log(std::abs(10.0)))*x198))) + 1.0) - 142.7664 - (21.360939999999999*x138*x204 + 11.28215*x138*std::log(14254.549999999999*x52 + 1.0) + 70.138370000000009*x139 - 4.7035149999999994*x142 - x203 - 203.11568)/(std::exp((x208)*std::log(std::abs(std::exp((-x205 - 9.3055640000000004)*std::log(std::abs(10.0)))*x198))) + 1.0))*std::log(std::abs(10.0))); + Real x209 = std::pow(10.0, 21.360939999999999*x138*x204 + 42.707410000000003*x138*x23 - 2.0273650000000001*x142 - x203 - 21467.790000000001*x52 - 1657.4099999999999*x52/(std::pow(std::pow(10.0, -x205 - 8.1313220000000008)*x198, x208) + 1.0) - 142.7664 - (21.360939999999999*x138*x204 + 11.28215*x138*std::log(14254.549999999999*x52 + 1.0) + 70.138370000000009*x139 - 4.7035149999999994*x142 - x203 - 203.11568)/(std::pow(std::pow(10.0, -x205 - 9.3055640000000004)*x198, x208) + 1.0)); Real x210 = 4.6331670000000003*x185; @@ -580,7 +580,7 @@ void rhs_specie(const burn_t& state, Real x213 = -2.0563129999999998*x206 + 0.58640729999999996*x207 + 0.82274429999999998; - Real x214 = std::exp((69.700860000000006*x138*x211 - 68.422430000000006*x139 + 43.20243*x140*x141 - x210 - 23705.700000000001*x52 - 2080.4099999999999*x52/(std::exp((x213)*std::log(std::abs(std::exp((-x212 - 13.656822)*std::log(std::abs(10.0)))*x198))) + 1.0) - 178.4239 - (69.700860000000006*x138*x211 + 19.734269999999999*x138*std::log(16780.950000000001*x52 + 1.0) - 14.509090000000008*x139 + 37.886913*x142 - x210 - 307.31920000000002)/(std::exp((x213)*std::log(std::abs(std::exp((-x212 - 14.82123)*std::log(std::abs(10.0)))*x198))) + 1.0))*std::log(std::abs(10.0))); + Real x214 = std::pow(10.0, 69.700860000000006*x138*x211 - 68.422430000000006*x139 + 43.20243*x140*x141 - x210 - 23705.700000000001*x52 - 2080.4099999999999*x52/(std::pow(std::pow(10.0, -x212 - 13.656822)*x198, x213) + 1.0) - 178.4239 - (69.700860000000006*x138*x211 + 19.734269999999999*x138*std::log(16780.950000000001*x52 + 1.0) - 14.509090000000008*x139 + 37.886913*x142 - x210 - 307.31920000000002)/(std::pow(std::pow(10.0, -x212 - 14.82123)*x198, x213) + 1.0)); Real x215 = X(10)*X(4); @@ -590,7 +590,7 @@ void rhs_specie(const burn_t& state, Real x218 = X(4)*x198*x5/((1.0 + 10000.0*std::exp(-600.0/(Tdust + 9.9999999999999993e-41)))*(7.9999999999999996e-6*x2 + x3 + 0.040000000000000001*std::sqrt(Tdust + x1))); - Real x219 = 1.1194685000500219e-20*std::exp((2.7400000000000002)*std::log(std::abs(x1)))*std::exp(-4740.0*x52) + 1.0e-99; + Real x219 = 1.1194685000500219e-20*std::pow(x1, 2.7400000000000002)*std::exp(-4740.0*x52) + 1.0e-99; Real x220 = X(10)*X(20); @@ -600,7 +600,7 @@ void rhs_specie(const burn_t& state, Real x223 = x219*x220 - x221*x222; - Real x224 = 1.9643512449766835e-11*std::exp((0.17000000000000001)*std::log(std::abs(x1)))*std::exp(-6400.0*x52) + 1.0e-99; + Real x224 = 1.9643512449766835e-11*std::pow(x1, 0.17000000000000001)*std::exp(-6400.0*x52) + 1.0e-99; Real x225 = X(10)*X(19); @@ -630,7 +630,7 @@ void rhs_specie(const burn_t& state, Real x238 = x235 - x236*x237; - Real x239 = 6.1739095063118665e-10*std::exp((0.40999999999999998)*std::log(std::abs(x1))) + 1.0e-99; + Real x239 = 6.1739095063118665e-10*std::pow(x1, 0.40999999999999998) + 1.0e-99; Real x240 = -x239; @@ -679,7 +679,7 @@ void rhs_specie(const burn_t& state, Real x257 = -x129 + x254*x255 + x255*x256; Real x258 = ((x1 <= 1167.4796423742259) ? ( - std::exp((5.8888600000000002*x139 + 7.1969200000000004*x142 + 2.2506900000000001*x185 - 56.473700000000001 - 2.1690299999999998*x180/((((x137)*(x137)))*(((x137)*(x137)))) + 0.31788699999999998*x181/((x137)*(x137)*(x137)*(x137)*(x137)))*std::log(std::abs(10))) + std::pow(10, 5.8888600000000002*x139 + 7.1969200000000004*x142 + 2.2506900000000001*x185 - 56.473700000000001 - 2.1690299999999998*x180/((((x137)*(x137)))*(((x137)*(x137)))) + 0.31788699999999998*x181/((x137)*(x137)*(x137)*(x137)*(x137))) ) : ( 0 @@ -705,7 +705,7 @@ void rhs_specie(const burn_t& state, Real x263 = x258*x259 + x259*x260 - x261*x262; - Real x264 = 3.5999999999999998e-16*std::exp((1.52)*std::log(std::abs(T)))*std::exp(-1740.0*x52) + 1.0e-99; + Real x264 = 3.5999999999999998e-16*std::pow(T, 1.52)*std::exp(-1740.0*x52) + 1.0e-99; Real x265 = X(10)*X(25); @@ -751,7 +751,7 @@ void rhs_specie(const burn_t& state, Real x278 = X(17)*X(25); Real x279 = ((x151) ? ( - 1.5557027645121538e-10*std::exp((-0.33900000000000002)*std::log(std::abs(x1)))*std::exp(-0.108*x52) + 1.5557027645121538e-10*std::pow(x1, -0.33900000000000002)*std::exp(-0.108*x52) ) : ( 0 @@ -781,7 +781,7 @@ void rhs_specie(const burn_t& state, Real x288 = X(25)*X(4); Real x289 = ((x1 <= 280.0) ? ( - 8.1009819130659227e-21*std::exp((2.7999999999999998)*std::log(std::abs(x1)))*std::exp(-1950.0*x52) + 8.1009819130659227e-21*std::pow(x1, 2.7999999999999998)*std::exp(-1950.0*x52) ) : ( 0 @@ -789,7 +789,7 @@ void rhs_specie(const burn_t& state, Real x290 = x285*x286 - x287*x288 - x288*x289 + x83*x84; - Real x291 = 8.648509591313017e-18*std::exp((-0.38)*std::log(std::abs(x1))) + 1.0e-99; + Real x291 = 8.648509591313017e-18*std::pow(x1, -0.38) + 1.0e-99; Real x292 = X(22)*X(4); @@ -881,7 +881,7 @@ void rhs_specie(const burn_t& state, Real x328 = X(17)*X(32); Real x329 = ((x1 > 1052.0) ? ( - 3.7991369449339675e-16*std::exp((1.54)*std::log(std::abs(x1)))*std::exp(613.0*x52) + 3.7991369449339675e-16*std::pow(x1, 1.54)*std::exp(613.0*x52) ) : ( 0 @@ -897,13 +897,13 @@ void rhs_specie(const burn_t& state, Real x334 = ((X(25))*(X(25))); - Real x335 = std::exp((1.1399999999999999)*std::log(std::abs(x1)))*std::exp(-50.0*x52); + Real x335 = std::pow(x1, 1.1399999999999999)*std::exp(-50.0*x52); Real x336 = x334*(2.4749459212614312e-15*x335 + 1.0e-99); Real x337 = ((X(22))*(X(22))); - Real x338 = std::exp((1.5800000000000001)*std::log(std::abs(x1))); + Real x338 = std::pow(x1, 1.5800000000000001); Real x339 = 1.0e-99 + 2.8e-11*std::exp(-26500.0*x52); @@ -1013,7 +1013,7 @@ Real rhs_eint(const burn_t& state, using namespace Rates; Real T = state.T; - Real x0 = 3.3470650363800003e-24*X(10) + 5.0186540981899997e-24*X(11) + 5.01956503638e-24*X(12) + 6.6902431600000005e-24*X(13) + 6.6911540981899994e-24*X(14) + 6.6920650363799998e-24*X(15) + 2.007528417094e-23*X(16) + 2.0076195109127999e-23*X(17) + 2.1749727627315999e-23*X(18) + 2.3423260145503998e-23*X(19) + 9.1093818800000008e-28*X(2) + 2.5096792663692001e-23*X(20) + 2.6767349207315999e-23*X(21) + 2.6768260145504001e-23*X(22) + 2.6770325181879997e-23*X(23) + 2.8440881725504001e-23*X(24) + 2.8441792663692003e-23*X(25) + 3.0114414243691998e-23*X(26) + 3.011532518188e-23*X(27) + 3.178794676188e-23*X(28) + 4.6843544316443998e-23*X(29) + 1.6726215800000001e-24*X(3) + 4.6844455254632e-23*X(30) + 5.3535609352819999e-23*X(31) + 5.3536520291008001e-23*X(32) + 7.3612715400136001e-23*X(33) + 1.6735325181900001e-24*X(4) + 1.6744434563800001e-24*X(5) + 3.3451215800000003e-24*X(6) + 3.3460325181899999e-24*X(7) + 3.3461540981899999e-24*X(8) + 3.3469434563800003e-24*X(9); + Real x0 = 3.3470650363800003e-24*X(10) + 5.0186540981899997e-24*X(11) + 5.01956503638e-24*X(12) + 6.6902431600000005e-24*X(13) + 6.6911540981899994e-24*X(14) + 6.6920650363799998e-24*X(15) + 2.007528417094e-23*X(16) + 2.0076195109127999e-23*X(17) + 2.1749727627315999e-23*X(18) + 2.3423260145503998e-23*X(19) + 9.1093818800000008e-28*X(2) + 2.5096792663692001e-23*X(20) + 2.6767349207315999e-23*X(21) + 2.6768260145504001e-23*X(22) + 2.6770325181879997e-23*X(23) + 2.8440881725504001e-23*X(24) + 2.8441792663692003e-23*X(25) + 3.0114414243691998e-23*X(26) + 3.011532518188e-23*X(27) + 3.178794676188e-23*X(28) + 4.6843544316443998e-23*X(29) + 1.6726215800000001e-24*X(3) + 4.6844455254632e-23*X(30) + 5.3535609352819999e-23*X(31) + 5.3536520291008001e-23*X(32) + 7.3612715400136001e-23*X(33) + 1.6735325181900001e-24*X(4) + 1.6744434563800001e-24*X(5) + 3.3451215800000003e-24*X(6) + 3.3460325181899999e-24*X(7) + 3.3461540981899999e-24*X(8) + 3.3469434563800003e-24*X(9); Real x1 = 1.0/x0; @@ -1073,18 +1073,18 @@ Real rhs_eint(const burn_t& state, 0.63095734448019325 ) : ((x9) ? ( - std::exp((-0.20000000000000001)*std::log(std::abs(T))) + std::pow(T, -0.20000000000000001) ) : ( - 0.81802668403695478*std::exp((-0.20000000000000001)*std::log(std::abs(x10))) + 0.81802668403695478*std::pow(x10, -0.20000000000000001) )))/(6.3095734448019361e-5*((x8) ? ( 5.011872336272722 ) : ((x9) ? ( - std::exp((0.69999999999999996)*std::log(std::abs(T))) + std::pow(T, 0.69999999999999996) ) : ( - 2.0198255098866533*std::exp((0.69999999999999996)*std::log(std::abs(x10))) + 2.0198255098866533*std::pow(x10, 0.69999999999999996) ))) + 1.0); Real x21 = x11*x17; @@ -1093,10 +1093,10 @@ Real rhs_eint(const burn_t& state, 0.67810976749343443 ) : ((x9) ? ( - std::exp((-0.16869999999999999)*std::log(std::abs(T))) + std::pow(T, -0.16869999999999999) ) : ( - 0.8441494880969781*std::exp((-0.16869999999999999)*std::log(std::abs(x10))) + 0.8441494880969781*std::pow(x10, -0.16869999999999999) ))); Real x23 = 1.0/x5; @@ -1111,7 +1111,7 @@ Real rhs_eint(const burn_t& state, Real x28 = 1.0/(1000000.0*x23*x26/(1.3999999999999999*X(10)*x26*std::exp(-12000.0/(x4 + 1200.0)) + 1.6000000000000001*X(4)*x26*std::exp(-160000.0/x27)) + 1.0); - Real x29 = std::exp((-0.25)*std::log(std::abs(x4))); + Real x29 = std::pow(x4, -0.25); Real x30 = std::log(T); @@ -1121,7 +1121,7 @@ Real rhs_eint(const burn_t& state, Real x33 = std::exp((-2)*std::log(std::abs(x31))); - Real x34 = 1.0/(std::exp((1.3*x32*(x30 - 9.2103403719761836) - 137.42519902360013*x33*((0.10857362047581294*x30 - 1)*(0.10857362047581294*x30 - 1)) - 4.8449999999999998)*std::log(std::abs(10.0)))*x25 + 1.0); + Real x34 = 1.0/(std::pow(10.0, 1.3*x32*(x30 - 9.2103403719761836) - 137.42519902360013*x33*((0.10857362047581294*x30 - 1)*(0.10857362047581294*x30 - 1)) - 4.8449999999999998)*x25 + 1.0); Real x35 = std::sqrt(M_PI); @@ -1192,12 +1192,12 @@ Real rhs_eint(const burn_t& state, Real x68 = x0 >= 9.9999999999999998e-13; Real x69 = ((x68) ? ( - std::exp((x67*std::log(((x66) ? ( + std::pow(10.0, x67*std::log(((x66) ? ( 0.5 ) : ( x0 - ))) + 2.1498900000000001)*std::log(std::abs(10.0))) + ))) + 2.1498900000000001) ) : ( 0.0 @@ -1211,7 +1211,7 @@ Real rhs_eint(const burn_t& state, Real x73 = x0 <= 9.9999999999999993e-41; - Real x74 = std::exp((2.1498900000000001 - 0.69317629274152892*x32)*std::log(std::abs(10.0))); + Real x74 = std::pow(10.0, 2.1498900000000001 - 0.69317629274152892*x32); Real x75 = x0*x35; @@ -1219,7 +1219,7 @@ Real rhs_eint(const burn_t& state, Real x77 = x74*x76; - Real x78 = std::exp((x67*std::log(x0) + 2.1498900000000001)*std::log(std::abs(10.0))); + Real x78 = std::pow(10.0, x67*std::log(x0) + 2.1498900000000001); Real x79 = x76*x78; @@ -1239,7 +1239,7 @@ Real rhs_eint(const burn_t& state, Real x87 = std::sqrt(x71*x86); - Real x88 = 4.985670872372847e-33*std::exp((3.7599999999999998)*std::log(std::abs(x4)))*std::exp(-2197000.0/((x4)*(x4)*(x4)))/(6.0142468035272636e-8*std::exp((2.1000000000000001)*std::log(std::abs(x4))) + 1.0) + 1.6e-18*std::exp(-11700.0*x6) + 6.7e-19*std::exp(-5860.0*x6) + 3.0e-24*std::exp(-510.0*x6); + Real x88 = 4.985670872372847e-33*std::pow(x4, 3.7599999999999998)*std::exp(-2197000.0/((x4)*(x4)*(x4)))/(6.0142468035272636e-8*std::pow(x4, 2.1000000000000001) + 1.0) + 1.6e-18*std::exp(-11700.0*x6) + 6.7e-19*std::exp(-5860.0*x6) + 3.0e-24*std::exp(-510.0*x6); Real x89 = x4 < 2000.0; @@ -1263,7 +1263,7 @@ Real rhs_eint(const burn_t& state, Real x99 = std::exp((8)*std::log(std::abs(x91))); - Real x100 = std::exp((5.0194035000000001*x32*x90 + 5627.2167698544854*x37*x95 + 86051.290034608537*x46*x96 + 9415777.8988952208*x47*x99 - 75.100986441619156*x92 - 1554.3387057364687*x94 - 428804.85473346239*x97 - 1662263.0320406025*x98 - 20.584225)*std::log(std::abs(10))); + Real x100 = std::pow(10, 5.0194035000000001*x32*x90 + 5627.2167698544854*x37*x95 + 86051.290034608537*x46*x96 + 9415777.8988952208*x47*x99 - 75.100986441619156*x92 - 1554.3387057364687*x94 - 428804.85473346239*x97 - 1662263.0320406025*x98 - 20.584225); Real x101 = x4 <= 10000.0; @@ -1279,53 +1279,53 @@ Real rhs_eint(const burn_t& state, x102 ))); - Real x104 = std::exp((25.0*x32)*std::log(std::abs(x4))); + Real x104 = std::pow(x4, 25.0*x32); - Real x105 = std::exp((-200.0 + 20000.0/((10.0 + 2.3538526683701997e+17/x104)*(1.6889118802245084e-48*x104 + 10.0)))*std::log(std::abs(10.0))); + Real x105 = std::pow(10.0, -200.0 + 20000.0/((10.0 + 2.3538526683701997e+17/x104)*(1.6889118802245084e-48*x104 + 10.0))); Real x106 = x37*x95; Real x107 = x46*x96; - Real x108 = std::exp((-339.5619991617852*x106 - 529.07725573213918*x107 + 2.0943374000000001*x32*x90 + 144.02112655888752*x39*x93 - 36.814414747418546*x92 - 23.962112000000001)*std::log(std::abs(10)))*X(10)*x105; + Real x108 = std::pow(10, -339.5619991617852*x106 - 529.07725573213918*x107 + 2.0943374000000001*x32*x90 + 144.02112655888752*x39*x93 - 36.814414747418546*x92 - 23.962112000000001)*X(10)*x105; Real x109 = x32*x90; - Real x110 = std::exp((-377.88183430702219*x106 + 3018.4974183098116*x107 + 2.1892372*x109 - 38.89917505778142*x92 + 95.70878894783884*x94 - 23.689236999999999)*std::log(std::abs(10)))*X(15); + Real x110 = std::pow(10, -377.88183430702219*x106 + 3018.4974183098116*x107 + 2.1892372*x109 - 38.89917505778142*x92 + 95.70878894783884*x94 - 23.689236999999999)*X(15); Real x111 = x4 > 10.0; Real x112 = x101 && x111; - Real x113 = std::exp((16.666666666666664*x32)*std::log(std::abs(x4))); + Real x113 = std::pow(x4, 16.666666666666664*x32); - Real x114 = std::exp((-200.0 + 20000.0/((10.0 + 785.77199422741614/x113)*(5.0592917094448065e-34*x113 + 10.0)))*std::log(std::abs(10.0))); + Real x114 = std::pow(10.0, -200.0 + 20000.0/((10.0 + 785.77199422741614/x113)*(5.0592917094448065e-34*x113 + 10.0))); Real x115 = 1.002560385050777e-22*X(15)*x114; - Real x116 = std::exp((-1161.2797752309887*x106 + 5059.6285287169567*x107 + 1.5714710999999999*x109 + 0.73442154540113413*x92 - 77.855706084264682*x94 - 22.089523)*std::log(std::abs(10)))*X(3); + Real x116 = std::pow(10, -1161.2797752309887*x106 + 5059.6285287169567*x107 + 1.5714710999999999*x109 + 0.73442154540113413*x92 - 77.855706084264682*x94 - 22.089523)*X(3); Real x117 = 1.1825091393820599e-21*X(3)*x114; - Real x118 = std::exp((45902.322591745004*x106 + 60522.293708798054*x107 + 37.383713*x109 + 2774.5177117396752*x92 + 16037.924047681272*x94 - 16.818342000000001)*std::log(std::abs(10)))*X(4); + Real x118 = std::pow(10, 45902.322591745004*x106 + 60522.293708798054*x107 + 37.383713*x109 + 2774.5177117396752*x92 + 16037.924047681272*x94 - 16.818342000000001)*X(4); Real x119 = x4 <= 100.0; - Real x120 = std::exp((-48562.751069188118*x106 - 66875.646562351845*x107 + 3.5692468000000002*x32*x90 - 540.77102118284597*x92 - 9179.8864335208946*x94 - 24.311209000000002)*std::log(std::abs(10)))*X(4); + Real x120 = std::pow(10, -48562.751069188118*x106 - 66875.646562351845*x107 + 3.5692468000000002*x32*x90 - 540.77102118284597*x92 - 9179.8864335208946*x94 - 24.311209000000002)*X(4); Real x121 = x4 <= 1000.0; - Real x122 = std::exp((-12547.661945180447*x106 + 24439.250555499191*x107 + 4.6450521*x109 - 177.55453097873294*x92 + 1956.911370108365*x94 - 24.311209000000002)*std::log(std::abs(10)))*X(4); + Real x122 = std::pow(10, -12547.661945180447*x106 + 24439.250555499191*x107 + 4.6450521*x109 - 177.55453097873294*x92 + 1956.911370108365*x94 - 24.311209000000002)*X(4); Real x123 = x4 <= 6000.0; - Real x124 = std::exp((17.997580222853362*x32)*std::log(std::abs(x4))); + Real x124 = std::pow(x4, 17.997580222853362*x32); - Real x125 = 1.8623144679125181e-22*std::exp((-200.0 + 20000.0/((10.0 + 2973.7534532281375/x124)*(1.3368457736780898e-34*x124 + 10.0)))*std::log(std::abs(10.0)))*X(4); + Real x125 = 1.8623144679125181e-22*std::pow(10.0, -200.0 + 20000.0/((10.0 + 2973.7534532281375/x124)*(1.3368457736780898e-34*x124 + 10.0)))*X(4); Real x126 = x47*x99; - Real x127 = std::exp((1672890.7229183144*x106 + 15471651.937466398*x107 + 16.815729999999999*x109 + 366063607.58415633*x126 + 4616.3011562659685*x92 + 113122.17137872758*x94 + 87115306.05744876*x97 + 273295393.17143697*x98 - 21.928795999999998)*std::log(std::abs(10)))*X(2); + Real x127 = std::pow(10, 1672890.7229183144*x106 + 15471651.937466398*x107 + 16.815729999999999*x109 + 366063607.58415633*x126 + 4616.3011562659685*x92 + 113122.17137872758*x94 + 87115306.05744876*x97 + 273295393.17143697*x98 - 21.928795999999998)*X(2); Real x128 = x4 <= 500.0; @@ -1333,7 +1333,7 @@ Real rhs_eint(const burn_t& state, Real x130 = x128 && x129; - Real x131 = std::exp((-10763.919849753534*x106 - 138531.11016116844*x107 + 1.6802758*x109 - 33025002.640084207*x126 + 44.525106942242758*x92 + 1331.8748828877385*x94 + 968783.44101153011*x97 + 4831859.3594864924*x98 - 22.921188999999998)*std::log(std::abs(10)))*X(2); + Real x131 = std::pow(10, -10763.919849753534*x106 - 138531.11016116844*x107 + 1.6802758*x109 - 33025002.640084207*x126 + 44.525106942242758*x92 + 1331.8748828877385*x94 + 968783.44101153011*x97 + 4831859.3594864924*x98 - 22.921188999999998)*X(2); Real x132 = x4 > 500.0; @@ -1846,7 +1846,7 @@ Real rhs_eint(const burn_t& state, Real x372 = std::max(0.0, X(4)); - Real x373 = x372*std::exp((0.070000000000000007)*std::log(std::abs(x4))); + Real x373 = x372*std::pow(x4, 0.070000000000000007); Real x374 = 5.6000000000000006e-6*x369*x371 + 1.1590975361199841e-9*x369*x373; @@ -1867,15 +1867,15 @@ Real rhs_eint(const burn_t& state, Real x382 = std::max(0.0, X(3)); Real x383 = x382*((x4 <= 2090.0) ? ( - 2.03e-11*std::exp((0.56000000000000005)*std::log(std::abs(x4))) + 2.03e-11*std::pow(x4, 0.56000000000000005) ) : ( - 3.43e-10*std::exp((0.19)*std::log(std::abs(x4))) + 3.43e-10*std::pow(x4, 0.19) )); - Real x384 = x372*std::exp((0.44)*std::log(std::abs(x4))); + Real x384 = x372*std::pow(x4, 0.44); - Real x385 = x370*std::exp((0.92600000000000005)*std::log(std::abs(x4))); + Real x385 = x370*std::pow(x4, 0.92600000000000005); Real x386 = x383 + 1.4500824124120476e-11*x384 + 1.08e-14*x385; @@ -1883,24 +1883,24 @@ Real rhs_eint(const burn_t& state, Real x388 = std::exp(-330.0*x6); - Real x389 = x370*std::exp((-0.025999999999999999)*std::log(std::abs(x4))); + Real x389 = x370*std::pow(x4, -0.025999999999999999); Real x390 = 9.7199999999999998e-11*x388*x389; - Real x391 = std::exp((0.80000000000000004)*std::log(std::abs(x4))); + Real x391 = std::pow(x4, 0.80000000000000004); Real x392 = x372*x391; Real x393 = 2.160222331098239e-13*x388*x392; Real x394 = x382*((x4 <= 511.0) ? ( - 6.1000000000000003e-13*std::exp((1.1000000000000001)*std::log(std::abs(x4))) + 6.1000000000000003e-13*std::pow(x4, 1.1000000000000001) ) : ((x4 <= 7510.0) ? ( - 2.1199999999999999e-12*std::exp((0.90000000000000002)*std::log(std::abs(x4))) + 2.1199999999999999e-12*std::pow(x4, 0.90000000000000002) ) : ( - 4.49e-10*std::exp((0.29999999999999999)*std::log(std::abs(x4))) + 4.49e-10*std::pow(x4, 0.29999999999999999) ))); Real x395 = 0.20000000000000001*x388*x394; @@ -1909,18 +1909,18 @@ Real rhs_eint(const burn_t& state, Real x397 = std::exp(-230.0*x6); - Real x398 = x370*std::exp((-0.074999999999999997)*std::log(std::abs(x4))); + Real x398 = x370*std::pow(x4, -0.074999999999999997); - Real x399 = x372*std::exp((0.67000000000000004)*std::log(std::abs(x4))); + Real x399 = x372*std::pow(x4, 0.67000000000000004); Real x400 = x382*((x4 <= 194.0) ? ( - 6.3800000000000002e-11*std::exp((0.40000000000000002)*std::log(std::abs(x4))) + 6.3800000000000002e-11*std::pow(x4, 0.40000000000000002) ) : ((x4 <= 3686.0) ? ( 7.7500000000000007e-12*x391 ) : ( - 2.6500000000000002e-10*std::exp((0.37)*std::log(std::abs(x4))) + 2.6500000000000002e-10*std::pow(x4, 0.37) ))); Real x401 = 3.0719999999999998e-10*x397*x398 + 2.5231268066741097e-12*x397*x399 + 0.59999999999999998*x397*x400; @@ -1935,7 +1935,7 @@ Real rhs_eint(const burn_t& state, Real x406 = X(22)/(x387*(-x390 - x393 - x395 + x404) + x401*(-x402 - x404) + x405*(x396 + x402)); - Real x407 = x372*std::exp((0.26000000000000001)*std::log(std::abs(x4))); + Real x407 = x372*std::pow(x4, 0.26000000000000001); Real x408 = std::max(0.0, 0.75*X(10)); @@ -1959,10 +1959,10 @@ Real rhs_eint(const burn_t& state, Real x415 = x4 <= 5000.0; Real x416 = x382*((x415) ? ( - std::exp((0.69999999999999996)*std::log(std::abs(x4)))*(1.7e-18*x27 - 2.2000000000000001e-14*x4 + 1.0e-10) + std::pow(x4, 0.69999999999999996)*(1.7e-18*x27 - 2.2000000000000001e-14*x4 + 1.0e-10) ) : ( - 9.1999999999999997e-9*std::exp((0.053499999999999999)*std::log(std::abs(x4))) + 9.1999999999999997e-9*std::pow(x4, 0.053499999999999999) )); Real x417 = x408*(2.8999999999999998e-10 - 1.8999999999999999e-10*x409) + x410*(1.8e-10*x411 - 2.5999999999999998e-10*x412 + 2.7e-10) + x414 + x416; @@ -1987,7 +1987,7 @@ Real rhs_eint(const burn_t& state, x4*(3.9000000000000001e-20*x27 - 5.9999999999999999e-16*x4 + 3.1000000000000001e-12) ) : ( - 2.2999999999999999e-9*std::exp((0.096500000000000002)*std::log(std::abs(x4))) + 2.2999999999999999e-9*std::pow(x4, 0.096500000000000002) )); Real x427 = 5.0*x419; @@ -2007,7 +2007,7 @@ Real rhs_eint(const burn_t& state, Real x432 = std::exp(-24.0*x6); - Real x433 = x372*std::exp((0.14000000000000001)*std::log(std::abs(x4))); + Real x433 = x372*std::pow(x4, 0.14000000000000001); Real x434 = std::exp(-0.01582278481012658*x4); @@ -2018,10 +2018,10 @@ Real rhs_eint(const burn_t& state, Real x437 = std::exp(-0.0045808520384791572*x4); Real x438 = x382*((x415) ? ( - std::exp((0.45000000000000001)*std::log(std::abs(x4)))*(1.8999999999999999e-18*x27 - 1.7999999999999999e-14*x4 + 9.6000000000000005e-11) + std::pow(x4, 0.45000000000000001)*(1.8999999999999999e-18*x27 - 1.7999999999999999e-14*x4 + 9.6000000000000005e-11) ) : ( - 8.9000000000000003e-10*std::exp((0.11700000000000001)*std::log(std::abs(x4))) + 8.9000000000000003e-10*std::pow(x4, 0.11700000000000001) )); Real x439 = 3.0*x432; @@ -2187,168 +2187,164 @@ Real rhs_eint(const burn_t& state, Real x517 = x450*x511; - Real x518 = x455*x456; + Real x518 = x450*x513; - Real x519 = x455*x458; + Real x519 = x455*x456; - Real x520 = x455*x460; + Real x520 = x455*x458; - Real x521 = x452*x508; + Real x521 = x455*x460; - Real x522 = x452*x509; + Real x522 = x452*x508; - Real x523 = x452*x510; + Real x523 = x452*x509; - Real x524 = x452*x512; + Real x524 = x452*x510; - Real x525 = x452*x513; + Real x525 = x452*x512; Real x526 = x455*x464; - Real x527 = x455*x467; + Real x527 = x455*x466; - Real x528 = x455*x469; + Real x528 = x455*x467; - Real x529 = x456*x508; + Real x529 = x456*x509; - Real x530 = x456*x509; + Real x530 = x456*x511; - Real x531 = x456*x511; + Real x531 = x457*x508; - Real x532 = x457*x508; + Real x532 = x457*x510; - Real x533 = x457*x509; + Real x533 = x458*x509; Real x534 = x459*x508; Real x535 = x464*x510; - Real x536 = x465*x508; + Real x536 = x466*x508; - Real x537 = x465*x511; + Real x537 = x466*x510; - Real x538 = x466*x508; + Real x538 = x467*x508; - Real x539 = x466*x510; + Real x539 = std::min(18.99981, x32*(x454 + std::log(std::max(9.9999999999999993e-41, X(27))))); - Real x540 = x467*x508; + Real x540 = ((x539)*(x539)); - Real x541 = std::min(18.99981, x32*(x454 + std::log(std::max(9.9999999999999993e-41, X(27))))); + Real x541 = ((x539)*(x539)*(x539)); - Real x542 = ((x541)*(x541)); + Real x542 = ((((x539)*(x539)))*(((x539)*(x539)))); - Real x543 = ((x541)*(x541)*(x541)); + Real x543 = ((x539)*(x539)*(x539)*(x539)*(x539)); - Real x544 = ((((x541)*(x541)))*(((x541)*(x541)))); + Real x544 = std::exp((6)*std::log(std::abs(x539))); - Real x545 = ((x541)*(x541)*(x541)*(x541)*(x541)); + Real x545 = ((x539)*(x539)*(x539)*(x539)*(x539)*(x539)*(x539)); - Real x546 = std::exp((6)*std::log(std::abs(x541))); + Real x546 = std::exp((8)*std::log(std::abs(x539))); - Real x547 = ((x541)*(x541)*(x541)*(x541)*(x541)*(x541)*(x541)); + Real x547 = std::min(13.99986, x449); - Real x548 = std::exp((8)*std::log(std::abs(x541))); + Real x548 = std::min(3.9999600000000002, x61); - Real x549 = std::min(13.99986, x449); + Real x549 = std::min(24.999750000000002, x32*(x454 + std::log(std::max(9.9999999999999993e-41, X(30))))); - Real x550 = std::min(3.9999600000000002, x61); + Real x550 = ((x547)*(x547)); - Real x551 = std::min(24.999750000000002, x32*(x454 + std::log(std::max(9.9999999999999993e-41, X(30))))); + Real x551 = ((x547)*(x547)*(x547)); - Real x552 = ((x549)*(x549)); + Real x552 = ((((x547)*(x547)))*(((x547)*(x547)))); - Real x553 = ((x549)*(x549)*(x549)); + Real x553 = ((x547)*(x547)*(x547)*(x547)*(x547)); - Real x554 = ((((x549)*(x549)))*(((x549)*(x549)))); + Real x554 = std::exp((6)*std::log(std::abs(x547))); - Real x555 = ((x549)*(x549)*(x549)*(x549)*(x549)); + Real x555 = ((x547)*(x547)*(x547)*(x547)*(x547)*(x547)*(x547)); - Real x556 = std::exp((6)*std::log(std::abs(x549))); + Real x556 = std::exp((8)*std::log(std::abs(x547))); - Real x557 = ((x549)*(x549)*(x549)*(x549)*(x549)*(x549)*(x549)); + Real x557 = ((x547)*(x547)*(x547)*(x547)*(x547)*(x547)*(x547)*(x547)*(x547)); - Real x558 = std::exp((8)*std::log(std::abs(x549))); + Real x558 = ((x548)*(x548)); - Real x559 = ((x549)*(x549)*(x549)*(x549)*(x549)*(x549)*(x549)*(x549)*(x549)); + Real x559 = ((x548)*(x548)*(x548)); - Real x560 = ((x550)*(x550)); + Real x560 = ((((x548)*(x548)))*(((x548)*(x548)))); - Real x561 = ((x550)*(x550)*(x550)); + Real x561 = ((x548)*(x548)*(x548)*(x548)*(x548)); - Real x562 = ((((x550)*(x550)))*(((x550)*(x550)))); + Real x562 = std::exp((6)*std::log(std::abs(x548))); - Real x563 = ((x550)*(x550)*(x550)*(x550)*(x550)); + Real x563 = ((x548)*(x548)*(x548)*(x548)*(x548)*(x548)*(x548)); - Real x564 = std::exp((6)*std::log(std::abs(x550))); + Real x564 = std::exp((8)*std::log(std::abs(x548))); - Real x565 = ((x550)*(x550)*(x550)*(x550)*(x550)*(x550)*(x550)); + Real x565 = ((x548)*(x548)*(x548)*(x548)*(x548)*(x548)*(x548)*(x548)*(x548)); - Real x566 = std::exp((8)*std::log(std::abs(x550))); + Real x566 = x553*x558; - Real x567 = ((x550)*(x550)*(x550)*(x550)*(x550)*(x550)*(x550)*(x550)*(x550)); + Real x567 = ((x549)*(x549)); - Real x568 = x555*x560; + Real x568 = ((x549)*(x549)*(x549)); - Real x569 = ((x551)*(x551)); + Real x569 = ((((x549)*(x549)))*(((x549)*(x549)))); - Real x570 = ((x551)*(x551)*(x551)); + Real x570 = ((x549)*(x549)*(x549)*(x549)*(x549)); - Real x571 = ((((x551)*(x551)))*(((x551)*(x551)))); + Real x571 = std::exp((6)*std::log(std::abs(x549))); - Real x572 = ((x551)*(x551)*(x551)*(x551)*(x551)); + Real x572 = ((x549)*(x549)*(x549)*(x549)*(x549)*(x549)*(x549)); - Real x573 = std::exp((6)*std::log(std::abs(x551))); + Real x573 = std::exp((8)*std::log(std::abs(x549))); - Real x574 = ((x551)*(x551)*(x551)*(x551)*(x551)*(x551)*(x551)); + Real x574 = ((x549)*(x549)*(x549)*(x549)*(x549)*(x549)*(x549)*(x549)*(x549)); - Real x575 = std::exp((8)*std::log(std::abs(x551))); + Real x575 = x549*x550; - Real x576 = ((x551)*(x551)*(x551)*(x551)*(x551)*(x551)*(x551)*(x551)*(x551)); + Real x576 = x549*x551; - Real x577 = x551*x552; + Real x577 = x550*x567; - Real x578 = x551*x553; + Real x578 = x558*x571; - Real x579 = x552*x569; + Real x579 = x561*x569; - Real x580 = x560*x573; + Real x580 = x548*x549; - Real x581 = x563*x571; + Real x581 = x547*x568; - Real x582 = x550*x551; + Real x582 = x547*x571; - Real x583 = x549*x570; + Real x583 = x547*x572; - Real x584 = x549*x573; + Real x584 = x548*x554; - Real x585 = x549*x574; + Real x585 = x547*x549; - Real x586 = x550*x556; + Real x586 = x550*x568; - Real x587 = x549*x551; + Real x587 = x551*x567; - Real x588 = x552*x570; + Real x588 = x551*x568; - Real x589 = x553*x569; + Real x589 = x553*x568; - Real x590 = x553*x570; + Real x590 = x547*x567; - Real x591 = x555*x570; + Real x591 = x559*x568; - Real x592 = x549*x569; + Real x592 = x549*x559; - Real x593 = x561*x570; + Real x593 = x552*x558; - Real x594 = x551*x561; + Real x594 = x550*x559; - Real x595 = x554*x560; + Real x595 = x559*x567; - Real x596 = x552*x561; - Real x597 = x561*x569; - - - return (x1*(-7.1777505408000004e-12*((X(10))*(X(10)))*(std::exp((x34)*std::log(std::abs(1.1800000000000001e-10*std::exp(-69500.0*x6))))*std::exp((1.0 - x34)*std::log(std::abs(8.1250000000000003e-8*std::exp((-0.5)*std::log(std::abs(T)))*(1.0 - std::exp(-6000.0*x6))*std::exp(-52000.0*x6)))) + 1.0e-99) + 7.1777505408000004e-12*X(10)*((X(4))*(X(4)))*x28*(2.5000000000000002e-32*x23 + 7.5000000000000001e-33*x29 + 1.0e-99) - 1.3854129535261706e-25*X(13)*x20 - 9.3799999999999993e-22*X(15)*x16*x18*std::exp(-285335.40000000002*x12) - 2.1299999999999999e-27*X(2)*x5*(4.0*X(13) + x2) - X(2)*x7*(4.3799999999999999e-10*std::exp((0.34999999999999998)*std::log(std::abs(x4)))*std::exp(-102000.0*x6) + 1.0e-99) - 5.6500000000000001e-36*X(2)*(-x3 + x4)*((((z + 1.0)*(z + 1.0)))*(((z + 1.0)*(z + 1.0)))) - 3.4635323838154264e-26*X(3)*x20 + 7.1777505408000004e-12*((X(4))*(X(4))*(X(4)))*x28*(2.0000000000000002e-31*x23 + 6.0000000000000001e-32*x29 + 1.0e-99) + 5.6556829037999995e-12*X(4)*X(5)*x28*((1.3500000000000001e-9*std::exp((0.098492999999999997)*std::log(std::abs(x4))) + 4.4350199999999998e-10*std::exp((0.55610000000000004)*std::log(std::abs(x4))) + 3.7408500000000004e-16*std::exp((2.1825999999999999)*std::log(std::abs(x4))))/(0.0061910000000000003*std::exp((1.0461)*std::log(std::abs(x4))) + 8.9711999999999997e-11*std::exp((3.0424000000000002)*std::log(std::abs(x4))) + 3.2575999999999999e-14*std::exp((3.7740999999999998)*std::log(std::abs(x4))) + 1.0) + 1.0e-99) + 1.75918975308e-21*X(4)*X(8)*x28 + 4.8065293799999994e-30*X(4)*user_dust2gas_ratio*x25*std::sqrt(x4)*(4.2000000000000002*x28 + 0.20000000000000001)/((1.0 + 10000.0*std::exp(-600.0/(Tdust + 9.9999999999999993e-41)))*(7.9999999999999996e-6*x27 + 0.002*x4 + 0.040000000000000001*std::sqrt(Tdust + x4) + 1.0)) - X(4)*x7*(std::exp((42.707410000000003*x32*x50 + 21.360939999999999*x32*x56 - x53 - 2.0273650000000001*x55 - 21467.790000000001*x6 - 1657.4099999999999*x6/(std::exp((x60)*std::log(std::abs(std::exp((-x57 - 8.1313220000000008)*std::log(std::abs(10.0)))*x25))) + 1.0) - 142.7664 - (21.360939999999999*x32*x56 + 11.28215*x32*std::log(14254.549999999999*x6 + 1.0) - x53 - 4.7035149999999994*x55 + 70.138370000000009*x61 - 203.11568)/(std::exp((x60)*std::log(std::abs(std::exp((-x57 - 9.3055640000000004)*std::log(std::abs(10.0)))*x25))) + 1.0))*std::log(std::abs(10.0))) + std::exp((69.700860000000006*x32*x63 + 43.20243*x33*x54 - 23705.700000000001*x6 - 2080.4099999999999*x6/(std::exp((x65)*std::log(std::abs(std::exp((-x64 - 13.656822)*std::log(std::abs(10.0)))*x25))) + 1.0) - 68.422430000000006*x61 - x62 - 178.4239 - (69.700860000000006*x32*x63 + 19.734269999999999*x32*std::log(16780.950000000001*x6 + 1.0) + 37.886913*x55 - 14.509090000000008*x61 - x62 - 307.31920000000002)/(std::exp((x65)*std::log(std::abs(std::exp((-x64 - 14.82123)*std::log(std::abs(10.0)))*x25))) + 1.0))*std::log(std::abs(10.0))) + 3.0000000000000001e-99) - 0.00022681492*x0*((((x4)*(x4)))*(((x4)*(x4))))*x69*((((x66 && x68 && x72 && x73 && x9) ? ( + return (x1*(-7.1777505408000004e-12*((X(10))*(X(10)))*(std::pow(1.1800000000000001e-10*std::exp(-69500.0*x6), x34)*std::pow(8.1250000000000003e-8*std::pow(T, -0.5)*(1.0 - std::exp(-6000.0*x6))*std::exp(-52000.0*x6), 1.0 - x34) + 1.0e-99) + 7.1777505408000004e-12*X(10)*((X(4))*(X(4)))*x28*(2.5000000000000002e-32*x23 + 7.5000000000000001e-33*x29 + 1.0e-99) - 1.3854129535261706e-25*X(13)*x20 - 9.3799999999999993e-22*X(15)*x16*x18*std::exp(-285335.40000000002*x12) - 2.1299999999999999e-27*X(2)*x5*(4.0*X(13) + x2) - X(2)*x7*(4.3799999999999999e-10*std::pow(x4, 0.34999999999999998)*std::exp(-102000.0*x6) + 1.0e-99) - 5.6500000000000001e-36*X(2)*(-x3 + x4)*((((z + 1.0)*(z + 1.0)))*(((z + 1.0)*(z + 1.0)))) - 3.4635323838154264e-26*X(3)*x20 + 7.1777505408000004e-12*((X(4))*(X(4))*(X(4)))*x28*(2.0000000000000002e-31*x23 + 6.0000000000000001e-32*x29 + 1.0e-99) + 5.6556829037999995e-12*X(4)*X(5)*x28*((1.3500000000000001e-9*std::pow(x4, 0.098492999999999997) + 4.4350199999999998e-10*std::pow(x4, 0.55610000000000004) + 3.7408500000000004e-16*std::pow(x4, 2.1825999999999999))/(0.0061910000000000003*std::pow(x4, 1.0461) + 8.9711999999999997e-11*std::pow(x4, 3.0424000000000002) + 3.2575999999999999e-14*std::pow(x4, 3.7740999999999998) + 1.0) + 1.0e-99) + 1.75918975308e-21*X(4)*X(8)*x28 + 4.8065293799999994e-30*X(4)*user_dust2gas_ratio*x25*std::sqrt(x4)*(4.2000000000000002*x28 + 0.20000000000000001)/((1.0 + 10000.0*std::exp(-600.0/(Tdust + 9.9999999999999993e-41)))*(7.9999999999999996e-6*x27 + 0.002*x4 + 0.040000000000000001*std::sqrt(Tdust + x4) + 1.0)) - X(4)*x7*(std::pow(10.0, 42.707410000000003*x32*x50 + 21.360939999999999*x32*x56 - x53 - 2.0273650000000001*x55 - 21467.790000000001*x6 - 1657.4099999999999*x6/(std::pow(std::pow(10.0, -x57 - 8.1313220000000008)*x25, x60) + 1.0) - 142.7664 - (21.360939999999999*x32*x56 + 11.28215*x32*std::log(14254.549999999999*x6 + 1.0) - x53 - 4.7035149999999994*x55 + 70.138370000000009*x61 - 203.11568)/(std::pow(std::pow(10.0, -x57 - 9.3055640000000004)*x25, x60) + 1.0)) + std::pow(10.0, 69.700860000000006*x32*x63 + 43.20243*x33*x54 - 23705.700000000001*x6 - 2080.4099999999999*x6/(std::pow(std::pow(10.0, -x64 - 13.656822)*x25, x65) + 1.0) - 68.422430000000006*x61 - x62 - 178.4239 - (69.700860000000006*x32*x63 + 19.734269999999999*x32*std::log(16780.950000000001*x6 + 1.0) + 37.886913*x55 - 14.509090000000008*x61 - x62 - 307.31920000000002)/(std::pow(std::pow(10.0, -x64 - 14.82123)*x25, x65) + 1.0)) + 3.0000000000000001e-99) - 0.00022681492*x0*((((x4)*(x4)))*(((x4)*(x4))))*x69*((((x66 && x68 && x72 && x73 && x9) ? ( 4.8339620236294848e-32/((x77 + 2.1986273043946046e-56)*(x77 + 2.1986273043946046e-56)) >= 1.0 ) : ( @@ -2502,26 +2498,26 @@ Real rhs_eint(const burn_t& state, 0.031622776601683791 ) : ((x9) ? ( - std::exp((-1.5)*std::log(std::abs(T))) + std::pow(T, -1.5) ) : ( - 0.22169506691164925*std::exp((-1.5)*std::log(std::abs(x10))) + 0.22169506691164925*std::pow(x10, -1.5) )))*std::exp(-470000.0*x12) - 1.5499999999999999e-26*x11*((x8) ? ( 2.3157944032250755 ) : ((x9) ? ( - std::exp((0.36470000000000002)*std::log(std::abs(T))) + std::pow(T, 0.36470000000000002) ) : ( - 1.4423428878129001*std::exp((0.36470000000000002)*std::log(std::abs(x10))) + 1.4423428878129001*std::pow(x10, 0.36470000000000002) ))) - 1.2700000000000001e-21*x16*x19*std::exp(-157809.10000000001*x12) - 4.9500000000000001e-22*x16*x21*std::exp(-631515.0*x12) - 7.4999999999999996e-19*x19*std::exp(-118348.0*x12) - 5.5399999999999998e-17*x21*((x8) ? ( 0.4008667176273028 ) : ((x9) ? ( - std::exp((-0.39700000000000002)*std::log(std::abs(T))) + std::pow(T, -0.39700000000000002) ) : ( - 0.67118683452634376*std::exp((-0.39700000000000002)*std::log(std::abs(x10))) + 0.67118683452634376*std::pow(x10, -0.39700000000000002) )))*std::exp(-473638.0*x12) - 5.0099999999999997e-27*x22*std::exp(-55338.0*x12) - 9.1000000000000001e-27*x22*std::exp(-13179.0*x12) - ((x4 < 2.0) ? ( 0 ) @@ -2826,14 +2822,14 @@ Real rhs_eint(const burn_t& state, ) : ( 0 - ))*std::min(1.0, std::exp((-0.45000000000000001)*std::log(std::abs(std::max(9.9999999999999993e-41, 1.2500000000000001e-10*X(0) + 1.2500000000000001e-10*X(1) + 1.2500000000000001e-10*X(10) + 1.2500000000000001e-10*X(11) + 1.2500000000000001e-10*X(12) + 1.2500000000000001e-10*X(13) + 1.2500000000000001e-10*X(14) + 1.2500000000000001e-10*X(15) + 1.2500000000000001e-10*X(16) + 1.2500000000000001e-10*X(17) + 1.2500000000000001e-10*X(18) + 1.2500000000000001e-10*X(19) + 1.2500000000000001e-10*X(2) + 1.2500000000000001e-10*X(20) + 1.2500000000000001e-10*X(21) + 1.2500000000000001e-10*X(22) + 1.2500000000000001e-10*X(23) + 1.2500000000000001e-10*X(24) + 1.2500000000000001e-10*X(25) + 1.2500000000000001e-10*X(26) + 1.2500000000000001e-10*X(27) + 1.2500000000000001e-10*X(28) + 1.2500000000000001e-10*X(29) + 1.2500000000000001e-10*X(3) + 1.2500000000000001e-10*X(30) + 1.2500000000000001e-10*X(31) + 1.2500000000000001e-10*X(32) + 1.2500000000000001e-10*X(33) + 1.2500000000000001e-10*X(4) + 1.2500000000000001e-10*X(5) + 1.2500000000000001e-10*X(6) + 1.2500000000000001e-10*X(7) + 1.2500000000000001e-10*X(8) + 1.2500000000000001e-10*X(9)))))) + ))*std::min(1.0, std::pow(std::max(9.9999999999999993e-41, 1.2500000000000001e-10*X(0) + 1.2500000000000001e-10*X(1) + 1.2500000000000001e-10*X(10) + 1.2500000000000001e-10*X(11) + 1.2500000000000001e-10*X(12) + 1.2500000000000001e-10*X(13) + 1.2500000000000001e-10*X(14) + 1.2500000000000001e-10*X(15) + 1.2500000000000001e-10*X(16) + 1.2500000000000001e-10*X(17) + 1.2500000000000001e-10*X(18) + 1.2500000000000001e-10*X(19) + 1.2500000000000001e-10*X(2) + 1.2500000000000001e-10*X(20) + 1.2500000000000001e-10*X(21) + 1.2500000000000001e-10*X(22) + 1.2500000000000001e-10*X(23) + 1.2500000000000001e-10*X(24) + 1.2500000000000001e-10*X(25) + 1.2500000000000001e-10*X(26) + 1.2500000000000001e-10*X(27) + 1.2500000000000001e-10*X(28) + 1.2500000000000001e-10*X(29) + 1.2500000000000001e-10*X(3) + 1.2500000000000001e-10*X(30) + 1.2500000000000001e-10*X(31) + 1.2500000000000001e-10*X(32) + 1.2500000000000001e-10*X(33) + 1.2500000000000001e-10*X(4) + 1.2500000000000001e-10*X(5) + 1.2500000000000001e-10*X(6) + 1.2500000000000001e-10*X(7) + 1.2500000000000001e-10*X(8) + 1.2500000000000001e-10*X(9)), -0.45000000000000001)) )) - ((x3 <= x4) ? ( - std::exp((0.92432999999999998*x32*x36 + 21.93385*x32*x42 + 0.77951999999999999*x33*x36*x42 + 0.54962*x33*x41 - 10.19097*x33*x45 + 0.11711000000000001*x36*x37*x44 - 0.54262999999999995*x36*x39*x45 - 0.0083499999999999998*x36*x43*x46 + 0.0027499999999999998*x37*x38 + 0.11864*x37*x40*x42 + 0.62343000000000004*x37*x41*x45 - 0.17333999999999999*x37*x43 - 0.0036600000000000001*x38*x42*x46 + 6.1920000000000003e-5*x38*x43*x47 - 0.00066631000000000004*x38*x44*x48 + 0.0025140000000000002*x38*x45*x49 - 0.076759999999999995*x39*x40 - 1.06447*x39*x41*x42 + 2.1990599999999998*x39*x44 - 0.001482*x40*x43*x48 + 0.017590000000000001*x40*x44*x49 - 0.073660000000000003*x40*x45*x46 + 0.0106*x41*x43*x49 - 0.13768*x41*x44*x46 - 42.567880000000002)*std::log(std::abs(10.0)))*X(12) + std::pow(10.0, 0.92432999999999998*x32*x36 + 21.93385*x32*x42 + 0.77951999999999999*x33*x36*x42 + 0.54962*x33*x41 - 10.19097*x33*x45 + 0.11711000000000001*x36*x37*x44 - 0.54262999999999995*x36*x39*x45 - 0.0083499999999999998*x36*x43*x46 + 0.0027499999999999998*x37*x38 + 0.11864*x37*x40*x42 + 0.62343000000000004*x37*x41*x45 - 0.17333999999999999*x37*x43 - 0.0036600000000000001*x38*x42*x46 + 6.1920000000000003e-5*x38*x43*x47 - 0.00066631000000000004*x38*x44*x48 + 0.0025140000000000002*x38*x45*x49 - 0.076759999999999995*x39*x40 - 1.06447*x39*x41*x42 + 2.1990599999999998*x39*x44 - 0.001482*x40*x43*x48 + 0.017590000000000001*x40*x44*x49 - 0.073660000000000003*x40*x45*x46 + 0.0106*x41*x43*x49 - 0.13768*x41*x44*x46 - 42.567880000000002)*X(12) ) : ( 0 )) - ((x368) ? ( - 16.0*x176*(std::exp((((x178) ? ( + 16.0*x176*(std::pow(10, ((x178) ? ( 21.754701117084799*x177 - 116.63708847273 ) : ((x179) ? ( @@ -3396,7 +3392,7 @@ Real rhs_eint(const burn_t& state, ) : ( -9.9999999999999997e+98 - )))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))*std::log(std::abs(10))) + std::exp((((x178) ? ( + ))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))) + std::pow(10, ((x178) ? ( 2.93326303855088*x177 - 34.858441143158402 ) : ((x179) ? ( @@ -3962,21 +3958,21 @@ Real rhs_eint(const burn_t& state, ) : ( -9.9999999999999997e+98 - ))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))*std::log(std::abs(10)))*Z)*((0.5*X(10) + 0.25*X(11) + 0.25*X(12) + 0.25*X(18) + 0.5*X(19) + 0.75*X(20) + X(23) + 0.25*X(24) + 0.25*X(25) + 0.5*X(26) + 0.5*X(27) + 0.75*X(28) + 0.25*X(3) + 0.25*X(4) + 0.25*X(5) + 0.5*X(8))*(0.5*X(10) + 0.25*X(11) + 0.25*X(12) + 0.25*X(18) + 0.5*X(19) + 0.75*X(20) + X(23) + 0.25*X(24) + 0.25*X(25) + 0.5*X(26) + 0.5*X(27) + 0.75*X(28) + 0.25*X(3) + 0.25*X(4) + 0.25*X(5) + 0.5*X(8))) + (1 - x176)*(((x451 || x453 || x455 < 9.0) ? ( + ))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))*Z)*((0.5*X(10) + 0.25*X(11) + 0.25*X(12) + 0.25*X(18) + 0.5*X(19) + 0.75*X(20) + X(23) + 0.25*X(24) + 0.25*X(25) + 0.5*X(26) + 0.5*X(27) + 0.75*X(28) + 0.25*X(3) + 0.25*X(4) + 0.25*X(5) + 0.5*X(8))*(0.5*X(10) + 0.25*X(11) + 0.25*X(12) + 0.25*X(18) + 0.5*X(19) + 0.75*X(20) + X(23) + 0.25*X(24) + 0.25*X(25) + 0.5*X(26) + 0.5*X(27) + 0.75*X(28) + 0.25*X(3) + 0.25*X(4) + 0.25*X(5) + 0.5*X(8))) + (1 - x176)*(((x451 || x453 || x455 < 9.0) ? ( 0.0 ) : ( - std::exp((-20.360212905100301*x450*x508 - 0.52797032562909996*x450*x510 - 0.0012516455574000001*x450*x512 + 2.45492893e-5*x450*x513 - 2.0491730000000001e-7*x450*x514 - 3.4429418092545001*x450*x521 - 0.55036731592760002*x450*x522 + 0.097706870306800003*x450*x523 + 0.0001978699223*x450*x524 - 2.1488524999999999e-6*x450*x525 - 100.1540946139759*x450*x526 - 0.1220236291332*x450*x527 - 0.0007440436186*x450*x528 + 0.049613747445099998*x450*x535 + 4.8520959300000001e-5*x450*x537 + 0.35426903937339999*x450*x538 - 4.3161844299999999e-5*x450*x539 - 0.060202510145399998*x450*x540 + 20.3649299915947*x450 - 75.299695359826899*x452*x455 - 0.0076851140205*x452*x511 - 3.209781e-7*x452*x514 + 54.085098837550902*x452*x515 - 0.0065739247922000001*x452*x517 + 19.184854295220699*x452*x518 + 0.1615293342819*x452*x519 - 5.6870301499999997e-5*x452*x520 + 0.00014410438979999999*x452*x531 + 0.00026475764900000002*x452*x534 - 47.559803051881197*x452 + 1.1940279999999999e-7*((x455)*(x455)*(x455)*(x455)*(x455)*(x455)*(x455)*(x455)*(x455)) - 3.0836667989033999*x455*x457 - 0.040347262524199998*x455*x459 - 6.5481524799999998e-5*x455*x461 + 8.6699919999999998e-7*x455*x462 + 47.556044042082497*x455*x465 - 86.020612066646194*x455*x466 - 22.2802823459682*x455*x468 - 0.13464931496829999*x455*x470 - 2.8556552777972999*x455*x474 - 0.002904776222*x455*x476 + 2.6345508e-6*x455*x478 + 12.8360066029064*x455*x481 + 0.0278872971314*x455*x484 - 0.044769368153399998*x455*x494 + 0.0048337398222999999*x455*x495 - 7.7849984599999994e-5*x455*x496 - 0.0096135394491999992*x455*x499 - 7.5822643719078*x455 - 0.043497398505999998*x456*x510 - 3.4867000999999999e-5*x456*x512 + 3.5333219999999999e-7*x456*x513 - 2.4301583262415001*x456*x521 + 0.1712810070345*x456*x522 - 0.0064759322805000002*x456*x523 - 1.6696918e-6*x456*x524 + 3.5133050137291*x456*x526 - 0.00070595390469999997*x456*x535 + 0.00077362074520000005*x456*x540 - 56.915400919787402*x456 + 0.0016854860852*x457*x510 - 3.3478540799999997e-5*x457*x511 + 1.697195e-7*x457*x512 + 0.28173719990879997*x457*x521 - 0.0137884415355*x457*x522 + 0.00029133323910000001*x457*x523 + 0.074096619680599995*x457*x526 + 0.00020615397360000001*x457*x527 + 3.6266469e-6*x457*x535 + 12.6066721363912*x457 - 0.032618157779500001*x458*x508 + 0.001902687729*x458*x509 - 6.0232547200000002e-5*x458*x510 + 7.8674779999999999e-7*x458*x511 - 0.0135100513841*x458*x521 + 0.00049797288530000004*x458*x522 - 6.7732749e-6*x458*x523 - 0.014159748673100001*x458*x526 - 8.5587269900000006e-5*x458*x536 - 2.7928933908445002*x458 - 4.87054635e-5*x459*x509 + 5.6505170000000002e-7*x459*x510 - 5.2183217e-6*x459*x522 + 0.00040007346859999998*x459*x526 + 0.34205790108779999*x459 - 4.61090426e-5*x460*x508 + 7.1519579999999996e-7*x460*x509 - 2.0907143999999999e-6*x460*x521 - 3.8957497000000001e-6*x460*x526 - 0.024213057736899999*x460 + 4.5347009999999999e-7*x461*x508 + 0.00095884092259999995*x461 - 1.7797473799999999e-5*x462 + 6.9368000000000003e-8*x463 - 14.1815605361866*x464*x508 - 2.7554623671076999*x464*x509 + 0.0075222770970999999*x464*x511 - 0.0004669134188*x464*x512 + 8.1984912999999997e-6*x464*x513 - 1.0905152389177999*x464*x516 + 0.0242714411203*x464*x530 + 0.013369313361399999*x464*x532 - 0.00049248667509999998*x464*x533 + 6.3431601000000002e-6*x464*x534 + 0.90563183927519997*x464 + 1.6463007811879999*x465*x509 - 0.16362513598780001*x465*x510 - 6.1688840300000002e-5*x465*x512 + 0.042057431356300003*x465*x516 + 0.0033260418458999999*x465*x518 + 0.00033784509460000001*x465*x530 + 33.569263775829498*x465 + 0.46079059024460001*x466*x509 + 6.0713992099999998e-5*x466*x511 - 1.8297653193047001*x466*x515 + 0.0034593370547*x466*x516 - 0.089280005816900002*x466*x518 + 1.11981398e-5*x466*x519 - 0.017098171269000002*x466*x529 + 5.7414321999999997e-5*x466*x530 + 0.00031254232230000002*x466*x532 - 17.3552073528923*x466 - 0.10105727851019999*x467*x509 - 0.0007373397229*x467*x510 - 0.00028850397480000001*x467*x516 + 0.086398191443600003*x467*x518 - 57.703059324335001*x467 - 0.1921360560176*x468*x508 + 0.0115453052291*x468*x509 - 0.0093329110050999996*x468*x518 + 47.789440447875599*x468 - 0.0196095956543*x469*x508 - 22.2858836676964*x469 + 5.3339578239099001*x470 - 0.51041682043059999*x471 + 14.164582630669701*x472 - 83.354317190196298*x473 - 1.6101414000000001e-6*x474*x511 + 15.0991587239831*x474 - 1.2258681523532999*x475 + 0.059113880622099997*x476 - 0.0017177123872999999*x477 + 2.5732575999999999e-5*x478 - 1.925221e-7*x479 + 14.8241641825581*x480*x508 - 0.0010489094365*x480*x511 + 6.7124765000000004e-6*x480*x512 + 61.439278554514303*x480 - 1.6009057790223*x481*x508 - 0.0024989369339000001*x481*x510 + 66.325985989249105*x481 - 7.8551891148808997*x482 - 12.244271306143601*x483 + 0.0051192934768000001*x484*x508 + 6.4011003727495996*x484 - 1.0865832613576001*x485 + 0.056776674990399997*x486 - 0.52947614193869996*x487*x508 + 8.3381011e-6*x487*x511 + 1.311635494626*x487 - 2.7623083129385999*x488 - 0.00039222307050000002*x489*x508 + 1.5185419100000001e-5*x489*x509 + 0.27348164083190002*x489 - 0.019100980547200001*x490 + 0.00074922383839999996*x491 - 1.22400281e-5*x492 + 0.049242766990500002*x493*x508 + 7.3490350000000004e-7*x493*x510 + 1.7166022505292999*x493 + 0.0011489194243000001*x494*x508 - 3.0581260599999998e-5*x494*x509 + 0.50309367074760003*x494 + 0.0013474467272000001*x495 + 0.0003300557909*x496 - 1.7934860499999999e-5*x497 - 3.6187820926523999*x498 - 0.090943285887500003*x499 - 0.0111295294489*x500 + 0.0002140404241*x501 + 1.8509211874314*x502 + 0.056834102826200002*x503 + 0.00043558306989999998*x504 - 0.56101466485609996*x505 - 0.0058434230991000003*x506 + 0.060192127702100001*x507 - 23.741334917810999*x508 + 10.4159294438888*x509 - 2.0869172860516998*x510 + 0.22605546648569999*x511 - 0.014638487079*x512 + 0.00057260663920000003*x513 - 1.25535688e-5*x514 + 23.837391969558599*x515 + 4.6242788943301996*x516 + 0.034087068909300001*x517 + 27.377460699605798*x518 + 0.44068107103589998*x519 + 0.0021941995751000001*x520 + 8.1723278893797993*x521 + 1.276214911337*x522 + 0.031193196973900002*x523 + 7.4001529100000004e-5*x524 + 1.16367671e-5*x525 + 39.764294983229803*x526 + 72.523921027240604*x527 + 3.0795133468160998*x528 - 5.8011752165720996*x529 + 0.66942669331799998*x530 + 0.0016319787938*x531 + 0.42168773563220002*x532 - 0.037291536482600002*x533 + 0.0016101104869*x534 + 0.094683646925100007*x535 + 21.686787075207999*x536 + 0.0046026070559999997*x537 - 18.068679356687898*x538 + 0.0024703528802999999*x539 + 3.3120227111728*x540 + 168.09888262271619)*std::log(std::abs(10)))*X(25)*x448 - )) + ((x451 || x453 || x541 < 9.0) ? ( + std::pow(10, -19.9926824317543*x450*x508 - 0.52101304223440004*x450*x510 - 0.0012371771563*x450*x512 - 2.025523e-7*x450*x514 - 99.329623204854599*x450*x526 - 0.15773801363449999*x450*x528 + 0.049654764947699999*x450*x535 + 0.34727291958810003*x450*x536 - 4.3902720399999997e-5*x450*x537 - 0.059327608505199997*x450*x538 + 20.0650035564842*x450 - 74.713459337991395*x452*x455 - 0.0065204020911999996*x452*x511 + 1.27871529e-5*x452*x513 - 3.3175559999999997e-7*x452*x514 - 0.53411509369950005*x452*x516 - 2.1390089999999998e-6*x452*x518 + 19.1923861231489*x452*x519 + 0.1616529475864*x452*x520 - 5.6705485e-5*x452*x521 + 0.17141782863419999*x452*x529 + 0.0001439867417*x452*x530 + 0.00029186501660000001*x452*x532 + 0.00049849519790000005*x452*x533 - 47.517481143431297*x452 + 1.191057e-7*((x455)*(x455)*(x455)*(x455)*(x455)*(x455)*(x455)*(x455)*(x455)) - 3.0330528171500002*x455*x457 - 0.040219489778600001*x455*x459 - 6.5397406699999996e-5*x455*x461 + 8.6616289999999995e-7*x455*x462 + 47.779027220289002*x455*x465 - 22.179350016916899*x455*x468 + 3.0661944361070002*x455*x469 - 0.13396404931969999*x455*x470 + 53.986875199917399*x455*x472 - 2.8571295162038002*x455*x474 - 0.0029105106853000002*x455*x476 + 2.6321730999999998e-6*x455*x478 + 12.300420924217301*x455*x481 - 0.044261169582100003*x455*x494 + 0.0048230242393000004*x455*x495 - 7.7770003900000002e-5*x455*x496 + 0.00020966390190000001*x455*x503 - 7.8949806449478004*x455 - 5.7277298688334*x456*x508 - 0.043223906760100003*x456*x510 - 3.4765150899999998e-5*x456*x512 + 3.5276910000000002e-7*x456*x513 - 2.4330613855768002*x456*x522 - 0.0064765491634*x456*x524 - 1.6670821999999999e-6*x456*x525 + 3.5303617069073998*x456*x526 - 0.00070624420869999996*x456*x535 - 55.772982130894199*x456 - 0.036865510067599999*x457*x509 - 3.3158269899999997e-5*x457*x511 + 1.670798e-7*x457*x512 + 0.28210704022389999*x457*x522 - 0.0138094104646*x457*x523 + 0.072221694161200001*x457*x526 - 0.0096793760875000003*x457*x527 + 3.6315933999999999e-6*x457*x535 + 12.4564129877271*x457 - 0.032344369329199998*x458*x508 - 5.9920136499999999e-5*x458*x510 + 7.8368329999999999e-7*x458*x511 - 0.0135249926194*x458*x522 - 6.7795023000000004e-6*x458*x524 - 0.0141013436985*x458*x526 - 2.7798956975681*x458 - 4.84912661e-5*x459*x509 + 5.6266539999999998e-7*x459*x510 + 0.00026507121149999999*x459*x522 + 0.00039906215160000002*x459*x526 + 0.34114118166250001*x459 - 4.59774543e-5*x460*x508 + 7.1354209999999996e-7*x460*x509 - 2.0937332000000001e-6*x460*x522 - 3.8862170999999999e-6*x460*x526 - 0.024161129759300001*x460 + 4.5234329999999999e-7*x461*x508 + 0.00095672488770000004*x461 - 1.7744158299999999e-5*x462 + 6.8746300000000002e-8*x463 - 14.5067682237128*x464*x508 - 2.7051669318233*x464*x509 + 0.0076047959769000003*x464*x511 - 0.00046768386289999998*x464*x512 + 8.1961122000000004e-6*x464*x513 - 1.0900116892503*x464*x516 - 0.0010505811818000001*x464*x517 + 0.024278624911799999*x464*x529 + 8.3404412999999993e-6*x464*x530 + 0.013406972941*x464*x531 + 1.51839181e-5*x464*x533 + 6.3477803999999997e-6*x464*x534 + 0.5975764242828*x464 + 21.656299270225201*x465*x508 + 1.6408874267053*x465*x509 - 0.16313400868620001*x465*x510 + 0.0045884667059000002*x465*x511 - 6.1550235000000001e-5*x465*x512 - 0.0009430205952*x465*x519 + 33.770910375285403*x465 + 0.45993780466700002*x466*x509 + 6.10066494e-5*x466*x511 - 1.6484720484398001*x466*x515 + 0.0035753188361000001*x466*x516 - 0.0887188773973*x466*x519 + 1.18590686e-5*x466*x520 + 5.71056327e-5*x466*x529 + 0.00031294391269999998*x466*x531 - 17.952175722342002*x466 - 0.10094611338950001*x467*x509 - 0.00073857723610000004*x467*x510 - 0.00029541020049999998*x467*x516 + 0.086355739615899998*x467*x519 - 58.190545272689903*x467 - 0.1934292691551*x468*x508 + 0.011545689602299999*x468*x509 + 0.031792738899999999*x468*x515 - 0.0093312353852000006*x468*x519 + 48.447216244679304*x468 - 0.019499399200999998*x469*x508 - 0.00092824833489999999*x469*x515 - 22.5429482843711*x469 + 5.3785945835476001*x470 - 0.51340012747920005*x471 - 3.5366521183169*x472*x508 + 0.096475646631500003*x472*x510 - 0.0065234438004999999*x472*x511 + 0.00019677491820000001*x472*x512 + 13.3445171763044*x472 - 82.661705464070906*x473 - 1.6155785e-6*x474*x511 + 15.0433737682245*x474 - 1.2233428327708*x475 - 5.2231879999999999e-6*x476*x509 + 0.059053482857800002*x476 - 0.0017181002083999999*x477 + 2.57966965e-5*x478 - 1.9363439999999999e-7*x479 + 14.7789280685707*x480*x508 + 6.7323494000000004e-6*x480*x512 + 61.052074661994297*x480 - 1.5729471883424999*x481*x508 + 0.041397466824299997*x481*x509 - 0.0024924218309000002*x481*x510 + 4.8513968500000001e-5*x481*x511 + 66.617992420320206*x481 - 7.9160224042617999*x482 - 12.2325618905237*x483 + 0.0050738190381999998*x484*x508 + 6.3942922953709003*x484 - 1.0845953667873001*x485 + 0.056584855567100002*x486 - 0.52970771916649995*x487*x508 + 0.46011881334179999*x487 - 0.00049289101450000002*x488*x509 - 2.7010241272834001*x488 - 0.00039286300609999998*x489*x508 + 0.27097298702860001*x489 - 0.019032918740700001*x490 + 0.00074815639190000001*x491 - 1.22336158e-5*x492 + 0.049179142615399997*x493*x508 + 0.00033973388480000003*x493*x509 + 7.3643890000000003e-7*x493*x510 + 2.2122323463722*x493 + 0.0011420736754*x494*x508 - 3.0554493800000001e-5*x494*x509 + 0.47664689686630002*x494 - 8.5497981599999998e-5*x495*x508 + 0.0020981266227000001*x495 + 0.00031715525929999999*x496 - 1.7836264199999999e-5*x497 - 0.017075233931400001*x498*x508 - 3.7983452612036999*x498 - 0.084331930384299994*x499 - 0.0112432241909*x500 + 0.00021497635970000001*x501 + 0.00077184455469999995*x502*x508 + 1.8908951729129999*x502 + 0.055918940364399997*x503 + 0.00044264301000000001*x504 - 0.56600117984999998*x505 - 0.0057888990350000004*x506 + 0.060458018568200003*x507 - 23.338685288072199*x508 + 10.277319996911199*x509 - 2.0652126800509998*x510 + 0.22418121793740001*x511 - 0.0145430957054*x512 + 0.00056976337950000004*x513 - 1.25080523e-5*x514 + 23.215967916352799*x515 + 4.554675353975*x516 + 0.033676237701199999*x517 + 2.42668793e-5*x518 + 26.887607489364399*x519 + 0.43769189354580001*x520 + 0.0021901593191999999*x521 + 7.9246235535301999*x522 + 1.367807670768*x523 + 0.016633908132799999*x524 + 2.33933078e-5*x525 + 40.107774411679202*x526 - 85.5850290476819*x527 + 72.171467367178494*x528 + 0.66354267554680002*x529 + 0.0016246752704*x530 + 0.41531456331650002*x531 + 0.0016694722784000001*x532 + 0.0018897404184*x533 + 0.0016025825331000001*x534 + 0.091638174636899994*x535 - 18.049219731801902*x536 + 0.0024669223464000001*x537 + 3.3145678538165999*x538 + 167.41038846894398)*X(25)*x448 + )) + ((x451 || x453 || x539 < 9.0) ? ( 0.0 ) : ( - std::exp((-15.6735419271706*x450*x541 + 8.5045766567269006*x450*x542 - 1.4972563957478*x450*x543 + 0.13047868949079999*x450*x544 - 0.0062773436733000003*x450*x545 + 0.00017407734850000001*x450*x546 - 2.9054825999999999e-6*x450*x547 + 2.7458599999999999e-8*x450*x548 - 9.8760835772856996*x450 + 5.4586553412771002*x452*x541 + 1.9760786678008*x452*x542 - 0.67472957414050005*x452*x543 + 0.17540631800670001*x452*x544 - 0.022834960298200001*x452*x545 + 0.0014418794812000001*x452*x546 - 4.4932579099999998e-5*x452*x547 + 5.6193980000000003e-7*x452*x548 + 4.5634379481376*x452 - 3.5528602795922*x456*x541 - 0.74231825885180003*x456*x542 + 0.21099124651350001*x456*x543 - 0.019394333998500001*x456*x544 + 0.00084390984799999997*x456*x545 - 1.7576673099999998e-5*x456*x546 + 1.400519e-7*x456*x547 + 20.9658337145544*x456 + 1.3758733773666001*x457*x541 - 0.053043503323299997*x457*x542 - 0.0076961297644000002*x457*x543 + 0.00082299662219999996*x457*x544 - 2.7446162600000001e-5*x457*x545 + 2.9510110000000001e-7*x457*x546 - 5.8665160212550997*x457 - 0.13655529547520001*x458*x541 + 0.006878461933*x458*x542 + 6.2950829000000007e-5*x458*x543 - 1.5859306199999999e-5*x458*x544 + 3.5343039999999998e-7*x458*x545 + 0.59960535355900002*x458 + 0.0071954965516000001*x459*x541 - 0.00025848591050000002*x459*x542 + 1.5560391e-6*x459*x543 + 5.4168300000000002e-8*x459*x544 - 0.0278172887885*x459 - 0.00024943263539999998*x460*x541 + 5.1790795999999998e-6*x460*x542 + 2.6661200000000001e-8*x460*x543 + 0.00013591247090000001*x460 + 5.2732847e-6*x461*x541 - 9.1609499999999997e-8*x461*x542 + 3.0908317500000001e-5*x461 - 3.8680900000000002e-8*x462*x541 + 1.494538e-7*x462 - 4.7120000000000001e-8*x463 + 1.3440244509784001*x464*x541 - 2.6323935986257001*x464*x542 - 0.63586010955950001*x464*x543 + 0.14674619922349999*x464*x544 - 0.0096970290444000004*x464*x545 + 0.00029708562290000002*x464*x546 - 3.6368400999999998e-6*x464*x547 + 2.5161884420602001*x464 - 1.0237576973172*x465*x541 + 6.1535122772876996*x465*x542 - 0.76291955278949997*x465*x543 + 0.0265998154454*x465*x544 - 0.00051001248800000001*x465*x545 + 4.4068406999999997e-6*x465*x546 + 2.4341312509119*x465 - 8.5623467147410999*x466*x541 - 0.2813319983316*x466*x542 + 0.1577743012293*x466*x543 - 0.0037305896661000001*x466*x544 + 3.8073638100000003e-5*x466*x545 - 0.53967782227049999*x466 + 3.3110368609715*x467*x541 - 0.50249921046150003*x467*x542 - 0.0115254870206*x467*x543 + 8.8841644699999995e-5*x467*x544 - 3.3247503247040999*x467 + 0.3220442657109*x468*x541 + 0.1148838792372*x468*x542 + 0.0006791871357*x468*x543 + 2.1452691279259999*x468 - 0.2464800834019*x469*x541 - 0.0100100109153*x469*x542 - 1.397474402394*x469 + 0.030497966722699998*x470*x541 + 0.38181813023950001*x470 - 0.039390464892000002*x471 - 17.595716520764402*x472*x541 + 4.3434614562203002*x472*x542 - 0.55339310618509996*x472*x543 + 0.042667300411499999*x472*x544 - 0.0020317800997999998*x472*x545 + 5.7698430099999997e-5*x472*x546 - 7.6956910000000005e-7*x472*x547 - 11.099775151034899*x472 + 1.6907992647648*x473*x541 - 0.90218066191159996*x473*x542 + 0.1063179056835*x473*x543 - 0.0055157020021999997*x473*x544 + 0.00013315511249999999*x473*x545 - 1.1633197e-6*x473*x546 + 18.1536876840762*x473 - 0.32374948483510002*x474*x541 + 0.091257030742499995*x474*x542 - 0.0074216902436000003*x474*x543 + 0.0002384532436*x474*x544 - 2.6686475000000001e-6*x474*x545 - 1.2096109992180999*x474 + 0.023144128920999998*x475*x541 - 0.0040382879026999996*x475*x542 + 0.00023255522560000001*x475*x543 - 4.0639182000000001e-6*x475*x544 + 0.051205700689499999*x475 - 0.00088091977179999997*x476*x541 + 5.9579273200000003e-5*x476*x542 - 1.7863086e-6*x476*x543 - 0.0035559076976*x476 + 2.83458978e-5*x477*x541 - 4.0861769999999998e-7*x477*x542 + 0.00068851345740000001*x477 - 3.6354329999999999e-7*x478*x541 - 5.5323727700000002e-5*x478 + 1.4514925000000001e-6*x479 + 5.4673964678251004*x480*x541 + 1.0660139132327*x480*x542 - 0.23292294092970001*x480*x543 + 0.015672395274900001*x480*x544 - 0.00050902837310000002*x480*x545 + 7.1523700999999999e-6*x480*x546 - 0.1628147174764*x480 - 15.963263559048499*x481*x541 + 1.8055455053091001*x481*x542 - 0.084817027749599996*x481*x543 + 0.0022681696148999998*x481*x544 - 3.0409404800000001e-5*x481*x545 + 18.501865221499301*x481 + 0.59703124368659999*x482*x541 - 0.1527626342914*x482*x542 + 0.0026212163114*x482*x543 + 1.4218711e-5*x482*x544 + 12.4406982350406*x482 + 0.2981554850184*x483*x541 + 0.018278996070200001*x483*x542 - 0.00035799879249999997*x483*x543 - 2.6159083748851*x483 - 0.094078833140799997*x484*x541 + 5.5556733700000001e-5*x484*x542 - 0.57547145134050004*x484 + 0.0053807608027000003*x485*x541 + 0.30498502561469998*x485 - 0.0301525351682*x486 + 3.9298735211797999*x487*x541 - 0.48456767114499999*x487*x542 + 0.028480890781399999*x487*x543 - 0.00077447721980000001*x487*x544 + 6.8596344999999997e-6*x487*x545 - 22.845731870181901*x487 - 0.34414960089390001*x488*x541 + 0.0282438466058*x488*x542 - 0.0010344421251*x488*x543 + 1.4952163600000001e-5*x488*x544 + 1.8754275566329*x488 + 0.0146909729356*x489*x541 - 0.00075693974790000005*x489*x542 + 1.28837281e-5*x489*x543 - 0.096485453374899999*x489 - 0.00027764600010000002*x490*x541 + 7.3709799999999997e-6*x490*x542 + 3.4704941499999997e-5*x490 + 2.1562743999999999e-6*x491*x541 + 0.00015311135120000001*x491 - 4.4968760000000003e-6*x492 + 0.65217956369329999*x493*x541 - 0.050484721527699998*x493*x542 + 0.0016061902897999999*x493*x543 - 1.23009568e-5*x493*x544 + 5.4283657019531004*x493 - 0.021749987409700001*x494*x541 + 0.0012715116938*x494*x542 - 3.2741474300000002e-5*x494*x543 - 0.063783476515700002*x494 + 0.0001901837358*x495*x541 - 9.4105437999999996e-6*x495*x542 + 0.013330079761100001*x495 - 5.5285418999999996e-6*x496*x541 - 0.00027119803169999999*x496 + 2.9662535000000002e-6*x497 - 0.0099310448876000004*x498*x541 + 0.000450542098*x498*x542 - 1.2332561399999999e-5*x498*x543 - 3.9627091042916001*x498 + 0.00049523192299999995*x499*x541 + 5.8076977699999998e-5*x499*x542 + 0.018026677275599999*x499 + 4.15462625e-5*x500*x541 - 0.0026947914996999999*x500 + 3.1513558400000001e-5*x501 - 0.00073381158360000004*x502*x541 - 0.00014651520950000001*x502*x542 + 1.2170419193973001*x502 - 0.00036387813799999998*x503*x541 + 0.0029760547547000001*x503 + 9.7334827099999994e-5*x504 + 0.0011345025923*x505*x541 - 0.20926143712279999*x505 - 2.6366922300000002e-5*x506 + 0.0136346787342*x507 - 3.5612899999999997e-8*((x541)*(x541)*(x541)*(x541)*(x541)*(x541)*(x541)*(x541)*(x541)) - 0.118672054406*x541 - 0.050008209197199997*x542 - 0.85938724395099997*x543 + 0.25027983039980001*x544 - 0.0340585153781*x545 + 0.002684927785*x546 - 0.0001252971097*x547 + 3.2291785999999998e-6*x548 - 10.717982853902399)*std::log(std::abs(10)))*X(27)*x448 - )) + ((x549 < -2.0 || x550 < 0.47712125 || x551 < -18.0) ? ( + std::pow(10, -15.6735419271706*x450*x539 + 8.5045766567269006*x450*x540 - 1.4972563957478*x450*x541 + 0.13047868949079999*x450*x542 - 0.0062773436733000003*x450*x543 + 0.00017407734850000001*x450*x544 - 2.9054825999999999e-6*x450*x545 + 2.7458599999999999e-8*x450*x546 - 9.8760835772856996*x450 + 5.4586553412771002*x452*x539 + 1.9760786678008*x452*x540 - 0.67472957414050005*x452*x541 + 0.17540631800670001*x452*x542 - 0.022834960298200001*x452*x543 + 0.0014418794812000001*x452*x544 - 4.4932579099999998e-5*x452*x545 + 5.6193980000000003e-7*x452*x546 + 4.5634379481376*x452 - 3.5528602795922*x456*x539 - 0.74231825885180003*x456*x540 + 0.21099124651350001*x456*x541 - 0.019394333998500001*x456*x542 + 0.00084390984799999997*x456*x543 - 1.7576673099999998e-5*x456*x544 + 1.400519e-7*x456*x545 + 20.9658337145544*x456 + 1.3758733773666001*x457*x539 - 0.053043503323299997*x457*x540 - 0.0076961297644000002*x457*x541 + 0.00082299662219999996*x457*x542 - 2.7446162600000001e-5*x457*x543 + 2.9510110000000001e-7*x457*x544 - 5.8665160212550997*x457 - 0.13655529547520001*x458*x539 + 0.006878461933*x458*x540 + 6.2950829000000007e-5*x458*x541 - 1.5859306199999999e-5*x458*x542 + 3.5343039999999998e-7*x458*x543 + 0.59960535355900002*x458 + 0.0071954965516000001*x459*x539 - 0.00025848591050000002*x459*x540 + 1.5560391e-6*x459*x541 + 5.4168300000000002e-8*x459*x542 - 0.0278172887885*x459 - 0.00024943263539999998*x460*x539 + 5.1790795999999998e-6*x460*x540 + 2.6661200000000001e-8*x460*x541 + 0.00013591247090000001*x460 + 5.2732847e-6*x461*x539 - 9.1609499999999997e-8*x461*x540 + 3.0908317500000001e-5*x461 - 3.8680900000000002e-8*x462*x539 + 1.494538e-7*x462 - 4.7120000000000001e-8*x463 + 1.3440244509784001*x464*x539 - 2.6323935986257001*x464*x540 - 0.63586010955950001*x464*x541 + 0.14674619922349999*x464*x542 - 0.0096970290444000004*x464*x543 + 0.00029708562290000002*x464*x544 - 3.6368400999999998e-6*x464*x545 + 2.5161884420602001*x464 - 1.0237576973172*x465*x539 + 6.1535122772876996*x465*x540 - 0.76291955278949997*x465*x541 + 0.0265998154454*x465*x542 - 0.00051001248800000001*x465*x543 + 4.4068406999999997e-6*x465*x544 + 2.4341312509119*x465 - 8.5623467147410999*x466*x539 - 0.2813319983316*x466*x540 + 0.1577743012293*x466*x541 - 0.0037305896661000001*x466*x542 + 3.8073638100000003e-5*x466*x543 - 0.53967782227049999*x466 + 3.3110368609715*x467*x539 - 0.50249921046150003*x467*x540 - 0.0115254870206*x467*x541 + 8.8841644699999995e-5*x467*x542 - 3.3247503247040999*x467 + 0.3220442657109*x468*x539 + 0.1148838792372*x468*x540 + 0.0006791871357*x468*x541 + 2.1452691279259999*x468 - 0.2464800834019*x469*x539 - 0.0100100109153*x469*x540 - 1.397474402394*x469 + 0.030497966722699998*x470*x539 + 0.38181813023950001*x470 - 0.039390464892000002*x471 - 17.595716520764402*x472*x539 + 4.3434614562203002*x472*x540 - 0.55339310618509996*x472*x541 + 0.042667300411499999*x472*x542 - 0.0020317800997999998*x472*x543 + 5.7698430099999997e-5*x472*x544 - 7.6956910000000005e-7*x472*x545 - 11.099775151034899*x472 + 1.6907992647648*x473*x539 - 0.90218066191159996*x473*x540 + 0.1063179056835*x473*x541 - 0.0055157020021999997*x473*x542 + 0.00013315511249999999*x473*x543 - 1.1633197e-6*x473*x544 + 18.1536876840762*x473 - 0.32374948483510002*x474*x539 + 0.091257030742499995*x474*x540 - 0.0074216902436000003*x474*x541 + 0.0002384532436*x474*x542 - 2.6686475000000001e-6*x474*x543 - 1.2096109992180999*x474 + 0.023144128920999998*x475*x539 - 0.0040382879026999996*x475*x540 + 0.00023255522560000001*x475*x541 - 4.0639182000000001e-6*x475*x542 + 0.051205700689499999*x475 - 0.00088091977179999997*x476*x539 + 5.9579273200000003e-5*x476*x540 - 1.7863086e-6*x476*x541 - 0.0035559076976*x476 + 2.83458978e-5*x477*x539 - 4.0861769999999998e-7*x477*x540 + 0.00068851345740000001*x477 - 3.6354329999999999e-7*x478*x539 - 5.5323727700000002e-5*x478 + 1.4514925000000001e-6*x479 + 5.4673964678251004*x480*x539 + 1.0660139132327*x480*x540 - 0.23292294092970001*x480*x541 + 0.015672395274900001*x480*x542 - 0.00050902837310000002*x480*x543 + 7.1523700999999999e-6*x480*x544 - 0.1628147174764*x480 - 15.963263559048499*x481*x539 + 1.8055455053091001*x481*x540 - 0.084817027749599996*x481*x541 + 0.0022681696148999998*x481*x542 - 3.0409404800000001e-5*x481*x543 + 18.501865221499301*x481 + 0.59703124368659999*x482*x539 - 0.1527626342914*x482*x540 + 0.0026212163114*x482*x541 + 1.4218711e-5*x482*x542 + 12.4406982350406*x482 + 0.2981554850184*x483*x539 + 0.018278996070200001*x483*x540 - 0.00035799879249999997*x483*x541 - 2.6159083748851*x483 - 0.094078833140799997*x484*x539 + 5.5556733700000001e-5*x484*x540 - 0.57547145134050004*x484 + 0.0053807608027000003*x485*x539 + 0.30498502561469998*x485 - 0.0301525351682*x486 + 3.9298735211797999*x487*x539 - 0.48456767114499999*x487*x540 + 0.028480890781399999*x487*x541 - 0.00077447721980000001*x487*x542 + 6.8596344999999997e-6*x487*x543 - 22.845731870181901*x487 - 0.34414960089390001*x488*x539 + 0.0282438466058*x488*x540 - 0.0010344421251*x488*x541 + 1.4952163600000001e-5*x488*x542 + 1.8754275566329*x488 + 0.0146909729356*x489*x539 - 0.00075693974790000005*x489*x540 + 1.28837281e-5*x489*x541 - 0.096485453374899999*x489 - 0.00027764600010000002*x490*x539 + 7.3709799999999997e-6*x490*x540 + 3.4704941499999997e-5*x490 + 2.1562743999999999e-6*x491*x539 + 0.00015311135120000001*x491 - 4.4968760000000003e-6*x492 + 0.65217956369329999*x493*x539 - 0.050484721527699998*x493*x540 + 0.0016061902897999999*x493*x541 - 1.23009568e-5*x493*x542 + 5.4283657019531004*x493 - 0.021749987409700001*x494*x539 + 0.0012715116938*x494*x540 - 3.2741474300000002e-5*x494*x541 - 0.063783476515700002*x494 + 0.0001901837358*x495*x539 - 9.4105437999999996e-6*x495*x540 + 0.013330079761100001*x495 - 5.5285418999999996e-6*x496*x539 - 0.00027119803169999999*x496 + 2.9662535000000002e-6*x497 - 0.0099310448876000004*x498*x539 + 0.000450542098*x498*x540 - 1.2332561399999999e-5*x498*x541 - 3.9627091042916001*x498 + 0.00049523192299999995*x499*x539 + 5.8076977699999998e-5*x499*x540 + 0.018026677275599999*x499 + 4.15462625e-5*x500*x539 - 0.0026947914996999999*x500 + 3.1513558400000001e-5*x501 - 0.00073381158360000004*x502*x539 - 0.00014651520950000001*x502*x540 + 1.2170419193973001*x502 - 0.00036387813799999998*x503*x539 + 0.0029760547547000001*x503 + 9.7334827099999994e-5*x504 + 0.0011345025923*x505*x539 - 0.20926143712279999*x505 - 2.6366922300000002e-5*x506 + 0.0136346787342*x507 - 3.5612899999999997e-8*((x539)*(x539)*(x539)*(x539)*(x539)*(x539)*(x539)*(x539)*(x539)) - 0.118672054406*x539 - 0.050008209197199997*x540 - 0.85938724395099997*x541 + 0.25027983039980001*x542 - 0.0340585153781*x543 + 0.002684927785*x544 - 0.0001252971097*x545 + 3.2291785999999998e-6*x546 - 10.717982853902399)*X(27)*x448 + )) + ((x547 < -2.0 || x548 < 0.47712125 || x549 < -18.0) ? ( 0.0 ) : ( - std::exp((-3.5355e-9*std::exp((10)*std::log(std::abs(x549))) + 0.0019497087522999999*x549*x550*x569 + 2.5033869e-6*x549*x550*x571 + 2.425819e-7*x549*x550*x572 + 2.74e-11*x549*x550*x575 - 0.25171915905909997*x549*x550 + 0.16154758856199999*x549*x551*x560 + 0.10319386725309999*x549*x551*x562 + 0.0091058356367999995*x549*x551*x564 + 7.0547189099999996e-5*x549*x551*x566 + 0.028910314039899999*x549*x551 + 0.00031554761819999999*x549*x560*x570 + 6.9799317000000001e-6*x549*x560*x571 - 8.3654700000000001e-8*x549*x560*x572 - 0.13936594077190001*x549*x560 + 0.0031773255127*x549*x561*x569 - 2.2029417000000002e-6*x549*x561*x571 + 3.32601e-8*x549*x561*x572 + 1.1211e-9*x549*x561*x573 + 0.40330227063039997*x549*x561 + 7.7272750100000005e-5*x549*x562*x570 + 4.3411179999999998e-7*x549*x562*x571 - 2.2911000000000001e-9*x549*x562*x572 + 0.053531392101600003*x549*x562 + 0.00074114505229999996*x549*x563*x569 - 0.4569416534031*x549*x563 + 7.4963150000000001e-7*x549*x564*x570 + 0.3477567512959*x549*x564 + 1.1028151899999999e-5*x549*x565*x569 - 0.11816708728189999*x549*x565 + 0.019397806516399999*x549*x566 - 0.0012520513382999999*x549*x567 + 0.00089752818999999998*x549*x569 + 8.7854088899999995e-5*x549*x570 - 1.36574957e-5*x549*x571 - 9.4252039999999997e-7*x549*x572 + 2.2287e-8*x549*x573 + 2.0822000000000001e-9*x549*x574 - 1.6900000000000001e-11*x549*x575 - 1.1e-12*x549*x576 - 5.0747599999999997e-8*x549*x581 - 0.092472151351600002*x549*x582 - 0.00022388895510000001*x549*x593 + 0.1098450355633*x549 + 0.0028626466126000002*std::exp((10)*std::log(std::abs(x550))) + 0.036943070713000001*x550*x551*x552 + 9.9116502200000001e-5*x550*x551*x555 + 1.5874459999999999e-7*x550*x551*x557 + 0.0208809736142*x550*x551 + 0.00074934921520000003*x550*x552*x569 - 2.1007040000000002e-6*x550*x552*x571 + 2.3848400000000001e-8*x550*x552*x572 + 1.7461e-9*x550*x552*x573 - 9.6999999999999995e-12*x550*x552*x574 + 0.12448499329609999*x550*x552 + 1.706318e-7*x550*x553*x571 - 2.0289999999999999e-9*x550*x553*x572 - 6.3199999999999999e-11*x550*x553*x573 + 0.0021851183530000001*x550*x553 + 9.2263300000000005e-6*x550*x554*x569 + 5.5200690000000001e-7*x550*x554*x570 - 7.8312e-9*x550*x554*x571 + 9.6999999999999995e-12*x550*x554*x572 + 0.0025381599434999998*x550*x554 - 2.020005e-7*x550*x555*x569 + 1.5340000000000001e-10*x550*x555*x571 - 0.0004386554592*x550*x555 + 9.9660000000000005e-10*x550*x556*x570 + 3.4105511399999997e-5*x550*x556 + 8.7899999999999996e-10*x550*x557*x569 - 5.5643800000000001e-8*x550*x557 - 2.6524680000000001e-7*x550*x558 + 1.29712e-8*x550*x559 + 0.12116895857649999*x550*x569 + 0.0042852113657000004*x550*x570 - 0.00031917813630000003*x550*x571 - 4.8222189999999996e-6*x550*x572 + 4.3293909999999998e-7*x550*x573 - 6.8428000000000004e-9*x550*x574 - 3.6380000000000003e-10*x550*x575 + 1.2200000000000001e-11*x550*x576 - 0.0001724768049*x550*x583 - 1.85014e-8*x550*x584 - 4.3599999999999999e-10*x550*x585 - 2.1681996300000001e-5*x550*x588 - 7.3379559800000005e-5*x550*x589 - 1.2760957e-6*x550*x590 - 3.5007499999999998e-8*x550*x591 + 2.3160339913389998*x550 - 9.0e-13*std::exp((10)*std::log(std::abs(x551))) + 0.0242280835853*x551*x552*x561 - 0.0073688190838000003*x551*x552*x562 + 0.00093238064600000004*x551*x552*x563 + 2.1883406499999999e-5*x551*x552*x564 - 1.1538743700000001e-5*x551*x552*x565 + 0.0027944244039999999*x551*x553*x560 + 0.00066004885729999997*x551*x553*x562 + 8.2064122999999992e-6*x551*x553*x564 + 0.00015823084750000001*x551*x554*x560 + 6.7963485999999997e-6*x551*x554*x561 - 4.6492789999999997e-6*x551*x554*x562 - 1.2901950000000001e-7*x551*x554*x563 + 0.00082868593670000002*x551*x554 + 1.112603e-6*x551*x555*x561 + 1.7737910000000001e-7*x551*x555*x562 - 5.5623862100000002e-5*x551*x555 + 9.0303320000000002e-7*x551*x556*x560 - 4.6753752e-6*x551*x556 - 1.2193400000000001e-8*x551*x557*x560 + 8.5708430000000003e-7*x551*x557 - 4.4505500000000001e-8*x551*x558 + 8.0859999999999998e-10*x551*x559 - 0.21663723691590001*x551*x560 + 0.1920547245942*x551*x561 + 0.11056734301350001*x551*x562 - 0.24899279055580001*x551*x563 + 0.15351117027309999*x551*x564 - 0.046268182456999998*x551*x565 + 0.0069910235136000001*x551*x566 - 0.00042396606349999998*x551*x567 - 2.02893876e-5*x551*x568 - 5.7396302000000004e-6*x551*x586 + 0.032905353185900002*x551 + 9.5462233999999995e-6*x552*x560*x570 + 2.5328600000000001e-7*x552*x560*x571 + 1.2e-9*x552*x560*x572 - 0.51226978719009997*x552*x560 + 0.00056809192320000001*x552*x561*x569 + 0.62437751308199996*x552*x561 + 1.7199399999999999e-8*x552*x562*x571 - 0.37407109166519997*x552*x562 + 4.1850417199999998e-5*x552*x563*x569 + 2.5775619999999998e-7*x552*x563*x570 + 0.11685657299130001*x552*x563 - 0.018009935602*x552*x564 + 0.00098798154369999992*x552*x565 + 2.0669974100000001e-5*x552*x566 + 1.3096201499999999e-5*x552*x570 + 1.7933759000000001e-6*x552*x571 - 2.9679800000000001e-8*x552*x572 - 4.7500000000000001e-10*x552*x573 + 6.1000000000000003e-12*x552*x574 - 2.9999999999999998e-13*x552*x575 - 2.0799999999999999e-11*x552*x580 - 1.7803739999999999e-7*x552*x593 + 0.063963953049699998*x552 - 1.1496e-9*x553*x560*x571 + 5.2320000000000001e-10*x553*x560*x572 + 0.017505307428100001*x553*x560 + 2.5435716e-6*x553*x561*x569 + 5.5091469999999995e-7*x553*x561*x570 - 1.7527e-9*x553*x561*x571 - 0.024079071690400002*x553*x561 + 0.018071886233*x553*x562 + 3.1931600000000003e-8*x553*x563*x569 - 0.0067334926452*x553*x563 + 0.0011495771008*x553*x564 - 7.1535481699999996e-5*x553*x565 + 9.6953757399999995e-5*x553*x569 + 4.3347996999999997e-6*x553*x570 + 3.7136899999999998e-8*x553*x571 + 1.5729999999999999e-10*x553*x572 - 2.5189999999999998e-10*x553*x573 - 9.0e-13*x553*x574 - 0.00042506434070000001*x553*x582 - 0.0020141106596*x553*x594 - 0.0090247676379999996*x553 + 2.7279999999999998e-10*x554*x560*x571 + 0.00013082126980000001*x554*x560 + 1.08747e-8*x554*x561*x570 - 9.8150547299999998e-5*x554*x561 + 1.7071000000000001e-8*x554*x562*x569 + 0.0002125030906*x554*x562 - 3.7537003600000001e-5*x554*x563 - 3.6726384000000002e-6*x554*x564 - 9.1657204000000005e-6*x554*x569 - 8.2257539999999999e-7*x554*x570 - 7.8846999999999995e-9*x554*x571 + 5.4019999999999995e-10*x554*x572 + 1.64e-11*x554*x573 - 0.00078838830499999998*x554*x582 - 1.14912e-7*x554*x597 - 0.0063418789521999998*x554 + 1.5290449999999999e-7*x555*x560*x569 - 1.9481793299999999e-5*x555*x561 + 1.7275312500000001e-5*x555*x562 + 2.6860743999999999e-6*x555*x563 + 9.0077900000000002e-8*x555*x569 + 3.12127e-8*x555*x570 - 1.026e-10*x555*x571 - 2.6600000000000001e-11*x555*x572 - 8.9480000000000001e-9*x555*x597 + 0.00074280899980000003*x555 - 1.8155e-9*x556*x560*x569 + 4.5213379300000001e-5*x556*x560 - 8.7344791e-6*x556*x561 - 1.3962928e-6*x556*x562 + 3.7880400000000001e-8*x556*x569 + 9.688e-10*x556*x570 + 1.1300000000000001e-11*x556*x571 - 6.0539900000000006e-8*x556*x594 - 2.2511410100000002e-5*x556 - 4.928726e-7*x557*x560 + 5.4556329999999998e-7*x557*x561 - 1.7335000000000001e-9*x557*x569 - 6.4400000000000005e-11*x557*x570 + 1.04782484e-5*x557 - 6.3537600000000001e-8*x558*x560 + 1.7300000000000001e-11*x558*x569 - 1.7361e-9*x558*x582 - 2.1144599000000001e-6*x558 + 1.4765149999999999e-7*x559 - 0.23243336161629999*x560*x569 - 0.0060047590352000002*x560*x570 + 0.00040085503570000001*x560*x571 + 7.2740530000000001e-6*x560*x572 - 1.2737000000000001e-9*x560*x574 - 1.7799999999999999e-11*x560*x575 - 0.041127044175299998*x560*x577 - 0.00070136213949999996*x560*x579 - 7.0284000000000003e-9*x560*x584 - 7.2e-12*x560*x585 - 5.9143245999999997e-6*x560*x589 - 1.1470484e-6*x560*x590 - 7.1189999999999999e-10*x560*x591 - 0.0035571581033999999*x560*x592 + 1.1160026329971999*x560 + 0.2431257553817*x561*x569 + 0.0041151083692999996*x561*x570 - 0.00027830134009999998*x561*x571 - 3.1120835000000001e-6*x561*x572 + 7.7286200000000006e-8*x561*x573 + 3.1390000000000002e-10*x561*x574 - 0.16738748127719999*x561*x587 - 0.55519169578280003*x561 - 0.15397014495149999*x562*x569 - 0.0016821148905000001*x562*x570 + 0.00010047376289999999*x562*x571 + 5.953574e-7*x562*x572 - 9.8590000000000008e-9*x562*x573 - 0.00021914885569999999*x562*x579 - 1.1690754000000001e-6*x562*x588 - 2.37214e-8*x562*x589 - 9.4609900000000006e-8*x562*x590 - 0.0020150238596*x562*x592 - 1.0529487515022*x562 + 0.0612481890093*x563*x569 + 0.0004202885447*x563*x570 - 4.5542200000000003e-8*x563*x572 - 0.0001096246207*x563*x578 - 1.27390263e-5*x563*x583 - 0.039115478906500002*x563*x587 + 0.13827178873939999*x563 - 0.014958490725999999*x564*x569 - 5.9086323300000002e-5*x564*x570 + 1.4636716e-6*x564*x571 - 3.3338918999999998e-6*x564*x579 - 0.0001417281462*x564*x592 + 1.0378186719538001*x564 + 0.0020487853523*x565*x569 + 3.5730467999999999e-6*x565*x570 - 0.0012093420751*x565*x587 - 0.83956762642930005*x565 - 0.0001203960609*x566*x569 + 0.28381726714819999*x566 - 0.045602120043099999*x567 - 0.00034135580390000001*x568 - 2.4955399999999999e-8*x569*x586 - 1.5002941e-6*x569*x595 - 0.030726740117700001*x569 - 4.8511999999999998e-8*x570*x595 - 0.0014092586279*x570 - 9.52185e-8*x571*x596 + 0.00019874633390000001*x571 - 1.8476999999999999e-9*x572*x596 + 4.6423740000000004e-6*x572 - 7.9429120000000001e-7*x573 - 1.0256300000000001e-8*x574 + 1.4564000000000001e-9*x575 + 7.7999999999999999e-12*x576 - 0.0142476868124*x577 - 0.0020219335184000001*x578 - 0.000655172836*x579 - 1.821314e-7*x580 - 1.88829785e-5*x581 - 27.169728828167202)*std::log(std::abs(10)))*X(30)*x448 + std::pow(10, -3.5355e-9*std::exp((10)*std::log(std::abs(x547))) + 0.0019497087522999999*x547*x548*x567 + 2.5033869e-6*x547*x548*x569 + 2.425819e-7*x547*x548*x570 + 2.74e-11*x547*x548*x573 - 0.25171915905909997*x547*x548 + 0.16154758856199999*x547*x549*x558 + 0.10319386725309999*x547*x549*x560 + 0.0091058356367999995*x547*x549*x562 + 7.0547189099999996e-5*x547*x549*x564 + 0.028910314039899999*x547*x549 + 0.00031554761819999999*x547*x558*x568 + 6.9799317000000001e-6*x547*x558*x569 - 8.3654700000000001e-8*x547*x558*x570 - 0.13936594077190001*x547*x558 + 0.0031773255127*x547*x559*x567 - 2.2029417000000002e-6*x547*x559*x569 + 3.32601e-8*x547*x559*x570 + 1.1211e-9*x547*x559*x571 + 0.40330227063039997*x547*x559 + 7.7272750100000005e-5*x547*x560*x568 + 4.3411179999999998e-7*x547*x560*x569 - 2.2911000000000001e-9*x547*x560*x570 + 0.053531392101600003*x547*x560 + 0.00074114505229999996*x547*x561*x567 - 0.4569416534031*x547*x561 + 7.4963150000000001e-7*x547*x562*x568 + 0.3477567512959*x547*x562 + 1.1028151899999999e-5*x547*x563*x567 - 0.11816708728189999*x547*x563 + 0.019397806516399999*x547*x564 - 0.0012520513382999999*x547*x565 + 0.00089752818999999998*x547*x567 + 8.7854088899999995e-5*x547*x568 - 1.36574957e-5*x547*x569 - 9.4252039999999997e-7*x547*x570 + 2.2287e-8*x547*x571 + 2.0822000000000001e-9*x547*x572 - 1.6900000000000001e-11*x547*x573 - 1.1e-12*x547*x574 - 5.0747599999999997e-8*x547*x579 - 0.092472151351600002*x547*x580 - 0.00022388895510000001*x547*x591 + 0.1098450355633*x547 + 0.0028626466126000002*std::exp((10)*std::log(std::abs(x548))) + 0.036943070713000001*x548*x549*x550 + 9.9116502200000001e-5*x548*x549*x553 + 1.5874459999999999e-7*x548*x549*x555 + 0.0208809736142*x548*x549 + 0.00074934921520000003*x548*x550*x567 - 2.1007040000000002e-6*x548*x550*x569 + 2.3848400000000001e-8*x548*x550*x570 + 1.7461e-9*x548*x550*x571 - 9.6999999999999995e-12*x548*x550*x572 + 0.12448499329609999*x548*x550 + 1.706318e-7*x548*x551*x569 - 2.0289999999999999e-9*x548*x551*x570 - 6.3199999999999999e-11*x548*x551*x571 + 0.0021851183530000001*x548*x551 + 9.2263300000000005e-6*x548*x552*x567 + 5.5200690000000001e-7*x548*x552*x568 - 7.8312e-9*x548*x552*x569 + 9.6999999999999995e-12*x548*x552*x570 + 0.0025381599434999998*x548*x552 - 2.020005e-7*x548*x553*x567 + 1.5340000000000001e-10*x548*x553*x569 - 0.0004386554592*x548*x553 + 9.9660000000000005e-10*x548*x554*x568 + 3.4105511399999997e-5*x548*x554 + 8.7899999999999996e-10*x548*x555*x567 - 5.5643800000000001e-8*x548*x555 - 2.6524680000000001e-7*x548*x556 + 1.29712e-8*x548*x557 + 0.12116895857649999*x548*x567 + 0.0042852113657000004*x548*x568 - 0.00031917813630000003*x548*x569 - 4.8222189999999996e-6*x548*x570 + 4.3293909999999998e-7*x548*x571 - 6.8428000000000004e-9*x548*x572 - 3.6380000000000003e-10*x548*x573 + 1.2200000000000001e-11*x548*x574 - 0.0001724768049*x548*x581 - 1.85014e-8*x548*x582 - 4.3599999999999999e-10*x548*x583 - 2.1681996300000001e-5*x548*x586 - 7.3379559800000005e-5*x548*x587 - 1.2760957e-6*x548*x588 - 3.5007499999999998e-8*x548*x589 + 2.3160339913389998*x548 - 9.0e-13*std::exp((10)*std::log(std::abs(x549))) + 0.0242280835853*x549*x550*x559 - 0.0073688190838000003*x549*x550*x560 + 0.00093238064600000004*x549*x550*x561 + 2.1883406499999999e-5*x549*x550*x562 - 1.1538743700000001e-5*x549*x550*x563 + 0.0027944244039999999*x549*x551*x558 + 0.00066004885729999997*x549*x551*x560 + 8.2064122999999992e-6*x549*x551*x562 + 0.00015823084750000001*x549*x552*x558 + 6.7963485999999997e-6*x549*x552*x559 - 4.6492789999999997e-6*x549*x552*x560 - 1.2901950000000001e-7*x549*x552*x561 + 0.00082868593670000002*x549*x552 + 1.112603e-6*x549*x553*x559 + 1.7737910000000001e-7*x549*x553*x560 - 5.5623862100000002e-5*x549*x553 + 9.0303320000000002e-7*x549*x554*x558 - 4.6753752e-6*x549*x554 - 1.2193400000000001e-8*x549*x555*x558 + 8.5708430000000003e-7*x549*x555 - 4.4505500000000001e-8*x549*x556 + 8.0859999999999998e-10*x549*x557 - 0.21663723691590001*x549*x558 + 0.1920547245942*x549*x559 + 0.11056734301350001*x549*x560 - 0.24899279055580001*x549*x561 + 0.15351117027309999*x549*x562 - 0.046268182456999998*x549*x563 + 0.0069910235136000001*x549*x564 - 0.00042396606349999998*x549*x565 - 2.02893876e-5*x549*x566 - 5.7396302000000004e-6*x549*x584 + 0.032905353185900002*x549 + 9.5462233999999995e-6*x550*x558*x568 + 2.5328600000000001e-7*x550*x558*x569 + 1.2e-9*x550*x558*x570 - 0.51226978719009997*x550*x558 + 0.00056809192320000001*x550*x559*x567 + 0.62437751308199996*x550*x559 + 1.7199399999999999e-8*x550*x560*x569 - 0.37407109166519997*x550*x560 + 4.1850417199999998e-5*x550*x561*x567 + 2.5775619999999998e-7*x550*x561*x568 + 0.11685657299130001*x550*x561 - 0.018009935602*x550*x562 + 0.00098798154369999992*x550*x563 + 2.0669974100000001e-5*x550*x564 + 1.3096201499999999e-5*x550*x568 + 1.7933759000000001e-6*x550*x569 - 2.9679800000000001e-8*x550*x570 - 4.7500000000000001e-10*x550*x571 + 6.1000000000000003e-12*x550*x572 - 2.9999999999999998e-13*x550*x573 - 2.0799999999999999e-11*x550*x578 - 1.7803739999999999e-7*x550*x591 + 0.063963953049699998*x550 - 1.1496e-9*x551*x558*x569 + 5.2320000000000001e-10*x551*x558*x570 + 0.017505307428100001*x551*x558 + 2.5435716e-6*x551*x559*x567 + 5.5091469999999995e-7*x551*x559*x568 - 1.7527e-9*x551*x559*x569 - 0.024079071690400002*x551*x559 + 0.018071886233*x551*x560 + 3.1931600000000003e-8*x551*x561*x567 - 0.0067334926452*x551*x561 + 0.0011495771008*x551*x562 - 7.1535481699999996e-5*x551*x563 + 9.6953757399999995e-5*x551*x567 + 4.3347996999999997e-6*x551*x568 + 3.7136899999999998e-8*x551*x569 + 1.5729999999999999e-10*x551*x570 - 2.5189999999999998e-10*x551*x571 - 9.0e-13*x551*x572 - 0.00042506434070000001*x551*x580 - 0.0020141106596*x551*x592 - 0.0090247676379999996*x551 + 2.7279999999999998e-10*x552*x558*x569 + 0.00013082126980000001*x552*x558 + 1.08747e-8*x552*x559*x568 - 9.8150547299999998e-5*x552*x559 + 1.7071000000000001e-8*x552*x560*x567 + 0.0002125030906*x552*x560 - 3.7537003600000001e-5*x552*x561 - 3.6726384000000002e-6*x552*x562 - 9.1657204000000005e-6*x552*x567 - 8.2257539999999999e-7*x552*x568 - 7.8846999999999995e-9*x552*x569 + 5.4019999999999995e-10*x552*x570 + 1.64e-11*x552*x571 - 0.00078838830499999998*x552*x580 - 1.14912e-7*x552*x595 - 0.0063418789521999998*x552 + 1.5290449999999999e-7*x553*x558*x567 - 1.9481793299999999e-5*x553*x559 + 1.7275312500000001e-5*x553*x560 + 2.6860743999999999e-6*x553*x561 + 9.0077900000000002e-8*x553*x567 + 3.12127e-8*x553*x568 - 1.026e-10*x553*x569 - 2.6600000000000001e-11*x553*x570 - 8.9480000000000001e-9*x553*x595 + 0.00074280899980000003*x553 - 1.8155e-9*x554*x558*x567 + 4.5213379300000001e-5*x554*x558 - 8.7344791e-6*x554*x559 - 1.3962928e-6*x554*x560 + 3.7880400000000001e-8*x554*x567 + 9.688e-10*x554*x568 + 1.1300000000000001e-11*x554*x569 - 6.0539900000000006e-8*x554*x592 - 2.2511410100000002e-5*x554 - 4.928726e-7*x555*x558 + 5.4556329999999998e-7*x555*x559 - 1.7335000000000001e-9*x555*x567 - 6.4400000000000005e-11*x555*x568 + 1.04782484e-5*x555 - 6.3537600000000001e-8*x556*x558 + 1.7300000000000001e-11*x556*x567 - 1.7361e-9*x556*x580 - 2.1144599000000001e-6*x556 + 1.4765149999999999e-7*x557 - 0.23243336161629999*x558*x567 - 0.0060047590352000002*x558*x568 + 0.00040085503570000001*x558*x569 + 7.2740530000000001e-6*x558*x570 - 1.2737000000000001e-9*x558*x572 - 1.7799999999999999e-11*x558*x573 - 0.041127044175299998*x558*x575 - 0.00070136213949999996*x558*x577 - 7.0284000000000003e-9*x558*x582 - 7.2e-12*x558*x583 - 5.9143245999999997e-6*x558*x587 - 1.1470484e-6*x558*x588 - 7.1189999999999999e-10*x558*x589 - 0.0035571581033999999*x558*x590 + 1.1160026329971999*x558 + 0.2431257553817*x559*x567 + 0.0041151083692999996*x559*x568 - 0.00027830134009999998*x559*x569 - 3.1120835000000001e-6*x559*x570 + 7.7286200000000006e-8*x559*x571 + 3.1390000000000002e-10*x559*x572 - 0.16738748127719999*x559*x585 - 0.55519169578280003*x559 - 0.15397014495149999*x560*x567 - 0.0016821148905000001*x560*x568 + 0.00010047376289999999*x560*x569 + 5.953574e-7*x560*x570 - 9.8590000000000008e-9*x560*x571 - 0.00021914885569999999*x560*x577 - 1.1690754000000001e-6*x560*x586 - 2.37214e-8*x560*x587 - 9.4609900000000006e-8*x560*x588 - 0.0020150238596*x560*x590 - 1.0529487515022*x560 + 0.0612481890093*x561*x567 + 0.0004202885447*x561*x568 - 4.5542200000000003e-8*x561*x570 - 0.0001096246207*x561*x576 - 1.27390263e-5*x561*x581 - 0.039115478906500002*x561*x585 + 0.13827178873939999*x561 - 0.014958490725999999*x562*x567 - 5.9086323300000002e-5*x562*x568 + 1.4636716e-6*x562*x569 - 3.3338918999999998e-6*x562*x577 - 0.0001417281462*x562*x590 + 1.0378186719538001*x562 + 0.0020487853523*x563*x567 + 3.5730467999999999e-6*x563*x568 - 0.0012093420751*x563*x585 - 0.83956762642930005*x563 - 0.0001203960609*x564*x567 + 0.28381726714819999*x564 - 0.045602120043099999*x565 - 0.00034135580390000001*x566 - 2.4955399999999999e-8*x567*x584 - 1.5002941e-6*x567*x593 - 0.030726740117700001*x567 - 4.8511999999999998e-8*x568*x593 - 0.0014092586279*x568 - 9.52185e-8*x569*x594 + 0.00019874633390000001*x569 - 1.8476999999999999e-9*x570*x594 + 4.6423740000000004e-6*x570 - 7.9429120000000001e-7*x571 - 1.0256300000000001e-8*x572 + 1.4564000000000001e-9*x573 + 7.7999999999999999e-12*x574 - 0.0142476868124*x575 - 0.0020219335184000001*x576 - 0.000655172836*x577 - 1.821314e-7*x578 - 1.88829785e-5*x579 - 27.169728828167202)*X(30)*x448 )) + ((x368) ? ( dustSemenov_cooling ) @@ -4009,6 +4005,183 @@ Real rhs_eint(const burn_t& state, )) + 0.00084373771595996178*x4*(1.3806479999999999e-16*X(0) + 1.3806479999999999e-16*X(1) + 1.3806479999999999e-16*X(10) + 1.3806479999999999e-16*X(11) + 1.3806479999999999e-16*X(12) + 1.3806479999999999e-16*X(13) + 1.3806479999999999e-16*X(14) + 1.3806479999999999e-16*X(15) + 1.3806479999999999e-16*X(16) + 1.3806479999999999e-16*X(17) + 1.3806479999999999e-16*X(18) + 1.3806479999999999e-16*X(19) + 1.3806479999999999e-16*X(2) + 1.3806479999999999e-16*X(20) + 1.3806479999999999e-16*X(21) + 1.3806479999999999e-16*X(22) + 1.3806479999999999e-16*X(23) + 1.3806479999999999e-16*X(24) + 1.3806479999999999e-16*X(25) + 1.3806479999999999e-16*X(26) + 1.3806479999999999e-16*X(27) + 1.3806479999999999e-16*X(28) + 1.3806479999999999e-16*X(29) + 1.3806479999999999e-16*X(3) + 1.3806479999999999e-16*X(30) + 1.3806479999999999e-16*X(31) + 1.3806479999999999e-16*X(32) + 1.3806479999999999e-16*X(33) + 1.3806479999999999e-16*X(4) + 1.3806479999999999e-16*X(5) + 1.3806479999999999e-16*X(6) + 1.3806479999999999e-16*X(7) + 1.3806479999999999e-16*X(8) + 1.3806479999999999e-16*X(9))/(std::sqrt(x1)*x35))); } +/* +void init_anytab2D(const std::string& filename, Array1D& x, Array1D& y, Array2D& z) { + + // Check if the sizes of x and z(1) match + // Check if the sizes of y and z(2) match + + std::cout << "Reading tables from " << filename << std::endl; + + // Open the file and check if it exists + std::ifstream file(filename); + if (!file.is_open()) { + throw std::runtime_error("ERROR: file " + filename + " not found!"); + } + + std::string line; + + // Skip comments and read until the first non-comment line + while (std::getline(file, line)) { + if (line[0] != '#') { + break; + } + } + + // Check if the first line is formatted correctly + if (line.find(',') == std::string::npos) { + throw std::runtime_error("ERROR: file " + filename + " should contain the number of rows and columns in the format 'RR, CC'"); + } + + // Process data line by line + for (int i = 0; i < x.size(); ++i) { + for (int j = 0; j < y.size(); ++j) { + if (!std::getline(file, line)) { + throw std::runtime_error("ERROR: Unexpected end of file while reading data."); + } + + std::istringstream iss(line); + Real x_val, y_val, z_val; + if (!(iss >> x_val >> y_val >> z_val)) { + throw std::runtime_error("ERROR: Insufficient data on line " + std::to_string(i * y.size() + j + 1)); + } + + x(i) = x_val; + y(j) = y_val; + z(i,j) = z_val; + } + + // Skip any remaining blanks in the current line + if (!std::getline(file, line)) { + throw std::runtime_error("ERROR: Unexpected end of file while skipping blanks."); + } + } +} +*/ + +AMREX_GPU_HOST_DEVICE AMREX_INLINE +std::pair compute_Semenov_Tdust(Real Tgas, const Array1D& composition, Real const user_dust2gas_ratio, Real const z, + Real krome_Semenov_Tdust, const Array1D& semenov_x, const Array1D& semenov_y, + const Array2D& semenov_z) { + + Real dustSemenov_cooling = 0.0; + + if (user_dust2gas_ratio == 0) { + return std::make_pair(krome_Semenov_Tdust, dustSemenov_cooling); + } + + Tgas = std::max(Tgas, 2.73 * (1.0 + z)); + + Array1D mass; + mass(0) = 0.0; + mass(1) = 0.0; + mass(2) = 9.10938188e-28; + mass(3) = 1.67262158e-24; + mass(4) = 1.673532518188e-24; + mass(5) = 1.6744434563759998e-24; + mass(6) = 3.34512158e-24; + mass(7) = 3.346032518188e-24; + mass(8) = 3.346154098188e-24; + mass(9) = 3.346943456376e-24; + mass(10) = 3.347065036376e-24; + mass(11) = 5.0186540981880006e-24; + mass(12) = 5.0195650363760004e-24; + mass(13) = 6.690243159999999e-24; + mass(14) = 6.691154098188e-24; + mass(15) = 6.692065036376e-24; + mass(16) = 2.007528417094e-23; + mass(17) = 2.0076195109128e-23; + mass(18) = 2.1749727627316e-23; + mass(19) = 2.3423260145503998e-23; + mass(20) = 2.5096792663692e-23; + mass(21) = 2.6767349207316e-23; + mass(22) = 2.6768260145504e-23; + mass(23) = 2.677032518188e-23; + mass(24) = 2.8440881725504e-23; + mass(25) = 2.8441792663692003e-23; + mass(26) = 3.0114414243692e-23; + mass(27) = 3.011532518188e-23; + mass(28) = 3.178794676188e-23; + mass(29) = 4.6843544316444e-23; + mass(30) = 4.6844455254632e-23; + mass(31) = 5.353560935282e-23; + mass(32) = 5.3536520291008e-23; + mass(33) = 7.3612715400136e-23; + + Real rhogas = 0.0; + Real sum_comp = 0.0; + for (int i = 0; i < NumSpec-1; ++i) { + rhogas += composition(i) * mass(i); + sum_comp += composition(i); + } + rhogas = std::max(rhogas, 1e-40); + + Real nH = 0.0; + nH = composition(3) + composition(4) + composition(5) + composition(10) * 2.0 + composition(8) * 2.0 + + composition(25) + composition(18) + composition(19) * 2.0 + composition(20) * 3.0 + + composition(27) * 2.0 + composition(23) * 4.0 + composition(24) + composition(26) * 2.0 + + composition(28) * 3.0 + composition(11) + composition(12); + + //std::cout << "compute semenov tdust nH: " << nH << std::endl; + + Real clipped_x = std::max(-18.0, std::min(std::log10(rhogas), -7.0)); + Real clipped_y = std::max(1e1, std::min(Tgas, 1e4)); + + //Real kappaP = interpolate_2d(semenov_x, semenov_y, semenov_z, clipped_x, clipped_y); + Real kappaP = std::pow((-1.58418734e+00*std::pow(std::log10(clipped_y),7) + 2.70026428e+01*std::pow(std::log10(clipped_y),6) - 1.90914865e+02*std::pow(std::log10(clipped_y),5) + \ + 7.24661740e+02*std::pow(std::log10(clipped_y),4) - 1.59299671e+03*std::pow(std::log10(clipped_y),3) + 2.02540166e+03*std::pow(std::log10(clipped_y),2) - \ + 1.37581206e+03*std::pow(std::log10(clipped_y),1) + 3.83747773e+02), 10); + + //std::cout << "compute semenov tdust kappaP: " << kappaP << std::endl; + + Real mu = rhogas / std::max(sum_comp, 1e-40) * (1.0/C::m_p); + Real ljeans = std::sqrt(M_PI * C::k_B * Tgas / rhogas / C::m_p / C::Gconst / mu); + + Real tau_d = rhogas * kappaP * user_dust2gas_ratio; + Real tau_g = 0.0; + Real tau = (tau_d + tau_g) * ljeans; + + Real besc = (tau < 1.0) ? 1.0 : std::pow(tau, -2); + + Real const alpha_gd = 3.2e-34; + Real const aR = 4 * C::sigma_SB / C::c_light; + + kappaP *= besc; + + Real intJRad = 0.0; + + Real A = rhogas * kappaP * user_dust2gas_ratio * aR * C::c_light; + Real B = std::pow(Tgas, 0.5) * std::pow(nH, 2) * alpha_gd * user_dust2gas_ratio; + Real phys_Tcmb = 2.73 * (1.0 + z); + Real C = -1.0 * (intJRad + A * std::pow(phys_Tcmb, 4) + std::pow(Tgas, 1.5) * std::pow(nH, 2) * alpha_gd * user_dust2gas_ratio); + + int iterr = 0; + Real Tdold = krome_Semenov_Tdust; + Real Tdnew = 0.0; + + while (true) { + Real fx = A * std::pow(Tdold, 4) + B * Tdold + C; + Real fdash_x = 4.0 * A * std::pow(Tdold, 3) + B; + + Tdnew = Tdold - fx / fdash_x; + Real rel_t = std::abs((Tdnew - Tdold) / Tdold); + Real abs_t = std::abs(Tdnew - Tdold); + Tdold = Tdnew; + ++iterr; + + if (abs_t < 1e-8 || rel_t < 1e-5) { + break; + } + + if (iterr > 1000) { + throw std::runtime_error("Maximum iterations reached in dust temperature NR-solver"); + } + } + + dustSemenov_cooling = A * std::pow(Tdnew, 4) - intJRad - A * std::pow(phys_Tcmb, 4); + + return std::make_pair(Tdnew, dustSemenov_cooling); +} AMREX_GPU_HOST_DEVICE AMREX_INLINE void actual_rhs (burn_t& state, Array1D& ydot) @@ -4021,9 +4194,19 @@ void actual_rhs (burn_t& state, Array1D& ydot) X(i) = state.xn[i]; } + Real krome_Semenov_Tdust = state.aux[0]; + Array1D semenov_x; + Array1D semenov_y; + Array2D semenov_z; + std::string filename = "Semenov_PlanckOpacity.dat"; + //init_anytab2D(filename, semenov_x, semenov_y, semenov_z); + // Call the compute_Semenov_Tdust function + auto result = compute_Semenov_Tdust(state.T, X, dust2gas_ratio, z, krome_Semenov_Tdust, semenov_x, semenov_y, semenov_z); // compute Tdust - Real Tdust = 2.73*(1.0 + z); - Real dustSemenov_cooling = 0.0; + Real Tdust = result.first; + // store Tdust in aux + state.aux[0] = result.first; + Real dustSemenov_cooling = result.second; // YDOTS @@ -4039,8 +4222,25 @@ void actual_rhs (burn_t& state, Array1D& ydot) ydot(net_ienuc) = edot; -} + /* + std::cout << "redshift " << redshift << ", " << metallicity << ", " << dust2gas_ratio << ", " << krome_Semenov_Tdust << ", " << dustSemenov_cooling << std::endl; + std::cout << "network X: "; + for (int i = 0; i < NumSpec; ++i) { + std::cout << "X" << i << ": " << X(i) << " "; + } + std::cout << std::endl; + std::cout << "Tgas: " << state.T << std::endl; + std::cout << "network Tdust: " << state.aux[0] << std::endl; + std::cout << "network ydots: "; + for (int i = 1; i < neqs; ++i) { + std::cout << "ydot" << i << ": " << ydot(i) << " "; + } + std::cout << std::endl; + std::cout << "network edot: " << edot << std::endl; + std::cout << " " << std::endl; + */ +} template