Rebaselined the problem. There were minor changes in the
solution. Moved the test to a csv comparison which can differentiate between large and small differences.
This commit is contained in:
parent
a5cb5b7f8a
commit
6a85d86b78
5 changed files with 190 additions and 99 deletions
|
|
@ -7,3 +7,6 @@ ct2ctml.log
|
|||
diff_test.out
|
||||
.depends
|
||||
*.d
|
||||
diff_csv.txt
|
||||
diff_out.txt
|
||||
table.csv
|
||||
|
|
|
|||
|
|
@ -46,6 +46,18 @@ int main(int argc, char **argv) {
|
|||
IndVar1[6] = 1.0E-2;
|
||||
int nj = 6;
|
||||
int ni = 7;
|
||||
FILE *FF = fopen("table.csv","w");
|
||||
int kk = gas->nSpecies();
|
||||
std::vector<double> Xmol(kk, 0.0);
|
||||
const std::vector<string> &snames = gas->speciesNames();
|
||||
fprintf(FF,"Temperature, Pressure,");
|
||||
for (int k = 0; k < kk; k++) {
|
||||
fprintf(FF, "%11s", snames[k].c_str());
|
||||
if (k < kk-1) {
|
||||
fprintf(FF,",");
|
||||
}
|
||||
}
|
||||
fprintf(FF,"\n");
|
||||
|
||||
for (int j=0; j<nj; j++) {
|
||||
for (int i=0; i<ni; i++) {
|
||||
|
|
@ -62,10 +74,23 @@ int main(int argc, char **argv) {
|
|||
tkelvin = gas->temperature();
|
||||
pres = gas->pressure();
|
||||
printf("Final T = %g, pres = %g atm\n", tkelvin, pres/OneAtm);
|
||||
|
||||
gas->getMoleFractions(DATA_PTR(Xmol));
|
||||
fprintf(FF,"%10.4g, %10.4g,", tkelvin, pres);
|
||||
for (int k = 0; k < kk; k++) {
|
||||
if (fabs(Xmol[k]) < 1.0E-130) {
|
||||
fprintf(FF," %10.4g", 0.0);
|
||||
} else {
|
||||
fprintf(FF," %10.4g", Xmol[k]);
|
||||
}
|
||||
if (k < kk-1) {
|
||||
fprintf(FF,",");
|
||||
}
|
||||
}
|
||||
fprintf(FF,"\n");
|
||||
}
|
||||
}
|
||||
delete gas;
|
||||
fclose(FF);
|
||||
}
|
||||
|
||||
catch (CanteraError) {
|
||||
|
|
|
|||
|
|
@ -89,7 +89,7 @@ Initial T = 207.284, pres = 5.8944e-07 atm
|
|||
N 9.55809e-114 4.63943e-114 -18.9446
|
||||
O 1.18094e-57 6.54771e-58 -20.4172
|
||||
NO 1.82556e-23 1.89829e-23 -39.3618
|
||||
NO+ 2.2952e-228 2.38659e-228 12.2863
|
||||
NO+ 2.2952e-228 2.3866e-228 12.2863
|
||||
Electron 2.2952e-228 4.33484e-233 -541.209
|
||||
N+ 0 0
|
||||
O+ 0 0
|
||||
|
|
@ -122,8 +122,8 @@ Initial T = 207.284, pres = 5.8944e-06 atm
|
|||
N 3.02253e-114 1.46712e-114 -17.7933
|
||||
O 3.73447e-58 2.07057e-58 -19.266
|
||||
NO 1.82556e-23 1.89829e-23 -37.0592
|
||||
NO+ 1.2588e-229 1.30893e-229 11.6857
|
||||
Electron 1.2588e-229 2.37745e-234 -541.81
|
||||
NO+ 1.25881e-229 1.30893e-229 11.6857
|
||||
Electron 1.25881e-229 2.37745e-234 -541.81
|
||||
N+ 0 0
|
||||
O+ 0 0
|
||||
N2+ 0 0
|
||||
|
|
@ -188,8 +188,8 @@ Initial T = 207.284, pres = 0.00058944 atm
|
|||
N 3.02253e-115 1.46712e-115 -15.4907
|
||||
O 3.73447e-59 2.07057e-59 -16.9634
|
||||
NO 1.82556e-23 1.89829e-23 -32.454
|
||||
NO+ 1.26648e-231 1.31691e-231 11.6917
|
||||
Electron 1.26648e-231 2.39194e-236 -541.804
|
||||
NO+ 1.2656e-231 1.316e-231 11.691
|
||||
Electron 1.2656e-231 2.39028e-236 -541.804
|
||||
N+ 0 0
|
||||
O+ 0 0
|
||||
N2+ 0 0
|
||||
|
|
@ -221,8 +221,8 @@ Initial T = 207.284, pres = 0.0058944 atm
|
|||
N 9.55809e-116 4.63943e-116 -14.3394
|
||||
O 1.18094e-59 6.54771e-60 -15.8121
|
||||
NO 1.82556e-23 1.89829e-23 -30.1515
|
||||
NO+ 2.40463e-232 2.50039e-232 12.3329
|
||||
Electron 2.40463e-232 4.54152e-237 -541.163
|
||||
NO+ 2.40459e-232 2.50035e-232 12.3329
|
||||
Electron 2.40459e-232 4.54145e-237 -541.163
|
||||
N+ 0 0
|
||||
O+ 0 0
|
||||
N2+ 0 0
|
||||
|
|
@ -254,12 +254,12 @@ Initial T = 414.143, pres = 1.17767e-08 atm
|
|||
N 2.06855e-53 1.00406e-53 -20.8551
|
||||
O 1.8737e-25 1.03887e-25 -22.3285
|
||||
NO 5.58575e-12 5.80828e-12 -43.1837
|
||||
NO+ 4.00599e-102 4.16552e-102 11.996
|
||||
Electron 4.00599e-102 7.56595e-107 -254.374
|
||||
N+ 7.16723e-207 3.47879e-207 34.3246
|
||||
O+ 3.0296e-168 1.6797e-168 32.8512
|
||||
N2+ 6.19807e-167 6.01689e-167 13.4694
|
||||
O2+ 3.22979e-125 3.58144e-125 10.5226
|
||||
NO+ 4.00605e-102 4.16558e-102 11.996
|
||||
Electron 4.00605e-102 7.56606e-107 -254.374
|
||||
N+ 7.16734e-207 3.47884e-207 34.3246
|
||||
O+ 3.02965e-168 1.67972e-168 32.8512
|
||||
N2+ 6.19817e-167 6.01698e-167 13.4694
|
||||
O2+ 3.22984e-125 3.58149e-125 10.5226
|
||||
|
||||
Final T = 414.143, pres = 1.17767e-08 atm
|
||||
Initial T = 414.143, pres = 1.17767e-07 atm
|
||||
|
|
@ -287,12 +287,12 @@ Initial T = 414.143, pres = 1.17767e-07 atm
|
|||
N 6.54134e-54 3.17512e-54 -19.7038
|
||||
O 5.92517e-26 3.28519e-26 -21.1772
|
||||
NO 5.58575e-12 5.80828e-12 -40.8811
|
||||
NO+ 5.13663e-103 5.34118e-103 12.2446
|
||||
Electron 5.13663e-103 9.70133e-108 -254.125
|
||||
N+ 2.90616e-208 1.41057e-208 33.4219
|
||||
O+ 1.22844e-169 6.81081e-170 31.9485
|
||||
N2+ 7.9474e-168 7.71507e-168 13.718
|
||||
O2+ 4.14136e-126 4.59225e-126 10.7712
|
||||
NO+ 5.13662e-103 5.34116e-103 12.2446
|
||||
Electron 5.13662e-103 9.7013e-108 -254.125
|
||||
N+ 2.90615e-208 1.41057e-208 33.4219
|
||||
O+ 1.22844e-169 6.81079e-170 31.9485
|
||||
N2+ 7.94737e-168 7.71505e-168 13.718
|
||||
O2+ 4.14134e-126 4.59223e-126 10.7712
|
||||
|
||||
Final T = 414.143, pres = 1.17767e-07 atm
|
||||
Initial T = 414.143, pres = 1.17767e-06 atm
|
||||
|
|
@ -320,12 +320,12 @@ Initial T = 414.143, pres = 1.17767e-06 atm
|
|||
N 2.06855e-54 1.00406e-54 -18.5525
|
||||
O 1.8737e-26 1.03887e-26 -20.026
|
||||
NO 5.58575e-12 5.80828e-12 -38.5785
|
||||
NO+ 2.82554e-104 2.93805e-104 11.6469
|
||||
Electron 2.82554e-104 5.33647e-109 -254.723
|
||||
N+ 5.05525e-210 2.45369e-210 31.6729
|
||||
O+ 2.13686e-171 1.18474e-171 30.1995
|
||||
N2+ 4.37167e-169 4.24387e-169 13.1203
|
||||
O2+ 2.27806e-127 2.52609e-127 10.1735
|
||||
NO+ 2.82555e-104 2.93806e-104 11.6469
|
||||
Electron 2.82555e-104 5.33649e-109 -254.723
|
||||
N+ 5.05526e-210 2.45369e-210 31.6729
|
||||
O+ 2.13687e-171 1.18474e-171 30.1995
|
||||
N2+ 4.37169e-169 4.24389e-169 13.1203
|
||||
O2+ 2.27807e-127 2.52609e-127 10.1735
|
||||
|
||||
Final T = 414.143, pres = 1.17767e-06 atm
|
||||
Initial T = 414.143, pres = 1.17767e-05 atm
|
||||
|
|
@ -353,12 +353,12 @@ Initial T = 414.143, pres = 1.17767e-05 atm
|
|||
N 6.54134e-55 3.17512e-55 -17.4013
|
||||
O 5.92517e-27 3.28519e-27 -18.8747
|
||||
NO 5.58575e-12 5.80828e-12 -36.2759
|
||||
NO+ 3.94044e-105 4.09735e-105 11.9795
|
||||
Electron 3.94044e-105 7.44214e-110 -254.39
|
||||
NO+ 3.94045e-105 4.09736e-105 11.9795
|
||||
Electron 3.94045e-105 7.44215e-110 -254.39
|
||||
N+ 2.22939e-211 1.08209e-211 30.8542
|
||||
O+ 9.42367e-173 5.22475e-173 29.3808
|
||||
N2+ 6.09665e-170 5.91843e-170 13.4529
|
||||
O2+ 3.17694e-128 3.52283e-128 10.5061
|
||||
O+ 9.42368e-173 5.22475e-173 29.3808
|
||||
N2+ 6.09666e-170 5.91843e-170 13.4529
|
||||
O2+ 3.17694e-128 3.52284e-128 10.5061
|
||||
|
||||
Final T = 414.143, pres = 1.17767e-05 atm
|
||||
Initial T = 414.143, pres = 0.000117767 atm
|
||||
|
|
@ -390,7 +390,7 @@ Initial T = 414.143, pres = 0.000117767 atm
|
|||
Electron 2.91676e-106 5.50875e-111 -254.691
|
||||
N+ 5.21845e-213 2.5329e-213 29.4021
|
||||
O+ 2.20585e-174 1.22298e-174 27.9287
|
||||
N2+ 4.51281e-171 4.38089e-171 13.1521
|
||||
N2+ 4.51281e-171 4.38088e-171 13.1521
|
||||
O2+ 2.35161e-129 2.60764e-129 10.2053
|
||||
|
||||
Final T = 414.143, pres = 0.000117767 atm
|
||||
|
|
@ -419,12 +419,12 @@ Initial T = 414.143, pres = 0.00117767 atm
|
|||
N 6.54134e-56 3.17512e-56 -15.0987
|
||||
O 5.92517e-28 3.28519e-28 -16.5721
|
||||
NO 5.58575e-12 5.80828e-12 -31.6707
|
||||
NO+ 3.85738e-107 4.01099e-107 11.9582
|
||||
Electron 3.85738e-107 7.28527e-112 -254.411
|
||||
N+ 2.1824e-214 1.05928e-214 28.5303
|
||||
O+ 9.22503e-176 5.11462e-176 27.0569
|
||||
N2+ 5.96814e-172 5.79367e-172 13.4316
|
||||
O2+ 3.10997e-130 3.44857e-130 10.4848
|
||||
NO+ 3.85557e-107 4.00911e-107 11.9577
|
||||
Electron 3.85557e-107 7.28185e-112 -254.412
|
||||
N+ 2.18137e-214 1.05878e-214 28.5298
|
||||
O+ 9.22071e-176 5.11222e-176 27.0564
|
||||
N2+ 5.96534e-172 5.79096e-172 13.4311
|
||||
O2+ 3.10852e-130 3.44696e-130 10.4843
|
||||
|
||||
Final T = 414.143, pres = 0.00117767 atm
|
||||
Initial T = 414.143, pres = 0.0117767 atm
|
||||
|
|
@ -452,12 +452,12 @@ Initial T = 414.143, pres = 0.0117767 atm
|
|||
N 2.06855e-56 1.00406e-56 -13.9474
|
||||
O 1.8737e-28 1.03887e-28 -15.4208
|
||||
NO 5.58575e-12 5.80828e-12 -29.3682
|
||||
NO+ 6.14689e-108 6.39166e-108 12.4242
|
||||
Electron 6.14689e-108 1.16094e-112 -253.946
|
||||
N+ 1.09976e-215 5.33793e-216 27.845
|
||||
O+ 4.64869e-177 2.57736e-177 26.3716
|
||||
N2+ 9.51046e-173 9.23244e-173 13.8976
|
||||
O2+ 4.95586e-131 5.49544e-131 10.9508
|
||||
NO+ 2.85055e-102 2.96406e-102 25.4712
|
||||
Electron 2.85055e-102 5.38371e-107 -240.898
|
||||
N+ 5.1e-210 2.47541e-210 40.892
|
||||
O+ 2.15578e-171 1.19522e-171 39.4186
|
||||
N2+ 4.41037e-167 4.28145e-167 26.9447
|
||||
O2+ 2.29823e-125 2.54845e-125 23.9978
|
||||
|
||||
Final T = 414.143, pres = 0.0117767 atm
|
||||
Initial T = 1153.16, pres = 3.27918e-08 atm
|
||||
|
|
@ -485,12 +485,12 @@ Initial T = 1153.16, pres = 3.27918e-08 atm
|
|||
N 2.86312e-15 1.38975e-15 -21.3563
|
||||
O 2.23425e-05 1.23879e-05 -22.8775
|
||||
NO 0.000134494 0.000139854 -44.2338
|
||||
NO+ 9.03856e-22 9.39859e-22 11.6755
|
||||
Electron 9.03856e-22 1.70709e-26 -69.7353
|
||||
N+ 1.59896e-54 7.76101e-55 34.553
|
||||
O+ 3.05445e-41 1.69349e-41 33.0319
|
||||
N2+ 1.21258e-44 1.17715e-44 13.1967
|
||||
O2+ 3.22543e-30 3.57664e-30 10.1544
|
||||
NO+ 7.02819e-19 7.30815e-19 18.3317
|
||||
Electron 7.02819e-19 1.3274e-23 -63.0791
|
||||
N+ 1.24332e-51 6.0348e-52 41.2092
|
||||
O+ 2.37508e-38 1.31683e-38 39.688
|
||||
N2+ 9.42879e-42 9.15326e-42 19.8529
|
||||
O2+ 2.50803e-27 2.78112e-27 16.8106
|
||||
|
||||
Final T = 1152.46, pres = 3.27723e-08 atm
|
||||
Initial T = 1152.46, pres = 3.27723e-07 atm
|
||||
|
|
@ -518,12 +518,12 @@ Initial T = 1152.46, pres = 3.27723e-07 atm
|
|||
N 9.11146e-16 4.42265e-16 -20.2051
|
||||
O 7.08903e-06 3.93051e-06 -21.7263
|
||||
NO 0.000134663 0.000140028 -41.9314
|
||||
NO+ 9.05883e-23 9.4196e-23 11.6643
|
||||
Electron 9.05883e-23 1.71091e-27 -69.7332
|
||||
N+ 5.12845e-56 2.48922e-56 33.3906
|
||||
O+ 9.75583e-43 5.40893e-43 31.8694
|
||||
N2+ 1.22376e-45 1.18799e-45 13.1855
|
||||
O2+ 3.24058e-31 3.59342e-31 10.1431
|
||||
NO+ 9.67397e-23 1.00592e-22 11.73
|
||||
Electron 9.67398e-23 1.82709e-27 -69.6675
|
||||
N+ 5.4767e-56 2.65826e-56 33.4563
|
||||
O+ 1.04183e-42 5.77622e-43 31.9351
|
||||
N2+ 1.30686e-45 1.26866e-45 13.2512
|
||||
O2+ 3.46064e-31 3.83743e-31 10.2088
|
||||
|
||||
Final T = 1152.61, pres = 3.27763e-07 atm
|
||||
Initial T = 1152.61, pres = 3.27763e-06 atm
|
||||
|
|
@ -551,12 +551,12 @@ Initial T = 1152.61, pres = 3.27763e-06 atm
|
|||
N 2.88709e-16 1.40138e-16 -19.0539
|
||||
O 2.24413e-06 1.24425e-06 -20.575
|
||||
NO 0.000134716 0.000140083 -39.6289
|
||||
NO+ 9.06539e-24 9.4264e-24 11.6607
|
||||
Electron 9.06539e-24 1.71214e-28 -69.7325
|
||||
N+ 1.62909e-57 7.90717e-58 32.2358
|
||||
O+ 3.09489e-44 1.7159e-44 30.7146
|
||||
N2+ 1.22735e-46 1.19147e-46 13.1819
|
||||
O2+ 3.24545e-32 3.5988e-32 10.1396
|
||||
NO+ 9.16628e-20 9.5313e-20 20.8821
|
||||
Electron 9.16628e-20 1.7312e-24 -60.5111
|
||||
N+ 1.64721e-53 7.99516e-54 41.4572
|
||||
O+ 3.12933e-40 1.73499e-40 39.936
|
||||
N2+ 1.24101e-42 1.20473e-42 22.4033
|
||||
O2+ 3.28157e-28 3.63885e-28 19.361
|
||||
|
||||
Final T = 1152.66, pres = 3.27775e-06 atm
|
||||
Initial T = 1152.66, pres = 3.27775e-05 atm
|
||||
|
|
@ -584,12 +584,12 @@ Initial T = 1152.66, pres = 3.27775e-05 atm
|
|||
N 9.1356e-17 4.43436e-17 -17.9026
|
||||
O 7.09896e-07 3.936e-07 -19.4238
|
||||
NO 0.000134733 0.000140101 -37.3264
|
||||
NO+ 9.06748e-25 9.42856e-25 11.6596
|
||||
Electron 9.06748e-25 1.71254e-29 -69.7322
|
||||
N+ 5.15899e-59 2.50404e-59 31.0834
|
||||
O+ 9.79676e-46 5.4316e-46 29.5622
|
||||
N2+ 1.22849e-47 1.19258e-47 13.1808
|
||||
O2+ 3.247e-33 3.60052e-33 10.1384
|
||||
NO+ 1.27093e-23 1.32154e-23 14.2998
|
||||
Electron 1.27093e-23 2.40034e-28 -67.092
|
||||
N+ 7.231e-58 3.50974e-58 33.7236
|
||||
O+ 1.37315e-44 7.61311e-45 32.2024
|
||||
N2+ 1.72189e-46 1.67155e-46 15.821
|
||||
O2+ 4.5511e-32 5.0466e-32 12.7787
|
||||
|
||||
Final T = 1152.67, pres = 3.27779e-05 atm
|
||||
Initial T = 1152.67, pres = 0.000327779 atm
|
||||
|
|
@ -598,7 +598,7 @@ Initial T = 1152.67, pres = 0.000327779 atm
|
|||
|
||||
temperature 1152.68 K
|
||||
pressure 33.2124 Pa
|
||||
density 1e-04 kg/m^3
|
||||
density 0.0001 kg/m^3
|
||||
mean mol. weight 28.8565 amu
|
||||
|
||||
1 kg 1 kmol
|
||||
|
|
@ -617,12 +617,12 @@ Initial T = 1152.67, pres = 0.000327779 atm
|
|||
N 2.88951e-17 1.40255e-17 -16.7513
|
||||
O 2.24513e-07 1.24481e-07 -18.2725
|
||||
NO 0.000134738 0.000140106 -35.0238
|
||||
NO+ 9.06814e-26 9.42925e-26 11.6592
|
||||
Electron 9.06814e-26 1.71266e-30 -69.7322
|
||||
N+ 1.63215e-60 7.92205e-61 29.9317
|
||||
O+ 3.099e-47 1.71817e-47 28.4106
|
||||
N2+ 1.22885e-48 1.19293e-48 13.1804
|
||||
O2+ 3.24749e-34 3.60106e-34 10.1381
|
||||
NO+ 5.59478e-26 5.81757e-26 11.1763
|
||||
Electron 5.59478e-26 1.05666e-30 -70.2151
|
||||
N+ 1.00699e-60 4.88768e-61 29.4488
|
||||
O+ 1.91199e-47 1.06006e-47 27.9276
|
||||
N2+ 7.58165e-49 7.36002e-49 12.6975
|
||||
O2+ 2.00361e-34 2.22175e-34 9.65516
|
||||
|
||||
Final T = 1152.68, pres = 0.00032778 atm
|
||||
Initial T = 1152.68, pres = 0.0032778 atm
|
||||
|
|
@ -650,12 +650,12 @@ Initial T = 1152.68, pres = 0.0032778 atm
|
|||
N 9.13803e-18 4.43553e-18 -15.6
|
||||
O 7.09996e-08 3.93655e-08 -17.1212
|
||||
NO 0.00013474 0.000140108 -32.7212
|
||||
NO+ 9.06835e-27 9.42946e-27 11.6591
|
||||
Electron 9.06835e-27 1.7127e-31 -69.7322
|
||||
N+ 5.16206e-62 2.50553e-62 28.7803
|
||||
O+ 9.80088e-49 5.43388e-49 27.2592
|
||||
N2+ 1.22896e-49 1.19304e-49 13.1803
|
||||
O2+ 3.24764e-35 3.60123e-35 10.138
|
||||
NO+ 2.9014e-21 3.01693e-21 24.335
|
||||
Electron 2.9014e-21 5.47974e-26 -57.0562
|
||||
N+ 1.65159e-56 8.01638e-57 41.4562
|
||||
O+ 3.13577e-43 1.73856e-43 39.9351
|
||||
N2+ 3.93204e-44 3.81709e-44 25.8562
|
||||
O2+ 1.03908e-29 1.1522e-29 22.8139
|
||||
|
||||
Final T = 1152.68, pres = 0.00327781 atm
|
||||
Initial T = 1152.68, pres = 0.0327781 atm
|
||||
|
|
@ -683,12 +683,12 @@ Initial T = 1152.68, pres = 0.0327781 atm
|
|||
N 2.88976e-18 1.40267e-18 -14.4487
|
||||
O 2.24523e-08 1.24486e-08 -15.9699
|
||||
NO 0.00013474 0.000140109 -30.4186
|
||||
NO+ 9.8979e-28 1.0292e-27 11.7466
|
||||
Electron 9.8979e-28 1.86937e-32 -69.6446
|
||||
N+ 1.78178e-63 8.6483e-64 27.7165
|
||||
O+ 3.38291e-50 1.87558e-50 26.1954
|
||||
N2+ 1.34142e-50 1.3022e-50 13.2678
|
||||
O2+ 3.54476e-36 3.93069e-36 10.2255
|
||||
NO+ 9.1752e-22 9.54057e-22 25.4863
|
||||
Electron 9.1752e-22 1.73288e-26 -55.9049
|
||||
N+ 1.65168e-57 8.01684e-58 41.4562
|
||||
O+ 3.13591e-44 1.73863e-44 39.935
|
||||
N2+ 1.24347e-44 1.20712e-44 27.0075
|
||||
O2+ 3.28593e-30 3.64369e-30 23.9652
|
||||
|
||||
Final T = 1152.68, pres = 0.0327781 atm
|
||||
Initial T = 3008.94, pres = 8.55633e-08 atm
|
||||
|
|
@ -1060,7 +1060,7 @@ Initial T = 3318.59, pres = 0.00114605 atm
|
|||
|
||||
temperature 3417.09 K
|
||||
pressure 119.127 Pa
|
||||
density 1e-04 kg/m^3
|
||||
density 0.0001 kg/m^3
|
||||
mean mol. weight 23.8496 amu
|
||||
|
||||
1 kg 1 kmol
|
||||
|
|
|
|||
|
|
@ -2,7 +2,7 @@
|
|||
#
|
||||
#
|
||||
temp_success="1"
|
||||
/bin/rm -f output.txt outputa.txt
|
||||
/bin/rm -f output.txt outputa.txt table.csv diff_out.txt diff_csv.txt
|
||||
testName="ChemEquil_ionizedGas"
|
||||
|
||||
#################################################################
|
||||
|
|
@ -19,14 +19,34 @@ then
|
|||
echo "ERROR: $testName: ionizedGasEquil returned with bad status, $retnStat, check output"
|
||||
fi
|
||||
|
||||
diff -w output.txt output_blessed.txt > diff_test.out
|
||||
retnStat=$?
|
||||
if [ $retnStat = "0" ]
|
||||
then
|
||||
echo "successful diff comparison on $testName test"
|
||||
temp_success="0"
|
||||
else
|
||||
echo "FAILED: unsuccessful diff comparison on $testName test"
|
||||
echo "FAILED" > csvCode.txt
|
||||
fi
|
||||
diff -w output.txt output_blessed.txt > diff_out.txt
|
||||
retnStat_txt=$?
|
||||
csvdiff -a 1.0E-50 table.csv table_blessed.csv > diff_csv.txt
|
||||
retnStat_csv=$?
|
||||
|
||||
eCode=1
|
||||
if test $retnStat_csv = "1"
|
||||
then
|
||||
eCode=0
|
||||
echo "Successful test comparison on " $testName
|
||||
if [ $retnStat_txt != "0" ]
|
||||
then
|
||||
echo " But, text output files have differences. See diff_out.txt"
|
||||
fi
|
||||
else
|
||||
echo "Unsuccessful test comparison on " $testName " test"
|
||||
if test $retnStat_csv != "1"
|
||||
then
|
||||
echo " csv files are different - see diff_csv.txt"
|
||||
fi
|
||||
if test $retnStat_txt != "0"
|
||||
then
|
||||
echo " And, text output files have differences. See diff_out.txt"
|
||||
fi
|
||||
fi
|
||||
#
|
||||
#
|
||||
exit 0
|
||||
|
||||
|
||||
|
||||
|
|
|
|||
43
test_problems/ChemEquil_ionizedGas/table_blessed.csv
Normal file
43
test_problems/ChemEquil_ionizedGas/table_blessed.csv
Normal file
|
|
@ -0,0 +1,43 @@
|
|||
Temperature, Pressure, N2, O2, N, O, NO, NO+, Electron, N+, O+, N2+, O2+
|
||||
207.3, 0.0005973, 0.7885, 0.2115, 9.558e-113, 1.181e-56, 1.826e-23, 0, 0, 0, 0, 0, 0
|
||||
207.3, 0.005973, 0.7885, 0.2115, 3.023e-113, 3.734e-57, 1.826e-23, 0, 0, 0, 0, 0, 0
|
||||
207.3, 0.05973, 0.7885, 0.2115, 9.558e-114, 1.181e-57, 1.826e-23, 0, 0, 0, 0, 0, 0
|
||||
207.3, 0.5973, 0.7885, 0.2115, 3.023e-114, 3.734e-58, 1.826e-23, 0, 0, 0, 0, 0, 0
|
||||
207.3, 5.973, 0.7885, 0.2115, 9.558e-115, 1.181e-58, 1.826e-23, 0, 0, 0, 0, 0, 0
|
||||
207.3, 59.73, 0.7885, 0.2115, 3.023e-115, 3.734e-59, 1.826e-23, 0, 0, 0, 0, 0, 0
|
||||
207.3, 597.3, 0.7885, 0.2115, 9.558e-116, 1.181e-59, 1.826e-23, 0, 0, 0, 0, 0, 0
|
||||
414.1, 0.001193, 0.7885, 0.2115, 2.069e-53, 1.874e-25, 5.586e-12, 4.006e-102, 4.006e-102, 0, 0, 0, 3.23e-125
|
||||
414.1, 0.01193, 0.7885, 0.2115, 6.541e-54, 5.925e-26, 5.586e-12, 5.137e-103, 5.137e-103, 0, 0, 0, 4.141e-126
|
||||
414.1, 0.1193, 0.7885, 0.2115, 2.069e-54, 1.874e-26, 5.586e-12, 2.826e-104, 2.826e-104, 0, 0, 0, 2.278e-127
|
||||
414.1, 1.193, 0.7885, 0.2115, 6.541e-55, 5.925e-27, 5.586e-12, 3.94e-105, 3.94e-105, 0, 0, 0, 3.177e-128
|
||||
414.1, 11.93, 0.7885, 0.2115, 2.069e-55, 1.874e-27, 5.586e-12, 2.917e-106, 2.917e-106, 0, 0, 0, 2.352e-129
|
||||
414.1, 119.3, 0.7885, 0.2115, 6.541e-56, 5.925e-28, 5.586e-12, 3.856e-107, 3.856e-107, 0, 0, 0, 3.109e-130
|
||||
414.1, 1193, 0.7885, 0.2115, 2.069e-56, 1.874e-28, 5.586e-12, 2.851e-102, 2.851e-102, 0, 0, 0, 2.298e-125
|
||||
1152, 0.003321, 0.7884, 0.2114, 2.863e-15, 2.234e-05, 0.0001345, 7.028e-19, 7.028e-19, 1.243e-51, 2.375e-38, 9.429e-42, 2.508e-27
|
||||
1153, 0.03321, 0.7884, 0.2115, 9.111e-16, 7.089e-06, 0.0001347, 9.674e-23, 9.674e-23, 5.477e-56, 1.042e-42, 1.307e-45, 3.461e-31
|
||||
1153, 0.3321, 0.7884, 0.2115, 2.887e-16, 2.244e-06, 0.0001347, 9.166e-20, 9.166e-20, 1.647e-53, 3.129e-40, 1.241e-42, 3.282e-28
|
||||
1153, 3.321, 0.7884, 0.2115, 9.136e-17, 7.099e-07, 0.0001347, 1.271e-23, 1.271e-23, 7.231e-58, 1.373e-44, 1.722e-46, 4.551e-32
|
||||
1153, 33.21, 0.7884, 0.2115, 2.89e-17, 2.245e-07, 0.0001347, 5.595e-26, 5.595e-26, 1.007e-60, 1.912e-47, 7.582e-49, 2.004e-34
|
||||
1153, 332.1, 0.7884, 0.2115, 9.138e-18, 7.1e-08, 0.0001347, 2.901e-21, 2.901e-21, 1.652e-56, 3.136e-43, 3.932e-44, 1.039e-29
|
||||
1153, 3321, 0.7884, 0.2115, 2.89e-18, 2.245e-08, 0.0001347, 9.175e-22, 9.175e-22, 1.652e-57, 3.136e-44, 1.243e-44, 3.286e-30
|
||||
1764, 0.005446, 0.7345, 0.1293, 7.047e-08, 0.1335, 0.00278, 5.61e-11, 5.61e-11, 1.279e-29, 2.263e-21, 1.383e-25, 1.154e-16
|
||||
1879, 0.05766, 0.7384, 0.1357, 1.608e-07, 0.1218, 0.004184, 1.472e-10, 1.472e-10, 4.308e-28, 2.092e-20, 3.125e-24, 6.605e-16
|
||||
2007, 0.6114, 0.7426, 0.1427, 3.489e-07, 0.1084, 0.006245, 3.68e-10, 3.68e-10, 1.256e-26, 1.732e-19, 6.362e-23, 3.535e-15
|
||||
2147, 6.489, 0.7472, 0.1506, 7.059e-07, 0.09302, 0.009208, 8.567e-10, 8.567e-10, 2.983e-25, 1.225e-18, 1.116e-21, 1.716e-14
|
||||
2299, 68.85, 0.752, 0.1591, 1.288e-06, 0.07561, 0.01332, 1.793e-09, 1.793e-09, 5.239e-24, 6.863e-18, 1.573e-20, 7.206e-14
|
||||
2457, 728.7, 0.7568, 0.1679, 2.009e-06, 0.05658, 0.01867, 3.186e-09, 3.186e-09, 5.818e-23, 2.697e-17, 1.587e-19, 2.417e-13
|
||||
2611, 7665, 0.7613, 0.1765, 2.472e-06, 0.0373, 0.02495, 4.397e-09, 4.398e-09, 3.215e-22, 6.159e-17, 9.61e-19, 5.738e-13
|
||||
2946, 0.01041, 0.6308, 1.393e-06, 0.02415, 0.3449, 0.0001029, 2.912e-06, 2.913e-06, 7.396e-12, 8.35e-10, 3.216e-12, 3.032e-12
|
||||
3080, 0.1085, 0.6361, 5.853e-06, 0.01774, 0.3459, 0.0002489, 3.28e-06, 3.281e-06, 6.298e-12, 8.265e-10, 4.524e-12, 9.624e-12
|
||||
3207, 1.126, 0.6411, 2.749e-05, 0.0117, 0.3466, 0.0006234, 3.384e-06, 3.384e-06, 3.787e-12, 6.566e-10, 4.893e-12, 2.838e-11
|
||||
3319, 11.61, 0.6451, 0.0001486, 0.006702, 0.3464, 0.001629, 3.11e-06, 3.11e-06, 1.453e-12, 3.921e-10, 3.762e-12, 7.586e-11
|
||||
3417, 119.1, 0.6477, 0.0008763, 0.003474, 0.3435, 0.004358, 2.623e-06, 2.623e-06, 4.057e-13, 1.906e-10, 2.273e-12, 1.873e-10
|
||||
3555, 1230, 0.6499, 0.00421, 0.002098, 0.3329, 0.01082, 2.485e-06, 2.485e-06, 1.882e-13, 1.254e-10, 2.028e-12, 4.968e-10
|
||||
3800, 1.297e+04, 0.6534, 0.01249, 0.001863, 0.3095, 0.02275, 3.185e-06, 3.186e-06, 3.119e-13, 1.785e-10, 4.815e-12, 1.621e-09
|
||||
3294, 0.0125, 0.5182, 1.58e-07, 0.1606, 0.3211, 4.646e-05, 1.403e-05, 1.408e-05, 4.794e-09, 5.129e-08, 4.044e-10, 1.185e-11
|
||||
3516, 0.1328, 0.5258, 5.179e-07, 0.1514, 0.3227, 0.0001044, 1.954e-05, 1.963e-05, 9.133e-09, 8.456e-08, 1.06e-09, 4.501e-11
|
||||
3767, 1.414, 0.5345, 1.729e-06, 0.1408, 0.3244, 0.0002361, 2.703e-05, 2.719e-05, 1.669e-08, 1.359e-07, 2.705e-09, 1.701e-10
|
||||
4050, 15.11, 0.5445, 5.91e-06, 0.1286, 0.3262, 0.0005385, 3.701e-05, 3.726e-05, 2.88e-08, 2.11e-07, 6.656e-09, 6.382e-10
|
||||
4372, 161.8, 0.556, 2.081e-05, 0.1146, 0.3281, 0.001241, 4.987e-05, 5.025e-05, 4.594e-08, 3.117e-07, 1.555e-08, 2.363e-09
|
||||
4734, 1736, 0.5688, 7.603e-05, 0.09874, 0.3294, 0.002891, 6.556e-05, 6.609e-05, 6.571e-08, 4.295e-07, 3.37e-08, 8.558e-09
|
||||
5140, 1.863e+04, 0.5827, 0.0002894, 0.08126, 0.3288, 0.006811, 8.306e-05, 8.378e-05, 8.126e-08, 5.379e-07, 6.601e-08, 2.99e-08
|
||||
|
Loading…
Add table
Reference in a new issue