Apply clang-format to source tree
[alexxy/gromacs.git] / src / programs / mdrun / tests / initialconstraints.cpp
index ab9f27f76d3abc3bfb975e7614fb5d93e0fc9d31..d4ae76ae5b8a2b81e3a7dae94334c8e19a9e24b5 100644 (file)
@@ -63,11 +63,12 @@ namespace
 {
 
 //! This type holds input integrators. Now it's holding names, but ei* enum values from md_enums.h could be used instead.
-using EnergyIntegratorType = const char *;
+using EnergyIntegratorType = const char*;
 
 //! Test fixture parametrized on integrators
-class InitialConstraintsTest : public gmx::test::MdrunTestFixture,
-                               public ::testing::WithParamInterface<EnergyIntegratorType>
+class InitialConstraintsTest :
+    public gmx::test::MdrunTestFixture,
+    public ::testing::WithParamInterface<EnergyIntegratorType>
 {
 };
 
@@ -78,19 +79,20 @@ TEST_P(InitialConstraintsTest, Works)
     auto              integrator = GetParam();
     const std::string integratorName(integrator);
     SCOPED_TRACE("Integrating with " + integratorName);
-    const std::string theMdpFile = formatString("nstcalcenergy           = 1\n"
-                                                "nstenergy               = 1\n"
-                                                "comm-mode               = linear\n"
-                                                "continuation            = no\n"
-                                                "constraints             = h-bonds\n"
-                                                "lincs_iter              = 2\n"
-                                                "verlet-buffer-tolerance = 1e-4\n"
-                                                "nsttcouple              = 1\n"    // for md-vv-avek
-                                                "nstpcouple              = 1\n"    // for md-vv-avek
-                                                "integrator              = %s\n"
-                                                "nsteps                  = %d\n"
-                                                "dt                      = %f\n",
-                                                integratorName.c_str(), nsteps, timestep);
+    const std::string theMdpFile = formatString(
+            "nstcalcenergy           = 1\n"
+            "nstenergy               = 1\n"
+            "comm-mode               = linear\n"
+            "continuation            = no\n"
+            "constraints             = h-bonds\n"
+            "lincs_iter              = 2\n"
+            "verlet-buffer-tolerance = 1e-4\n"
+            "nsttcouple              = 1\n" // for md-vv-avek
+            "nstpcouple              = 1\n" // for md-vv-avek
+            "integrator              = %s\n"
+            "nsteps                  = %d\n"
+            "dt                      = %f\n",
+            integratorName.c_str(), nsteps, timestep);
 
     runner_.useStringAsMdpFile(theMdpFile);
 
@@ -101,14 +103,15 @@ TEST_P(InitialConstraintsTest, Works)
     runner_.edrFileName_ = fileManager_.getTemporaryFilePath(inputFile + ".edr");
     ASSERT_EQ(0, runner_.callMdrun());
 
-    auto energyReader = openEnergyFileToReadTerms(runner_.edrFileName_, {"Total Energy", "Kinetic En."});
-    real totalEnergy  = 0.0, prevTotalEnergy = 0.0;
-    auto tolerance    = ulpTolerance(0); // The real value is set below from starting kinetic energy
+    auto energyReader =
+            openEnergyFileToReadTerms(runner_.edrFileName_, { "Total Energy", "Kinetic En." });
+    real totalEnergy = 0.0, prevTotalEnergy = 0.0;
+    auto tolerance = ulpTolerance(0); // The real value is set below from starting kinetic energy
     for (int i = 0; i <= nsteps; i++)
     {
         EnergyFrame frame = energyReader->frame();
-        prevTotalEnergy = totalEnergy;
-        totalEnergy     = frame.at("Total Energy");
+        prevTotalEnergy   = totalEnergy;
+        totalEnergy       = frame.at("Total Energy");
         if (i == 0)
         {
             // We set the tolerance for total energy based on magnitude of kinetic energy.
@@ -124,10 +127,10 @@ TEST_P(InitialConstraintsTest, Works)
 }
 
 //! Integrators with energy conservation to test
-const EnergyIntegratorType c_integratorsToTest [] = {"md", "md-vv", "md-vv-avek"};
+const EnergyIntegratorType c_integratorsToTest[] = { "md", "md-vv", "md-vv-avek" };
 
 INSTANTIATE_TEST_CASE_P(Checking, InitialConstraintsTest, ::testing::ValuesIn(c_integratorsToTest));
 
-}  // namespace
-}  // namespace test
-}  // namespace gmx
+} // namespace
+} // namespace test
+} // namespace gmx