research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

[Test] Make reaction rate comparisons more granular

Browse source
This commit is contained in:
Ray Speth 2016-10-15 11:44:11 -04:00
parent e7ac87918d
commit 50d0f9913f
1 changed files with 19 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

26
test/kinetics/kineticsFromScratch.cpp
View file

@ -420,24 +420,36 @@ public:
vector_fp k(kin.nReactions()), k_ref(kin_ref.nReactions());
vector_fp w(kin.nTotalSpecies()), w_ref(kin_ref.nTotalSpecies());
kin.getCreationRates(w.data());
kin_ref.getCreationRates(w_ref.data());
for (size_t i = 0; i < kin.nTotalSpecies(); i++) {
size_t iref = p_ref.speciesIndex(p.speciesName(i));
EXPECT_DOUBLE_EQ(w_ref[iref], w[i]) << "sp = " << p.speciesName(i) << "; N = " << N;
}
kin.getFwdRateConstants(k.data());
kin_ref.getFwdRateConstants(k_ref.data());
for (size_t i = 0; i < kin.nReactions(); i++) {
EXPECT_DOUBLE_EQ(k_ref[i], k[i]) << "i = " << i << "; N = " << N;
}
kin.getFwdRatesOfProgress(k.data());
kin_ref.getFwdRatesOfProgress(k_ref.data());
for (size_t i = 0; i < kin.nReactions(); i++) {
EXPECT_DOUBLE_EQ(k_ref[i], k[i]) << "i = " << i << "; N = " << N;
}
kin.getRevRateConstants(k.data());
kin_ref.getRevRateConstants(k_ref.data());
for (size_t i = 0; i < kin.nReactions(); i++) {
EXPECT_DOUBLE_EQ(k_ref[i], k[i]) << "i = " << i << "; N = " << N;
}
kin.getRevRatesOfProgress(k.data());
kin_ref.getRevRatesOfProgress(k_ref.data());
for (size_t i = 0; i < kin.nReactions(); i++) {
EXPECT_DOUBLE_EQ(k_ref[i], k[i]) << "i = " << i << "; N = " << N;
}
kin.getCreationRates(w.data());
kin_ref.getCreationRates(w_ref.data());
for (size_t i = 0; i < kin.nTotalSpecies(); i++) {
size_t iref = p_ref.speciesIndex(p.speciesName(i));
EXPECT_DOUBLE_EQ(w_ref[iref], w[i]) << "sp = " << p.speciesName(i) << "; N = " << N;
}
}
};