2 * This file is part of the GROMACS molecular simulation package.
4 * Copyright (c) 2012,2013, by the GROMACS development team, led by
5 * David van der Spoel, Berk Hess, Erik Lindahl, and including many
6 * others, as listed in the AUTHORS file in the top-level source
7 * directory and at http://www.gromacs.org.
9 * GROMACS is free software; you can redistribute it and/or
10 * modify it under the terms of the GNU Lesser General Public License
11 * as published by the Free Software Foundation; either version 2.1
12 * of the License, or (at your option) any later version.
14 * GROMACS is distributed in the hope that it will be useful,
15 * but WITHOUT ANY WARRANTY; without even the implied warranty of
16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
17 * Lesser General Public License for more details.
19 * You should have received a copy of the GNU Lesser General Public
20 * License along with GROMACS; if not, see
21 * http://www.gnu.org/licenses, or write to the Free Software Foundation,
22 * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
24 * If you want to redistribute modifications to GROMACS, please
25 * consider that scientific software is very special. Version
26 * control is crucial - bugs must be traceable. We will be happy to
27 * consider code for inclusion in the official distribution, but
28 * derived work must not be called official GROMACS. Details are found
29 * in the README & COPYING files - if they are missing, get the
30 * official version at http://www.gromacs.org.
32 * To help us fund GROMACS development, we humbly ask that you cite
33 * the research papers on the package. Check out http://www.gromacs.org.
37 * Tests for functionality of analysis data histogram modules.
39 * These tests check that classes in histogram.h compute histograms correctly
40 * with simple input data. Also different ways of initializing the histograms
42 * Checking is done using gmx::test::AnalysisDataTestFixture and reference
43 * data. Also the input data is written to the reference data to catch
44 * out-of-date reference.
46 * \author Teemu Murtola <teemu.murtola@gmail.com>
47 * \ingroup module_analysisdata
49 #include <gtest/gtest.h>
51 #include "gromacs/analysisdata/analysisdata.h"
52 #include "gromacs/analysisdata/modules/histogram.h"
54 #include "testutils/datatest.h"
55 #include "testutils/testasserts.h"
57 using gmx::test::AnalysisDataTestInput;
62 /********************************************************************
63 * Tests for gmx::AnalysisHistogramSettings.
65 * These tests check that gmx::AnalysisHistogramSettings objects can be
66 * initialized from various types of values, and that the bin positions are
67 * computed correctly based on the input values.
70 TEST(AnalysisHistogramSettingsTest, InitializesFromBins)
72 gmx::AnalysisHistogramSettings settings(
73 gmx::histogramFromBins(1.0, 5, 0.5));
74 EXPECT_FLOAT_EQ(1.0, settings.firstEdge());
75 EXPECT_EQ(5, settings.binCount());
76 EXPECT_FLOAT_EQ(0.5, settings.binWidth());
77 EXPECT_FLOAT_EQ(3.5, settings.lastEdge());
81 TEST(AnalysisHistogramSettingsTest, InitializesFromBinsWithIntegerBins)
83 gmx::AnalysisHistogramSettings settings(
84 gmx::histogramFromBins(1.0, 5, 0.5).integerBins());
85 EXPECT_FLOAT_EQ(0.75, settings.firstEdge());
86 EXPECT_EQ(5, settings.binCount());
87 EXPECT_FLOAT_EQ(0.5, settings.binWidth());
88 EXPECT_FLOAT_EQ(3.25, settings.lastEdge());
92 TEST(AnalysisHistogramSettingsTest, InitializesFromRangeWithBinCount)
94 gmx::AnalysisHistogramSettings settings(
95 gmx::histogramFromRange(1.0, 4.0).binCount(6));
96 EXPECT_FLOAT_EQ(1.0, settings.firstEdge());
97 EXPECT_FLOAT_EQ(4.0, settings.lastEdge());
98 EXPECT_EQ(6, settings.binCount());
99 EXPECT_FLOAT_EQ(0.5, settings.binWidth());
103 TEST(AnalysisHistogramSettingsTest, InitializesFromRangeWithBinWidth)
105 gmx::AnalysisHistogramSettings settings(
106 gmx::histogramFromRange(1.0, 4.0).binWidth(0.5));
107 EXPECT_FLOAT_EQ(1.0, settings.firstEdge());
108 EXPECT_FLOAT_EQ(4.0, settings.lastEdge());
109 EXPECT_FLOAT_EQ(0.5, settings.binWidth());
110 EXPECT_EQ(6, settings.binCount());
114 TEST(AnalysisHistogramSettingsTest, InitializesFromRangeWithBinCountAndIntegerBins)
116 gmx::AnalysisHistogramSettings settings(
117 gmx::histogramFromRange(1.0, 4.0).binCount(7).integerBins());
118 EXPECT_FLOAT_EQ(0.75, settings.firstEdge());
119 EXPECT_FLOAT_EQ(4.25, settings.lastEdge());
120 EXPECT_EQ(7, settings.binCount());
121 EXPECT_FLOAT_EQ(0.5, settings.binWidth());
125 TEST(AnalysisHistogramSettingsTest, InitializesFromRangeWithBinWidthAndIntegerBins)
127 gmx::AnalysisHistogramSettings settings(
128 gmx::histogramFromRange(1.0, 4.0).binWidth(0.5).integerBins());
129 EXPECT_FLOAT_EQ(0.75, settings.firstEdge());
130 EXPECT_FLOAT_EQ(4.25, settings.lastEdge());
131 EXPECT_FLOAT_EQ(0.5, settings.binWidth());
132 EXPECT_EQ(7, settings.binCount());
136 TEST(AnalysisHistogramSettingsTest, InitializesFromRangeWithRoundedRange)
138 gmx::AnalysisHistogramSettings settings(
139 gmx::histogramFromRange(1.2, 3.8).binWidth(0.5).roundRange());
140 EXPECT_FLOAT_EQ(1.0, settings.firstEdge());
141 EXPECT_FLOAT_EQ(4.0, settings.lastEdge());
142 EXPECT_FLOAT_EQ(0.5, settings.binWidth());
143 EXPECT_EQ(6, settings.binCount());
147 /********************************************************************
148 * Tests for gmx::AnalysisDataSimpleHistogramModule.
151 //! Test fixture for gmx::AnalysisDataSimpleHistogramModule.
152 typedef gmx::test::AnalysisDataTestFixture SimpleHistogramModuleTest;
154 // Input data for gmx::AnalysisDataSimpleHistogramModule tests.
155 class SimpleInputData
158 static const AnalysisDataTestInput &get()
160 static SimpleInputData singleton;
161 return singleton.data_;
164 SimpleInputData() : data_(1, true)
166 using gmx::test::AnalysisDataTestInputFrame;
167 data_.setColumnCount(0, 1);
168 AnalysisDataTestInputFrame &frame1 = data_.addFrame(1.0);
169 frame1.addPointSetWithValues(0, 0, 0.7);
170 frame1.addPointSetWithValues(0, 0, 1.1);
171 frame1.addPointSetWithValues(0, 0, 2.3);
172 frame1.addPointSetWithValues(0, 0, 2.9);
173 AnalysisDataTestInputFrame &frame2 = data_.addFrame(2.0);
174 frame2.addPointSetWithValues(0, 0, 1.3);
175 frame2.addPointSetWithValues(0, 0, 2.2);
176 AnalysisDataTestInputFrame &frame3 = data_.addFrame(3.0);
177 frame3.addPointSetWithValues(0, 0, 3.3);
178 frame3.addPointSetWithValues(0, 0, 1.2);
179 frame3.addPointSetWithValues(0, 0, 1.3);
183 AnalysisDataTestInput data_;
186 TEST_F(SimpleHistogramModuleTest, ComputesCorrectly)
188 const AnalysisDataTestInput &input = SimpleInputData::get();
189 gmx::AnalysisData data;
190 ASSERT_NO_THROW_GMX(setupDataObject(input, &data));
192 gmx::AnalysisDataSimpleHistogramModulePointer module(
193 new gmx::AnalysisDataSimpleHistogramModule(
194 gmx::histogramFromRange(1.0, 3.0).binCount(4)));
195 data.addModule(module);
197 ASSERT_NO_THROW_GMX(addStaticCheckerModule(input, &data));
198 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("InputData", &data));
199 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("Histogram", module.get()));
200 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("HistogramAverage",
201 &module->averager()));
202 ASSERT_NO_THROW_GMX(presentAllData(input, &data));
203 ASSERT_NO_THROW_GMX(module->averager().done());
207 TEST_F(SimpleHistogramModuleTest, ComputesCorrectlyWithAll)
209 const AnalysisDataTestInput &input = SimpleInputData::get();
210 gmx::AnalysisData data;
211 ASSERT_NO_THROW_GMX(setupDataObject(input, &data));
213 gmx::AnalysisDataSimpleHistogramModulePointer module(
214 new gmx::AnalysisDataSimpleHistogramModule(
215 gmx::histogramFromRange(1.0, 3.0).binCount(4).includeAll()));
216 data.addModule(module);
218 ASSERT_NO_THROW_GMX(addStaticCheckerModule(input, &data));
219 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("InputData", &data));
220 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("Histogram", module.get()));
221 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("HistogramAverage",
222 &module->averager()));
223 ASSERT_NO_THROW_GMX(presentAllData(input, &data));
224 ASSERT_NO_THROW_GMX(module->averager().done());
228 /********************************************************************
229 * Tests for gmx::AnalysisDataWeightedHistogramModule.
232 //! Test fixture for gmx::AnalysisDataWeightedHistogramModule.
233 typedef gmx::test::AnalysisDataTestFixture WeightedHistogramModuleTest;
235 // Input data for both weighted histogram and bin average module tests.
236 class WeightedInputData
239 static const AnalysisDataTestInput &get()
241 static WeightedInputData singleton;
242 return singleton.data_;
245 WeightedInputData() : data_(1, true)
247 using gmx::test::AnalysisDataTestInputFrame;
248 data_.setColumnCount(0, 2);
249 AnalysisDataTestInputFrame &frame1 = data_.addFrame(1.0);
250 frame1.addPointSetWithValues(0, 0, 0.7, 0.5);
251 frame1.addPointSetWithValues(0, 0, 1.1, 1.0);
252 frame1.addPointSetWithValues(0, 0, 2.3, 1.0);
253 frame1.addPointSetWithValues(0, 0, 2.9, 2.0);
254 AnalysisDataTestInputFrame &frame2 = data_.addFrame(2.0);
255 frame2.addPointSetWithValues(0, 0, 1.3, 1.0);
256 frame2.addPointSetWithValues(0, 0, 2.2, 3.0);
257 AnalysisDataTestInputFrame &frame3 = data_.addFrame(3.0);
258 frame3.addPointSetWithValues(0, 0, 3.3, 0.5);
259 frame3.addPointSetWithValues(0, 0, 1.2, 2.0);
260 frame3.addPointSetWithValues(0, 0, 1.3, 1.0);
264 AnalysisDataTestInput data_;
267 TEST_F(WeightedHistogramModuleTest, ComputesCorrectly)
269 const AnalysisDataTestInput &input = WeightedInputData::get();
270 gmx::AnalysisData data;
271 ASSERT_NO_THROW_GMX(setupDataObject(input, &data));
273 gmx::AnalysisDataWeightedHistogramModulePointer module(
274 new gmx::AnalysisDataWeightedHistogramModule(
275 gmx::histogramFromRange(1.0, 3.0).binCount(4)));
276 data.addModule(module);
278 ASSERT_NO_THROW_GMX(addStaticCheckerModule(input, &data));
279 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("InputData", &data));
280 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("Histogram", module.get()));
281 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("HistogramAverage",
282 &module->averager()));
283 ASSERT_NO_THROW_GMX(presentAllData(input, &data));
284 ASSERT_NO_THROW_GMX(module->averager().done());
288 TEST_F(WeightedHistogramModuleTest, ComputesCorrectlyWithAll)
290 const AnalysisDataTestInput &input = WeightedInputData::get();
291 gmx::AnalysisData data;
292 ASSERT_NO_THROW_GMX(setupDataObject(input, &data));
294 gmx::AnalysisDataWeightedHistogramModulePointer module(
295 new gmx::AnalysisDataWeightedHistogramModule(
296 gmx::histogramFromRange(1.0, 3.0).binCount(4).includeAll()));
297 data.addModule(module);
299 ASSERT_NO_THROW_GMX(addStaticCheckerModule(input, &data));
300 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("InputData", &data));
301 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("Histogram", module.get()));
302 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("HistogramAverage",
303 &module->averager()));
304 ASSERT_NO_THROW_GMX(presentAllData(input, &data));
305 ASSERT_NO_THROW_GMX(module->averager().done());
309 /********************************************************************
310 * Tests for gmx::AnalysisDataBinAverageModule.
313 //! Test fixture for gmx::AnalysisDataBinAverageModule.
314 typedef gmx::test::AnalysisDataTestFixture BinAverageModuleTest;
316 TEST_F(BinAverageModuleTest, ComputesCorrectly)
318 const AnalysisDataTestInput &input = WeightedInputData::get();
319 gmx::AnalysisData data;
320 ASSERT_NO_THROW_GMX(setupDataObject(input, &data));
322 gmx::AnalysisDataBinAverageModulePointer module(
323 new gmx::AnalysisDataBinAverageModule(
324 gmx::histogramFromRange(1.0, 3.0).binCount(4)));
325 data.addModule(module);
327 ASSERT_NO_THROW_GMX(addStaticCheckerModule(input, &data));
328 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("InputData", &data));
329 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("HistogramAverage", module.get()));
330 ASSERT_NO_THROW_GMX(presentAllData(input, &data));
334 TEST_F(BinAverageModuleTest, ComputesCorrectlyWithAll)
336 const AnalysisDataTestInput &input = WeightedInputData::get();
337 gmx::AnalysisData data;
338 ASSERT_NO_THROW_GMX(setupDataObject(input, &data));
340 gmx::AnalysisDataBinAverageModulePointer module(
341 new gmx::AnalysisDataBinAverageModule(
342 gmx::histogramFromRange(1.0, 3.0).binCount(4).includeAll()));
343 data.addModule(module);
345 ASSERT_NO_THROW_GMX(addStaticCheckerModule(input, &data));
346 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("InputData", &data));
347 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("HistogramAverage", module.get()));
348 ASSERT_NO_THROW_GMX(presentAllData(input, &data));
352 /********************************************************************
353 * Tests for gmx::AbstractAverageHistogram.
355 * This class derives from gmx::AbstractAnalysisArrayData, and is tested using
356 * corresponding facilities in gmx::test::AnalysisDataTestFixture.
359 //! Test fixture for gmx::AbstractAverageHistogram.
360 typedef gmx::test::AnalysisDataTestFixture AbstractAverageHistogramTest;
362 // Input data for gmx::AbstractAverageHistogram tests.
363 class AverageInputData
366 static const AnalysisDataTestInput &get()
368 static AverageInputData singleton;
369 return singleton.data_;
372 AverageInputData() : data_(1, false)
374 data_.setColumnCount(0, 2);
375 data_.addFrameWithValues(1.0, 2.0, 1.0);
376 data_.addFrameWithValues(1.5, 1.0, 1.0);
377 data_.addFrameWithValues(2.0, 3.0, 2.0);
378 data_.addFrameWithValues(2.5, 4.0, 2.0);
379 data_.addFrameWithValues(3.0, 2.0, 1.0);
380 data_.addFrameWithValues(3.5, 0.0, 3.0);
381 data_.addFrameWithValues(4.0, 1.0, 3.0);
385 AnalysisDataTestInput data_;
389 * Mock object for testing gmx::AbstractAverageHistogram.
391 * Exposes necessary methods from gmx::AbstractAverageHistogram to use with
392 * gmx::test::AnalysisDataTestFixture::setupArrayData().
394 class MockAverageHistogram : public gmx::AbstractAverageHistogram
397 MockAverageHistogram() {}
398 //! Creates a histogram module with defined bin parameters.
399 explicit MockAverageHistogram(const gmx::AnalysisHistogramSettings &settings)
400 : AbstractAverageHistogram(settings)
404 using AbstractAverageHistogram::init;
405 using AbstractAverageHistogram::setColumnCount;
406 using AbstractAverageHistogram::setRowCount;
407 using AbstractAverageHistogram::allocateValues;
408 using AbstractAverageHistogram::setValue;
412 TEST_F(AbstractAverageHistogramTest, ClonesCorrectly)
414 const AnalysisDataTestInput &input = AverageInputData::get();
415 MockAverageHistogram data(
416 gmx::histogramFromBins(1.0, input.frameCount(), 0.5).integerBins());
417 setupArrayData(input, &data);
419 ASSERT_NO_THROW_GMX(addStaticCheckerModule(input, &data));
420 gmx::AverageHistogramPointer copy(data.clone());
421 ASSERT_NO_THROW_GMX(addStaticCheckerModule(input, copy.get()));
422 ASSERT_NO_THROW_GMX(copy->done());
423 ASSERT_NO_THROW_GMX(data.done());
424 gmx::AverageHistogramPointer copy2(data.clone());
425 ASSERT_NO_THROW_GMX(addStaticCheckerModule(input, copy2.get()));
426 ASSERT_NO_THROW_GMX(copy2->done());
430 TEST_F(AbstractAverageHistogramTest, ResamplesAtDoubleBinWidth)
432 const AnalysisDataTestInput &input = AverageInputData::get();
433 MockAverageHistogram data(
434 gmx::histogramFromBins(1.0, input.frameCount(), 0.5).integerBins());
435 setupArrayData(input, &data);
437 ASSERT_NO_THROW_GMX(addStaticCheckerModule(input, &data));
438 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("InputData", &data));
439 gmx::AverageHistogramPointer resampled(data.resampleDoubleBinWidth(false));
440 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("ResampledHistogram", resampled.get()));
441 ASSERT_NO_THROW_GMX(data.done());
442 ASSERT_NO_THROW_GMX(resampled->done());
446 TEST_F(AbstractAverageHistogramTest, ResamplesAtDoubleBinWidthWithIntegerBins)
448 const AnalysisDataTestInput &input = AverageInputData::get();
449 MockAverageHistogram data(
450 gmx::histogramFromBins(1.0, input.frameCount(), 0.5).integerBins());
451 setupArrayData(input, &data);
453 ASSERT_NO_THROW_GMX(addStaticCheckerModule(input, &data));
454 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("InputData", &data));
455 gmx::AverageHistogramPointer resampled(data.resampleDoubleBinWidth(true));
456 ASSERT_NO_THROW_GMX(addReferenceCheckerModule("ResampledHistogram", resampled.get()));
457 ASSERT_NO_THROW_GMX(data.done());
458 ASSERT_NO_THROW_GMX(resampled->done());