src/programs/mdrun/tests/energycomparison.h

   1 /*
   2  * This file is part of the GROMACS molecular simulation package.
   3  *
   4  * Copyright (c) 2018,2019, by the GROMACS development team, led by
   5  * Mark Abraham, David van der Spoel, Berk Hess, and Erik Lindahl,
   6  * and including many others, as listed in the AUTHORS file in the
   7  * top-level source directory and at http://www.gromacs.org.
   8  *
   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.
  13  *
  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.
  18  *
  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.
  23  *
  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.
  31  *
  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.
  34  */
  35
  36 /*! \internal \file
  37  * \brief Declares function for comparing energy frames.
  38  *
  39  * \author Mark Abraham <mark.j.abraham@gmail.com>
  40  * \ingroup module_mdrun_integration_tests
  41  */
  42 #ifndef GMX_PROGRAMS_MDRUN_TESTS_ENERGYCOMPARISON_H
  43 #define GMX_PROGRAMS_MDRUN_TESTS_ENERGYCOMPARISON_H
  44
  45 #include <string>
  46 #include <unordered_map>
  47
  48 #include "testutils/testasserts.h"
  49
  50 namespace gmx
  51 {
  52
  53 class EnergyFrame;
  54
  55 namespace test
  56 {
  57
  58 class TestReferenceChecker;
  59
  60 //! Convenience type
  61 using EnergyTermsToCompare = std::unordered_map<std::string, FloatingPointTolerance>;
  62
  63 /*! \internal
  64  * \brief Function object to compare all energy terms of reference
  65  * with all matching terms from test within the given tolerances. */
  66 class EnergyComparison
  67 {
  68     public:
  69         //! Defaults for energy comparisons
  70         static const EnergyTermsToCompare s_defaultEnergyTermsToCompare;
  71         //! Constructor
  72         EnergyComparison(const EnergyTermsToCompare &energyTermsToCompare);
  73         /*! \brief Return the names of energies that will be compared
  74          *
  75          * This function can be used to provide an input for
  76          * openEnergyFileToReadTerms().
  77          *
  78          * \todo This returns a copy of the keys, which is convenient, but
  79          * inefficient. Alternatively, this could return a view of the keys
  80          * from a range rather than a container, but there's no implementation
  81          * of that in C++11 at the moment. */
  82         std::vector<std::string> getEnergyNames() const;
  83         /*! \brief Compare \c reference with \c test within \c
  84          * energyTermsToCompare_
  85          *
  86          * Ignore any key found in either \c reference or \c test that is not
  87          * found in the other. For all keys found in both frames, compare the
  88          * values with EXPECT_REAL_EQ_TOL and the given tolerance for that
  89          * key. */
  90         void operator()(const EnergyFrame &reference, const EnergyFrame &test) const;
  91
  92         //! Energy terms to match with given tolerances.
  93         EnergyTermsToCompare energyTermsToCompare_;
  94 };
  95
  96 /*! \brief Check a subset of the energies found in an energy file
  97  * against reference data.
  98  *
  99  * Opens the energy file, loops over all frames, matching the
 100  * indicated energies against refdata at the given tolerance.
 101  *
 102  * \param[in]  energyFilename        The name of an energy file.
 103  * \param[in]  energyTermsToCompare  Set of energies to match at given tolerances.
 104  * \param[in]  checker               Root checker for reference data.
 105  *
 106  * \todo This is quite similar to the functionality used in PmeTest,
 107  * and we should consider reducing the duplication.
 108  */
 109 void
 110 checkEnergiesAgainstReferenceData(const std::string          &energyFilename,
 111                                   const EnergyTermsToCompare &energyTermsToCompare,
 112                                   TestReferenceChecker       *checker);
 113
 114 }  // namespace test
 115 }  // namespace gmx
 116
 117 #endif