Pass gmx::ForceFlags to CPU nbnxm dispatch code

Also added one last missing flag, ForceFlags.computeDhdl.

Change-Id: Ibcfea7f3975685f2219c5a5e3c8e45c716be1789

Pass the new gmx::ForceFlags to the nbnxm module

- Changed in NB kernel dispatch
- NB GPU transfer launch and wait

Change-Id: Idd2738797ddcdd372e90cdfcc066b056a29d8de2

Remove incorrect check with rf-zero

Change-Id: Ia6ae62bc128c38470d42fd1ab919215ce13bc6ec

Domain decomposition and PME load balancing for modular simulator

This change introduces two infrastructure elements responsible for
domain decomposition and PME load balancing, respectively. These
encapsulate function calls which are important for performance, but
outside the scope of this effort. They rely on legacy data structures
for the state (both) and the topology (domdec).

The elements do not implement the ISimulatorElement interface, as
the Simulator is calling them explicitly between task queue population
steps. This allows elements to receive the new topology before
deciding what functionality they need to run.

This commit is part of the commit chain introducing the new modular
simulator. Please see docs/doxygen/lib/modularsimulator.md for details
on the chosen approach. As the elements of the new simulator cannot all
be introduced in one commit, it might be worth to view Iaae1e205 to see
a working prototype of the approach.

Change-Id: I1be444270e79cf1391f5a228c8ce3a9934d92701

Relative entropy as density similarity measure

Adds a new density similarity measure that uses the relative entropy of
two densities to compare them to one another. Also adds the gradient of
this measure.

refs #2282

Change-Id: I3be18ec3fad3c2551ce7c9fe896dc257d5789cdf

Add gmx convert-trj

Adds a new module gmx convert-trj, aimed at providing the minimal
functionality to convert different GROMACS supported trajectory formats
into each other, as well as supporting selections to chose atoms for
writing to disk.

The tool is based on the OutputManager and OutputAdapter framework for
writing new trajectory files and for setting meta information in the
coordinate datastructures.

Refs #2115, #2396

Change-Id: I1413721faa1383a9d3b8f14cdbea2bcc3089273f

Extract Cartesian data from DDRankSetup

The data related to the Cartesian MPI communication setup is now
extracted from DDRankSetup to CartesianRankSetup.

This change is only refactoring.

Change-Id: Ifa05c6fecad7ad0e80d0b64ee7f0b1df5f65b327

Rename DDSetup to DDGridSetup

Moved the definition of DDGridSetup and declarations for
omdec_setup.cpp to a new file called domdec_setup.h.
Also added documentation.

Change-Id: Ic14416e8acfa3b832ceff4bf529abdaf8ceb3728

Add option to normalize densities prior to density fitting

Adds an mdp option that makes sure that the sum of voxel values
for both input densities is unity.

refs #2282

Change-Id: I9d3c63d7548bca72a083cdafa6818ded4fcdd400

Clean up PME rank variables

Clarified the naming of several PME rank count variables and added
a boolean to DDRankSetup that tells whether we use separate PME ranks.
Fixes several misuses of npmenodes, all of which did not cause issues
in pratice.

Todo: Remove npmenodes from t_commrec.

Change-Id: I03b4547c975e3ffa354b53b35abcd91b4f3f6b26

Adding every N-steps option to densityfitting

Adds an option to the density fitting module to apply forces only every
N steps.

refs #2282

Change-Id: I8d264734e1d9a82ff9717d270552e38a28b34729

Refactor init_domain_decomposition()

The creation of both the whole (Cartesian) communicator and the PP
and PME communicators is now extracted to separate functions.
The first part is done before initializing gmx_domdec_t, so the role
of nodes is know before gmx_domdec_t is created.

This change is only refactoring.

Change-Id: I6d02b8c3af55570353331fa05f998a7b6de79572

Make pme_loadbal_do() independent of t_state

Similar to I28b58df4 and other recent patches, this changes
pme_loadbal_do() to take the const box and a const view of the
positions as input instead of the full t_state object. This makes
subsequent changes to the t_state object easier.

Change-Id: I6ba3459fa1de80c0ed8bcf0310f4565be78d3d97

ForceElement and ShellFCElement for the modular simulator

These elements encapsulate the calls to do_force and relax_shell_flexcon,
respectively. Rewriting calls to force routines is not within the scope
of this effort, so this allows us not to worry about it.

This commit is part of the commit chain introducing the new modular
simulator. Please see docs/doxygen/lib/modularsimulator.md for details
on the chosen approach. As the elements of the new simulator cannot all
be introduced in one commit, it might be worth to view Iaae1e205 to see
a working prototype of the approach.

Change-Id: I0923ddfc9f8f477b8920d212bcbdcf641075f80b

Adding energy output field for density fitting

Added a new energy output field for energies that stem from the density
fitting code.

refs #2282

Change-Id: I78c64213f4b958f00a938f0a67c586c2f4e9033e

ComputeGlobalsElement for the modular simulator

This element encapsulates the calls to `compute_globals`. A new
approach to fit the global computation calls into the client approach has
been proposed in I44a7193c, but a complete rewriting is currently outside
the scope of this effort. This element therefore aims at offering an
interface to the legacy implementation which is compatible with the new
simulator approach.

The element comes in 3 (templated) flavors: the leap-frog case, the first
call during a velocity-verlet integrator, and the second call during a
velocity-verlet integrator.

This commit is part of the commit chain introducing the new modular
simulator. Please see docs/doxygen/lib/modularsimulator.md for details
on the chosen approach. As the elements of the new simulator cannot all
be introduced in one commit, it might be worth to view Iaae1e205 to see
a working prototype of the approach.

Change-Id: I35f310daed6355dc9ea6b574d855929ec2fee6fe

Clean up commrec uniformly

This was leaking commrec from gmx help during the building of
command-line completions, which the ASAN build sometimes complained
about.

Change-Id: I4ce53ba13187ea71102f5a93118978525da7061a

Add gpuupdate and gpucomm OnDemand tokens

These allow triggering the GPU update-constraints and the GPU direct
communcation development features (both PP-PME and halo exchange),
respectively.

Also added matrix configs for both with configs identical to those in
the gpubufferops for the sake of simplicity.

Added an OpenCL config to ensure the additons of uninmpled features on
that codepath don't break things.

Change-Id: Iaf308d9b55ff10d8e43b655c2e1708eb05a5b43c

Periodic boundary correction for density fitting

Corrects the periodic boundaries for density fittinng
such that always the periodic image closest to the center
of the density is chosen.

ref #2282

Change-Id: I46d9546d7dbb659e03070b4ba6ffd401eb7368c4

Extract getting of DDRankSetup

This change is only refactoring.

Change-Id: I9a7f8ce16fe81f0016b33139933336fcb49b48b4

Add PP grid info to DDRankSetup

This enables many functions to use only DDRankSetup instead of
the whole gmx_domdec_t.
This change is only refactoring.

Change-Id: I6467481e3e1fa27c035c70d34591a6db0eeae90d

Introduce topology class for modular simulator

The topology object owns the local topology and holds a pointer to the global
topology owned by the ISimulator.

The local topology is only infrequently changed if domain decomposition is
on, and never otherwise. The topology therefore offers elements to register
as ITopologyClients. If they do so, they get an updated local topology whenever
it is changed, and can rely that their copy is valid until the next update. The
domain decomposition element is defined as friend class to be able to update
the local topology when needed.

This commit is part of the commit chain introducing the new modular
simulator. Please see docs/doxygen/lib/modularsimulator.md for details
on the chosen approach. As the elements of the new simulator cannot all
be introduced in one commit, it might be worth to view Iaae1e205 to see
a working prototype of the approach.

Change-Id: I89e4006809ac74785205c549781a90b2d8233688

Update mdrun test comparison infrastructure

This change is refactoring and extension, no use is made of any new
capabilities.

Referred to position coordinates as coordinates, in line with
agreement elsewhere. Also referred to energy terms consistently,
rather than energy fields or tolerances.

Extended the range of mdp option parameters able to be tested
configurably with the contents of the simulation database.

Made energy and trajectory comparison machinery easier to use.

Trajectory comparison is now more configurable. Sometimes the test
runner will know that it is OK to compare the data in a test frame
with a reference frame even though one of them has data that the other
does not have, or will differ, etc. So, it is now possible to compare
e.g. a test trajectory that has a frame with velocities and not forces
with a reference trajectory whose frame has both velocities and
forces, and have the comparison be valid or not according to whether
the test runner expected the mismatch in contents.

Fixed issue with trajectory frame comparison assuming that coordinates
would always be present.

More use of std::string.

Change-Id: I6926380e271344e7bd1a7c6594cc604d7a68e67a

Rename MicroState to StatePropagatorData for modular simulator

The name of the t_state replacement is changed from MicroState
to StatePropagatorData. Also, the write and read access to x/v/f
data is renamed to xxxView() and constXxxView() for additional
clarity and to avoid confusion with file writing operations.

This also removes the constraint element from the list of friend
classes, as there is no need for direct access to the state
anymore thanks to I28b58df4.

Change-Id: I1b5b966aafbfa6fe5c44c36416ebb25f38505b58

Extract DDRankSetup from gmx_domdec_comm_t

Change-Id: I5b16d624a375478928546788d1f8319e9d796100

Add reference density file option to densityfitting

Add option to density fitting to read reference density for
density fitting code.

refs #2282

Change-Id: I598b679a3c1c63b913a104d7ccff8e1643b3946d

Adapt inheritance intendation to uncrustify

As discussed in I35f310da and some other commits along the modular
simulator chain, uncrustify insists on aligning the first inheritance
with the `final` keyword, as in

template <ComputeGlobalsAlgorithm algorithm>
class ComputeGlobalsElement final :
    public                  ISimulatorElement,
    public IEnergySignallerClient,
    public ITrajectorySignallerClient,
    public ITopologyHolderClient

It was therefore decided to align all inheritances to the same level,
as in

template <ComputeGlobalsAlgorithm algorithm>
class ComputeGlobalsElement final :
    public                  ISimulatorElement,
    public                  IEnergySignallerClient,
    public                  ITrajectorySignallerClient,
    public                  ITopologyHolderClient

This formatting will be reevaluated when uncrustify is replaced by
clang-format.

This commit applies this decision to the classes previously introduced
along the commit chain. It also addresses two typos in trajectoryelement.h.

Change-Id: I03b9fd142d822f7a7523e351f08a4044bff2cff4

Make do_constrain_first() independent of t_state

This changes the do_constrain_first() function to directly take the
required state data (natoms, x, v, box, lambda[efptBONDED]) as input
instead of a pointer to the full t_state object. This makes subsequent
changes to the t_state object easier.

This also adds some trivial const qualifier for the box matrix in the
constraining functions.

Change-Id: I28b58df45481549cd93076334b0778e23e228154

Split off getting of the DD setup

Change-Id: I19e956b04a14df75e91b721954bba52abe49d1d8

Collect settings for DD in DDSettings

This change is only refactoring.

Change-Id: Icb32f8f5211ff3033c4d5c70afbfd24d4a8ffbf3

Extract getting of DDSystemInfo

This change is only refactoring.

Change-Id: Ie7127f0fd8219108de1385978a781a3c15fc6239

Simplify handling of DD bonded distances

To simplify and clarify the DD setup code, we now always always store
the systemInfo.minCutoffForMultiBody and use a separate flag to tell
if we should increase the cut-off distance for bonded communication.
There is a minor behavioral change in that with large domains and
bonded communication filtering or DLB, the bonded cut-off is now
5% of the bonded cut-off longer as the margin is now included.
This has a negligible effect on performance in all cases.

Change-Id: Id409353c517181ac56e8d3f1f36c22c705aa8077

Move constraint and bonded filtering info into DDSystemInfo

Change-Id: Ic3cc4b5da309ab3f8e719739a5a162cd2c2d62bf

Bump ARM HPC toolchain in post-submit to v19.3

Change-Id: Id27642e016ee47888c7c3a7dac426309288aff9e

Fix double overload in LocalAtomSet code when gmx::index is int

In cases where the int type matches the gmx::index (aka ptrdiff) type
two overloaded constructors in LocalAtomSetData became the same, as well
as the overloaded add method. This is now avoided by an enable_if
clause.

Change-Id: I139c7feff141928f71c71d6f6db3ee1bffe325b7

Add back a clang-tidy exception

Previous commit (658819) removed a clang-tidy exception that was
needed with clang 8, but the SIMD + GPU path have not been enabled
yet with clang-8 tidy. To be able to enable the pending CI configs
this change adds back the nolint exception.

Change-Id: I6be4856c64161a1a5f57849d13b4cdea335cc3bb

Read mrc file using file name

Provide a higher level class for reading an mrc file given file name and
convert the information read from the density file into coordinate
transformation and three-dimensional density data.

Change-Id: I30ac5c044c5433043224be288340b794455f8fc9

Move update grouping to DDSystemInfo

Change-Id: Ida28b87a4a0ff004d25ca4d6f0f78f30a19536dd

Inline functions of calculations of FE interaction types.

The functions are so far only templated on the real
requirements of the softcore that is used.

This is one step towards templating the calculations
for SIMD.

Refs #2997.

Change-Id: I3fd119dce30f95eba4f8cd6f139f99260acf0e22

Change DD pme force buffer to std::vector

Change-Id: I8cf00535fab713911e3dbb24a4daed3eeb0b46ee

Initialize all members of gmx_domdec_t

Change-Id: I78958c9c0e72de77d6226685401b9a5ff19434aa

Centralize more MD/force shedule flags

This change creates an MdScheduleWorkload top-level class to store
information related to per-step, per search-phase, or per run schedule
information and moving the previous PpForceWorkload call into it.

This change also migrates force flags bitmask into a new class
stored in the MdWorkload class.

Change-Id: I142a76d9515e82b16b9a4100ffb00469d7f3d3c2

Extract DDSystemInfo from gmx_domdec_comm_t

This is only refactoring.
This is preparation for separating the DD grid choice from
the initialization of gmx_domdec_t.

Change-Id: I83d8ecdd1a4d21d0a791c463d092f66d0b28df51

Unify manual and auto DD grid choices

This change is only refactoring, except that now an assert is added
for consistency of the auto grid choice with the cell size limit.

Change-Id: I36306495881826de809e3d048d0ced52543674e2

Remove gmx_domdec_t from dd_choose_grid()

This change is only refactoring.

Change-Id: I2d658739da1f0ff185ab89a5efe2361da0c7fa82

Extact UnitCellInfo from gmx_domdec_t

This required adding an explicit constructor for gmx_domdec_t,
which is currently incomplete.

Change-Id: Iba39fec36945465bffaf97d499c56501500ecb3b

Python wrapping code for gmxapi mdrun bindings.

Expose the mdrun bindings to the Python interface.

Updates workflow.py and context.py from 0.0.7.4 tagged release of external project.

Refs #2893
Refs #2912
Refs #3027

Change-Id: I44a66646eceb191a2ec755169671af326f22526a

Introduce energy element for modular simulator

The energy element owns energy related data (energies for output,
tensors for virial and pressure, mu_tot vector), and holds pointers
to energy-related data owned by ISimulator (enerd, ekind). It is
responsible for calculating energies at steps, and writing energies
to trajectory files. It offers interface access to add data to the
virials, but also offers raw pointers, as not a lot of functionality
across GROMACS relies on getting raw pointers to the tensors /
mu_tot vector.

This commit is part of the commit chain introducing the new modular
simulator. Please see docs/doxygen/lib/modularsimulator.md for details
on the chosen approach. As the elements of the new simulator cannot all
be introduced in one commit, it might be worth to view Iaae1e205 to see
a working prototype of the approach.

Improved handling of MSVC warning 4244 about possible lossy conversion.

Change-Id: Iae8d666e1c38b5a817d7ef2018adb13567d4760b

Simplify energy-file comparisons in tests

Now that we use newer infrastructure, we can use lambdas.

Also made FramePairManager and ContinuationFramePairManager
more consistent.

Change-Id: I8ccc9dd3134e5e188c3ea2b8b6ab335b4f267b10

Add Parrinello-Rahman pressure coupling to CUDA Leap-Frog integrator

Parrinello-Rahman isotropic pressure coupling scheme is now
added to CUDA version of the Leap-Frog integrator. The
Leap-Frog tests are updated to check the new addition.

Refs. #2887, #2888.

Change-Id: Icf42667621c16a994e68baf5158ea4abac387928

Move initialization of clFFT

Gave ClfftInitializer the responsibility for mutual exclusion, which
means the initialization is now convenient to do alongside other
PME-on-GPU initialization tasks. This simplifies the code.

Removed mention of lazy initialization, which was not implemented at
all.

Refs #2535

Change-Id: I429767b059ddc3b4f16f2f4a8830b54ed1f49fab

Only use fshift when needed

The shift force buffer in ForceWithShiftForces is now empty when
virial computation is not requested. This will simplify refactoring
of the force buffer management.

Change-Id: I6c8e543766703b038c178a0e95ccb536f367d5e1

Add checkpointing for MdModules

Adds the possiblity to write and read from checkpoint files for MdModules.

refs #2282

Change-Id: Id3ab57d1c9a40b6e1fb0193430e5b66ffa87d787

Fix CPU detection and SIMD on modern Mac/cmake combos

With recent versions of cmake on Mac, the CPU detection
program failed to compile since it was very fragile (buggy)
with respect to how compiler flags were defined - and that
in turn will mean losing all SIMD support.
This will fix the build, and it also adds a more vocal warning
if the CPU detection code cannot be compiled at all.

Change-Id: Idcdc8b59a4252b3faf8572311d2ca63a4d2a1bc0

Constructing local atom sets for density fitting

Constructs local atom sets for density fitting during simulation
setup.

refs #2282

Change-Id: I0300376a6818f3ad8bb6474b09a9b0f1cbb46eae

Fix build matrices

Onee agent has a non-SIMD version of fftw that CMake warns us about. As
fftw is well tested elsewhere, we can just avoid testing it here.

icc 19 (only) has a codegen issue for pre-avx2 SIMD paths that mean
some exp SIMD unit tests fail. While we are diagnosing the issue,
this build configuration is changed to something that will pass.

Change-Id: I0e39cc245563b2afc3aba2b76858b17373ee91a3

LocalAtomSet construction with vectors of gmx::index

LocalAtomSet construction is updated to follow the transition to use
gmx::index whenever objects, like atoms, are of type gmx::index in
an array.

Keeping the old int based constructors for legacy compability.

Change-Id: Ic0973f5260a8de9e87bf5391040e14e533012d88

Import base code for Python-driven simulation from gmxapi 0.0.7.4.

Copy workflow.py and simulation.py from
https://github.com/kassonlab/gmxapi/releases/tag/v0.0.7.4
to unexposed submodules in gmxapi.simulation.

Follow-up change will update these sources and expose functionality
for continued development in the GROMACS repository.

Change-Id: I5f62808371b325de3286f13ee34b87f3a231b67d

Changed some variables from SCReal to real in FE kernel.

This means that there is a reduced precision in some of the calculations
when using sc-r-power 48 if GROMACS is compiled with single precision.

This is to achieve the behaviour intended by
commit a7c8e19121.

Change-Id: Iff754dfad51d1e468d96c368859db0a48ba51ebb

Implementing the density fitting force provider

refs #2282

Change-Id: Ib65679a82168f18992920f4529596db2494cccc3

Some quotes I bumped into last week

Change-Id: I071d73bbdf9dca016d5d229baaad9db770509067

More quotes

Change-Id: Ia4f11f76be7871ee304167c78366ef2602533b5f

Introduce MicroState for modular simulator

The MicroState replaces the t_state structure. It owns the position, velocity
and force vectors, as well as the box tensor. It also keeps a backup of the
previous position and box. It is therefore having more data than the pure
statistical-physical micro state definition. TODO: Decide if we should give
it another name to avoid misunderstandings!

While the modular approach aims at having as much data as possible localized,
the data in this class are so universally used in molecular simulations that it
is reasonable to keep them in a central object. Conveniently, this also is all
data commonly used in trr files.

Note that the MicroState can be converted to and from the legacy t_state
object. This is useful when dealing with functionality which has not yet
been adapted to use the new data approach - currently, domain decomposition,
PME load balancing, and the initial constraining are using this.

This commit is part of the commit chain introducing the new modular
simulator. Please see docs/doxygen/lib/modularsimulator.md for details
on the chosen approach. As the elements of the new simulator cannot all
be introduced in one commit, it might be worth to view Iaae1e205 to see
a working prototype of the approach.

Change-Id: Id1a22892c2dcf2f30127a9fd780b89d4c860640f

Introduce signallers for modular integrator

In the current do_md() function, code flow is governed by booleans
set on every iteration of the simulator loop (e.g. bNS, do_ene,
bCalcEner, etc). In a modular simulator, we want to avoid common,
quasi-global variables. This change proposes an approach to replace
these booleans by a signaller-client system. An object ("signaller")
is thereby responsible to decide when specific events are happening,
and is informing its clients about it. This approach is implemented
for five types of events, namely neighbor-searching steps, logging
steps, energy-calculation steps (energy, free energy and virial
calculation), the last step, and trajectory (x/v/f and energy)
writing steps.

These signallers are designed to run at the beginning of scheduling
for a specific step.

If signallers depend on each others, they need to be ran in their
order of dependency. To illustrate this, consider the case of energy
calculation. This signaller can run by itself, checking if the
current step is an energy-calculation step. It can also be triggered
by the trajectory signaller, since the energy signaller is a client -
another way to say that we should calculate energies if we want to
write to the energy trajectory in this step. While we could imagine
that the trajectory signaller triggers a callback which in turn
informs the energy clients about the next step, this becomes shaky
when a signaller depends on more than one other signaller. Running
the signallers in their order of dependence, on the other hand, is
stable as long as there are no cyclic dependencies (which should
never happen). In our example, we therefore require the trajectory
signaller to run before the energy signaller. The builder, which
also registers signallers which each others, is responsible to keep
the signallers in order. This is trivially done by adding the
signallers to a vector right after building them.

The TrajectoryElement is both a simulator element (as it needs to run
periodically to output state or energy information to file) and a
signaller (as it needs to warn the state / energy objects that a
writing step is coming up). The TrajectoryElement is built around the
gmx_mdoutf object, which it owns. It can therefore not easily be
split into a element and a signaller part, or into a state writing and a
energy writing part. During the simulator run phase, it is
calling its trajectory clients (which do not necessarily need
to be identical with the signaller clients), passing them a valid
output pointer and letting them write to trajectory. Unlike the
legacy implementation, the trajectory element itself knows nothing
about the data that is written to file - it is only responsible
to inform clients about trajectory steps, and providing a valid
file pointer to the objects that need to write to trajectory.

This commit is part of the commit chain introducing the new modular
simulator. Please see docs/doxygen/lib/modularsimulator.md for details
on the chosen approach. As the elements of the new simulator cannot all
be introduced in one commit, it might be worth to view Iaae1e205 to see
a working prototype of the approach.

Change-Id: Ic58a18a374ff0ec5a47d0853c514e439ba1a45a3

Fix return value

While testing on gitlab, clang-tidy complained about the return value
of a MPI detection function.

Change-Id: I09b4af9ad22d450bf48d360397e378eb10ff1610

Have selections available for mdrun only

Due to the inclusion of the densityfitting code, selections need to be
available also when only building mdrun.

Fixes post-submit build failures.

Change-Id: I53e63fb5ca9f678dc3fcfcafa180367dd6c62c20

Fix integration issues

Post-submit is warning about CUDA, and an automatic rebase missed out
on the change to remove gmx_install_headers().

Change-Id: I64d6fd149a4fa0d72cbd7e7d210b325f9f66a112

Introduce the main modular simulator interfaces and documentation

This change introduces the two main interfaces for the modulator
simulator approach. It also introduces a new module which will contain
all files introduced for the modular simulator, making differentiation
between legacy code and new implementations easier. It further starts
the documentation of the new approach.

The ISimulatorElement interface is for any element that wants to
take part in the simulator loop. It requires to define a
scheduleRun method, which allows the element to decide whether
it needs to run at a specific time step. If it does, it adds one or more
function pointers to the run queue (registration function is passed
as argument, allowing flexibility in the exact implementation of the
run queue and the number of registered function pointers). The
Simulator can hence pre-compute all run function needed in one or
more future simulation step, allowing it to then run through the
function pointers without the need for branching

The ISignaller interface is used to spread information about upcoming
steps. Signallers are guaranteed to be ran before the elements
get queried about running at a specific step. This allows, for
example, to signal to the element that a neighbor-searching step
or a logging step is about to happen. Elements can register to
signallers they are interested in being notified by.

This commit is part of the commit chain introducing the new modular
simulator. Please see docs/doxygen/lib/modularsimulator.md for details
on the chosen approach. As the elements of the new simulator cannot all
be introduced in one commit, it might be worth to view Iaae1e205 to see
a working prototype of the approach.

Change-Id: Iefc9a0c729223cc5da0b500f0b410027c14d1f68

Change tolerance for normal mode tests

Testing on gitlab showed that the comparisson for the normal mode
tests was too tight, so I relaxed the test tolerances a bit.

Change-Id: I2c0443d8deb9aa2e49559fd75f6b5fe6953d7912

Improved usability of grompp by rejecting more invalid input

Improved test coverage and code documentation.

Also renamed some test cases so people understand the intent better.

Change-Id: I06f704d2d5516888ef8c84e78b68fa8c53fa7948

Add clang-tidy nightly test with SIMD + OpenCL

This also triggers the offline OpenCL kernel compilation.

Fixed some remaining clang-tidy issues found in that configuration.

Change-Id: Ia43404aebdc690f75cf0778391b1395e8770ecaf

Move nonbonded GPU wait timing into the GPU module

Extended the wallcycle module to allow only incrementing a counter --
used to simplify the polling GPU wait that would otherwise require a
regular start/stop to accomplish the same.

Also eliminating some TODOs.

Change-Id: I058c7213ab03d3fb5b438db58f69180049771454

Add sample_restraint Python tests to `check` target.

When building as part of the GROMACS project (not as a standalone
project), add a CMake testing target to invoke pytest for the
sample_restraint package. Include this target in the GROMACS project
`check` target. Make some minor CMake infrastructure changes to help
the sample_restraint tests find the gmxapi Python package.

Change-Id: I753951c1bb5f5a9c6caa20d4d632eba214895a1f

Add integration tests for exact restarts

These tests demonstrates the extent to which mdrun checkpoint restarts
reproduce the same run that would have taken place without the
restart.

I've been working on these, and the bugs they exposed, for a few
years, but the code has been fixed for a few years now.

The tests don't run with OpenCL because they have caused driver out of
memory issues.

Refs #1137, #1793, #1882, #1883

Change-Id: I8bc441d945f13158bbe10f097e772ea87cc6a559

Add density fitting parameters

Adds the mdp parameters that are required for running a
density fitting simulation.

refs #2282

Change-Id: Id22ceee0f0036450173522bd5e85af06542a2c39

Access to GaussTransform3D data

Return a non-constant view to the GaussTransform3D to allow
sum reduction of Gauss transformation data via MPI. The alternative
to have the GaussTransform3D take a commrec to perform the
sum reduction itself leads to a cyclic dependency between
math and mdtypes modules.

refs #2282

Change-Id: I8b786965c3562ed84fe4e97fa0f68cf89ea0f03a

Add index group option to density fitting

Adds the option to select a specific group for density fitting, stores
the indices in the tpr file and restores them after reading.

Adapted the densityfittingoptions test to test fixture to better reflect
the increased required complexity for testing the storage of indices.

refs #2282

Change-Id: I3e33508bbcadf968b0ddf978a73e7f6cf5753b5b

Remove unused forward declarations.

A lot of stale forward declarations were cluttering up the public
interface for the Restraint module. The remaining forward declaration
should be revisited as the public / installed API is revised.

Change-Id: I1a4d8a913ade93158ad54ac99ca9da1ced66f347

Remove section checking api level

Removes the check for api levels in check-source to allow modernization
without having to account for the legacy api layout.

Refs #988

Change-Id: I12cd3f6765bc57801ff4dd81583b7836dc4f18fb

Make FEP rerun tests cover the intended lambda space

These tests are meant to look at whether a rerun can recompute
quantites for a range of lambda values, but something went wrong at
some point so that only state 3 (bonded-lambda 0.5) has been being
tested.

TODO make a redmine that grompp should reject lines that aren't
recognized as either comments or valid input.

Change-Id: I1eb405aa55bfb77c07d9a7e9c48803380c96c1df

Remove installed headers from CMake

Removed the gmx_install_header sections from CMake files, as well as the
CMake code used to add and check them.

Refs #2382, #2139, #988

Change-Id: I4525ea14d2967f83d940300daeb2ade08717ed5d

Add temperature coupling to CUDA Leap-Frog integrator

Velocity rescaling added to the integrator. None, one,
or multiple rescaling groups are supported. Tests that
include temperature coupling are added.

NoseHoover temperature control is not implemented.

Refs. #2887, #2888.

Change-Id: I1e6850eeb74de829554260fba6a6e6c1c63ceb46

Rearrange exposure of gmxapicompat simulation input tools.

* Add gmxapi/gmxapicompat.h installed header. Move the most benign API
  aspects to this header, including exceptions. Retain some extra
  details in gmxapi/compat/tpr.h and gmxapi/compat/mdparams.h.
* Move the provisional typing enum to gmxapi namespace.
* Change usages of several types to unique_ptr handles to allow opaque
  pointers in gmxapi/gmxapicompat.h

Change-Id: I9c0ae50250b6bc1a1ad3e1885356a5ae718c610e

Converted a CUDA config to use clang host-only compilation

This is distinct from the clang_cuda configuration because that one
uses clang also for device compilation. This one doesn't use the
modified build system, and will be easier to update when NVIDIA
support newer clang versions on Linux.

Refs #2899

Change-Id: I82bdf700c5553fcf367809209481ab4e27b037a2

Move gmxapicompat TPR tools to libgmxapi.

Make a few adjustments to satisfy clang-tidy and other warnings.

Change-Id: Iaede422d6726d0ca301436806b496b0bb64f2fdf

Fix pairlist bug

Introduced recently in 658819431c0dabd while fixing issues for
upgrading clang-tidy version

Change-Id: I4c0dd4212e93cccf11f2ce60fbb500923eb136e9

Reduce artificial header dependencies in restraint plugin machinery.

Change-Id: I8c9c6170d89e8395a1e9a82d03067283d4999154

Activation of density fitting and move of module notifications

The densityfitting code is activated by adding it to the MDModules.

Moved the MdModuleNotifications to mdrunutility to break cyclic module
dependencies between mdrun and applied_forces via MdModules and
DensityFittingModuleInfo.

Observe that this patch changes the default .mdp output. Updated
reference data accordingly. Note that the changes in .mdp output
lead to "check mdp file differences" warning in the regression tests.

refs #2283

Change-Id: I93c0dd85285cac3f9e94d2dfc7c7fcc434aa81db

Make use of recent changes to prepareKernelArguments(...)

Recent change allows to avoid assigning pointers to intermediate
variables before passing them to prepareKernelArguments(...)
functions in CUDA builds.

Change-Id: I3b70d7663cbd687e9f22e0422e7baa3fcdfc1bf3

Index groups for MdModules from grompp

Make index groups available for mdrun modules during preprocessing.
Notify mdrun modules with a DefaultIndexGroupsAndNames class that
translates index group names into atom indices.

refs #2282

Change-Id: Idc1bc3a5840e7b7ec3c82904416887b02417c92c

clang-tidy 8 fixes for the GPU host code-path

Tested with OpenCL and CUDA (with nvcc as device compiler).

Change-Id: I6c7f31efdf89aa863e6d831b1a5c41bfc2ef332a

Further fixes for clang-tidy clang version 8

These eliminate warnings with the standard recommended flags.

Change-Id: I541f1bfa918489219c5e02d8b0b92619e7c3f293

Misceleneous improvements to the Leap Frog tests

1. The contained with runnerss is made static so it will be
   initialized only once for the test case.
2. The member initializer is now only use constructor arguments,
   and hence does not impose requirements on the class fields
   ordering.
3. mdAtoms_.cFREEZE is set to nullptr explicetely, to make the
   code less sesitive to possible further changes in intialization
   of t_mdatoms.

Change-Id: I6e1db5cb3d5d81028480457989ce2e2b387018f4

Removes checker for installed files

As the plan is to remove the current handling of installed header files,
this removes the checker parts from the doxygen check-source script for
this attribute.

Refs #2382, #2139, #988

Change-Id: I76b519f39a5c793c9f1ea8c1eb5eebb39b4a9352

Fix clang-tidy warnings in the OpenCL build

Change-Id: I4cc6a1402d068431feab4523911a7b60e0e66d59

Fix clang tidy warnings in CUDA build

Tested with clang 7 on a build with clang for both host and device code.

Change-Id: Ie2be427ae63f90f4eec4e79cc85da3ab13fcfe4a

Update clang-tidy to clang version 8

Some new warnings are well intentioned, but not practical to turn on
in GROMACS given the limited suppression capabilities.

Rationales for suppressions of new warnings can be found in
src/.clang-tidy

Change-Id: I857d3790450e33a669c6a68ad846cae3783875e7

Merge branch release-2019

Change-Id: I71b540586e32f54f3536c5a59eb99d89ad0b3639

Improve clang-analyzer build time

The tests for selection and utility modules take a long time to
link, so these are moved earlier in the CMake so that linking
will start sooner.

Change-Id: I2fb50515d8ffe5f93e0bcf13c6a780cb5725aa3e

C++ code and Python bindings for TPR parameter read/write.

Core functionality to allow TPR files to be rewritten with altered
parameters, ported from https://github.com/kassonlab/gmxapi v0.0.7.4
with modifications to naming and coding style. (Refer to patch history.)

A lot of the C++ code is not intended to be a long term solution,
but demonstrates the use cases that will need to be addressed
as modules become able to self-describe their inputs and outputs.

However, this change is necessary to support near-term data portability
between gmxapi operations for preparing and wrangling chains or
ensembles of simulations.

Refs #2993

Change-Id: I54677c861dfb19c9f34b11d2c30456e6ee5dbe8d