@@ -129,27 +129,26 @@ TEST( testKineticReactions, computeReactionRatesTest_carbonateSystem )
129129 3.76e-1 , // HCO3-
130130 3.87e-2 , // Ca+2
131131 3.21e-2 , // SO4-2
132- 1.89 , // Cl-
132+ 1.89 , // Cl-
133133 1.65e-2 , // Mg+2
134- 1.09 // Na+1
134+ 1.09 // Na+1
135135 };
136136
137- double const expectedReactionRates[11 ] = { -9.9632648e-01 , -1.41376e-01 , -3.7599919536e-01 , -1.41376e-01 , -1.45512e-02 , -1.455119956008e-02 , -1.24227e-03 , -7.3143e-02 , -1.3824027e-01 ,
138- -5.2965000e-04 , -3.4989e-02 };
137+ double const expectedReactionRates[11 ] = { -0.001424736 ,-12610.7392 ,-0.175591624 ,-473.6096 ,-269197.19999999984 ,-0.00012441275624000003 ,-18012.914999999986 ,-1.56526019999999e6 ,-346983.07769999903 ,-14247.58499999999 ,-2.316271799999999e6 };
139138 double const expectedReactionRatesDerivatives[11 ][18 ] =
140139 {
141- { 3.67352e+13 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 9.77e-03 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0. 0 },
142- { 0.0 , 4.37e-07 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , - 3.76e-01 , - 3.76e-01 , 0.0 , 0.0 , 0.0 , 0.0 , 0. 0 },
143- { 0.0 , 0.0 , 8.0464e+09 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 2.14e-06 , - 1.00e+00 , 0.0 , 0.0 , 0.0 , 0.0 , 0. 0 },
144- { 0.0 , 0.0 , 0.0 , 1.70e-04 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , - 3.76e-01 , - 3.76e-01 , 0.0 , 0.0 , 0.0 , 0.0 , 0. 0 },
145- { 0.0 , 0.0 , 0.0 , 0.0 , 8.13e-02 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , - 3.87e-02 , - 3.76e-01 , 0.0 , 0.0 , 0.0 , 0. 0 },
146- { 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 4.3992e+06 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 1.17e-09 , - 3.87e-02 , - 3.76e-01 , 0.0 , 0.0 , 0.0 , 0. 0 },
147- { 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 6.92e-03 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , - 3.21e-02 , - 3.87e-02 , 0.0 , 0.0 , 0. 0 },
148- { 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 4.68e+00 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , - 1.89e+00 , 0.0 , - 3.87e-02 , 0.0 , 0. 0 },
149- { 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 3.98e+00 , 0.0 , 0.0 , 0.0 , 0.0 , - 3.5721e+00 , 0.0 , - 1.46286e-01 , 0.0 , 0. 0 },
150- { 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 3.72e-03 , 0.0 , 0.0 , 0.0 , 0.0 , - 1.65e-02 , 0.0 , - 3.21e-02 , 0. 0 },
151- { 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 1.51e-01 , 0.0 , 0.0 , 0.0 , - 1.09e+00 , 0.0 , 0.0 , - 3.21e-02 }
152- };
140+ { 5.264e10 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0.000014 , 0 , 0 , 0 , 0 , 0 , 0 },
141+ { 0 , 0.039 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ,- 33539.2 ,- 33539.2 , 0 , 0 , 0 , 0 , 0 },
142+ { 0 , 0 , 3.76e9 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 . e - 6 ,- 0.467 , 0 , 0 , 0 , 0 , 0 },
143+ { 0 , 0 , 0 , 0.57 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ,- 1259.6 ,- 1259.6 , 0 , 0 , 0 , 0 , 0 },
144+ { 0 , 0 , 0 , 0 , 1.5e6 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ,- 715950 .,- 6.956e6 , 0 , 0 , 0 , 0 },
145+ { 0 , 0 , 0 , 0 , 0 , 37600 ., 0 , 0 , 0 , 0 , 0 , 1 . e - 11 ,- 0.000330885 ,- 0.0032148000000000003 , 0 , 0 , 0 , 0 },
146+ { 0 , 0 , 0 , 0 , 0 , 0 , 100000 ., 0 , 0 , 0 , 0 , 0 , 0 ,- 465449.99999999994 ,- 561150 ., 0 , 0 , 0 },
147+ { 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 . e8 , 0 , 0 , 0 , 0 , 0 ,- 4.0446e7 , 0 ,- 828180 ., 0 , 0 },
148+ { 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 . e7 , 0 , 0 , 0 , 0 ,- 8.965971e6 , 0 ,- 367177.86 , 0 , 0 },
149+ { 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 100000 ., 0 , 0 , 0 , 0 ,- 443850 ., 0 ,- 863489.9999999999 , 0 },
150+ { 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 . e7 , 0 , 0 , 0 ,- 7.2158e7 , 0 , 0 ,- 2.12502e6 }
151+ } ;
153152
154153 computeReactionRatesTest< double , false >( carbonateSystem,
155154 initialSpeciesConcentration,
0 commit comments