12345qiupeng
/
gtsam
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

update test

release/4.3a0
Varun Agrawal 2022-12-24 09:08:47 +05:30
parent ff8a58671d
commit 0be2a679c0
1 changed files with 16 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
gtsam/hybrid/tests/testHybridBayesNet.cpp
View File

@ -356,7 +356,7 @@ TEST(HybridBayesNet, Sampling) {
size_t num_samples = 1000; size_t num_samples = 1000;
for (size_t i = 0; i < num_samples; i++) { for (size_t i = 0; i < num_samples; i++) {
// Sample // Sample
HybridValues sample = bn->sample(&gen, noise_model); HybridValues sample = bn->sample(&gen);
discrete_samples.push_back(sample.discrete()[M(0)]); discrete_samples.push_back(sample.discrete()[M(0)]);
@ -366,18 +366,23 @@ TEST(HybridBayesNet, Sampling) {
average_continuous += sample.continuous(); average_continuous += sample.continuous();
} }
} }
double discrete_sum =
std::accumulate(discrete_samples.begin(), discrete_samples.end(),
decltype(discrete_samples)::value_type(0));
// regression for specific RNG seed EXPECT_LONGS_EQUAL(2, average_continuous.size());
EXPECT_DOUBLES_EQUAL(0.477, discrete_sum / num_samples, 1e-9); EXPECT_LONGS_EQUAL(num_samples, discrete_samples.size());
VectorValues expected; // Regressions don't work across platforms :-(
expected.insert({X(0), Vector1(-0.0131207162712)}); // // regression for specific RNG seed
expected.insert({X(1), Vector1(-0.499026377568)}); // double discrete_sum =
// regression for specific RNG seed // std::accumulate(discrete_samples.begin(), discrete_samples.end(),
EXPECT(assert_equal(expected, average_continuous.scale(1.0 / num_samples))); // decltype(discrete_samples)::value_type(0));
// EXPECT_DOUBLES_EQUAL(0.477, discrete_sum / num_samples, 1e-9);
// VectorValues expected;
// expected.insert({X(0), Vector1(-0.0131207162712)});
// expected.insert({X(1), Vector1(-0.499026377568)});
// // regression for specific RNG seed
// EXPECT(assert_equal(expected, average_continuous.scale(1.0 /
// num_samples)));
} }
/* ************************************************************************* */ /* ************************************************************************* */