Merge branch release-5-1

Several changes from release-5-1 have been fixed separately
in master, such changes are omitted here.

Change-Id: Ib03a2091ad6f76e785f7a6cc9ea91493d246da3b

Corrected value of electric conversion factor f in the manual

- As suggested by Christopher Neale
- Made a new latex command for the value so that it can appear
  consistently throughout all uses in the .tex manual
- changed the value to f = 138.935 458(9), which is value and
  standard deviation calculated from NIST 2010 CODATA, as also
  used by GROMACS in units.h (full value calculated from
  numbers in units.h is f = 138.935 457 839 ... )
- If I calculated correctly, the pre-June-2015 value of the std-dev
  was correct, just the last two digits of the value itself
  were swapped in some places 58(9) <-> 85(9)
- I believe commit 1ee3b15530a5 accidentally overwrote the std-dev
  with digits from the value itself

Change-Id: I293768b657a81af0a4dc0865db8c4a1b0eebd4ad

Spelling fixes

Thanks to: Nicholas Breen / DebiChem
https://anonscm.debian.org/viewvc/debichem?view=revision&revision=6293

Change-Id: I5bfecc4cc623276b65cd2dd6e455eff697d4a227

Do not create empty man/man7 directory

Thanks to: Nicholas Breen / DebiChem

Change-Id: I42d0bf14ea774c623ccf9af687ac54bbe634da4c

AVX512 transposeScatterIncr/DecrU with load/store

Also remove one more mask

Change-Id: I0d27eb42eead92d2f50725b2f3831af7e9ee229e

Fix hwloc narrowing warning on 32-bit architectures

hwloc always uses a 64-bit integer for cache size, while
size_t might be 32 bits on a 32-bit architecture. Since we
do not expect any 32-bit architectures to come with more
than 4GB of cache, we simply cast it explicitly to size_t.

Change-Id: I35b064de60a6b20dcf9ab6af751b4721037e271c

Mark more variables as advanced

CMake 3.x introduced a few more non-advanced CUDA cache variables that
most users do not care about, so this commit hides them.

Refs #1764

Change-Id: I926b03ece080d8f8255e37dc5471275d48956d49

Fix compile error for IBM VMX on ppc64

gcc-5.1 does not provide vec_mul(), only vec_madd().
This has already been fixed in the master branch with
the new implementation.

Fixes #1812.

Change-Id: I02c691b16bc0827cfa799aa99fd4fc7b84045e3f

Fix SIMD detection on new AMD AVX CPUs w/o fma

We earlier assumed that all AMD CPUs had fma
support and could use AVX_128_FMA, but with this
change we properly detect it and revert to
AVX_256 otherwise.

Fixes #1906 for release-5.1.

Change-Id: I803a6bfa75c5069b73023688a477cfca03c1ddcd

Correct nrdf for 1D/2D systems

With COMM removal, grompp would subtract degrees of freedom also for
VCM groups with fully frozen dimensions, i.e. 1D/2D systems.
Also fixed division by zero for groups with #DOF=0 with VV.

Fixes #1923.

Change-Id: I0ba2535df0495947d9bbb6ee1ef29f519635c221

Fix unit test segfaults with gcc-6.0

Work around a bug in gmock-1.7 by adding the flag
-fno-delete-null-pointer-checks to for the gmock
and testutil sources.

Fixes #1911.

Change-Id: I3e2d3fde99898191437d9eed986659f6f6fb056a

Fix incorrect return type for scalar version of simd andNot

The return type for the integer version was incorrectly set
to float. Fixed by properly adding separate integer, float
and double versions. This has never been used in any higher
level code yet, so it won't have caused any problems.

Refs #1930.

Change-Id: I84a4dfe51cf898664631ec2b0313ecf8e11a4a3a

Added references to ProtSqueeze and g_membed to the user guide

Change-Id: Ia05a3710eb170173bae9257993f05d2377c02db6

Fix coolquotes

The unseeded generator produced always the same number.

Change-Id: I7ce678e61b75aa789e7535bfedb0021c1f2117e7

Revert to treating -1 as "generate a seed" in .mdp

Commit 2d0247f6 changed the grompp behaviour while refactoring the rng
seed setup, without mentioning it in the docs or commit message.
Several of the tools also had their functionality changed, but in most
cases there was no previous documentation and now there is
some. However, we can't change grompp behaviour without catering for
all the scripts whose behaviour would change.

This change restores the previous functionality for grompp.

Change-Id: Id77098e1f03ebf96207f8c9908b2dada28c59e6a

Fix includes in readpull.cpp

Change-Id: Iaf055fc196a17dd59795c6579caa0afc210a0c82

Add option for external pull potential

An module other than the pull code can now provide the potential
and force for a pull coordinate, when the pull coordinate type is set
to external-potential. It is the responsibility of the module to call
a function that applies the scalar force to the atom in the groups
involved in th pull coordinate. A registering function has to be
called by the provider module at setup. This function checks the
consistency of the pull setup and checks a provider string the user
provided against the string of the providing module.
This functionality is intended to be used by the AWH module.

Removed the set_pull_coord_reference_value() function. Instead of
externally modifying the reference value, an external pull potential
should be used.

Change-Id: I231fe43fe3a2148a7a77d979b0fdb234d3de5c9b

Add test and improved docs for custom checkSequence

I tried to use this, but ran into problems if I passed a non-unique id
that was not NULL to the compound checker. Passing NULL is the
behaviour for the implementation of checkSequenceArray, and that seems
to be often appropriate here, too.

Change-Id: I9b02e4632085ad995ef67bbb52959f1533d56c34

Revert behavior for handling SETTLE errors

SETTLE errors were changed to fatal in 5.1.2. This causes problems with
minimization of some systems, since a SETTLE error is likely to occur at the
beginning of the minimization. This commit reverts that such that
it only issues a warning as before.

Fixes #1915.

Change-Id: I43ea1878cb61aecf5ae2c5ebcc67fcbda27ba50f

refactor CUDA texture-based table initialization

Also makes conditional the allocation of non-bonded parameter table and
related CUDA texture, skipping it when combination rule kernels are
used.

Change-Id: I6f7e25ff3906d8f92da22d2a227292bcb9e9fc74

avoid CPU oversubscription with GPUs

If the user specifies the number of OpenMP threads on the command line N
such that N*#GPU > #hardware threads, mdrun will determine the number of
thread-MPI ranks to start based on the number of GPUs without taking
into account the total number of threads that such a run would use. This
will result in CPU oversubscription.

This change will reduce the number of tMPI threads started to prevent
oversubscription.

Fixes #1338

Change-Id: I6bceaa2230686fee2e410caf8ae97c43847a5c28

Fix multiple tMPI ranks per OpenCL device

The OpenCL context and program objects were stored in the gpu_info
struct which was assumed to be a constant per compute host and therefore
shared across the tMPI ranks. Hence, gpu_info was initialized once
and a single pointer pointing to the data used by all ranks.
This led to the OpenCL context and program objects of different ranks
sharing a single device get overwritten/corrupted by one another.

Notes:
- MPI still segfaults in clCreateContext() with multiple ranks per node
  both with and without GPU sharing, so no changes on that front.
- The AMD OpenCL runtime overhead with all hw threads used is quite
  significant; as a short-term solution we should consider avoiding
  using HT by launching less threads (and/or warning the user).

Refs #1804

Change-Id: I7c6c53a3e6a049ce727ae65ddf0978f436c04579

LJ combination rule kernels for OpenCL

The current implementation enables combination rules for both AMD and
NVIDIA OpenCL (also ports the changes to the "nowarp" test/CPU kernel).

Like in the CUDA implementation, all kernels support it, but only for
plain cut-off are combination rules used.

Notes:
- On AMD tested on Hawaii, Fiji, Spectre and Oland devices;
  combination rules in all cases improve performance, although combined
  with the i-prefetching, the improvement is typically only ~10%.
- On NVIDIA tested on Kepler and Maxwell; in most cases the combination
  rule kernels are fastest.
  However, with certain inputs these kernels are 25% slower on Maxwell
  (e.g. pure water box, cut-off LJ, pot shift), but not on Kepler.
  This is likely a compiler mis-optimization, so we'll just leave the
  defaults the same as AMD.

Change-Id: I05396e000cdf93c1d872729e6b477192af152495

Re-enable i-atom type local mem prefetch in OpenCL

For reasons unknown this has been disabled in the original OpenCL
implementation. However, it turns out that prefetching does have
substantial performance benefits, especially on AMD (>10%) and in some
cases on NVIDIA too (although not on Maxwell).

This change re-enables prefetching code-path and turns it on
for AMD devices. For NVIDIA the decision will be revisited later.

The GMX_OCL_ENABLE_I_PREFETCH/GMX_OCL_DISABLE_I_PREFETCH environment
variables allow testing prefetching with future architectures/compilers.

Change-Id: I8324d62d3d78e0a1577dd3125edf059d3b311c2f

Correct mdrun tMPI (non-)parallel error message

Fixes #1931

Change-Id: Ifad46c7f62099a2cd80d70ccbe46bf3f2b5751e0

Make includesorter.py work with Doxygen 1.8.10+

Newer Doxygen versions write relative paths to the XML output, so make
the filtering that limits what to parse work with relative paths, using
XML entries that should give a relative path with all Doxygen versions.

Change-Id: I6172c625cd7e5293764b9a4dfeacf1c7de10ce78

Describe AMD APPSDK 3.0 issues in the guide

No better solution has been found to the issue, so we can only warn
users and provide advice for a workaround.

Fixes #1921

Change-Id: I4f6c8efa72e02d39ccb59977a2093802c0590563

Fix error reporting for documentation build

Make explicit file handling in the documentation build work, now that
the current working directory is managed separately.  In some cases,
this could cause missing some warnings (probably mostly from
check-source) and marking a build incorrectly as successful.  However,
the main effect was that reporting of the unsuccessful reason back to
Gerrit was less than ideal.

Fix issues that had been introduced while the checking was not working.

This can possibly be improved further by adding methods in releng for
wrapping all this file access, but that can wait for later, once there
is a bit of time to design a good, stable API for that.

Change-Id: I2f1dc927a2134365dba15a6b7dcec1a754b58433

Clean up vsite code

No functionality changes. Added some const correctness. Reduced the
scope of a lot of variables.

Change-Id: Ia6caba6686ecf9348223932f7cc812e601d0109c

More support for release workflow

- Add a matrix used for the deployment tests for tarballs.
- Make gromacs.py support builds with static linking, make it
  possible to run the tests (including regression tests) with 'make
  check', and do an install with a normal + mdrun-only build.
- Adapt the website build script to releng changes, and do not run all
  the Doxygen targets twice for the website build.

Some TODOs in gromacs.py identify what would still be necessary to get
to the level the Test_Tarballs_* jobs currently do, but that can also be
improved separately.  Now this, together with the releng changes,
produces something useful reasonably robustly.

Change-Id: Idca8974c5c9e2ee20e04f0a309e0f306a9c615d6

Ignore empty GMXLIB

If the GMXLIB env.var. was set, but empty, it triggered an assert
(at least) in pdb2gmx. Now this case is silently handled as if
the variable wasn't set at all.

Part of #1928.

Change-Id: I51646ec659d119b20b7cb429f9292bffa5775aaa

Add real/SIMD template for LJ-14 pairs

Added a real/SIMD templated code-path for plain LJ 1-4 pair
interactions without free-energy, virial and energy.
We should add virial and energy, but the virial requires a different
treatment of the shift forces, that we should introduce at the same
time for angles and dihedrals.
The gives a few factors speed improvement. The main improvement comes
from simplified analytical LJ instead of tables; SIMD helps a bit.

Change-Id: Ifdefd875e32f6c5b7d609c3bd66a864fe77fc13e

Removed srenew that was making array smaller incorrectly.

An array y was first allocated to 6 entries, then reduced
to 3 entries using srenew, but after that the indices
3, 4 and 5 were used anyway. Removing the srenew fixes the probem.

Fixes #1920

Change-Id: I3f35f0af3f56d33435bf54a9c7cda6273fb3f05a

Avoid assertion failure in pme_loadbal_do

When ir->init_step % ir->nstlist != 0 an assertion would fail in
pme_loadbal_do. If it does not fail every time, the assertion can safely
be replaced by a conditional return, which is now done in this case.

Fixes #1870.

Change-Id: Ie586022c87d3993f2de5a6fdf9210ae56aa7855b

Make gmx vanhove work without PBC

Change-Id: I37d59fc7edf6d6da36778ac3f342c71a876e2872

Clean-up code for swapcoords.

Fixes one of the TODOs left from change gerrit.gromacs.org/#/c/5660/

Change-Id: If9f443447f2cfd7fb757a386447a51a9f958dd41

LJ combination rule kernels for CUDA

Implemented LJ combination rule parameter lookup in the CUDA kernels
and enabled it for plain LJ cut-off, as was already present in the
SIMD kernels.

Change-Id: Ic1222e9433eb21e8bd1fb81658ebbbfe42c1d2c2

Moved pull initialization back to mdrunner

Commit 29943fe5 moved the init_pull and init_rot calls from
mdrunner() to do_md(), but both are needed in energy minimization
as well and init_pull should be called before init_constraints().
Added an assertion to init_constraints() for pull initialization.

Fixes #1924.

Change-Id: I9420940f772afc2be93416619f464a6ef7472372

Add more const correctness

This helps reduce the size of a future change to how we handle reading
and writing conformation files.

In particular, the way we write symmetric low-level I/O routines that
do both reading or writing according to the value of a parameter
requires that parameters for values to be written are passed as
non-const at some level of the call stack. In several cases, this
change moves that point lower down, so that routines whose job is
writing complex data to files take const parameters. Improving this
helps understand what is going on when e.g. an analysis tool reads in
coordinates, perhaps modifies them, and then writes them. The down
side is that we need to use some const_cast as we get close to the
low-level routines.

Change-Id: I5c22d8681fa5f247f302409bf67b062ebc9fe766

Fix warning from GCC 6

These omega variables could not actually be used when uninitialized
but the code is anyway improved by this change.

Change-Id: Ic038ecf1f5ae1ff1f43bbc77b6de023d94b020b2

Improved some mdrun integration tests

Made some better naming.

Avoided using the whole of OPLS/AA just to get water parameters, which
may have been unstable when calling grompp many times in a binary, for
reasons unknown.

Removed warning suppressions for unsupported icc compiler versions.

Change-Id: I5107fca64f6a73e3c0b85146f28347453ded705a

Warn about OpenMPI 1.8.6

Refs #1897

Change-Id: I529e59f095bb38a35569da57d395658e038e3f8f

Fix non-normalization of FFT autocorr

The function many_auto_correl incorrectly divided by the number of
data points, leading to incorrect absolute values of autocorrelation
functions, as e.g. out by gmx analyze -ac -nonormalize.
Also simplified the FFT size calculation.

Added a test for -nonormalize functionality.

Fixes #1914.

Change-Id: If0961c91d6e0a78f88e4390ce55b2a56143dd5e5

Fix typo in Fourier dihedral equation

This now matches the version in the preceding subsection, and seems
consistent with thoses equations and the code in convparm.c.

Change-Id: Ied7b74d07463f36ad6b1aadfb1c5f4a5b9e776a7

Added pull coordinate geometry angle-axis

The new geometry is described in the docs.
Some checks in readpull.cpp where reorganized since adding new
geometries made some old logic a bit convoluted.

Change-Id: I310891740d65e3c47948a06726e7a9e7c974bc44

Added pull coordinate geometry dihedral (angle)

How to use the new geometry is explained in the docs.

Change-Id: Ic5fc106df923857e5e94409583e3793dcd7f9fb8

Added pull coordinate geometry angle.

A new subsection was added to the docs explaining the new geometry.

Change-Id: Ie3ab77f6f2e9f35b3afe1a99414a9c0b2bafb991

Fixed division by zero in polarize.

Replaced dr2*gmx::invsqrt(dr2) by std::sqrt(dr2) in polarize.
invsqrt does not work for values that can be zero.

Change-Id: I1ae75bb64c35ac51afd237b2138cbacc86a425dc

Clarified license text in the .tex manual

This is the remainder of patch https://gerrit.gromacs.org/#/c/5659/
(which is to be abandoned), which was not already fixed
by https://gerrit.gromacs.org/#/c/5655/

Change-Id: I170f63d94a0ab0c4258f8a66b8c1a524aaf975ef

Various minor spelling fixes - messages and code comments only.

Subsequent revision: Also update copyright years as necessary.

Change-Id: I04ad9fd4948ea6c8598262cdd4b1fb87e5348248

Fix membed initialization

In commit 29943fe59dca471b2  the membed initialization was copied over to do_md,
but in do_md, state_global instead of state (which is still a NULL pointer
at that point) needs to be passed to init_membed.

Change-Id: Ic43e66f054f0e9b139ff9b402eb7cf4d08d20731

Fixed OpenCL e/fshift clear kernel first invocation

During refactoring the energy/shift force buffer clearing OpenCL kernel
call was moved before the kernels get initialized. This could never
work, but the lacking error checking in the code and (apparently) lax
runtimes allowed it to go unnoticed for some time.

Change-Id: I4bc13bf5e4e33c0797d6f08ad3cd7185d548af96

Merge branch release-5-0 into release-5-1

Change-Id: Ic214521c998b898175a3c9317534008faccef280

Add tests for empty strings in refdata

This is apparently fragile functionality when trying to rewrite the code
with something else than libxml2.

Change-Id: Iff6a3e76745e43212fc67a6481341d6413e72873

Fix tabulated wall potential.

The dispersion and repulsion of the tabulated wall potential were a factor
too small by a factor 6 and 12 respectively.

Fixes #1912.

Change-Id: I56c4b81e7f40f5d5dc9e0f2a184cf4a1bb01e3e5

Remove restriction on dihdraltypes order

Wildcards are allowed for dihedraltypes by grompp. But it was assumed,
with only code and no user documentation, that wildcard types come
after fully specified types. If this order was switched, grompp would
silently choose the wildcard matches. Now the match with the most real
atom type matches is chosen.
Added a paragraph to the bonded parameter section in the manual on
dihedral wildcard matching.

Fixes #1901.

Change-Id: I7e44ff19e4d069d1b186ea470438a47d48f1a72d

Support replacing solvent in insert-molecules

Make it possible to specify the solvent (or other set of atoms) with
-replace (as a selection) for gmx insert-molecules, and make the tool
replace residues from this set with the inserted molecules, instead of
not inserting there.  It is assumed that the solvent consists of
single-residue molecules, since molecule information would require a tpr
input, which might not be commonly available when preparing the system.

Add a basic test for the functionality.

Change-Id: I3a60cdcec9a1675a116c13d83465b07a69b3388d

Update use of isnan in linearalgebra

In C++11, cmath now provides std::isnan, so we should use it.

Change-Id: Id3c9f3103722c2ae94c93f09cea74a27852bcbc2

Find hwloc quietly if it's already been found

Change-Id: I9463decfd16731ddb8da4e6badb679bd53bf8900

Split off Fermi kernel from CUDA kernel

There are now two CUDA nbnxn kernel files, one for Fermi and one for
all architectures after Fermi. Compute capability 2.x does not support
many features used in the newer kernel, so this change removes many
preprocessing macro from the code and simplifies the shmem math.
This reorganization is even more useful with the coming combination
rule kernels that add more macros and shmem branching.
Also removed the ci_offset variable and replaced &~4 by &3.

Change-Id: I6fa50970215ad32e4d3e4430af04b403db83abc2

Fix MPI_IN_PLACE test

Failed for Intel-MPI because cmake adds libmpicxx and it depends
on libstdc++.

Change-Id: I7b4e485502779ce923fe511ed99ad224166b9aa1

Add infrastructure for tests to compare mdruns

These can now compare energies and forces within a given tolerance,
though nothing in this commit actually uses the new infrastructure.

Added and updated input files for various kinds of argon, water,
peptide and lipid boxes.

Added some functionality so that scoped_ptr can work with the legacy
APIs for reading files and managing data structures.

Change-Id: I463f7902568ec778871e0ae15898f240ff51d3ad

Merge branch release-5-0 into release-5-1

Change-Id: Icb8d7bc8d5639a824fe5f45e72eede3ec7b9d3c9

Improve AVX512

Reduce usage of masking when not necessary.
Fixed FPE in debug mode.

Change-Id: I88e6324535471669faf1c30642128242d5f54e7d

Require 2015 version for MSVC

C++11 support is best in this compiler in its latest version, which is
an acceptable compromise for this platform.

This is not good for a CUDA build, which won't officially support MSVC
2015 until CUDA 8, which is unlikely to be released before GROMACS
2016. Thus, there is likely to be a few months where CUDA-enabled
GROMACS 2016 cannot be built by a supported MSVC host compiler.

MSVC 2015 adds warnings for illegal implicit narrowing of double to
float, when used in a brace initializer. In some cases
* we intend the interpretation as real, which is now explicit
* we can just use double
* we can suppress the warning
* in some test code, it is more convenient to initialize as double and copy

Change-Id: Ic6b2f9165b6f1aaa3dc59ce05cd6ffb3abe8861c

Fix Ewald 3DC dipole with Verlet scheme

With the Verlet scheme, the dipole for the Ewald 3DC correction was
summed twice over the PP-ranks at non-NS steps, leads to a correction
that was a factor #PP-ranks too high.

Fixes #1916.

Change-Id: Iefda70ea42f1b6a4c6bd4cbc7c0982114282c15c

Random engines & distributions as proper C++11 classes

This change implements the ThreeFry2x64 random engine with flexible
number of encryption rounds and internal counter bits. The class is
compatible with the C++11 random number generators, and the GROMACS
tabulated normal distribution has likewise been turned into a
random distribution compatible with C++11, meaning they can be used in
almost any combination with the standard library distributions.
- The ThreeFry2x64 implementation uses John Salmon's idea of a template-
  selected internal counter so a number of bits are reserved to generate
  an arbitrary random stream. This makes it possible to use ThreeFry as
  a normal random engine, and even in counter mode it is possible to
  draw an arbitrary amount of random numbers before restarting counters.
- Both accurate (20-round) and fast (13-round) versions are available.
- There is a gmx::DefaultRandomEngine when we don't care about details.
- gmx::GammaDistribution has been added to work around bugs in
  libstdc++-4.4.7 headers, and to avoid getting different results
  for libstdc++ vs. libc++.
- Custom Uniform, normal, and exponential distributions have been added
  to make all results reproducible across platforms since stdlibc++ and
  libc++ do not use the same generating algorithms.
- Code using random numbers has been updated, but no changes have been
  made to turn random seeds into 64bits yet.
- The selection nbsearch unit test was a bit fragile and very sensitive
  to the coordinate specific values; this has been fixed so it should
  be resilient no matter what RNG is used in the future.

Change-Id: I47a04d03e2f264e1a6ef0aa0a2174cb464ed9af7

Support for portable hardware locality (hwloc)

Added CMake support to detect hwloc, and extended the
HardwareTopology class to read information from hwloc when present,
with the CpuInfo class as fallback. All available information is
dumped to the log file at the beginning of the run. Tested to
work on x86 (Intel, AMD, both single- and dual socket), Power8 and
32-bit Arm, and hwloc versions back to 1.5 (for older hardware).

Change-Id: Id98819254bdfdcb5725f1aee625ecb65dc40f752

Add SETTLE unit tests

Now that there is a decent seam for how to set up and run SETTLE, we
can have some unit tests that will support making further changes.

Also taught the include sorter that tuple is a C++11 header that we
now use.

Change-Id: I02c2ee19df0322e7c36dd7e65bb992d76f63a5da

Converted csettle() to use real/SimdReal template

The csettle() routine for constraining coordinates has been converted
to call a template function that is instantiated with either real or
SimdReal. On Haswell this makes settle a factor 5 faster.
Added reorganized indices to the SETTLE data structure, these are set
by a new function settle_set_constraints analogous to LINCS.
Reorganized the SETTLE initialization.
The settle fatal error no longer report the atom index of (one of the)
problematic water molecule, but the pdb dumps are more useful anyhow.

Change-Id: I61d81ad8a0add6fe234f8c7b5b44dc8c7084ace9

Fix length of SIMD4 boolean vector in reference build

While it will be correct for default builds, this length
should always be 4 even if the user changes the float
SIMD width.

Change-Id: Iea76382fa8c4819f39dc90c7177ad167e9c36ae1

Minor corrections to license file COPYING

Change-Id: I31424b0919c4229069700d097550d0caa7a64039

Use SIMD transpose scatter in bondeds

The angle and dihedral SIMD functions now use the SIMD transpose
scatter functions for force reduction. This change gives a massive
performance improvement for bondeds, mainly because the dihedral
force update did a lot of vector operations without SIMD that are
now fully replaced by SIMD operations.

Change-Id: Id08e6c83d4c9943d790bfe2a40c70fa4697077af

Use rvec4 for listed force buffer

To allow aligned 4-wide SIMD force reduction, changed the force
buffer passed to ifunc from rvec to rvec4.
Because the listed and normal force buffers now have a different
format, the listed force buffer reduction over threads also needs
to reduce the buffer of thread 0 and is thus now also used without
OpenMP threading. This extra cost is more than outweighed by the
faster listed force updates.
The current force updates in bonded.cpp will be replaced by SIMD
transposeScatter in a separate commit.

Change-Id: I9a9027d983efaf0e919634a55bf2c906aad1fa90

Remove epsfac from GPU kernel inner-loops

Multiplying the i-charge instead of the j-charge with epsfac in
the GPU kernels removes a flop form the inner-loop. On Maxwell with
CUDA the gain is 2-3% in the force kernels and 3-4% in the energy
kernels (which is more than a flop, probably a register less is used).
There is a single additional division in the energy kernel, but the
gain more than compensates for this.

Change-Id: I6924b32f3f61d3b7bbe532f95361b7fefb577609

Minor code reordering in GPU kernels

Updating bCalcFshift just before use instead at the top of the kernel
improves performance by 1-2% on CUDA. This also improves readability.
Making specialized (no)shift kernels will only add 1% gain.
Also updated the OpenCL kernels for consistency and readability
(the perfromance impact is negligible with current hardware/compiler).

Change-Id: I309f90ad61e5815726d55254e2cd38d5e4e7662d

Fix leftover load() calls in the simd util headers

Commit dc2488c8 left behind in a number of headers load() calls. This
change replaces these wiht simdLoad() fixing compilation on ARM NEON,
IBM VSX, VMX, and QPX.

Fixes #1909

Change-Id: I6b5e70f738380ea95d862c207ba45bb48fe2c805

Fix compilation with Simd4Float, without Simd4Double

Commit 971335dc added Simd4Real in nbnxn_search.cpp when Simd4Float
was supported, but Simd4Double not, e.g. SSE in double precision.

Fixes #1910.

Change-Id: I92464aa9e63c08901e272c7b9337bfe8c741b0aa

Separated test of CDATA end tag handling from others

Also documented what we're trying to achieve with such a test, since
it has twice been not obvious to me while simply reading the code.

Also clarified where it is and is not appropriate to test how the XML
code handles '\r'.

Change-Id: Ic17af40e2a64477a20b168ade666af588e7a3b41

Replace cppcheck suppression in checkpoint code

We used to have this suppression, and somehow this didn't get flagged
by Jenkins before the recent cppcheck update to support 1.72.

Change-Id: I1805937210d6cb6b46a6daa4d108b8d763c49267

Move init of specialized dynamical algorithms to do_md

These algorithms only make sense with dynamical integrators, so they
should be initialized in do_md.

Noted TODO for essential dynamics, which has a non-trivial interaction
with the constraint code, and no test coverage.

Noted TODO for swapcoords, which has no clean-up code.

Change-Id: I584ad264ef1b0ecfcb4fb95bf7b595a02d668114

Separate table construction

Construction of tables for the group scheme, pair interactions and
dispersion correction are now separated. The resulting tables are
never re-used for something else. This uses slightly more memory, but
makes the logic rather more simple. Some of the tables are now held by
reference by their owners, rather than by value, which might improve
cache locality a little.

With this change, we can implement the table support for the Verlet
scheme without getting involved with the group-scheme code, and will
have an easier time removing the group scheme.

Refs #1666, #1852

Change-Id: I8ca608f0e41b02723e6080b80b04d9e7ff048900

Merge "Merge branch release-5-1"

Fix OpenCL table interpolations of Ewald correction

The interpolation function used in the tabulated Ewald kernels was
implemented incorrectly. This did not affect results because by default
anaytical Ewald is used.

Change-Id: I4449ad13c06d9511725be9b19eff98c8d647f622

Merge branch release-5-1

Change-Id: I8935d1351656d6b0888e923c8316fce5d5c27cde

Make both ftp and http links for source download

Apparently some people aren't allowed to use ftp.

Change-Id: Ia9270a2692ca3e95569f1de938f1ba0a8fab1a5f

Fix excessive list splitting

Due to a possible integer overflow, the pair list splitting code could
end up over-splitting pair lists and causing large performance
degradation. Due to the larger processor count, runs using AMD GPUs,
using 100k+ simulation systems are more prone to suffer from the issue.

Fixes #1904

Change-Id: I29139ec80aa75c78fa93de0858f7c60cdae88d5b

Fix compiler warnings with ICC

Fixes two compiler remarks:

remark #111: statement is unreachable (nbnxn_search.cpp)

remark #181: argument of type "int" is incompatible with format
             "%3u", expecting argument of type "unsigned int"
             (atoms.cpp and espio.cpp)

Change-Id: Id118e9061c8634e06e5be99600810e74e17e23a6

Fix copyright

Fix COPYING to say that GROMACS is also covered under later versions
of LGPL.

Change-Id: I96e9608509d0b11d03813a094f82ec857ae13861

Small update to reference manual

This line was so long that the last word was not even on the page.
Added a space so that proper text wrapping is now possible.

Change-Id: I074174a500039b47ef93896fa742f24d6f65baac

Updated reference manual text about PME tuning

Also removed some junk that has long been commented out.

Change-Id: I0d1d962bd1d41fb35131c64b2c5d77edf36c8791

reintroduce supercl interaction mask constant

5822b64 removed the constant and placed the long expression directly
into the conditional. This made the code slightly less readable,
so this change restores the constant, but at function scope.

Additionally, the change is adopted in OpenCL kernels too to keep the
GPU kernels in sync.

Change-Id: I405b5ff73019b750e5304ff925f12edefc051573

change macros to constants in GPU code

Following the CPU-side changes, further changes and de-capitalization
was done.
Assertion for CUDA warp size to match the cluster split was also added.

Change-Id: Ib4d6c2f1f5680ccaeb7b4065d10ac8ffced334fb

Renamed nbnxn static const ints

Renamed all static const integers in the nbnxn code from Hungarian
to CamelCase.
Also replaced two functions in nbnxn_atomdata.c by a template.

Change-Id: Ie465c6252767f3042054452d2b85759a14f7d7e9

Same fix for XeonPhi (double precision)

Breakege was caused by I646b8b37f14e498f1a9a3a6e0d7a45761951ee0e

Change-Id: I7ed633fc1a02e4fff08a52f8ab509a0392de60fa
Signed-off-by: Alexey Shvetsov <alexxy@omrb.pnpi.spb.ru>

Merge branch release-5-0 into release-5-1

Change-Id: I4268a1ced2a874dc2844decf6c2853590230ad30

Changed norm(...) > 0 condition to norm2(...) > 0 in init_ekindata()

Fixes issue #1856 on K computer

Change-Id: Ib4f301f17124077bb6cc2aa6b955d01ccdcaec1b

Replaced pull-print-com1,2 mdp option with pull-print-com.

Changes were made to the pull output order and naming.

Change-Id: Ic1a5ef460c0979b8ea72a11e96adf51b4bf3f6de

Improved/enlarged Fig. 4.9 in manual for better readability

- Made the whole figure somewhat larger
- Made the axis text larger
- Added/updated the units:
  - kJ mole-1 -> kJ mol-1
  - added degree unit
- Changed plot type from .eps to .pdf

Change-Id: Ibc77c661ed89d9c814d3e2da8490f8b472c4c828

Make GROMACS compile again on XeonPhi

Breakege was caused by I646b8b37f14e498f1a9a3a6e0d7a45761951ee0e

Change-Id: Ia73fbd80cf1366cfa7eb6599cc3943f39e195179
Signed-off-by: Alexey Shvetsov <alexxy@omrb.pnpi.spb.ru>