7a8f2e737a3c29228b2ad12ddf9a7bfaee26e3d3
[alexxy/gromacs.git] / src / gromacs / gmxpreprocess / tests / solvate.cpp
1 /*
2  * This file is part of the GROMACS molecular simulation package.
3  *
4  * Copyright (c) 2013,2014,2015,2017,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 /*! \internal \file
36  * \brief
37  * Tests for solvation.
38  *
39  * \author Mark Abraham <mark.j.abraham@gmail.com>
40  * \author Teemu Murtola <teemu.murtola@gmail.com>
41  */
42
43 #include "gmxpre.h"
44
45 #include "gromacs/gmxpreprocess/solvate.h"
46
47 #include "gromacs/utility/futil.h"
48 #include "gromacs/utility/textreader.h"
49
50 #include "testutils/cmdlinetest.h"
51 #include "testutils/conftest.h"
52 #include "testutils/refdata.h"
53 #include "testutils/testfilemanager.h"
54 #include "testutils/textblockmatchers.h"
55
56 namespace
57 {
58
59 using gmx::test::CommandLine;
60 using gmx::test::ConfMatch;
61 using gmx::test::ExactTextMatch;
62
63 class SolvateTest : public gmx::test::CommandLineTestBase
64 {
65 public:
66     SolvateTest() { setOutputFile("-o", "out.gro", ConfMatch()); }
67
68     void runTest(const CommandLine& args)
69     {
70         CommandLine& cmdline = commandLine();
71         cmdline.merge(args);
72
73         ASSERT_EQ(0, gmx_solvate(cmdline.argc(), cmdline.argv()));
74         checkOutputFiles();
75     }
76 };
77
78 TEST_F(SolvateTest, cs_box_Works)
79 {
80     // use default solvent box (-cs without argument)
81     const char* const cmdline[] = { "solvate", "-cs", "-box", "1.1" };
82     runTest(CommandLine(cmdline));
83 }
84
85 TEST_F(SolvateTest, cs_cp_Works)
86 {
87     // use default solvent box (-cs without argument)
88     const char* const cmdline[] = { "solvate", "-cs" };
89     setInputFile("-cp", "spc-and-methanol.gro");
90     runTest(CommandLine(cmdline));
91 }
92
93 TEST_F(SolvateTest, cs_cp_p_Works)
94 {
95     // use default solvent box (-cs without argument)
96     const char* const cmdline[] = { "solvate", "-cs" };
97     setInputFile("-cp", "spc-and-methanol.gro");
98     setModifiableInputFile("-p", "spc-and-methanol.top");
99
100     runTest(CommandLine(cmdline));
101 }
102
103 TEST_F(SolvateTest, shell_Works)
104 {
105     // use default solvent box (-cs without argument)
106     const char* const cmdline[] = { "solvate", "-cs" };
107     setInputFile("-cp", "spc-and-methanol.gro");
108     commandLine().addOption("-shell", 1.0);
109
110     runTest(CommandLine(cmdline));
111 }
112
113 TEST_F(SolvateTest, update_Topology_Works)
114 {
115     // use solvent box with 2 solvents, check that topology has been updated
116     const char* const cmdline[] = { "solvate" };
117     setInputFile("-cs", "mixed_solvent.gro");
118     setInputFile("-cp", "simple.gro");
119     setInputAndOutputFile("-p", "simple.top", ExactTextMatch());
120
121     runTest(CommandLine(cmdline));
122 }
123
124 } // namespace