From 4812f9af7c790d2c0dee3f5f1601f58a1e1f8097 Mon Sep 17 00:00:00 2001 From: Piyush Sharda Date: Tue, 17 Sep 2024 17:21:14 +0200 Subject: [PATCH 1/7] handle ices correctly --- EOS/metal_chem/_parameters | 4 +- networks/metal_chem/actual_rhs.H | 978 +++++++++++++++---------------- 2 files changed, 491 insertions(+), 491 deletions(-) diff --git a/EOS/metal_chem/_parameters b/EOS/metal_chem/_parameters index 8fa847fd8..238fad957 100644 --- a/EOS/metal_chem/_parameters +++ b/EOS/metal_chem/_parameters @@ -3,11 +3,11 @@ eos_gamma_default real 1.4 # define the specie names, and their masses and gammas -species_1_name string "co_ice" +species_1_name string "co_total" species_1_gamma real 5./3. species_1_mass real 0.0 -species_2_name string "h2o_ice" +species_2_name string "h2o_total" species_2_gamma real 5./3. species_2_mass real 0.0 diff --git a/networks/metal_chem/actual_rhs.H b/networks/metal_chem/actual_rhs.H index 3e96e81cd..bf64a1e17 100644 --- a/networks/metal_chem/actual_rhs.H +++ b/networks/metal_chem/actual_rhs.H @@ -32,968 +32,968 @@ void rhs_specie(const burn_t& state, Real T = state.T; - Real x0 = 2.0e-99*X(0); + Real x0 = std::max(T, 2.73*z + 2.73); - Real x1 = std::max(T, 2.73*z + 2.73); + Real x1 = ((x0)*(x0)); - Real x2 = ((x1)*(x1)); + Real x2 = 0.002*x0 + 1.0; - Real x3 = 0.002*x1 + 1.0; + Real x3 = std::sqrt(x0); - Real x4 = std::sqrt(x1); + Real x4 = user_dust2gas_ratio*x3; - Real x5 = user_dust2gas_ratio*x4; + Real x5 = x4*(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))/(8.0000000000000013e-6*x1 + x2 + 0.040000000000000008*std::sqrt(x0 + 1)); - Real x6 = x5*(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))/(8.0000000000000013e-6*x2 + x3 + 0.040000000000000008*std::sqrt(x1 + 1)); + Real x6 = 1195783.1580871048*x5 + 2.0e-99; - Real x7 = 1195783.1580871048*x6 + 2.0e-99; + Real x7 = X(16)*X(32); - Real x8 = 2.0e-99*X(1); + Real x8 = 4.5299999999999999e-10*x7; - Real x9 = 1491378.0618531078*x6 + 2.0e-99; + Real x9 = std::sqrt(x0); - Real x10 = 1.0/x4; + Real x10 = 1.0/x0; - 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 x11 = 1.6974097914175e-15*x9*std::exp(-7550.0*x10) + 1.0e-99; - Real x12 = X(13)*X(2); + Real x12 = X(10)*X(33); - Real x13 = -x11*x12; + Real x13 = x11*x12; - Real x14 = std::pow(x1, -0.69999999999999996); + Real x14 = 7.5e-10*X(29)*X(4); - Real x15 = X(2)*X(31); + Real x15 = 1.0e-99 + 2.8e-11*std::exp(-26500.0*x10); - Real x16 = x15*(-1.0568649366703787e-5*x14 - 1.0e-99); + Real x16 = X(22)*X(33); - Real x17 = 7.1999999999999996e-8*x10 + 1.0e-99; + Real x17 = 1.0e-99 + 4.1999999999999999e-12*std::exp(-24000.0*x10); - Real x18 = X(11)*X(2); + Real x18 = X(30)*X(32); - Real x19 = -x17*x18; + Real x19 = x15*x16 - x17*x18; - Real x20 = 9.7355109978586414e-8*std::pow(x1, -0.47999999999999998) + 1.0e-99; + Real x20 = 1.0e-13*X(25)*X(30); - Real x21 = X(2)*X(24); + Real x21 = 1.0e-99 + 2.54e-10*std::exp(-13300.0*x10); - Real x22 = -x20*x21; + Real x22 = X(33)*X(4); - Real x23 = std::log(x1); + Real x23 = -x20 + x21*x22; - Real x24 = 0.10684732509875319*x23 - 1; + Real x24 = std::pow(x0, -0.34000000000000002); - Real x25 = ((x24)*(x24)); + Real x25 = ((x0 <= 1052.0) ? ( + 3.2682787072746432e-10*x24 +) +: ( + 0 +)) + 1.0e-99; - Real x26 = ((x24)*(x24)*(x24)); + Real x26 = X(17)*X(32); - Real x27 = ((((x24)*(x24)))*(((x24)*(x24)))); + Real x27 = ((x0 > 1052.0) ? ( + 3.7991369449339675e-16*std::pow(x0, 1.54)*std::exp(613.0*x10) +) +: ( + 0 +)) + 1.0e-99; - Real x28 = ((x24)*(x24)*(x24)*(x24)*(x24)); + Real x28 = x25*x26 + x26*x27; - Real x29 = std::exp((6)*std::log(std::abs(x24))); + Real x29 = x0 <= 10.0; - Real x30 = ((x24)*(x24)*(x24)*(x24)*(x24)*(x24)*(x24)); + Real x30 = ((x29) ? ( + 7.0509999999999997e-11 +) +: ( + 0 +)) + 1.0e-99; - Real x31 = std::exp((8)*std::log(std::abs(x24))); + Real x31 = X(17)*X(25); - 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 x32 = x0 > 10.0; - Real x33 = 4.3524079114767552e-117*x32 + 1.0e-99; + Real x33 = ((x32) ? ( + 1.5557027645121538e-10*std::pow(x0, -0.33900000000000002)*std::exp(-0.108*x10) +) +: ( + 0 +)) + 1.0e-99; - Real x34 = X(15)*X(2); + Real x34 = 1.0e-99 + 1.02e-10*std::exp(-914.0*x10); - Real x35 = -x33*x34; + Real x35 = X(18)*X(22); - Real x36 = std::pow(x1, -0.75); + Real x36 = x30*x31 + x31*x33 + x34*x35; - Real x37 = 2.5950363272655348e-10*x36 + 1.0e-99; + Real x37 = -X(30)*x6 + x13 + x14 + x19 + x23 + x28 + x36 + x8; - Real x38 = 2*X(2)*X(6); + Real x38 = 1491378.0618531078*x5 + 2.0e-99; - Real x39 = -x37*x38; + Real x39 = ((X(25))*(X(25))); - Real x40 = 6.3350403285625223e-6*std::pow(x1, -0.55000000000000004) + 1.0e-99; + Real x40 = std::pow(x0, 1.1399999999999999)*std::exp(-50.0*x10); - Real x41 = X(2)*X(29); + Real x41 = x39*(2.4749459212614312e-15*x40 + 1.0e-99); - Real x42 = -x40*x41; + Real x42 = ((x29) ? ( + 4.2020000000000002e-8 +) +: ( + 0 +)) + 1.0e-99; - 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 x43 = X(27)*X(3); - Real x44 = 3.7903999274394518e-18*x43 + 1.0e-99; + Real x44 = 1.0/x9; - Real x45 = X(2)*X(5); + Real x45 = ((x32) ? ( + 1.2730571616977901e-7*x44 + 1.767e-9 +) +: ( + 0 +)) + 1.0e-99; - Real x46 = -x44*x45; + Real x46 = -x42*x43 - x43*x45; - Real x47 = 1.3300135414628029e-18*std::pow(x1, 0.94999999999999996)*std::exp(-0.00010729613733905579*x1) + 1.0e-99; + Real x47 = 3.5999999999999998e-16*std::pow(T, 1.52)*std::exp(-1740.0*x10) + 1.0e-99; - Real x48 = X(2)*X(7); + Real x48 = X(10)*X(25); - Real x49 = -x47*x48; + Real x49 = 1.8706148721743872e-6*x44 + 1.0e-99; - Real x50 = std::pow(x1, 0.25); + Real x50 = X(2)*X(28); - Real x51 = 1.0/x50; + Real x51 = 6.3600000000000002e-14*x0*std::exp(-9610.0*x10) + 1.0e-99; - Real x52 = 1.0/x1; + Real x52 = X(27)*X(4); - Real x53 = 130666.72639118577*x52; + Real x53 = x47*x48 + x49*x50 - x51*x52; - Real x54 = x51*std::exp(-x53)/(x53 + 0.193); + Real x54 = -X(27)*x38 + x13 + x41 + x46 + x53; - Real x55 = 1.3023623244766063e-6*x54 + 1.0e-99; + Real x55 = 1.0/x3; - Real x56 = X(17)*X(2); + Real x56 = 5.7884371785482823e-10*x55*std::pow(0.00060040841663220993*x3 + 1.0, -1.7524)*std::pow(0.32668576019240059*x3 + 1.0, -0.24759999999999999) + 1.0e-99; - Real x57 = -x55*x56; + Real x57 = X(13)*X(2); - Real x58 = std::pow(x1, -0.34000000000000002); + Real x58 = -x56*x57; - Real x59 = 157821.26811013557*x52; + Real x59 = std::pow(x0, -0.69999999999999996); - Real x60 = x58*std::exp(-x59)/(x59 + 0.072999999999999995); + Real x60 = X(2)*X(31); - Real x61 = 2.1012161411986507e-6*x60 + 1.0e-99; + Real x61 = x60*(-1.0568649366703787e-5*x59 - 1.0e-99); - Real x62 = X(2)*X(22); + Real x62 = 7.1999999999999996e-8*x55 + 1.0e-99; - Real x63 = -x61*x62; + Real x63 = X(11)*X(2); - Real x64 = std::pow(x1, 0.34999999999999998)*std::exp(-102000.0*x52); + Real x64 = -x62*x63; - Real x65 = 4.3799999999999999e-10*x64 + 1.0e-99; + Real x65 = 9.7355109978586414e-8*std::pow(x0, -0.47999999999999998) + 1.0e-99; - Real x66 = X(10)*X(2); + Real x66 = X(2)*X(24); - 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 x67 = -x65*x66; - Real x68 = X(14)*X(2); + Real x68 = std::log(x0); - 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 x69 = 0.10684732509875319*x68 - 1; - Real x70 = X(2)*X(4); + Real x70 = ((x69)*(x69)); - Real x71 = ((x1 <= 1160.0) ? ( - 1.4643482606109061e-16*std::pow(x1, 1.78186) -) -: ( - 0 -)); + Real x71 = ((x69)*(x69)*(x69)); + + Real x72 = ((((x69)*(x69)))*(((x69)*(x69)))); + + Real x73 = ((x69)*(x69)*(x69)*(x69)*(x69)); + + Real x74 = std::exp((6)*std::log(std::abs(x69))); + + Real x75 = ((x69)*(x69)*(x69)*(x69)*(x69)*(x69)*(x69)); + + Real x76 = std::exp((8)*std::log(std::abs(x69))); + + Real x77 = std::pow(x0, 23.915965629999999)*std::exp(-941.91483008144996*x70 + 2506.9866529060901*x71 - 4361.9927099007555*x72 + 4879.7345146260486*x73 - 3366.4639698826941*x74 + 1300.3028484326148*x75 - 214.82451513312137*x76); + + Real x78 = 4.3524079114767552e-117*x77 + 1.0e-99; + + Real x79 = X(15)*X(2); + + Real x80 = -x78*x79; + + Real x81 = std::pow(x0, -0.75); + + Real x82 = 2.5950363272655348e-10*x81 + 1.0e-99; + + Real x83 = 2*X(2)*X(6); + + Real x84 = -x82*x83; + + Real x85 = 6.3350403285625223e-6*std::pow(x0, -0.55000000000000004) + 1.0e-99; + + Real x86 = X(2)*X(29); + + Real x87 = -x85*x86; + + Real x88 = std::pow(x0, 2.360852208681)*std::exp(-24.766609674457612*x70 + 13.307984239358756*x71 - 258.18559308467116*x72 + 846.15238706523724*x73 - 1113.0879095147111*x74 + 671.95094388835207*x75 - 154.90262957142161*x76); + + Real x89 = 3.7903999274394518e-18*x88 + 1.0e-99; + + Real x90 = X(2)*X(5); + + Real x91 = -x89*x90; + + Real x92 = 1.3300135414628029e-18*std::pow(x0, 0.94999999999999996)*std::exp(-0.00010729613733905579*x0) + 1.0e-99; + + Real x93 = X(2)*X(7); + + Real x94 = -x92*x93; + + Real x95 = std::pow(x0, 0.25); + + Real x96 = 1.0/x95; + + Real x97 = 130666.72639118577*x10; + + Real x98 = x96*std::exp(-x97)/(x97 + 0.193); + + Real x99 = 1.3023623244766063e-6*x98 + 1.0e-99; + + Real x100 = X(17)*X(2); + + Real x101 = -x100*x99; + + Real x102 = 157821.26811013557*x10; + + Real x103 = x24*std::exp(-x102)/(x102 + 0.072999999999999995); + + Real x104 = 2.1012161411986507e-6*x103 + 1.0e-99; + + Real x105 = X(2)*X(22); + + Real x106 = -x104*x105; + + Real x107 = std::pow(x0, 0.34999999999999998)*std::exp(-102000.0*x10); + + Real x108 = 4.3799999999999999e-10*x107 + 1.0e-99; + + Real x109 = X(10)*X(2); + + Real x110 = std::pow(x0, 43.933476326349997)*std::exp(-1618.789587733125*x70 + 3854.4033653120221*x71 - 5902.1601240760483*x72 + 5825.9326359379538*x73 - 3578.1439181805955*x74 + 1242.7294446825149*x75 - 186.35635455381879*x76); - Real x72 = x71 + 1.0e-99; + Real x111 = X(14)*X(2); - Real x73 = X(4)*X(5); + Real x112 = std::pow(x0, 13.536555999999999)*std::exp(-502.72883252679094*x70 + 1281.477767828706*x71 - 2207.4643501257692*x72 + 2500.8077583366976*x73 - 1768.8867461266502*x74 + 704.19926629500367*x75 - 120.0438480494693*x76); - Real x74 = ((x24)*(x24)*(x24)*(x24)*(x24)*(x24)*(x24)*(x24)*(x24)); + Real x113 = X(2)*X(4); - Real x75 = ((x1 > 1160.0) ? ( - 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) + Real x114 = ((x0 <= 1160.0) ? ( + 1.4643482606109061e-16*std::pow(x0, 1.78186) ) : ( 0 )); - Real x76 = x75 + 1.0e-99; + Real x115 = x114 + 1.0e-99; - Real x77 = std::sqrt(x1); + Real x116 = X(4)*X(5); - Real x78 = 1.0/x77; + Real x117 = ((x69)*(x69)*(x69)*(x69)*(x69)*(x69)*(x69)*(x69)*(x69)); - Real x79 = 4.4686910835277033e-6*x78 + 1.0e-99; + Real x118 = ((x0 > 1160.0) ? ( + 3.3178155742407601e-14*std::pow(x0, 1.139449335841631)*std::exp(-44.454280878123605*x117 - 12.447178055372776*x70 + 6.9391784778399117*x71 - 10.993097527150175*x72 + 14.449862906216714*x73 + 58.228375789703179*x74 - 162.59852239006702*x75 + 144.55426734953477*x76) +) +: ( + 0 +)); - Real x80 = X(2)*X(28); + Real x119 = x118 + 1.0e-99; - Real x81 = 1.8706148721743872e-6*x78 + 1.0e-99; + Real x120 = 4.4686910835277033e-6*x44 + 1.0e-99; - Real x82 = -x79*x80 - x80*x81; + Real x121 = -x120*x50 - x49*x50; - Real x83 = 1.4895636945092344e-6*x78 + 1.0e-99; + Real x122 = 1.4895636945092344e-6*x44 + 1.0e-99; - Real x84 = X(2)*X(26); + Real x123 = X(2)*X(26); - Real x85 = 6.7549981495186206e-7*x78 + 1.0e-99; + Real x124 = 6.7549981495186206e-7*x44 + 1.0e-99; - Real x86 = -x83*x84 - x84*x85; + Real x125 = -x122*x123 - x123*x124; - Real x87 = x1 <= 10000.0; + Real x126 = x0 <= 10000.0; - Real x88 = ((x87) ? ( - -5.5279999999999998e-28*((x1)*(x1)*(x1)*(x1)*(x1)) + 3.3467999999999999e-23*((((x1)*(x1)))*(((x1)*(x1)))) - 7.5474000000000004e-19*((x1)*(x1)*(x1)) - 2.3088e-11*x1 + 7.3427999999999993e-15*x2 + 4.2277999999999996e-8 + Real x127 = ((x126) ? ( + -5.5279999999999998e-28*((x0)*(x0)*(x0)*(x0)*(x0)) + 3.3467999999999999e-23*((((x0)*(x0)))*(((x0)*(x0)))) - 7.5474000000000004e-19*((x0)*(x0)*(x0)) - 2.3088e-11*x0 + 7.3427999999999993e-15*x1 + 4.2277999999999996e-8 ) : ( 0 )); - Real x89 = X(2)*X(8); + Real x128 = X(2)*X(8); - Real x90 = 1.0e-8*std::pow(x1, -0.40000000000000002) + 1.0e-99; + Real x129 = 1.0e-8*std::pow(x0, -0.40000000000000002) + 1.0e-99; - Real x91 = X(3)*X(5); + Real x130 = X(3)*X(5); - Real x92 = x89*(-x88 - 1.0e-99) + x90*x91; + Real x131 = x128*(-x127 - 1.0e-99) + x129*x130; - Real x93 = ((x1 <= 400.0) ? ( - 1.2999999999999999e-10*std::pow(x1, -0.64000000000000001) + Real x132 = ((x0 <= 400.0) ? ( + 1.2999999999999999e-10*std::pow(x0, -0.64000000000000001) ) : ( 0 )) + 1.0e-99; - Real x94 = X(2)*X(21); + Real x133 = X(2)*X(21); - Real x95 = std::pow(x1, -1.5); + Real x134 = std::pow(x0, -1.5); - Real x96 = ((x1 > 400.0) ? ( - 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) + Real x135 = ((x0 > 400.0) ? ( + 1.41e-10*std::pow(x0, -0.66000000000000003) + 0.00073999999999999999*x134*(1.0 + 0.062*std::exp(-145000.0*x10))*std::exp(-175000.0*x10) ) : ( 0 )) + 1.0e-99; - Real x97 = -x93*x94 - x94*x96; + Real x136 = -x132*x133 - x133*x135; - Real x98 = 1.4981088130721367e-10*std::pow(x1, -0.63529999999999998); + Real x137 = 1.4981088130721367e-10*std::pow(x0, -0.63529999999999998); - Real x99 = x1 > 2.73; + Real x138 = x0 > 2.73; - Real x100 = ((x99 && x1 <= 9280.0) ? ( - x98 + Real x139 = ((x138 && x0 <= 9280.0) ? ( + x137 ) : ( 0 )) + 1.0e-99; - Real x101 = x1 <= 100000000.0; + Real x140 = x0 <= 100000000.0; - Real x102 = ((x101 && x1 > 9280.0) ? ( - 1250086.112245841*x95*(1.5400000000000001e-9 + 4.6200000000000001e-10*std::exp(-93988.701501924661*x52))*std::exp(-469943.50750964211*x52) + x98 + Real x141 = ((x140 && x0 > 9280.0) ? ( + 1250086.112245841*x134*(1.5400000000000001e-9 + 4.6200000000000001e-10*std::exp(-93988.701501924661*x10))*std::exp(-469943.50750964211*x10) + x137 ) : ( 0 )) + 1.0e-99; - Real x103 = 3.8571873359681448e-209*x67 + 1.0e-99; + Real x142 = 3.8571873359681448e-209*x110 + 1.0e-99; - Real x104 = -x100*x68 - x102*x68 - x103*x68; + Real x143 = -x111*x139 - x111*x141 - x111*x142; - Real x105 = 2.6534040307116387e-9*std::pow(x1, -0.10000000000000001) + 1.0e-99; + Real x144 = 2.6534040307116387e-9*std::pow(x0, -0.10000000000000001) + 1.0e-99; - Real x106 = X(4)*X(9); + Real x145 = X(4)*X(9); - Real x107 = X(5)*X(7); + Real x146 = X(5)*X(7); - Real x108 = x105*x106 + x105*x107; + Real x147 = x144*x145 + x144*x146; - Real x109 = ((x1 <= 7950.0) ? ( - 1.4308352583277889e-10*std::pow(x1, -0.59999999999999998) + Real x148 = ((x0 <= 7950.0) ? ( + 1.4308352583277889e-10*std::pow(x0, -0.59999999999999998) ) : ( 0 )) + 1.0e-99; - Real x110 = X(16)*X(2); + Real x149 = X(16)*X(2); - Real x111 = ((x1 > 21140.0) ? ( - 0.00023813936486848612*std::pow(x1, -1.3700000000000001)*std::exp(-115786.2*x52) + Real x150 = ((x0 > 21140.0) ? ( + 0.00023813936486848612*std::pow(x0, -1.3700000000000001)*std::exp(-115786.2*x10) ) : ( 0 )) + 1.0e-99; - Real x112 = ((x1 <= 21140.0 && x1 > 7950.0) ? ( - 8.3535906003219973e-24*std::pow(x1, 2.4900000000000002)*std::exp(21845.599999999999*x52) + Real x151 = ((x0 <= 21140.0 && x0 > 7950.0) ? ( + 8.3535906003219973e-24*std::pow(x0, 2.4900000000000002)*std::exp(21845.599999999999*x10) ) : ( 0 )) + 1.0e-99; - Real x113 = -x109*x110 - x110*x111 - x110*x112; + Real x152 = -x148*x149 - x149*x150 - x149*x151; - Real x114 = ((x99 && x1 <= 5500.0) ? ( - x98 + Real x153 = ((x138 && x0 <= 5500.0) ? ( + x137 ) : ( 0 )) + 1.0e-99; - Real x115 = X(2)*X(3); + Real x154 = X(2)*X(3); - Real x116 = ((x101 && x1 > 5500.0) ? ( - 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) + Real x155 = ((x140 && x0 > 5500.0) ? ( + 3.2867337024382733e-10*std::pow(x0, -0.7241125657826851)*std::exp(-1.6921001126637108*x117 - 1.7746868094247411*x70 - 1.9518356165136789*x71 - 2.4649195146505537*x72 - 1.0207737270119371*x73 + 3.3530579587656565*x74 + 3.6203127646377788*x75 - 1.0930705283186734*x76) ) : ( 0 )) + 1.0e-99; - Real x117 = -x114*x115 - x115*x116; - - Real x118 = 35.5*std::pow(x1, -2.2799999999999998)*std::exp(-46707.0*x52) + 1.0e-99; + Real x156 = -x153*x154 - x154*x155; - 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 x157 = 35.5*std::pow(x0, -2.2799999999999998)*std::exp(-46707.0*x10) + 1.0e-99; - Real x120 = -x118*x66 + x119*x73 - x65*x66; + Real x158 = (1.3500000000000001e-9*std::pow(x0, 0.098492999999999997) + 4.4350199999999998e-10*std::pow(x0, 0.55610000000000004) + 3.7408500000000004e-16*std::pow(x0, 2.1825999999999999))/(0.0061910000000000003*std::pow(x0, 1.0461) + 8.9711999999999997e-11*std::pow(x0, 3.0424000000000002) + 3.2575999999999999e-14*std::pow(x0, 3.7740999999999998) + 1.0) + 1.0e-99; - Real x121 = 5.9082438637265071e-70*x69 + 1.0e-99; + Real x159 = -x108*x109 - x109*x157 + x116*x158; - Real x122 = 1.4000000000000001e-18*std::pow(x1, 0.92800000000000005)*std::exp(-6.1728395061728397e-5*x1) + 1.0e-99; + Real x160 = 5.9082438637265071e-70*x112 + 1.0e-99; - Real x123 = -x121*x70 - x122*x70; + Real x161 = 1.4000000000000001e-18*std::pow(x0, 0.92800000000000005)*std::exp(-6.1728395061728397e-5*x0) + 1.0e-99; - Real x124 = 7.9674337148168363e-7*x78 + 1.0e-99; + Real x162 = -x113*x160 - x113*x161; - Real x125 = X(3)*X(9); + Real x163 = 7.9674337148168363e-7*x44 + 1.0e-99; - Real x126 = -x124*x125; + Real x164 = X(3)*X(9); - Real x127 = 2.0000000000000001e-9*X(3)*X(32); + Real x165 = -x163*x164; - Real x128 = -x127; + Real x166 = 2.0000000000000001e-9*X(3)*X(32); - Real x129 = 7.5e-10*X(29)*X(4); + Real x167 = -x166; - Real x130 = 2.8833736969617052e-16*x50 + 1.0e-99; + Real x168 = 2.8833736969617052e-16*x95 + 1.0e-99; - Real x131 = X(14)*X(4); + Real x169 = X(14)*X(4); - Real x132 = ((x1 > 10000.0) ? ( - 4.0000000000000003e-37*std::pow(x1, 4.7400000000000002) + Real x170 = ((x0 > 10000.0) ? ( + 4.0000000000000003e-37*std::pow(x0, 4.7400000000000002) ) : ( 0 )) + 1.0e-99; - Real x133 = X(15)*X(3); + Real x171 = X(15)*X(3); - Real x134 = ((x87) ? ( - 1.26e-9*x36*std::exp(-127500.0*x52) + Real x172 = ((x126) ? ( + 1.26e-9*x81*std::exp(-127500.0*x10) ) : ( 0 )) + 1.0e-99; - Real x135 = x130*x131 - x132*x133 - x133*x134; - - Real x136 = 6.3999999999999996e-10*X(11)*X(4); - - Real x137 = M_LN10; - - Real x138 = 1.0/x137; - - Real x139 = x138*x23; + Real x173 = x168*x169 - x170*x171 - x171*x172; - Real x140 = ((x23)*(x23)); + Real x174 = 6.3999999999999996e-10*X(11)*X(4); - Real x141 = std::exp((-2)*std::log(std::abs(x137))); + Real x175 = M_LN10; - Real x142 = x140*x141; + Real x176 = 1.0/x175; - Real x143 = 8.4600000000000008e-10*x139 - 1.3700000000000002e-10*x142 + 4.1700000000000001e-10; + Real x177 = x176*x68; - Real x144 = X(10)*X(6); + Real x178 = ((x68)*(x68)); - Real x145 = 1.0e-99 + 1.0000000000000001e-9*std::exp(-457.0*x52); + Real x179 = std::exp((-2)*std::log(std::abs(x175))); - Real x146 = X(12)*X(3); + Real x180 = x178*x179; - Real x147 = x136 + x143*x144 - x145*x146; + Real x181 = 8.4600000000000008e-10*x177 - 1.3700000000000002e-10*x180 + 4.1700000000000001e-10; - Real x148 = x1 <= 10.0; + Real x182 = X(10)*X(6); - Real x149 = ((x148) ? ( - 4.2020000000000002e-8 -) -: ( - 0 -)) + 1.0e-99; + Real x183 = 1.0e-99 + 1.0000000000000001e-9*std::exp(-457.0*x10); - Real x150 = X(27)*X(3); + Real x184 = X(12)*X(3); - Real x151 = x1 > 10.0; + Real x185 = x174 + x181*x182 - x183*x184; - Real x152 = ((x151) ? ( - 1.2730571616977901e-7*x78 + 1.767e-9 -) -: ( - 0 -)) + 1.0e-99; - - Real x153 = -x149*x150 - x150*x152; - - Real x154 = ((x151 && x1 <= 100000.0) ? ( - -7.7700000000000002e-13*x1 + 2.96e-6*x10 + 2.5000000000000002e-10*x4 - 1.73e-9 + Real x186 = ((x32 && x0 <= 100000.0) ? ( + -7.7700000000000002e-13*x0 + 2.5000000000000002e-10*x3 + 2.96e-6*x55 - 1.73e-9 ) : ( 0 )); - Real x155 = -x90*x91 + x91*(-x154 - 1.0e-99); + Real x187 = -x129*x130 + x130*(-x186 - 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 x188 = std::pow(10.0, -0.12690000000000001*amrex::Math::powi<-3>(x175)*((x68)*(x68)*(x68)) + 1.1180000000000001*amrex::Math::powi<-2>(x175)*((x68)*(x68)) - 1.5229999999999999*x177 - 19.379999999999999) + 1.0e-99; - Real x157 = -x156; + Real x189 = -x188; - Real x158 = X(3)*X(7); + Real x190 = X(3)*X(7); - Real x159 = x1 > 50.0; + Real x191 = x0 > 50.0; - Real x160 = ((x159) ? ( - 2.0000000000000001e-10*std::pow(x1, 0.40200000000000002)*std::exp(-37.100000000000001*x52) - 3.3099999999999998e-17*std::pow(x1, 1.48) + Real x192 = ((x191) ? ( + 2.0000000000000001e-10*std::pow(x0, 0.40200000000000002)*std::exp(-37.100000000000001*x10) - 3.3099999999999998e-17*std::pow(x0, 1.48) ) : ( 0 )) + 1.0e-99; - Real x161 = ((x159) ? ( - 2.0299999999999998e-9*std::pow(x1, -0.33200000000000002) + 2.0600000000000001e-10*std::pow(x1, 0.39600000000000002)*std::exp(-33.0*x52) + Real x193 = ((x191) ? ( + 2.0299999999999998e-9*std::pow(x0, -0.33200000000000002) + 2.0600000000000001e-10*std::pow(x0, 0.39600000000000002)*std::exp(-33.0*x10) ) : ( 0 )) + 1.0e-99; - Real x162 = X(4)*X(6); + Real x194 = X(4)*X(6); - Real x163 = x157*x158 - x158*x160 + x161*x162; + Real x195 = x189*x190 - x190*x192 + x193*x194; - Real x164 = 8.7882738061352424e-22*std::pow(x1, 1.96)*std::exp(-170000.0*x52) + 1.0e-99; + Real x196 = 8.7882738061352424e-22*std::pow(x0, 1.96)*std::exp(-170000.0*x10) + 1.0e-99; - Real x165 = X(16)*X(4); + Real x197 = X(16)*X(4); - Real x166 = 3.8999999999999998e-16*std::pow(x1, 0.21299999999999999) + 1.0e-99; + Real x198 = 3.8999999999999998e-16*std::pow(x0, 0.21299999999999999) + 1.0e-99; - Real x167 = X(17)*X(3); + Real x199 = X(17)*X(3); - Real x168 = x164*x165 - x166*x167; + Real x200 = x196*x197 - x198*x199; - Real x169 = ((x148) ? ( + Real x201 = ((x29) ? ( 3.7450000000000001e-8 ) : ( 0 )) + 1.0e-99; - Real x170 = X(25)*X(3); + Real x202 = X(25)*X(3); - Real x171 = ((x151) ? ( - 1.1352943524561152e-7*x78 + 1.55e-9 + Real x203 = ((x32) ? ( + 1.1352943524561152e-7*x44 + 1.55e-9 ) : ( 0 )) + 1.0e-99; - Real x172 = -x169*x170 - x170*x171; + Real x204 = -x201*x202 - x202*x203; - Real x173 = 7.5399999999999998e-10*std::pow(x1, -0.45800000000000002) + 4.9899999999999997e-11*std::pow(x1, 0.40500000000000003) + 1.0e-99; + Real x205 = 7.5399999999999998e-10*std::pow(x0, -0.45800000000000002) + 4.9899999999999997e-11*std::pow(x0, 0.40500000000000003) + 1.0e-99; - Real x174 = X(21)*X(4); + Real x206 = X(21)*X(4); - 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 x207 = (4.0000000000000001e-10*std::pow(x0, 0.0066899999999999998) + 1.0799999999999999e-11*std::pow(x0, 0.51700000000000002))*std::exp(-227.0*x10) + 1.0e-99; - Real x176 = X(22)*X(3); + Real x208 = X(22)*X(3); - Real x177 = x173*x174 - x175*x176; + Real x209 = x205*x206 - x207*x208; - Real x178 = 6.0e-10*X(4)*X(8); + Real x210 = 6.0e-10*X(4)*X(8); - Real x179 = ((x23)*(x23)*(x23)); + Real x211 = ((x68)*(x68)*(x68)); - Real x180 = ((((x23)*(x23)))*(((x23)*(x23)))); + Real x212 = ((((x68)*(x68)))*(((x68)*(x68)))); - Real x181 = ((x23)*(x23)*(x23)*(x23)*(x23)); + Real x213 = ((x68)*(x68)*(x68)*(x68)*(x68)); - Real x182 = ((x1 <= 30000.0 && x1 > 100.0) ? ( - (-1.4491368e-7*x140 + 3.4172804999999998e-8*x179 - 4.7813727999999997e-9*x180 + 3.9731542e-10*x181 + 3.5311931999999998e-13*((x23)*(x23)*(x23)*(x23)*(x23)*(x23)*(x23)) - 1.8171411000000001e-11*std::exp((6)*std::log(std::abs(x23))) + 3.3735381999999997e-7*x23 - 3.3232183000000002e-7)*std::exp(-21237.150000000001*x52) + Real x214 = ((x0 <= 30000.0 && x0 > 100.0) ? ( + (-1.4491368e-7*x178 + 3.4172804999999998e-8*x211 - 4.7813727999999997e-9*x212 + 3.9731542e-10*x213 + 3.5311931999999998e-13*((x68)*(x68)*(x68)*(x68)*(x68)*(x68)*(x68)) - 1.8171411000000001e-11*std::exp((6)*std::log(std::abs(x68))) + 3.3735381999999997e-7*x68 - 3.3232183000000002e-7)*std::exp(-21237.150000000001*x10) ) : ( 0 )) + 1.0e-99; - Real x183 = X(10)*X(3); + Real x215 = X(10)*X(3); - Real x184 = x178 - x182*x183; + Real x216 = x210 - x214*x215; - Real x185 = x179/((x137)*(x137)*(x137)); + Real x217 = x211/((x175)*(x175)*(x175)); - Real x186 = ((x1 > 30.0) ? ( - std::pow(10, -3.194*x139 + 1.786*x140*x141 - 0.2072*x185 - 18.199999999999999) + Real x218 = ((x0 > 30.0) ? ( + std::pow(10, -3.194*x177 + 1.786*x178*x179 - 0.2072*x217 - 18.199999999999999) ) : ( 0 )) + 1.0e-99; - Real x187 = X(3)*X(4); + Real x219 = X(3)*X(4); - Real x188 = std::pow(x1, -0.14999999999999999); + Real x220 = std::pow(x0, -0.14999999999999999); - Real x189 = ((x1 <= 30.0) ? ( - 3.4977396723747635e-20*x188 + Real x221 = ((x0 <= 30.0) ? ( + 3.4977396723747635e-20*x220 ) : ( 0 )) + 1.0e-99; - Real x190 = -x186*x187 - x187*x189; - - Real x191 = 6.3999999999999996e-10*X(10)*X(26); - - Real x192 = -x136; - - Real x193 = 1.0000000000000001e-17*X(17); + Real x222 = -x218*x219 - x219*x221; - Real x194 = X(4)*x193; + Real x223 = 6.3999999999999996e-10*X(10)*X(26); - Real x195 = -x194; + Real x224 = -x174; - Real x196 = ((X(10))*(X(10))); + Real x225 = 1.0000000000000001e-17*X(17); - Real x197 = std::log(T); + Real x226 = X(4)*x225; - 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 x227 = -x226; - 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 x228 = ((X(10))*(X(10))); - 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 x229 = std::log(T); - Real x201 = 2*x200 + 2.0e-99; + Real x230 = 2.0*X(1) + 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) + 3.0*X(28) + X(3) + X(4) + X(5) + 2.0*X(8); - Real x202 = ((X(4))*(X(4))*(X(4))); + Real x231 = 1.0/(std::pow(10.0, 1.3*x176*(x229 - 9.2103403719761836) - 137.42519902360013*x179*((0.10857362047581294*x229 - 1)*(0.10857362047581294*x229 - 1)) - 4.8449999999999998)*x230 + 1.0); - Real x203 = 0.25820969999999999*x185; + Real x232 = std::pow(1.1800000000000001e-10*std::exp(-69500.0*x10), x231)*std::pow(8.1250000000000003e-8*std::pow(T, -0.5)*(1.0 - std::exp(-6000.0*x10))*std::exp(-52000.0*x10), 1.0 - x231); - Real x204 = std::log(27535.310000000001*x52 + 1.0); + Real x233 = 2*x232 + 2.0e-99; - Real x205 = -2.4640089999999999*x138*x23 + 0.19859550000000001*x142 + 743.05999999999995*x52; + Real x234 = ((X(4))*(X(4))*(X(4))); - Real x206 = std::exp(-0.0022727272727272726*x1); + Real x235 = 0.25820969999999999*x217; - Real x207 = std::exp(-0.00054054054054054055*x1); + Real x236 = std::log(27535.310000000001*x10 + 1.0); - Real x208 = 2.9375070000000001*x206 + 0.23588480000000001*x207 + 0.75022860000000002; + Real x237 = 743.05999999999995*x10 - 2.4640089999999999*x176*x68 + 0.19859550000000001*x180; - 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 x238 = std::exp(-0.0022727272727272726*x0); - Real x210 = 4.6331670000000003*x185; + Real x239 = std::exp(-0.00054054054054054055*x0); - Real x211 = std::log(40870.379999999997*x52 + 1.0); + Real x240 = 2.9375070000000001*x238 + 0.23588480000000001*x239 + 0.75022860000000002; - Real x212 = -4.8909149999999997*x138*x23 + 0.47490300000000002*x142 - 133.82830000000001*x52; + Real x241 = std::pow(10.0, -21467.790000000001*x10 - 1657.4099999999999*x10/(std::pow(std::pow(10.0, -x237 - 8.1313220000000008)*x230, x240) + 1.0) + 21.360939999999999*x176*x236 + 42.707410000000003*x176*x68 - 2.0273650000000001*x180 - x235 - 142.7664 - (21.360939999999999*x176*x236 + 11.28215*x176*std::log(14254.549999999999*x10 + 1.0) + 70.138370000000009*x177 - 4.7035149999999994*x180 - x235 - 203.11568)/(std::pow(std::pow(10.0, -x237 - 9.3055640000000004)*x230, x240) + 1.0)); - Real x213 = -2.0563129999999998*x206 + 0.58640729999999996*x207 + 0.82274429999999998; + Real x242 = 4.6331670000000003*x217; - 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 x243 = std::log(40870.379999999997*x10 + 1.0); - Real x215 = X(10)*X(4); + Real x244 = -133.82830000000001*x10 - 4.8909149999999997*x176*x68 + 0.47490300000000002*x180; - Real x216 = 1.5e-32*x51 + 5.0000000000000004e-32*x78 + 2.0e-99; + Real x245 = -2.0563129999999998*x238 + 0.58640729999999996*x239 + 0.82274429999999998; - Real x217 = X(10)*((X(4))*(X(4))); + Real x246 = std::pow(10.0, -23705.700000000001*x10 - 2080.4099999999999*x10/(std::pow(std::pow(10.0, -x244 - 13.656822)*x230, x245) + 1.0) + 69.700860000000006*x176*x243 - 68.422430000000006*x177 + 43.20243*x178*x179 - x242 - 178.4239 - (69.700860000000006*x176*x243 + 19.734269999999999*x176*std::log(16780.950000000001*x10 + 1.0) - 14.509090000000008*x177 + 37.886913*x180 - x242 - 307.31920000000002)/(std::pow(std::pow(10.0, -x244 - 14.82123)*x230, x245) + 1.0)); - 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 x247 = X(10)*X(4); - Real x219 = 1.1194685000500219e-20*std::pow(x1, 2.7400000000000002)*std::exp(-4740.0*x52) + 1.0e-99; + Real x248 = 5.0000000000000004e-32*x44 + 1.5e-32*x96 + 2.0e-99; - Real x220 = X(10)*X(20); + Real x249 = X(10)*((X(4))*(X(4))); - Real x221 = 2.1555555555555562e-20*((x1)*(x1)*(x1))*std::exp(-4045.0*x52) + 1.0e-99; + Real x250 = X(4)*x230*x4/((1.0 + 10000.0*std::exp(-600.0/(Tdust + 9.9999999999999993e-41)))*(7.9999999999999996e-6*x1 + x2 + 0.040000000000000001*std::sqrt(Tdust + x0))); - Real x222 = X(23)*X(4); + Real x251 = 1.1194685000500219e-20*std::pow(x0, 2.7400000000000002)*std::exp(-4740.0*x10) + 1.0e-99; - Real x223 = x219*x220 - x221*x222; + Real x252 = X(10)*X(20); - Real x224 = 1.9643512449766835e-11*std::pow(x1, 0.17000000000000001)*std::exp(-6400.0*x52) + 1.0e-99; + Real x253 = 2.1555555555555562e-20*((x0)*(x0)*(x0))*std::exp(-4045.0*x10) + 1.0e-99; - Real x225 = X(10)*X(19); + Real x254 = X(23)*X(4); - Real x226 = 1.0e-99 + 1.0e-10*std::exp(-7600.0*x52); + Real x255 = x251*x252 - x253*x254; - Real x227 = X(20)*X(4); + Real x256 = 1.9643512449766835e-11*std::pow(x0, 0.17000000000000001)*std::exp(-6400.0*x10) + 1.0e-99; - Real x228 = x224*x225 - x226*x227; + Real x257 = X(10)*X(19); - Real x229 = 2.1999999999999999e-10*X(19)*X(4); + Real x258 = 1.0e-99 + 1.0e-10*std::exp(-7600.0*x10); - Real x230 = 1.0e-99 + 5.4599999999999998e-10*std::exp(-1943.0*x52); + Real x259 = X(20)*X(4); - Real x231 = X(10)*X(18); + Real x260 = x256*x257 - x258*x259; - Real x232 = -x229 + x230*x231; + Real x261 = 2.1999999999999999e-10*X(19)*X(4); - Real x233 = 1.69e-9*X(10)*X(21); + Real x262 = 1.0e-99 + 5.4599999999999998e-10*std::exp(-1943.0*x10); - Real x234 = x169*x170 + x170*x171 + x233; + Real x263 = X(10)*X(18); - Real x235 = 1.0e-13*X(25)*X(30); + Real x264 = -x261 + x262*x263; - Real x236 = 1.0e-99 + 2.54e-10*std::exp(-13300.0*x52); + Real x265 = 1.69e-9*X(10)*X(21); - Real x237 = X(33)*X(4); + Real x266 = x201*x202 + x202*x203 + x265; - Real x238 = x235 - x236*x237; + Real x267 = x20 - x21*x22; - Real x239 = 6.1739095063118665e-10*std::pow(x1, 0.40999999999999998) + 1.0e-99; + Real x268 = 6.1739095063118665e-10*std::pow(x0, 0.40999999999999998) + 1.0e-99; - Real x240 = -x239; + Real x269 = -x268; - Real x241 = -x105; + Real x270 = -x144; - Real x242 = x106*x240 + x106*x241 + x107*x239; + Real x271 = x145*x269 + x145*x270 + x146*x268; - Real x243 = 1.0e-99 + 6.6399999999999998e-10*std::exp(-11700.0*x52); + Real x272 = 1.0e-99 + 6.6399999999999998e-10*std::exp(-11700.0*x10); - Real x244 = X(10)*X(17); + Real x273 = X(10)*X(17); - Real x245 = 1.0e-99 + 1.3100000000000001e-10*std::exp(-80.0*x52); + Real x274 = 1.0e-99 + 1.3100000000000001e-10*std::exp(-80.0*x10); - Real x246 = X(18)*X(4); + Real x275 = X(18)*X(4); - Real x247 = x243*x244 - x245*x246; + Real x276 = x272*x273 - x274*x275; - Real x248 = 1.01e-9*X(10)*X(24); + Real x277 = 1.01e-9*X(10)*X(24); - Real x249 = x149*x150 + x150*x152 + x248; + Real x278 = x277 + x42*x43 + x43*x45; - Real x250 = x157*x162 + x158*x160 - x161*x162; + Real x279 = x189*x194 + x190*x192 - x193*x194; - Real x251 = -x178 + x182*x183; + Real x280 = -x210 + x214*x215; - Real x252 = -x173*x174 + x175*x176; + Real x281 = -x205*x206 + x207*x208; - Real x253 = -x130*x131 + x132*x133 + x133*x134; + Real x282 = -x168*x169 + x170*x171 + x171*x172; - Real x254 = ((x148) ? ( + Real x283 = ((x29) ? ( 1.371e-8 ) : ( 0 )) + 1.0e-99; - Real x255 = X(16)*X(25); + Real x284 = X(16)*X(25); - Real x256 = ((x151) ? ( - 4.1551773299893819e-8*x78 + 5.673e-10 + Real x285 = ((x32) ? ( + 4.1551773299893819e-8*x44 + 5.673e-10 ) : ( 0 )) + 1.0e-99; - Real x257 = -x129 + x254*x255 + x255*x256; + Real x286 = -x14 + x283*x284 + x284*x285; - Real x258 = ((x1 <= 1167.4796423742259) ? ( - 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))) + Real x287 = ((x0 <= 1167.4796423742259) ? ( + std::pow(10, 5.8888600000000002*x177 + 7.1969200000000004*x180 + 2.2506900000000001*x217 - 56.473700000000001 - 2.1690299999999998*x212/((((x175)*(x175)))*(((x175)*(x175)))) + 0.31788699999999998*x213/((x175)*(x175)*(x175)*(x175)*(x175))) ) : ( 0 )) + 1.0e-99; - Real x259 = X(10)*X(7); + Real x288 = X(10)*X(7); - Real x260 = ((x1 > 1167.4796423742259) ? ( - 3.1699999999999999e-10*std::exp(-5207.0*x52) + Real x289 = ((x0 > 1167.4796423742259) ? ( + 3.1699999999999999e-10*std::exp(-5207.0*x10) ) : ( 0 )) + 1.0e-99; - Real x261 = ((x1 > 200.0) ? ( - 5.25e-11*std::exp(-4430.0*x52 + 173900.0/x2) + Real x290 = ((x0 > 200.0) ? ( + 5.25e-11*std::exp(-4430.0*x10 + 173900.0/x1) ) : ( 0 )) + 1.0e-99; - Real x262 = X(12)*X(4); - - Real x263 = x258*x259 + x259*x260 - x261*x262; - - Real x264 = 3.5999999999999998e-16*std::pow(T, 1.52)*std::exp(-1740.0*x52) + 1.0e-99; - - Real x265 = X(10)*X(25); - - Real x266 = 6.3600000000000002e-14*x1*std::exp(-9610.0*x52) + 1.0e-99; + Real x291 = X(12)*X(4); - Real x267 = X(27)*X(4); + Real x292 = x287*x288 + x288*x289 - x290*x291; - Real x268 = x264*x265 - x266*x267 + x80*x81; - - Real x269 = ((x1 <= 150.0) ? ( + Real x293 = ((x0 <= 150.0) ? ( 4.9970000000000003e-11 ) : ( 0 )) + 1.0e-99; - Real x270 = X(22)*X(25); + Real x294 = X(22)*X(25); - Real x271 = ((x1 > 150.0) ? ( - 2.4000000000000001e-11*std::exp(110.0*x52) + Real x295 = ((x0 > 150.0) ? ( + 2.4000000000000001e-11*std::exp(110.0*x10) ) : ( 0 )) + 1.0e-99; - Real x272 = 1.0e-99 + 3.1319999999999996e-10*std::exp(-8156.0*x52); + Real x296 = 1.0e-99 + 3.1319999999999996e-10*std::exp(-8156.0*x10); - Real x273 = X(32)*X(4); + Real x297 = X(32)*X(4); - Real x274 = x269*x270 + x270*x271 - x272*x273; + Real x298 = x293*x294 + x294*x295 - x296*x297; - Real x275 = x118*x66 - x119*x73 - x72*x73 - x73*x76; + Real x299 = x109*x157 - x115*x116 - x116*x119 - x116*x158; - Real x276 = -x164*x165 + x166*x167; + Real x300 = -x196*x197 + x198*x199; - Real x277 = ((x148) ? ( - 7.0509999999999997e-11 -) -: ( - 0 -)) + 1.0e-99; + Real x301 = 1.0e-25*X(4)*X(7); - Real x278 = X(17)*X(25); + Real x302 = x163*x164 - x301 + x62*x63; - Real x279 = ((x151) ? ( - 1.5557027645121538e-10*std::pow(x1, -0.33900000000000002)*std::exp(-0.108*x52) -) -: ( - 0 -)) + 1.0e-99; - - Real x280 = 1.0e-99 + 1.02e-10*std::exp(-914.0*x52); - - Real x281 = X(18)*X(22); - - Real x282 = x277*x278 + x278*x279 + x280*x281; + Real x303 = 1.0e-99 + 1.46e-12*std::exp(-9650.0*x10); - Real x283 = 1.0e-25*X(4)*X(7); + Real x304 = X(10)*X(22); - Real x284 = x124*x125 + x17*x18 - x283; - - Real x285 = 1.0e-99 + 1.46e-12*std::exp(-9650.0*x52); - - Real x286 = X(10)*X(22); - - Real x287 = ((x1 > 280.0) ? ( + Real x305 = ((x0 > 280.0) ? ( 5.4499999999999998e-17 ) : ( 0 )) + 1.0e-99; - Real x288 = X(25)*X(4); + Real x306 = X(25)*X(4); - Real x289 = ((x1 <= 280.0) ? ( - 8.1009819130659227e-21*std::pow(x1, 2.7999999999999998)*std::exp(-1950.0*x52) + Real x307 = ((x0 <= 280.0) ? ( + 8.1009819130659227e-21*std::pow(x0, 2.7999999999999998)*std::exp(-1950.0*x10) ) : ( 0 )) + 1.0e-99; - Real x290 = x285*x286 - x287*x288 - x288*x289 + x83*x84; - - Real x291 = 8.648509591313017e-18*std::pow(x1, -0.38) + 1.0e-99; + Real x308 = x122*x123 + x303*x304 - x305*x306 - x306*x307; - Real x292 = X(22)*X(4); + Real x309 = 8.648509591313017e-18*std::pow(x0, -0.38) + 1.0e-99; - Real x293 = x20*x21 - x291*x292; + Real x310 = X(22)*X(4); - Real x294 = 4.9999999999999996e-6*x10 + 1.0e-99; + Real x311 = -x309*x310 + x65*x66; - Real x295 = X(5)*X(8); + Real x312 = 4.9999999999999996e-6*x55 + 1.0e-99; - Real x296 = x202*(6.0000000000000001e-32*x51 + 2.0000000000000002e-31*x78 + 1.0e-99) + x215*(-x209 - x214 - 3.0000000000000001e-99) + x294*x295; + Real x313 = X(5)*X(8); - Real x297 = -x294*x295; + Real x314 = x234*(2.0000000000000002e-31*x44 + 6.0000000000000001e-32*x96 + 1.0e-99) + x247*(-x241 - x246 - 3.0000000000000001e-99) + x312*x313; - Real x298 = x106*x239 + x107*x240 + x107*x241; + Real x315 = -x312*x313; - Real x299 = X(6)*X(9); + Real x316 = x145*x268 + x146*x269 + x146*x270; - Real x300 = x299*(-9.8726896031426014e-7*x78 - 1.0e-99); + Real x317 = X(6)*X(9); - Real x301 = -x143*x144 + x145*x146; + Real x318 = x317*(-9.8726896031426014e-7*x44 - 1.0e-99); - Real x302 = -x258*x259 - x259*x260 + x261*x262; + Real x319 = -x181*x182 + x183*x184; - Real x303 = -x248; + Real x320 = -x287*x288 - x288*x289 + x290*x291; - Real x304 = 1.6974097914175e-15*x77*std::exp(-7550.0*x52) + 1.0e-99; + Real x321 = -x277; - Real x305 = X(10)*X(33); + Real x322 = -x11*x12; - Real x306 = -x304*x305; + Real x323 = -x223; - Real x307 = -x191; + Real x324 = -x265; - Real x308 = -x233; + Real x325 = -x251*x252 + x253*x254; - Real x309 = -x219*x220 + x221*x222; + Real x326 = -x256*x257 + x258*x259; - Real x310 = -x224*x225 + x226*x227; + Real x327 = x261 - x262*x263; - Real x311 = x229 - x230*x231; + Real x328 = X(10)*x225; - Real x312 = X(10)*x193; + Real x329 = -x272*x273 + x274*x275 - x328; - Real x313 = -x243*x244 + x245*x246 - x312; + Real x330 = -x47*x48 + x51*x52; - Real x314 = -x264*x265 + x266*x267; + Real x331 = x123*x124 - x303*x304 + x305*x306 + x306*x307; - Real x315 = -x285*x286 + x287*x288 + x288*x289 + x84*x85; + Real x332 = X(17)*X(31); - Real x316 = X(17)*X(31); + Real x333 = 5.2000000000000001e-11*x332; - Real x317 = 5.2000000000000001e-11*x316; + Real x334 = x333 - 7.9500000000000005e-10*x7; - Real x318 = X(16)*X(32); + Real x335 = -x283*x284 - x284*x285; - Real x319 = x317 - 7.9500000000000005e-10*x318; - - Real x320 = -x254*x255 - x255*x256; - - Real x321 = ((x1 <= 300.0) ? ( + Real x336 = ((x0 <= 300.0) ? ( 2.5000000000000002e-18 ) : ( 0 )) + 1.0e-99; - Real x322 = X(16)*X(22); - - Real x323 = ((x1 > 300.0) ? ( - 7.3875058242066447e-18*x188*std::exp(-68.0*x52) -) -: ( - 0 -)) + 1.0e-99; - - Real x324 = -x321*x322 - x322*x323; + Real x337 = X(16)*X(22); - Real x325 = -1.04e-10*x316; - - Real x326 = x40*x41; - - Real x327 = ((x1 <= 1052.0) ? ( - 3.2682787072746432e-10*x58 + Real x338 = ((x0 > 300.0) ? ( + 7.3875058242066447e-18*x220*std::exp(-68.0*x10) ) : ( 0 )) + 1.0e-99; - Real x328 = X(17)*X(32); - - Real x329 = ((x1 > 1052.0) ? ( - 3.7991369449339675e-16*std::pow(x1, 1.54)*std::exp(613.0*x52) -) -: ( - 0 -)) + 1.0e-99; - - Real x330 = -x327*x328 - x328*x329; - - Real x331 = -x277*x278 - x278*x279; - - Real x332 = -x280*x281; - - Real x333 = 4.5299999999999999e-10*x318; - - Real x334 = ((X(25))*(X(25))); - - 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::pow(x1, 1.5800000000000001); - - Real x339 = 1.0e-99 + 2.8e-11*std::exp(-26500.0*x52); + Real x339 = -x336*x337 - x337*x338; - Real x340 = X(22)*X(33); + Real x340 = -1.04e-10*x332; - Real x341 = 1.0e-99 + 4.1999999999999999e-12*std::exp(-24000.0*x52); + Real x341 = x85*x86; - Real x342 = X(30)*X(32); + Real x342 = -x25*x26 - x26*x27; - Real x343 = -x339*x340 + x341*x342; + Real x343 = -x30*x31 - x31*x33; - Real x344 = x317 + 3.4200000000000001e-10*x318; + Real x344 = -x34*x35; - Real x345 = x327*x328 + x328*x329; + Real x345 = ((X(22))*(X(22))); - Real x346 = -x269*x270 - x270*x271 + x272*x273; + Real x346 = std::pow(x0, 1.5800000000000001); - Real x347 = -x235 + x236*x237; + Real x347 = -x15*x16 + x17*x18; - Real x348 = x304*x305; + Real x348 = x333 + 3.4200000000000001e-10*x7; - Real x349 = x339*x340 - x341*x342; + Real x349 = -x293*x294 - x294*x295 + x296*x297; - ydot(1) = X(30)*x7 - x0; + ydot(1) = X(30)*x6 + x37; - ydot(2) = X(27)*x9 - x8; + ydot(2) = X(27)*x38 + x54; - ydot(3) = x104 + x108 + x113 + x117 + x120 + x123 + x13 + x16 + x19 + x22 + x34*(8.7048158229535104e-117*x32 + 2.0e-99) + x35 + x39 + x42 + x45*(7.5807998548789035e-18*x43 + 2.0e-99) + x46 + x49 + x56*(2.6047246489532126e-6*x54 + 2.0e-99) + x57 + x62*(4.2024322823973014e-6*x60 + 2.0e-99) + x63 + x65*x66 + x68*(7.7143746719362896e-209*x67 + 2.0e-99) + x70*(1.1816487727453014e-69*x69 + 2.0e-99) + x72*x73 + x73*x76 + x82 + x86 + x92 + x97; + ydot(3) = x100*(2.6047246489532126e-6*x98 + 2.0e-99) + x101 + x105*(4.2024322823973014e-6*x103 + 2.0e-99) + x106 + x108*x109 + x111*(7.7143746719362896e-209*x110 + 2.0e-99) + x113*(1.1816487727453014e-69*x112 + 2.0e-99) + x115*x116 + x116*x119 + x121 + x125 + x131 + x136 + x143 + x147 + x152 + x156 + x159 + x162 + x58 + x61 + x64 + x67 + x79*(8.7048158229535104e-117*x77 + 2.0e-99) + x80 + x84 + x87 + x90*(7.5807998548789035e-18*x88 + 2.0e-99) + x91 + x94; - ydot(4) = x117 + x121*x70 + x126 + x128 + x129 + x135 + x147 + x153 + x155 + x163 + x168 + x172 + x177 + x184 + x190; + ydot(4) = x113*x160 + x14 + x156 + x165 + x167 + x173 + x185 + x187 + x195 + x200 + x204 + x209 + x216 + x222 + x46; - ydot(5) = x114*x115 + x115*x116 + x123 + x127 + x190 + x191 + x192 + x195 + x196*x201 + x202*(-1.8e-31*x51 - 6.0000000000000005e-31*x78 - 3.0000000000000001e-99) + x215*(3*x209 + 3*x214 + 9.0000000000000006e-99) - x216*x217 - 5.9999999999999997e-18*x218 + x223 + x228 + x232 + x234 + x238 + x242 + x247 + x249 + x250 + x251 + x252 + x253 + x257 + x263 + x268 + x274 + x275 + x276 + x282 + x284 + x290 + x293 + x296 + x44*x45 + x66*(8.7599999999999997e-10*x64 + 2.0e-99) + x73*(2*x71 + 2.0e-99) + x73*(2*x75 + 2.0e-99) + x80*(8.9373821670554065e-6*x78 + 2.0e-99) + x89*(2*x88 + 2.0e-99) + x91*(2*x154 + 2.0e-99); + ydot(5) = x109*(8.7599999999999997e-10*x107 + 2.0e-99) + x116*(2*x114 + 2.0e-99) + x116*(2*x118 + 2.0e-99) + x128*(2*x127 + 2.0e-99) + x130*(2*x186 + 2.0e-99) + x153*x154 + x154*x155 + x162 + x166 + x222 + x223 + x224 + x227 + x228*x233 + x234*(-6.0000000000000005e-31*x44 - 1.8e-31*x96 - 3.0000000000000001e-99) + x247*(3*x241 + 3*x246 + 9.0000000000000006e-99) - x248*x249 - 5.9999999999999997e-18*x250 + x255 + x260 + x264 + x266 + x267 + x271 + x276 + x278 + x279 + x280 + x281 + x282 + x286 + x292 + x298 + x299 + x300 + x302 + x308 + x311 + x314 + x36 + x50*(8.9373821670554065e-6*x44 + 2.0e-99) + x53 + x89*x90; - ydot(6) = x122*x70 + x155 + x275 + x297 + x298 + x46; + ydot(6) = x113*x161 + x187 + x299 + x315 + x316 + x91; - ydot(7) = x250 + x300 + x301 + x39; + ydot(7) = x279 + x318 + x319 + x84; - ydot(8) = x163 + x284 + x298 + x299*(1.9745379206285203e-6*x78 + 2.0e-99) + x302 + x37*x38 + x49; + ydot(8) = x195 + x302 + x316 + x317*(1.9745379206285203e-6*x44 + 2.0e-99) + x320 + x82*x83 + x94; - ydot(9) = x186*x187 + x187*x189 + x251 + x297 + x92; + ydot(9) = x131 + x218*x219 + x219*x221 + x280 + x315; - ydot(10) = x126 + x242 + x300 + x47*x48; + ydot(10) = x165 + x271 + x318 + x92*x93; - ydot(11) = x120 + x184 - x196*x201 + x196*(x200 + 1.0e-99) + x216*x217 + x217*(-7.5000000000000001e-33*x51 - 2.5000000000000002e-32*x78 - 1.0e-99) + 2.9999999999999998e-18*x218 + x296 + x301 + x302 + x303 + x306 + x307 + x308 + x309 + x310 + x311 + x313 + x314 + x315; + ydot(11) = x159 + x216 - x228*x233 + x228*(x232 + 1.0e-99) + x248*x249 + x249*(-2.5000000000000002e-32*x44 - 7.5000000000000001e-33*x96 - 1.0e-99) + 2.9999999999999998e-18*x250 + x314 + x319 + x320 + x321 + x322 + x323 + x324 + x325 + x326 + x327 + x329 + x330 + x331; - ydot(12) = x156*x158 + x156*x162 + x19 + x192; + ydot(12) = x188*x190 + x188*x194 + x224 + x64; - ydot(13) = x108 + x147 + x263 + x283; + ydot(13) = x147 + x185 + x292 + x301; - ydot(14) = x103*x68 + x13; + ydot(14) = x111*x142 + x58; - ydot(15) = x104 + x11*x12 + x253 + x33*x34; + ydot(15) = x143 + x282 + x56*x57 + x78*x79; - ydot(16) = x100*x68 + x102*x68 + x135 + x35; + ydot(16) = x111*x139 + x111*x141 + x173 + x80; - ydot(17) = x113 + x276 + x319 + x320 + x324 + x55*x56; + ydot(17) = x100*x99 + x152 + x300 + x334 + x335 + x339; - ydot(18) = x109*x110 + x110*x111 + x110*x112 + x168 + x195 + x313 + x325 + x326 + x330 + x331 + x57; + ydot(18) = x101 + x148*x149 + x149*x150 + x149*x151 + x200 + x227 + x329 + x340 + x341 + x342 + x343; - ydot(19) = x194 + x247 + x311 + x332; + ydot(19) = x226 + x276 + x327 + x344; - ydot(20) = x232 + x310 + x312; + ydot(20) = x264 + x326 + x328; - ydot(21) = x228 + x309; + ydot(21) = x260 + x325; - ydot(22) = x252 + x308 + x333 + x61*x62 + x97; + ydot(22) = x104*x105 + x136 + x281 + x324 + x8; - ydot(23) = x15*(2.1137298733407573e-5*x14 + 2.0e-99) + x177 + x293 + x315 + x324 + x326 + x332 + x336 + x337*(-1.1950181274729681e-23*x338 - 2.0e-99) + x343 + x344 + x345 + x346 + x63 + x93*x94 + x94*x96; + ydot(23) = x106 + x132*x133 + x133*x135 + x209 + x28 + x311 + x331 + x339 + x341 + x344 + x345*(-1.1950181274729681e-23*x346 - 2.0e-99) + x347 + x348 + x349 + x41 + x60*(2.1137298733407573e-5*x59 + 2.0e-99); - ydot(24) = x223; + ydot(24) = x255; - ydot(25) = x22 + x234 + x303; + ydot(25) = x266 + x321 + x67; - ydot(26) = x172 + x290 + x291*x292 + x314 + x320 + x331 + x334*(-4.9498918425228625e-15*x335 - 2.0e-99) + x346 + x347 + x79*x80; + ydot(26) = x120*x50 + x204 + x23 + x308 + x309*x310 + x330 + x335 + x343 + x349 + x39*(-4.9498918425228625e-15*x40 - 2.0e-99); - ydot(27) = x249 + x307 + x86; + ydot(27) = x125 + x278 + x323; - ydot(28) = -X(27)*x9 + x153 + x268 + x336 + x348 + x8; + ydot(28) = 2.0e-99*X(1) + x54; - ydot(29) = x191 + x82; + ydot(29) = x121 + x223; - ydot(30) = x257 + x321*x322 + x322*x323 + x344 + x42; + ydot(30) = x286 + x336*x337 + x337*x338 + x348 + x87; - ydot(31) = -X(30)*x7 + x0 + x129 + x282 + x333 + x345 + x347 + x348 + x349; + ydot(31) = 2.0e-99*X(0) + x37; - ydot(32) = x127 + x16 + x325; + ydot(32) = x166 + x340 + x61; - ydot(33) = x128 + x274 + x319 + x330 + x337*(5.9750906373648406e-24*x338 + 1.0e-99) + x349; + ydot(33) = x167 + x19 + x298 + x334 + x342 + x345*(5.9750906373648406e-24*x346 + 1.0e-99); - ydot(34) = x238 + x306 + x343; + ydot(34) = x267 + x322 + x347; @@ -1103,7 +1103,7 @@ Real rhs_eint(const burn_t& state, Real x24 = X(11) + X(12) + X(18) + X(24) + X(25) + X(5); - Real x25 = 2.0*X(10) + 2.0*X(19) + 3.0*X(20) + 4.0*X(23) + 2.0*X(26) + 2.0*X(27) + 3.0*X(28) + X(3) + X(4) + 2.0*X(8) + x24; + Real x25 = 2.0*X(1) + 2.0*X(10) + 2.0*X(19) + 3.0*X(20) + 4.0*X(23) + 2.0*X(26) + 3.0*X(28) + X(3) + X(4) + 2.0*X(8) + x24; Real x26 = 1.0/x25; @@ -3962,7 +3962,7 @@ Real rhs_eint(const burn_t& state, ) : ( -9.9999999999999997e+98 - ))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))*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(1) + 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.75*X(28) + 0.25*X(3) + 0.25*X(4) + 0.25*X(5) + 0.5*X(8))*(0.5*X(1) + 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.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 ) : ( @@ -4117,8 +4117,8 @@ std::pair compute_Semenov_Tdust(Real Tgas, const Array1D Date: Tue, 17 Sep 2024 17:38:15 +0200 Subject: [PATCH 2/7] update reference solutions --- .../reference_solution_1.out | 74 ++++++++-------- .../reference_solution_1_z10.out | 78 ++++++++--------- .../reference_solution_1e-1.out | 74 ++++++++-------- .../reference_solution_1e-2.out | 74 ++++++++-------- .../reference_solution_1e-3.out | 74 ++++++++-------- .../reference_solution_1e-4.out | 74 ++++++++-------- .../reference_solution_1e-5.out | 76 ++++++++--------- .../reference_solution_1e-6.out | 84 +++++++++---------- 8 files changed, 304 insertions(+), 304 deletions(-) diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1.out b/unit_test/burn_cell_metal_chem/reference_solution_1.out index f81ada569..27a1ffb23 100644 --- a/unit_test/burn_cell_metal_chem/reference_solution_1.out +++ b/unit_test/burn_cell_metal_chem/reference_solution_1.out @@ -1,45 +1,45 @@ -T final = 2033.7719 +T final = 2033.663091 Tdust initial = 2.73 -Tdust final = 2033.610864 +Tdust final = 2033.629887 Eint initial = 3.061675486e+10 -Eint final = 1.764013907e+11 +Eint final = 1.763898902e+11 rho initial = 2.203783613e-25 -rho final = 1.969175663e-06 -numdens final = 5.200225029e+17 +rho final = 1.969712207e-06 +numdens final = 5.200218306e+17 ------------------------------------ New number densities: 8.329619835e+13 -2.373088458e+14 -2.040789842 -6.529533105e-08 -1.944958889e+15 -2.554410631e-05 -6.055749336e-13 -9860014349 -1.817920945e-10 -1.199955649e-10 -4.481098213e+17 -2.444960273e-20 -8.975750696e+12 -2.630091902e-09 -8.82330545e-21 +2.373088567e+14 +2.040909947 +1.074117805e-06 +1.943614057e+15 +2.553528356e-05 +9.961696151e-12 +9853431998 +2.99089071e-09 +1.199571059e-10 +4.481104939e+17 +4.021635872e-19 +8.975754713e+12 +1.746215894e-09 +5.862794549e-21 6.963813204e+16 -9.871889761e-21 -6.739096013e-14 -2.597905176e-13 -5.714209204e-11 -1.703374921e-08 -2.53711863e-40 -3.383022442e-15 -1.998329514e-07 -3.16125345e-35 -5.387126028e-11 -1.825600604e-31 -1.766483831e-07 -2.040750381 -1.559962404e-34 -8.563324351e-29 -6.493533745e-05 -2.06061654e-32 -2.973172266e-55 +1.092755289e-19 +4.6280631e-13 +1.784798056e-12 +3.928263544e-10 +1.171832259e-07 +5.089162611e-39 +4.125246338e-15 +1.375692779e-06 +5.196274139e-34 +5.382850015e-11 +3.0339037e-30 +1.766616999e-07 +2.040849206 +1.406529162e-33 +7.715200326e-28 +8.519569305e-05 +1.858338724e-31 +2.678201778e-54 AMReX (22.07-589-g9b0d63e40412) finalized diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1_z10.out b/unit_test/burn_cell_metal_chem/reference_solution_1_z10.out index e485c3b44..9158256e4 100644 --- a/unit_test/burn_cell_metal_chem/reference_solution_1_z10.out +++ b/unit_test/burn_cell_metal_chem/reference_solution_1_z10.out @@ -1,47 +1,47 @@ -Tgas = 31.53250536 -Tgas = 31.52704064 -T final = 2033.846636 +Tgas = 31.53251506 +Tgas = 31.5270498 +T final = 2033.744028 Tdust initial = 30.03 -Tdust final = 2033.692617 +Tdust final = 2033.708884 Eint initial = 3.061675486e+10 -Eint final = 1.764079063e+11 +Eint final = 1.763969463e+11 rho initial = 2.203783613e-25 -rho final = 1.969175663e-06 -numdens final = 5.20022965e+17 +rho final = 1.969712208e-06 +numdens final = 5.200223311e+17 ------------------------------------ -New number densities: +New number densities: 8.329620029e+13 2.373088525e+14 -2.226446088 -0.000297491388 -1.945883073e+15 -2.787567024e-05 -2.759062723e-09 -9864501378 -8.281880969e-07 -1.309455569e-10 -4.481093592e+17 -1.114010181e-16 -8.97571489e+12 -1.601389437e-09 -5.858032102e-21 +2.22753933 +0.001917923831 +1.944614328e+15 +2.787872492e-05 +1.778752712e-08 +9858318342 +5.339966185e-06 +1.309633359e-10 +4.481099941e+17 +7.18142559e-16 +8.975743083e+12 +6.009135232e-10 +2.200891197e-21 6.963813204e+16 -1.415072594e-20 -1.996517884e-16 -7.694480337e-16 -1.691689432e-13 -5.040356928e-11 -5.482635383e-37 -1.604555857e-15 -5.910343021e-10 -1.715097579e-31 -6.415084542e-11 -1.001202566e-27 -2.10230157e-07 -2.226165405 -7.279438235e-38 -3.99809758e-32 -1.023291941e-05 -9.614519081e-36 -1.652332248e-58 +2.032448461e-18 +4.174154493e-16 +1.609285669e-15 +3.540277018e-13 +1.055530723e-10 +2.763957945e-36 +1.254725284e-15 +1.238523826e-09 +1.105069858e-30 +6.41122004e-11 +6.456231456e-27 +2.102757459e-07 +2.225642417 +2.245978095e-38 +1.232689556e-32 +1.509127756e-06 +2.967131355e-36 +5.094305713e-59 AMReX (22.07-875-gb05f69fe3748) finalized diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1e-1.out b/unit_test/burn_cell_metal_chem/reference_solution_1e-1.out index 5681bb98d..3a930088c 100644 --- a/unit_test/burn_cell_metal_chem/reference_solution_1e-1.out +++ b/unit_test/burn_cell_metal_chem/reference_solution_1e-1.out @@ -1,45 +1,45 @@ -T final = 2580.849045 +T final = 2580.84748 Tdust initial = 2.73 -Tdust final = 2580.903736 +Tdust final = 2580.896055 Eint initial = 3.074857008e+10 -Eint final = 2.252782311e+11 +Eint final = 2.25277789e+11 rho initial = 2.193512909e-25 -rho final = 1.989589323e-06 -numdens final = 5.416474951e+17 +rho final = 1.98964343e-06 +numdens final = 5.416473705e+17 ------------------------------------ New number densities: 8.412917854e+12 -2.39681895e+13 -2.207979967 -6.56484998e-06 -3.496528463e+16 -0.0001134569489 -5.932464766e-11 -1.640674045e+11 -9.26261325e-09 -4.948497742e-10 -4.363062397e+17 -3.520507679e-18 -8.912257632e+12 -4.109464144e-09 -6.788802573e-22 +2.396818951e+13 +2.208070971 +4.917822618e-05 +3.496505964e+16 +0.0001134618309 +4.44409566e-10 +1.640663651e+11 +6.938765994e-08 +4.948711095e-10 +4.363063401e+17 +2.637260563e-17 +8.912258695e+12 +9.51889453e-09 +1.613163821e-21 7.033451336e+16 -1.460565197e-19 -4.226130576e-11 -2.962304509e-11 -4.321097966e-10 -6.410593413e-09 -2.847622487e-35 -3.523568564e-12 -4.116188787e-09 -1.579040656e-30 -2.800477243e-08 -2.517542637e-28 -2.481297579e-06 -2.208086624 -1.820442242e-35 -2.493711987e-27 -2.172346472e-07 -1.03303842e-33 -1.308968162e-52 +4.313617275e-21 +5.114046531e-11 +3.584704032e-11 +5.229021966e-10 +7.757603344e-09 +2.133310512e-34 +3.523760341e-12 +4.981125586e-09 +1.182951876e-29 +2.800643178e-08 +1.886056751e-27 +2.481465418e-06 +2.208135 +1.690486453e-35 +2.334074296e-27 +1.667014038e-07 +9.596752642e-34 +1.22525698e-52 AMReX (22.07-589-g9b0d63e40412) finalized diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1e-2.out b/unit_test/burn_cell_metal_chem/reference_solution_1e-2.out index e6a2063ee..27b0abf28 100644 --- a/unit_test/burn_cell_metal_chem/reference_solution_1e-2.out +++ b/unit_test/burn_cell_metal_chem/reference_solution_1e-2.out @@ -1,45 +1,45 @@ -T final = 2811.207219 +T final = 2810.374696 Tdust initial = 2.73 -Tdust final = 2811.228114 +Tdust final = 2811.399316 Eint initial = 3.076181952e+10 -Eint final = 2.478723199e+11 +Eint final = 2.477868218e+11 rho initial = 2.192485838e-25 -rho final = 1.989661731e-06 -numdens final = 5.670189018e+17 +rho final = 1.989668195e-06 +numdens final = 5.668962568e+17 ------------------------------------ New number densities: -8.412914755e+11 +8.412914754e+11 2.396807474e+12 -2.236235976 -6.405982886e-05 -8.572287708e+16 -0.0001611585155 -5.332298561e-10 -4.037425745e+11 -6.855532471e-08 -7.073005674e-10 -4.109491978e+17 -3.761054318e-17 -8.671678263e+12 -1.456049659e-10 -1.650308794e-21 +2.237388464 +6.169419758e-07 +8.547697113e+16 +0.0001610647153 +5.137205748e-12 +4.02556912e+11 +6.608419452e-10 +7.068343694e-10 +4.110724588e+17 +3.621113956e-19 +8.672863935e+12 +6.173898911e-08 +3.888219294e-21 7.033451336e+16 -2.473602851e-09 -0.000502828669 -0.000195826651 -0.001167253794 -0.006497881953 -1.865521147e-32 -2.352999816e-10 -0.001602637809 -5.537533357e-28 -9.719529679e-07 -2.822461066e-26 -2.69412668e-05 -2.131195609 -4.275851205e-23 -2.534923404e-13 -0.1051373945 -1.86318788e-21 -8.869168988e-21 +2.542238142e-09 +0.0004989372543 +0.0001946900676 +0.001163928336 +0.00650045374 +1.786797422e-34 +2.340063332e-10 +0.001608372808 +5.313130942e-30 +9.68531334e-07 +2.718197375e-28 +2.694990265e-05 +2.129003659 +4.392878947e-23 +-4.560006446e-12 +0.1085451266 +2.537886844e-21 +-4.677302387e-19 AMReX (22.07-589-g9b0d63e40412) finalized diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1e-3.out b/unit_test/burn_cell_metal_chem/reference_solution_1e-3.out index da53a37ac..5387e4fc2 100644 --- a/unit_test/burn_cell_metal_chem/reference_solution_1e-3.out +++ b/unit_test/burn_cell_metal_chem/reference_solution_1e-3.out @@ -1,45 +1,45 @@ -T final = 2831.090224 +T final = 2831.019478 Tdust initial = 2.73 -Tdust final = 2831.289176 +Tdust final = 2831.686588 Eint initial = 3.076314515e+10 -Eint final = 2.499224717e+11 +Eint final = 2.499151417e+11 rho initial = 2.192383131e-25 -rho final = 1.989668738e-06 -numdens final = 5.700271327e+17 +rho final = 1.98966925e-06 +numdens final = 5.700162122e+17 ------------------------------------ New number densities: -8.412915268e+10 +8.412915269e+10 2.396823174e+11 -4.033710506 -6.777594776e-05 -9.174097417e+16 -0.0002983374384 -5.593606493e-10 -4.328369028e+11 -7.097532349e-08 -1.311908161e-09 -4.079422459e+17 -4.006694765e-17 -8.64258412e+12 -2.245209778e-14 -2.215435849e-21 +4.033631317 +7.268898769e-05 +9.171911853e+16 +0.0002983039113 +5.999266876e-10 +4.327309304e+11 +7.612608637e-08 +1.311751089e-09 +4.07953181e+17 +4.297040753e-17 +8.642689996e+12 +2.354952341e-10 +2.398645645e-21 7.033451336e+16 -1.304995286e-08 -8.203048039e-05 -3.050552029e-05 -0.0001694822058 -0.0008735603933 -5.43479155e-31 -6.483863045e-09 -0.0001998289953 -1.534470988e-26 -2.532177639e-05 -7.160676166e-25 -0.0006407136702 -0.4462317747 -2.224187897e-22 -8.810161969e-12 -3.587709209 -9.496345982e-21 -1.037846593e-34 +1.310908423e-08 +8.203187782e-05 +3.051098989e-05 +0.0001695545589 +0.0008741705667 +5.825557904e-31 +6.480280989e-09 +0.0002000217592 +1.645050698e-26 +2.53122014e-05 +7.67908607e-25 +0.0006406765401 +0.4462423571 +2.229330167e-22 +-6.520895357e-12 +3.587614486 +-2.492556219e-18 +-4.418583788e-19 AMReX (22.07-589-g9b0d63e40412) finalized diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1e-4.out b/unit_test/burn_cell_metal_chem/reference_solution_1e-4.out index 02948a346..9b01e0cfc 100644 --- a/unit_test/burn_cell_metal_chem/reference_solution_1e-4.out +++ b/unit_test/burn_cell_metal_chem/reference_solution_1e-4.out @@ -1,45 +1,45 @@ -T final = 2884.322233 +T final = 2884.408445 Tdust initial = 2.73 -Tdust final = 2885.222529 +Tdust final = 2885.384517 Eint initial = 3.076327772e+10 -Eint final = 2.555031196e+11 +Eint final = 2.555122564e+11 rho initial = 2.19237286e-25 -rho final = 1.989670334e-06 -numdens final = 5.787722231e+17 +rho final = 1.98966997e-06 +numdens final = 5.787870732e+17 ------------------------------------ New number densities: -8412915325 +8412915330 2.396818542e+10 -5.962439927 -0.0005137258492 -1.092306741e+17 -0.0004716327583 -4.142389741e-09 -5.183837544e+11 -5.083763652e-07 -2.087320476e-09 -3.991979278e+17 -3.093116661e-16 -8.557037312e+12 -4.381451187e-08 -3.7384033e-20 +5.962823224 +0.001071375051 +1.092606255e+17 +0.000471712832 +8.638627485e-09 +5.185317229e+11 +1.060123122e-06 +2.087700204e-09 +3.991828266e+17 +6.450887085e-16 +8.556889344e+12 +2.412076997e-07 +9.717358293e-20 7.033451336e+16 -1.431469116e-07 -3.125096897e-05 -1.030321161e-05 -4.764610267e-05 -0.000200901188 -1.246503729e-28 -1.967791718e-07 -3.775902673e-05 -3.167115752e-24 -0.000678373265 -1.175791827e-22 -0.01353231898 -0.03355686142 -1.176079223e-22 -5.66867762e-11 -5.928840231 -4.761295975e-21 -1.381019684e-32 +1.420660126e-07 +3.101380105e-05 +1.022306037e-05 +4.726165595e-05 +0.0001992165272 +2.601230644e-28 +1.969050239e-07 +3.743069657e-05 +6.608038109e-24 +0.0006786635515 +2.452346741e-22 +0.01353299323 +0.03354856481 +1.166801e-22 +5.625847973e-11 +5.928673305 +4.620188648e-21 +1.370612271e-32 AMReX (22.07-589-g9b0d63e40412) finalized diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1e-5.out b/unit_test/burn_cell_metal_chem/reference_solution_1e-5.out index 69a2fca45..8304a7591 100644 --- a/unit_test/burn_cell_metal_chem/reference_solution_1e-5.out +++ b/unit_test/burn_cell_metal_chem/reference_solution_1e-5.out @@ -1,45 +1,45 @@ -T final = 2770.253479 +T final = 2770.260502 Tdust initial = 2.73 -Tdust final = 2770.396907 +Tdust final = 2771.031412 Eint initial = 3.076329098e+10 -Eint final = 2.437024468e+11 +Eint final = 2.437031597e+11 rho initial = 2.192371833e-25 -rho final = 1.989670057e-06 -numdens final = 5.612534538e+17 +rho final = 1.989669807e-06 +numdens final = 5.61254438e+17 ------------------------------------ New number densities: -199898176.6 -3038212036 -2.259991418 -0.09252421839 -7.419371216e+16 -0.000154168177 -7.834842422e-07 -3.484514097e+11 -0.0001036533126 -6.744074981e-10 -4.167161489e+17 -5.354504797e-14 -8.726969668e+12 -3.849333801e-12 -1.421503406e-18 +199897170.8 +3038213042 +2.203629741 +0.0003491632575 +7.419566308e+16 +0.0001503269971 +2.956664728e-09 +3.484607261e+11 +3.91158392e-07 +6.576046293e-10 +4.167151823e+17 +2.020661145e-16 +8.726960346e+12 +5.643506785e-09 +5.771775307e-21 7.033451336e+16 -0.00281312011 -21789165.94 -9352664.009 -64650614.55 -423270591.2 -1.537707812e-26 -1.341686703e-07 -122330252.1 -5.051058737e-22 -0.0006197903089 -3.102910261e-20 -0.02081899296 -2.164703811 -4.251138744e-23 -0.00765019574 -1.387355782e-34 -1.697403972e-27 -3.012881382e-24 +1.063497408e-05 +21790854.1 +9353212.954 +64652603.8 +423271205.5 +5.758672328e-29 +1.331455347e-07 +122327033.3 +1.891574487e-24 +0.0006150522586 +1.161974271e-22 +0.02065917213 +2.203419865 +1.59480483e-25 +0.007594234715 +-3.519279316e-31 +-3.283812312e-21 +2.967855969e-24 AMReX (22.07-589-g9b0d63e40412) finalized diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1e-6.out b/unit_test/burn_cell_metal_chem/reference_solution_1e-6.out index 372f204db..e9152d6c8 100644 --- a/unit_test/burn_cell_metal_chem/reference_solution_1e-6.out +++ b/unit_test/burn_cell_metal_chem/reference_solution_1e-6.out @@ -1,45 +1,45 @@ -T final = 2812.050392290817 +T final = 2813.763186 Tdust initial = 2.73 -Tdust final = 2812.235837368595 -Eint initial = 30763292303.77594 -Eint final = 247958162899.1352 -rho initial = 2.192371730563381e-25 -rho final = 1.989669757431512e-06 -numdens final = 5.671429626664702e+17 +Tdust final = 2813.526575 +Eint initial = 3.07632923e+10 +Eint final = 2.481342021e+11 +rho initial = 2.192371731e-25 +rho final = 1.98967001e-06 +numdens final = 5.673973163e+17 ------------------------------------ -New number densities: -6438308.410300142 -317322562.9538858 -2.486641265170912 -2.274499630590923e-09 -8.597266896217706e+16 -0.0001794039501637353 -1.8925886867936e-14 -404944748318.3216 -2.431869387308584e-12 -7.874340701282926e-10 -4.108267045200611e+17 -1.335778922651524e-21 -8670476376287.059 -3.501410855685503e-14 -6.130614268934865e-26 -7.033451336163859e+16 -2.268083235873295e-12 -3948863.726870825 -1534860.916182474 -9121355.520456472 -50610347.48632588 -4.088062980895315e-28 -0.1452240345710818 -12442512.37027289 -1.216206468005654e-23 -601.0948939083418 -6.175692426101071e-22 -16596.84230665624 -2.486820667629221 -2.85628902206954e-26 -32951.98449995967 -1.632493779270958e-29 -1.471009580362413e-15 -1.014011717966889e-11 +New number densities: +6472049.128 +317338389.7 +2.473688688 +2.437815373e-09 +8.648123216e+16 +0.0001788715759 +2.026991226e-14 +4.073976485e+11 +2.601564064e-12 +7.85222469e-10 +4.10572495e+17 +1.432546932e-21 +8.668023467e+12 +5.262395419e-09 +3.603849018e-22 +7.033451336e+16 +2.490210408e-12 +4012171.391 +1553236.62 +9174440.11 +50566325.55 +4.2605529e-28 +0.1412494562 +12350866.07 +1.262582542e-23 +582.0669493 +6.362171187e-22 +15945.26769 +2.473867548 +3.01732562e-26 +32833.01421 +1.610264744e-29 +1.375039614e-15 +9.699753847e-12 AMReX (22.07-589-g9b0d63e40412) finalized From 445133212f8d82335e87ef4afc7b2af6be2325aa Mon Sep 17 00:00:00 2001 From: psharda Date: Tue, 17 Sep 2024 18:02:04 +0200 Subject: [PATCH 3/7] update ref file lines --- .github/workflows/burn_cell_metal_chem.yml | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/.github/workflows/burn_cell_metal_chem.yml b/.github/workflows/burn_cell_metal_chem.yml index 03386c464..7f2f18fa2 100644 --- a/.github/workflows/burn_cell_metal_chem.yml +++ b/.github/workflows/burn_cell_metal_chem.yml @@ -39,7 +39,7 @@ jobs: cd unit_test/burn_cell_metal_chem declare -A line_numbers_map=( - ["Z=1"]="4461 4463 4465 4467 4468 4471 4472 4475 4481 4486" + ["Z=1"]="4450 4452 4454 4456 4457 4460 4461 4464 4470 4476" ["Z=1_z10"]="1613 1615 4456 4458 4460 4463 4470 4476 4481" ["Z=1e-1"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" ["Z=1e-2"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" @@ -92,7 +92,7 @@ jobs: # Adjust the line number for the reference file if [[ "$Z" == "Z=1" ]]; then - reference_line_number=$((line_number - 4460)) + reference_line_number=$((line_number - 4449)) elif [[ "$Z" == "Z=1_z10" ]]; then reference_line_number=${ref_line_number_z10[$index]} else From 7b6a6af48240cf0193bf6cafdec46823aa364888 Mon Sep 17 00:00:00 2001 From: psharda Date: Tue, 17 Sep 2024 18:06:22 +0200 Subject: [PATCH 4/7] update ref file lines --- .github/workflows/burn_cell_metal_chem.yml | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/.github/workflows/burn_cell_metal_chem.yml b/.github/workflows/burn_cell_metal_chem.yml index 7f2f18fa2..78645fc11 100644 --- a/.github/workflows/burn_cell_metal_chem.yml +++ b/.github/workflows/burn_cell_metal_chem.yml @@ -39,7 +39,7 @@ jobs: cd unit_test/burn_cell_metal_chem declare -A line_numbers_map=( - ["Z=1"]="4450 4452 4454 4456 4457 4460 4461 4464 4470 4476" + ["Z=1"]="4450 4452 4454 4456 4457 4460 4461 4464 4470 4475" ["Z=1_z10"]="1613 1615 4456 4458 4460 4463 4470 4476 4481" ["Z=1e-1"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" ["Z=1e-2"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" From ceb00c0359bc5295aaf608d3415edd1e2a890f64 Mon Sep 17 00:00:00 2001 From: psharda Date: Tue, 17 Sep 2024 18:14:15 +0200 Subject: [PATCH 5/7] update ref file lines --- .github/workflows/burn_cell_metal_chem.yml | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/.github/workflows/burn_cell_metal_chem.yml b/.github/workflows/burn_cell_metal_chem.yml index 78645fc11..b23999ff5 100644 --- a/.github/workflows/burn_cell_metal_chem.yml +++ b/.github/workflows/burn_cell_metal_chem.yml @@ -46,7 +46,7 @@ jobs: ["Z=1e-3"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" ["Z=1e-4"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" ["Z=1e-5"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" - ["Z=1e-6"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" + ["Z=1e-6"]="4441 4453 4455 4457 4458 4461 4462 4465 4471 4476" ) declare -A ref_map=( @@ -95,6 +95,8 @@ jobs: reference_line_number=$((line_number - 4449)) elif [[ "$Z" == "Z=1_z10" ]]; then reference_line_number=${ref_line_number_z10[$index]} + elif [[ "$Z" == "Z=1e-6" ]]; then + reference_line_number=$((line_number - 4440)) else reference_line_number=$((line_number - 4437)) fi From d0665117139a18fdf533375616cc0e5c6787fa53 Mon Sep 17 00:00:00 2001 From: psharda Date: Tue, 17 Sep 2024 22:52:23 +0200 Subject: [PATCH 6/7] updated ref soln --- .github/workflows/burn_cell_metal_chem.yml | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/.github/workflows/burn_cell_metal_chem.yml b/.github/workflows/burn_cell_metal_chem.yml index b23999ff5..193447266 100644 --- a/.github/workflows/burn_cell_metal_chem.yml +++ b/.github/workflows/burn_cell_metal_chem.yml @@ -46,7 +46,7 @@ jobs: ["Z=1e-3"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" ["Z=1e-4"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" ["Z=1e-5"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" - ["Z=1e-6"]="4441 4453 4455 4457 4458 4461 4462 4465 4471 4476" + ["Z=1e-6"]="4451 4453 4455 4457 4458 4461 4462 4465 4471 4476" ) declare -A ref_map=( From a33d5418ccb5b9b212176434ed9e503a7e4abd54 Mon Sep 17 00:00:00 2001 From: psharda Date: Tue, 17 Sep 2024 23:21:53 +0200 Subject: [PATCH 7/7] updated ref soln --- .github/workflows/burn_cell_metal_chem.yml | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/.github/workflows/burn_cell_metal_chem.yml b/.github/workflows/burn_cell_metal_chem.yml index 193447266..0631ca629 100644 --- a/.github/workflows/burn_cell_metal_chem.yml +++ b/.github/workflows/burn_cell_metal_chem.yml @@ -96,7 +96,7 @@ jobs: elif [[ "$Z" == "Z=1_z10" ]]; then reference_line_number=${ref_line_number_z10[$index]} elif [[ "$Z" == "Z=1e-6" ]]; then - reference_line_number=$((line_number - 4440)) + reference_line_number=$((line_number - 4450)) else reference_line_number=$((line_number - 4437)) fi