2 * This file is part of the GROMACS molecular simulation package.
4 * Copyright (c) 2016,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.
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.
35 /*! \libinternal \file
37 * Declares gmx::MDModules.
39 * \author Teemu Murtola <teemu.murtola@gmail.com>
41 * \ingroup module_mdrun
43 #ifndef GMX_MDRUN_MDMODULES_H
44 #define GMX_MDRUN_MDMODULES_H
46 #include "gromacs/utility/classhelpers.h"
48 struct ForceProviders
;
55 class KeyValueTreeObjectBuilder
;
56 class KeyValueTreeObject
;
57 class IKeyValueTreeErrorHandler
;
58 class IKeyValueTreeTransformRules
;
59 class IMDOutputProvider
;
60 class KeyValueTreeObject
;
63 /*! \libinternal \brief
64 * Manages the collection of all modules used for mdrun.
66 * This class acts as a central place for constructing modules for mdrun
67 * and wiring up dependencies between them. This class should be the only
68 * place that contains the full list of modules, although in the future, some
69 * code (e.g., in tools) may benefit from the ability to only create one or a
70 * few modules and use them.
72 * The general idea is that each module takes care of its own data rather than
73 * mdrun having to know about all the details of each type of force calculation.
74 * Initially this is applied for simple things like electric field calculations
75 * but later more complex forces will be supported too.
77 * Currently, where the set of modules needs to be accessed, either a pointer
78 * to MDModules is passed around, or an instance of IMDOutputProvider or
79 * ForceProviders returned from MDModules. These objects returned from
80 * MDModules call the corresponding methods in the relevant modules.
81 * In the future, some additional logic may need to be introduced at
82 * the call sites that can also influence the signature of the methods,
83 * similar to what ForceProviders already does for force computation.
85 * The assignOptionsToModules() and adjustInputrecBasedOnModules() methods of
86 * this class also take responsibility for wiring up the options (and their
87 * defaults) for each module.
90 * \ingroup module_mdrunutility
99 * Initializes a transform from mdp values to sectioned options.
101 * \see IMdpOptionProvider::initMdpTransform()
103 * Initializes the combined transform from all modules.
105 void initMdpTransform(IKeyValueTreeTransformRules
*rules
);
107 /*! \brief Initializes a builder of flat mdp-style key-value pairs
108 * suitable for output.
110 * If used as input to initMdpTransform(), the key-value pairs
111 * resulting from this function would leave the module
112 * settings unchanged.
114 * Once the transition from mdp to key-value input is
115 * complete, this method will probably not exist.
117 void buildMdpOutput(KeyValueTreeObjectBuilder
*builder
);
120 * Sets input parameters from `params` for each module.
122 * \param[in] params Contains keys and values from user
123 * input (and defaults) to configure modules that have
124 * registered options with those keys.
125 * \param[out] errorHandler Called to report errors.
127 void assignOptionsToModules(const KeyValueTreeObject
¶ms
,
128 IKeyValueTreeErrorHandler
*errorHandler
);
131 * Normalizes inputrec parameters to match current code version.
133 * This orders the parameters in `ir->param` to match the current code
134 * and adds any missing defaults. It also throws an error if the
135 * inputrec contains parameters that are not recognized by any module.
137 void adjustInputrecBasedOnModules(t_inputrec
*ir
);
140 * Returns an interface for initializing and finalizing output for modules.
142 IMDOutputProvider
*outputProvider();
144 * Returns an object for computing forces from the modules.
146 ForceProviders
*initForceProviders();
149 * \brief Add a module to the container.
151 * An object may be added by a client to the bound MD Modules at run time.
152 * Both the client and the MDModules object may need to extend the life
153 * of the provided object. However, the MDModules container guarantees
154 * to extend the life of a provided object for as long as its consumers
155 * may attempt to use its the interfaces accessible through IMDModule
158 * \param module implements some sort of modular functionality for MD.
160 * \note: There is not yet a way to add a IMDModule object between
161 * creation of the MDModules container and the execution of the various
162 * initialization protocols it supports.
165 * Adding a module at an arbitrary point in the MDModules life breaks
166 * some assumptions in the protocol of the other member functions. If
167 * MDModules should not change after some point, we should move this
168 * to a builder class.
170 void add(std::shared_ptr
<gmx::IMDModule
> module
);
175 PrivateImplPointer
<Impl
> impl_
;