Documentation house-keeping in file swapcoords.c

Modified Doxygen documentation of functions such that they now have
\brief descriptions. This way, Teemu's automated Doxygen documentation
checker won't complain about many missing brief descriptions in that
file. I am not making this patch dependent on the documentation
checker (gerrit issue #3213) because the brief descriptions should be
there anyway.

Change-Id: I4fc2c9de17b5ce8e9abdcd8fc51d002a245fbf2a

Interactive Molecular Dynamics (IMD)

IMD allows to interact with and to monitor a running molecular dynamics
simulation. The protocol goes back to 2001 ("A system for interactive
molecular dynamics simulation", JE Stone, J Gullingsrud, K Schulten,
P Grayson, in: ACM symposium on interactive 3D graphics, Ed. JF Hughes
and CH Sequin, pp. 191--194, ACM SIGGRAPH). The user can watch the
running simulation (e.g. using VMD) and optionally interact with
it by pulling atoms or residues with a mouse or a force-feedback
device.
Communitcation between GROMACS and VMD is achieved via TCP sockets
and thus enables controlling an mdrun locally or one running on a
remote cluster. Every N steps, mdrun receives the applied forces from
the VMD client and sends the new positions to VMD.
Other features:
- correct PBC treatment, molecules of a (parallel) simulation are made
  whole (with respect to the configuration found in the .tpr file)
- in the .mdp file, one can define an IMD group (including the protein
  but not the water for example is useful). Only the coordinates of
  atoms belonging to this group are then transferred between mdrun and
  VMD. This can be used to reduce the performance impact to an almost
  negligible level.
- adds only two single-line function calls in the main MD loop
- and mdrun test fixture checks whether grompp and mdrun understand
  the IMD options

Change-Id: I235e07e204f2fb77f05c2f06a14b37efca5e70ea

Remove gmxpreprocess from include path

One step closer to the final goal of only having src/ there.
This wasn't even necessary except for a few rare cases.
Also cleaned up the cppcheck rules slightly.

Change-Id: Ifc65fbb48699d0e4bb5c55b90a2c41b0255d353d

Keep CMake quiet about non-detection of optional stuff

Repeated calls to cmake would otherwise warn every time that these
optional things were not detected. That lack is only rarely a problem
that needs solving, so they should not be emitting status messages all
the time.

Change-Id: Iad6f0bdf72ba7848bcc819db6937b70b37900148

Merge branch release-4-6

Change-Id: Ie4431fb0002cd36799c8e664dcf4e096e66bc0ea

Fixed bug in parallel v/f constraining

Constraining v or f with 3 or more decomposition domains in one
or more dimensions could lead to modification of communicated
v and f components by the box size for inter charge-group constraints.
Fixes #1462

Change-Id: Idece9d2d0d8f48e65a654d5c2892fbe1ff836ba0

CUDA cut-off kernels now shift exclusion energies

This is to make the GPU kernels consistent with the CPU nbnxn
kernels, which for eeltype=cut-off modifier=potential-shift
effectively do reaction-field with epsilon_rf=1.
Also implemented shifting of the Coulomb potential for the group
cut-off scheme for non-excluded pairs only.
The manual now explains these details.

Also fixed the generic group kernel with an exact cutoff for
either Coulomb or VdW. I think this was supposed to not be supported
but neither grompp nor mdrun checked for this. The fix is trivial.

Change-Id: I48bff73587e43338162f90fa7d526e1909ce5ad1

Fix print_date_and_time

The old call signature and lack of documentation of this function made
it easy to misuse, and (e.g.) mdrun was reporting the same initial and
final time stamp. Fixed signature and documented.

Change-Id: Ib4d9a07bc6e1165affa3bfdf7cd3b54c615d4765

Fix dependencies on config.h for typedefs.h

The include in simple.h hid other dependencies and was only
used for gmx_inline and gmx_restrict. Defining those
in just one place simplifies it.

Also fix missing includes in public headers.
All public headers should be able to be compiled by themselfes
and shouldn't rely on other headers aready being included.

Also fixes that math module headers were installed to maths.

thread-mpi headers and commrec.h are not installed anymore.
The installed headers have no dependency on GMX_LIB_MPI and
GMX_THREAD_MPI anymore.

Change-Id: Ie1f26913739d3e42d653a51fa2c2d71e6ea3ef7a

Revert "Added new amber99sb-ildn-bsc0 force field."

This reverts commit 3f7bdff15b41f7e0ffbfd0bf86b9c677ae51db77.

This is a temporary measure, pending testing of this variant
that is comparable with that for other AMBER forcefields, and/or
resolution of questions raised on gmx-developers.

Also removed line from the manual that was generated when this force
field was in the repo.

Change-Id: I1c2269d2ed554b827fb183d54f1db789a5762a8f

Renaming in PME-routines

Fixes naming issues in PME-routines by renaming variables
and functions to be more general, and not using names specific
to electrostatics

Fixes #1459

Change-Id: I7abb624b9f3d31e247404007dbe5d9c4e475687e

Bugfix in LJ-PME without electrostatics

Simulations using cut-off for electrostatics, and
PME for LJ resulted in segfault since pme.c tried
to access non-existing grids.

Change-Id: Ic32c8cd2eb6afb3e49265d6c9a7da81418b18dca

Removed more stuff related to gmx.ff and gmx2.ff

Code and documentation no longer refers to these force fields. There
are corresponding changes to examples in the manual.

Updated list of force fields in the manual, plus minor correction
to description of CHARMM27 there and in the forcefield.doc.

It would be nice to leave gmx2?.ff/forcefield.itp intact with an error
message, but this means gmx pdb2gmx continues to show these force
fields. Instead, I made a README to note that these force fields have
gone.

Change-Id: I10480a0fad169493d13a743c43e8b2321128155d

Remove interactive help from core selection parser

The help is now handled in the loop that processes user input, and it
never gets into the parser itself.  This removes several special cases
from the parser and simplifies it.

Add a few string utilities and tests for them to make it easy to handle
the help commands.

Change-Id: Ida58bccd3b9ced4b7060485cc2f064b73796c996

Improve interactive selection input

- Print the static index groups always at the start of the prompt.
- Add a hint for that an empty line provides information about the
  available index groups.
- Print the whole prompt again after an empty input line.
- Print only the selections from the current interactive prompt in the
  output.  This makes the listing clearer for commands like
    gmx sasa -surface -output
  when providing the latter set of selections interactively.
- Print the number of selections still required after the above list.
- Work better (on Unix) if stdin is from a pipe (using the same
  mechanism that was there in 4.5 and 4.6).

Change-Id: I38152fb4fccd142aeee318739ab3782fa6a775d8

Cleaned up endif/else-s in CMakeLists.txt and *.cmake files

Removed FindBLAS/LAPACK.cmake files which come anyway with
CMake > 2.8.8. Included properly the licenses in some of the
BlueGene*.cmake files.

Change-Id: I69d68566865facf0af759526a579f32c5c444464

House-keeping for MTTK

Leapfrog + MTTK silently produced a .tpr that is probably broken;
certainly the documentation only states support for Velocity-verlet
integrators.

Added note about deprecation and planned removal of MTTK + constraints.

Refs #1292

Change-Id: Iec2cf0dd866242735ce04a954e585b2461f6e701

Merge "Merge branch release-4-6"

Merge branch release-4-6

Conflicts:
src/gromacs/pulling/pull_rotation.c
Adopted the bugfixes from release-4-6

Uncrustified incoming code

Change-Id: I7d90bb49ef83750b9607f677e97a94f2b5793eeb

Fix perturbation to Verlet scheme when not perturbing

A .top that has B-state parameters for non-bonded interactions
should run the same with free-energy = no as a .top that has no
B-state parameters. Renamed variable to better reflect its use,
and used it to correctly trigger the FEP code path.

Change-Id: I47656b29ef4ca2462af0fb587e0424d0eab3d5c4

Update manual analysis section

- Reference tools as 'gmx <tool>' instead of the old names.
- Expand the section on selections now that there actually are some
  tools that can use the full power.
- Rewrite sections that were referencing obsolete tools (g_bond and
  g_sgangle).
- Remove references to tools that aren't there.

Change-Id: Iee99f69e2d4effa585d98930ef1d78c1367e2ec0

Use a constant instead of sizeof(int) for SSE2 immediate constant.

Apparently gcc-4.8.2 was not able to evaluate sizeof(int) to the
constant 4 at compile time, but thought it was either a 64-bit
value or a non-const value. Since the register sizes are fixed
we simply hardcode the value 4 bytes instead.

Change-Id: If739b733727b032eb7e5b7fa0a75e7cd51989a0b

Avoid cross product with zero vector in rotational pulling.

Fixes #1431 (rotation/flex-t regression test failing on BG/Q)

In do_flex_lowlevel() we checked (by mistake!) for xj-xcn being
zero, although we need to check for yj0-ycn being zero, since
we use yj0-ycn in a cross product in the following lines of code.
I now also replaced the direct check (0 == norm(...)) by checking
what gmx_numzero(norm(...)) returns. The latter replacement
was also applied in the do_flex2_lowlevel() routine. Note that there
the check for small xj-xcn was and is actually correct.

Change-Id: I972b6d67a81e30f297db286cd2224f66753a20aa

The van der Waals radius must equal the Coulomb radius with Verlet

This patch fixes an issue that can occur when using g_tune_pme
with a Verlet pair-list .tpr input file. If the .tpr itself has
large cutoffs (e.g. 1.2 nm) and one asks g_tune_pme to scale
_down_ the Coulomb radius, the van der Waals radius is not
scaled down with the Coulomb radius (only upscaling worked).
One ends up with an unusable .tpr file because rVdW != rCoul.
This patch ensures that van der Waals and Coulomb radii
are always equal with Verlet pair-lists. Fixes #1460. Thanks
to Joao Rodrigues for reporting the issue!

Change-Id: I5ef30e71a35cd83838040057e16e52e09ea82e9a

Install enums.h

Related to #1452

Change-Id: I94db8d0f9661b7398e5924f4adefbfdcc82da483

Fix MKL linking

Setting the same for CMAKE_STATIC_LINKER_FLAGS broke because that
variable got applied to calls to ar for linking libgmock.a. Probably
not important because Intel probably encourages shared linking for
MKL.

Change-Id: I6816d8ba318c704a0efad7659f5afecaeb48b092

Avoid array-bounds warnings with gcc 4.8.2

The asserts fix RelWithAssert build and for Relase
build the warning is suppressed.

Change-Id: I0a7e6576014bd32648a1e2f08d88e1f2eef3335e

Removed .itp files from top level directory.

There is no need for these files, the water models
should be included from the corresponding .ff directory.
Since many users are likely to still use those files,
an error message is printed.

Also re-added ions.itp and treated it similarly.

Also update an integration test that was using a one of the top-level
water .itp files that is now removed.

Change-Id: I599e203329340587afa7fc5efced485d09b9ce1c

Simplify gmxTestLibXml2

Change-Id: Ia652633f529cc1cf6ba09e29159025ff762efad2

PME load balancing now works with LJ-PME

PME load balancing only works with LJ-PME when PME is also used
for electrostatics.

Change-Id: I7145b242fd5dd00690e6b789b9cc6f3e147d430b

Fix misc. Doxygen issues

- Fix some missing brief descriptions.
- Fix some function definitions so that Doxygen correctly pairs them
  with the corresponding declarations.
- Fix some documentation visibility issues.  Not all these are strictly
  required, in particular in the SIMD parts, where \cond is used to
  exclude full classes.  But the presence of \cond is
  difficult/expensive to test automatically.
- Also some misc.  other fixes.
- Document some of these issues better on the Doxygen page.
- Exclude parser.h from Doxygen, since it doesn't contain any meaningful
  comments, and confuses Doxygen at least with the XML output.

Also, explicitly list input file extensions for Doxygen.
This excludes some random Python files from the documentation.  At least
for now, we probably want to keep the documentation such that it applies
mainly to the C/C++ code, and not bother how to also document Python
code in a Doxygen-compatible way.
CUDA files are still left out, as they were previously, since they
currently produce several warnings.

Most of these were found by the current work-in-progress checker script.
After these, the checker mainly complains about a lot of missing brief
descriptions in thread-MPI headers and in swapcoords.c.

Change-Id: Ib8a492fde6327ce3d814a23d887b05323bd6d334

Don't disable MSVC warnings in config.h

All other warning flags are passed as compiler flags, thus this
was inconsistent. Also didn't work correctly if source file
doesn't include config.h.

Also makes ICC warning flags consistent between Linux and
Windows.

Change-Id: Ic52f6eb3c2dc5751b8aac1bba836ad5ff7d7c3cf

Fix CheckCXXCompilerFlag

ICC on Windows was always failing the CheckCXXCompilerFlag.
This was caused by the warning:
LINK : warning LNK4224: /INCREMENTAL:YES is no longer supported;  ignored
being matched by the compiler fail_regex for aix because of a missing
escale for the semicolon. A single escape isn't sufficient. Using double
escape ("\\\;") doesn't produce this problem, but I'm not sure whether it
always matches correctly thus instead a "." is used (works because it is
a regex). This fix doesn't address either the fact that the deprecated
linker flag is used nor does it address the problem that if incorrect
flags are somehow added (either by cmake or the user), this macro
fails without a diagnosible error message or fatal error (to notice it).

Change-Id: Ia45d6cd4fbf4f464bcb03e0d945085987a97aaad

Bump required CUDA version to 4.0

- Added warning with CUDA 4.0;
- Fixed some typos in comments;
- Fixed the comment referring to GPU code generation with CUDA v5.0+.

Fixes #1382

Change-Id: Iee182ace765022fc6a179531d42e965594cff104

Replace inline with gmx_inline

We always try to avoid defining std names in our headers.
Thus we shouldn't define inline for those C compilers (=MSCC)
which don't have it.

Change-Id: Ia2023435c916892cdaf2ce97905b55ec10cd1a08

Fix dependency on swapcoords enums

state.h needs these enums, but it seems wrong to introduce
a header dependency on the whole enums.h, when so few places need
these enums at all.

Change-Id: Ia8f080d5c0dc73b19a54f5ff77fbf7ef76f16139

Remove some unnecessary declarations

- topexcl.h was an outdated copy of gpp_nextnb.h, and wasn't used
  anywhere.  Rename the source file to gpp_nextnb.c to match the header.
- Remove one declaration from gpp_nextnb.h, which wasn't used outside
  the source file.  Make it static in the source file.
- Remove duplicate declarations of (init|done)_top().
- Remove unused create_visibility() from src/programs/view/.

Change-Id: I87982e5581aff8caca3dc753492e01ee457ab694

Remove KERNEL_WIDTH

Wasn't used anywhere

Change-Id: I579297d71225ba3c0360f8419856b0f5e72a4b5a

Removed the deprecated gmx.ff and gmx2.ff.

There are still a few places in the manual that need to be fixed.

Change-Id: I5c265470c1b73cb755918c5a000208da3d1e7b22

Fix formatting in install guide

Change-Id: I28833dee21220aacbf06ccc7d06f8f5e6a328352

Fix CMake syntax error

Change-Id: I8a39e50a741a23da3640fe7ccbcd83c844bf0d69

Use gmx_int64_t for gmx_off_t

We anyhow require 64-bit file support and 64-bit integer
support. Thus the defintion of gmx_off_t can be dont
trivially.

Partial fix for #1452

Change-Id: I6ae70454cfa8ace3ca84a047149b6f4648e84b00

Fixed exit with FEP with PME nodes without LJ-PME

With FEP runs with separate PME nodes without LJ-PME, mdrun would
exit with a gmx_incons.
Reorganized the PME some flags and conditionals for clarity.

Change-Id: I04513a2d84682617b36c0c8ad4afe632ce3f7ee1

Martini CG - Combined angle-torsion potentials

- Potential important for Martini
- Most Important diffence: bondfree.c
- Rest of files: definition of types
- Different MD groups want to have it
- Avoid the numerical instabilities
of having three beads in a row and
therefore an undefined dihedral
angle. This might not happen easily
at the atomistic level, but for coarse-grain
models this is a real pain.
- Paper at
http://md.chem.rug.nl/~monica/Bulacu_Marrink_2013.pdf

Adopted new tpx_version scheme.

Paper published.
Change-Id: I8ed470c464d843eb0218c87aeea3a44c3a1f64f8

Six new quotes from classic literature.

Change-Id: Id881baea8d6eeedc70dec3975a4c3e20e788285a

Add RelWithAssert build type

Works mostly like Release just without -DNDEBUG.

Also generalized some machinery for managing compiler/linker flags so
that we can have an easier time later when we reconsider how to manage
all of our build types.

Also minor corrective changes to the behaviour of other build types.

Change-Id: I70283cb6069cb772c287c0aef329d95b8eacd987

Remove dependency on TNG in trxio.h

Partial fix for #1452

Change-Id: I59474ea08d40362fa9d10745a5fb1d4e12254b68

Report MIC SIMD correctly

Avoid cmake warning because SIMD cannot be detected. And
avoid that version information prints "None". This is
needed because MIC doesn't use GMX_SIMD so that the same
setting (__MIC__) can be used both for native and (future)
offload.

Also added installguide section.

Change-Id: I4af3b13678aa5c3c69c60bda499c3691ac83dc53

Merge "Merge release-4-6 into master"

Buffer estimate now works with rvd-switch=0

With rvdw-switch=0, the Verlet buffer estimation code would divide
by 0. The estimate for vdw-switch functions now use the exact 3rd
derivative of the potential at the cut-off iso estimating the 2nd.
grompp now prints a note with rvdw-switch < 0.5*rvdw.
Updated estimate formula in the manual and added a missing factor
of 2 (the code did use the factor of 2).

Change-Id: Ia259d703515f69b0ed3390e2f7be671e0a15bd5f

Check zlib can actually be linked

Only use zlib compression in tng if zlib can really be linked.
This is equivalent to ee26e1264 for libxml2.

It is slightly simpler than the approach in gmxTestLibXml2
and CheckLibraryExists could be used there too.
This does not detect potential linking errors with shared
libraries similar to those solved in #740 for FFTW.

Fixes #1435

Change-Id: I5009e50b6caf810a12e71ecd83bde2edb21a2b32

Avoid segfault when comparing tpr

Segfault was possible if number of idef didn't match

Change-Id: I62315e3340eec8acf14c01d14917392bc63e46e0

Add enum and docs for using tpx_version

Documented procedure for having hassle-free feature branches and
ensuring merges into main branch cannot silently cause problems.
Added an enum to make the process work smoothly.

Change-Id: I69b87f9bac0b51784e70f8a618808ef086951343

Merge release-4-6 into master

Conflicts:
admin/installguide/installguide.tex

3643d2a Add fatal error for Andersen+constraints+DD ignored
because already in master.

Change-Id: I6b3346503036e5b31bcc929daa5bcc9a04dde3bc

Fix calling SIMD group kernels

The #define of GMX_SIMD*_OR_HIGHER takes place in simd.h, which got
changed late in the evolution of the recent SIMD-module patch, in
order to reduce unwanted module dependencies. However, the explicit
dependency of the group kernel setup code on the SIMD module was not
added. (It seems that the dependency on the old SIMD macros.h was been
satisfied transitively.) In effect, this turned off group-scheme SIMD,
but the silent fallback to C kernels meant there was no red flag until
we got around to measuring performance.

Change-Id: I78842cb191f6188fe4188800399efd2ed99e5543

Fix incorrect grid cell size in g_sas -nopbc

Fixes #1445

Change-Id: I798fc8fe96608633f26d9a3500f83f39b44af008

Move default library path search to program context

The default library path search is now implemented in
CommandLineProgramContext, and accessible through the generic interface.
This provides the following benefits:
 - All the magic for finding stuff from the binary path is in a single
   place, under the commandline/ module, where it belongs.
 - Code using the library path is now in C++, fixing some potential
   buffer overruns and removing the need for some utility code.
 - Programs using the Gromacs library can customize how the library data
   files are found.
 - PATH_SEPARATOR is no longer exposed outside path.cpp, encapsulating
   this detail.
 - The library path is only searched once during the execution of the
   binary, reducing the amount of file I/O done if multiple data files
   are accessed.
 - The library path is printed in the startup header, making it explicit
   where they are read (helpful in particular if there are problems with
   it, or the user may have multiple installations).  The GMXLIB
   override doesn't currently affect this printout; maybe it should, but
   I don't want to put the environment variable (and the support for
   multiple directories) into the program context.

Change-Id: I4d099215ef7e1f269d4c7a42989835d980bc1566

Add fatal error for Andersen+constraints+DD

This combination produced a temperature that was 6.5 degrees higher
than the same with one domain, for a 1ns PME lysozyme in water
simulation.

Change-Id: I9f80276c47de955a5053bcabb6fe7c9bfdceaf0e

switch ld_seed from int to gmx_int64_t

* allows more than 4,294,967,295 different seeds for ld
* tpx version bumped, used new style for magic version constants
* changed printf format strings to conform to new type

Change-Id: Id661ca331a100dc942ac978cc1b50cfdade09731

Added CUDA LJ-PME nbnxn kernels

This change implements CUDA non-bonded kernels for LJ-PME introduced
in the Verlet scheme with 99029d.

The CUDA kernels implement geometric as well as Lorentz-Berthelot (LB)
combinations rules (unlike the CPU SIMD) mostly because even though PME
is very slow with LB, it is still beneficial to let the user offload the
non-bondeds to a GPU and potentially bump up the cut-off to further
reduce the CPU PME load.

Note that as now we have 120 kernels compiled for up to four different
target architectures, the nbnxn_cuda module takes a very long time to
build and can become the bottleneck during compilation. We will deal
with this later.

Change-Id: I819b59a8948da0c8492eac6a43d4a7fb6dc98354

Bumped TNG to latest version.

Change-Id: Id3ca6c18c1f51640997829bd8ad251a7ca91d04f

Improved CUDA non-bonded kernel performance

Some old tweak which was supposed to improve performance had in fact
the opposite effect. Removing this tweak and with it eliminating
shared memory bank conflicts it caused improved performance by up
to 2.5% in the force-only CUDA kernel.

Change-Id: I7fcb24defed2c68627457522c39805afc83b3276

Residual SIMD Doxygen improvements

- Add a link to the overview documentation from the module page.
- Remove some trailing spaces.

Some issues noted, but not fixed here:
- \nosubgrouping would make the file pages look more consistent, if it
  worked for files...
- Some of the SIMD routines are not included in any named group, making
  them appear separately at the end of the documentation.
- The list of members on the SIMD module page is exceedingly long,
  making it hard to find things there.  Would require some thinking on
  how to best solve this.

Change-Id: I628f321172a3c3e8f4874e166c48e12763d1490e

Misc. Doxygen configuration changes

- Adjust Doxygen preprocessor behavior such that it actually parses
  #ifdefs in the code.  Otherwise, some parts of the SIMD code seemed to
  get ignored.
- Extract static functions on all levels of documentation for
  consistency (and we may end up with some static inline members in
  public headers; currently at least vec.h qualifies).
- Remove unnecessary excludes now that VERBATIM_HEADERS=NO already
  excludes the unwanted files.
- Some other minor cleanup.
- Move member groups higher up in the layouts, since they generally are
  the most interesting part.

Change-Id: Iea04772628716ed12f44cd4465c2027f4670bfae

Convert gmx-sas to the C++ framework

An example of converting a non-trivial analysis tool to the new
framework.  Most functionality is intact, but the command-line interface
is changed a bit.  Changes in behavior:
 - Writing a position restraint file is not implemented.  Would not be
   very difficult to add.
 - Solvation free energy estimates are written into a separate file
   instead of the main output file, and only if requested with an -odg
   option.
 - The calculation group is specified as a selection with -surface.
   Can be a dynamic selection.  The total area of this group is always
   written out.
 - There can be 0..* extra output groups, specified as selections with
   -output.  All can be dynamic, and must always be a subset of the
   calculation group.
 - There is no longer a separate concept of hydrophilic and hydrophobic
   surface area.  The user can trivially calculate such things by
   providing two different output groups.  The old hydrophobicity
   definition is easily expressed as "... and charge {-0.2 to 0.2}".
 - tpr file is no longer required for most of the output.
 - Output precision is different from the old tool.  The current
   precision (fixed to three decimal places) should be sufficient for
   most uses, but it could be considered to use %g for formatting like
   the old tool if really necessary.

Future TODOs:
 - Consider if the legends on the plots could be improved.
 - Add unit tests for the actual surface area calculation routine.

Part of #665.

Change-Id: Iee60c13b927b3b63b6e218164cd961971f2f3fce

Add support for non-uniform X axis in array data

It is now possible to construct AnalysisArrayData objects that have a
non-uniform X axis.  AnalysisDataAverageModule can use this to provide a
more meaningful X axis in the output plot in case the values correspond
to, e.g., a set of atoms.

Change-Id: I578f6944d4ff51c315d068b011b705be2a4be00a

Improve support for data sets in C++ tools

- If a data set does not set the number of columns, the framework now
  leaves it uninitialized (does not call startData() etc.), and the tool
  can check this when processing the frame.
- Make it possible to specify an output file in a tool unit test without
  actually testing the contents of the file.  This supports testing
  tools that only produce data in response to certain output options
  being set.
- Add support for specifying a comparison tolerance for a data set in a
  unit test.

Change-Id: I49794c93336c22c07f1a745d28d1b4f6f856b25b

CUDA kernels for switched LJ

This change adds CUDA GPU kernels for force and potential switching,
implemented in the Verlet scheme for CPUs in f3acaf.

Change-Id: I73a83d352924379d9005104fa59f3478d3e684b0

Improve floating point comparisons in unit tests

Improve basic floating point comparison support:
- Replace EXPECT_NEAR_TOL with a more flexible alternative that allows
  specifying various tolerances.
- Some Doxygen improvements in testasserts.h.

Additionally, make the reference data framework use the new approach:
- Add a tolerance setting to the reference data framework, and make the
  default tolerance a lot stricter.
- Store the reference data with enough precision to allow lossless
  roundtrips.  This makes them more difficult to inspect visually with
  the XSLT files, but that can be addressed separately.
- Run FFT tests with GMX_FFT_FLAG_CONSERVATIVE to make the results with
  FFTW have less variance.  The 3D FFT tests were already using this,
  but others weren't.  Set the tolerance for the FFT tests based on a
  bit of testing of what works.  It should now pass even without the
  conservative flag, at least most of the time.  The tolerance is still
  lower than what it was earlier even for single, and a lot stricter for
  double.
- Regenerated all reference values from a double-precision build.
  For the FFT tests, used fftpack without optimizations.

Change-Id: I187dd9be20fd17cc8b3e779dcd9dbad34043bb00

Remove Cygwin+MingW from installguide

MingW was never supported (#448). And Cygwin is currently
not working. Until/Unless it is working the guide shouldn't
claim support.

Change-Id: I6558b95bc47b80f7d97a7960587e83e959682990

Free-energy now works with the Verlet scheme

Implemented perturbed non-bonded interactions with the Verlet
by masking them from the nbnxn pair lists and storing them in an
old-fashioned group neighbor list. This way we keep a single
free energy kernel.
The free energy kernel uses OpenMP atomic for reductions.

Also fixed grompp setting nstcalcenergy=0 with nstenergy=0.

Change-Id: I5a6a7e84b46e06250d141e2e08cb3a11077cddab

Fix -Wmaybe-uninitialized warnings in some files

Also clean-up and clarification of clean_vsite_dihs.

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

Misc. Doxygen build system improvements

- Check Doxygen version when running one of the Doxygen targets, and
  produce a note that states the expected version if there is a
  mismatch.
- Check the result of the Doxygen process and fail the build if the
  program returns an error code.
- Add a tag to all warnings that identify the type of documentation that
  produces them.  This should make it easier to understand the Jenkins
  warnings in the case not all documentation types produce a certain
  warning.

Change-Id: I030c0b2c69958fe41f0b964399cb1c2e1b250e77

Added message to warn that group kernels are deprecated.

Nothing will change for now, but both grompp and mdrun now
warn that we will remove the support for group kernels in
a future version when all interaction forms are supported by
the verlet cutoff scheme (tentatively a year from now).

Change-Id: Id00a67e44e60818cb0e6f1897d6939f53ea01540

Updating and improving Andersen thermostat documentation

Added mdp options for andersen thermostat, and updated to
reflect changes in 5.0, plus add a citation for one of the assertions.

Change-Id: Ia713a4f2dd39dd5224952be334d8953b6613a66e

ICC aliasing optimization flags

GCC already uses -fstrict-aliasing (-ansi-alias). Thus it is safe
to also use this for optimization in ICC. const/non-const aliasing
shouldn't be used anywhere either and should be safe too.

Change-Id: I6a1b10c3d00bc2473a25a15e26934f61c046157d

Use reduce in pme

Use the reduce added to the simd module. Also fixes that
aligned store was used for possible unaligned stack buffer.

Change-Id: I5e08b0a0998fe5371148d1af9668fd6f2ecd5934

Redesigned SIMD module and unit tests.

Second part of the commit (after renaming), which contains
the actual functionality changes. This new version
implements a complete interface layer that
is architecture-agnostic, while each architecture has a
separate implementation file with both single, double and
simd4 versions (that can be used simultaneously). simd.h
contains a number of defines that describe the capabilities
of the instruction set, and there is a new documentation
module for Doxygen. This will be used in a later patch
for modularized verlet kernels. With that, we hope to remove
all architecture-specific SIMD code from the rest of Gromacs.
All SIMD math functions have been redesigned so they work even
with instruction sets that do not support integers (and even
on sets that do not support logical operations), and
accuracy has been improved for double precision sincos() by
removing the table implementation. To try to reduce
the size of this relatively large patch I have kept a few
header files (in particular the math files in gromacs/simd)
to avoid touching all the group kernels. With this new kernel
module, 256-bit AVX2 SIMD acceleration will now automatically be
enabled for the verlet kernels. Group kernels will use AVX_256
in this case. Also incorporates changes from Teemu to make
static and gmx_inline library functions appear correctly in
the Doxygen documentation.
Relocation of nbnxn SIMD setup from nb_verlet.h to nbnxn_simd.h
in mdlib by Berk. Now the nbnxn SIMD setup is completely internal.
Replaced the pr4 functions in the nbnxn kernels by simd4.
The nbnxn kernel selection is now nearly architecture agnostic.
Also enabled FMA again for pmecorr SIMD functions in double.

Change-Id: I643da75f346f120500682bcc4bcc1333a635db70

Replace all mdrun rngs with cycle based rng

The stateful random rumber generator (rng) used previously doesn't
produce reproducible results in parallel for sd/bd and doesn't
produce reproducible results for continutation for replica exchange.
The rng state has been removed from the checkpoint file.

Fixes #995

Change-Id: Id2a5d064cf363c54db3c16a0675cfeba553feeaa

Fix GMX_PRIu32 and GMX_PRIu64 for MSVC

Change-Id: Ib7490c66687d98a135338d51de41bfda35b8938d

Add ThreadSanitizer build type

This might be useful for helping detect misuse of thread
operations.

It does run quite slowly, so use in an automatic per-patch Jenkins
build is probably not a good idea. In any case, it should be run
on a real workload, as well as toy examples.

As discussed in the CMake file that manages the new build type, using
TSan requires that CMake can demand that atomic operations are not
used. Accordingly, the atomic.h functionality and the CMake atomics
testing that uses it have some changes so all the ways atomic.h is
used work correctly. Documentation of how and why things work is
expanded.

A benign race condition with the (unused) return value of
mdrunner_start_threads is fixed.

Change-Id: I0da01c0bedeec2e49d4456d757c54fbe58c56ef4

Fixes seperate PME nodes without LJ-PME

The very recent LJ-PME nbnxn patch dereferenced a NULL pointer
with seperate PME nodes without LJ-PME. Also LJ parameters were
communicated when not needed.

Change-Id: I40bd0c8f96294b54599afee6021139b71cc8c860

Keep compilers happy on BG/Q

Change-Id: I99439b002f95145504c359fddc61dc701431d8f0

Clean up options for domain decomposition

Documented mdrun options better. Documented the environment variables.
Renamed the environment variables to better describe their function.
Renamed the helper function because not all of these environment
variables relate to numbers of steps.

Change-Id: I12dba40f8cdf489ea23c3f4dd8810efde54a4aaf

Implemented LJ-PME nbnxn kernels

Implemented nbnxn SIMD and C kernels for LJ-PME. CUDA kernels will
follow later.
The nbnxn LJ-PME kernels subtract the full mesh contribution, using
the mesh C6 coefficients (usually geometrically combined). This gives
much higher accuracy for non-geometric force-fields.
The group kernels will soon be changed to use the same setup.

Change-Id: Ibbc6eef1c62a97ff6b2e3013f19594946581e44e

Fixed bug with LJ pot-switch in 4xN kernels

The recent commit of LJ switch functions in the nbnxn kernels
contained a bug in the 4xN kernels which caused the force on some
pairs not to be switched.

Change-Id: Ib4303223924e92d7e84641071a4a56955de70274

Permit warning-free use of Andersen thermostat

Simulations using Andersen thermostat don't need to remove center of
mass motion, because that is intrinsic to the algorithm. Using it with
nstcomm > 1 generated a warning that center-of-mass removal is
unnecessary. Using it with nstcomm == 0 generated a warning that
ice-cube artifacts might be generated. Resolved by suppressing the
latter warning for both Andersen thermostats.

Change-Id: I6b2c1594cabd81964b3ba1ebc5dd61e0f1debb5e

Suppress icc warning for Boost 1.55 use of #pragma once

Change-Id: Ib187cf34ba2c72a2c57fc55cfd26154508e9543b

Removed deprecated Encad FF.

Change-Id: Ib7f5c3d9920307a808dde9ef13b0f042827eb8ca

Fix load_lj_pair_params2 for MIC

Change-Id: I181bf5c6190d19ee2cb7bcf629c1b81e849e6373

Updated Boost to latest 1.55.0

Shouldn't have any disadvantage because boost still has support
for e.g. GCC 3.3. Newest version should work better with newer
compilers (e.g. no warnings).

Reduced the number of changes in boost to two:
- Virtual dtor fix is not required anymore
- Removed #error which was aimed at producing compiler errors
  for removed dependencies. Because this internal boost is only
  meant to be used by Gromacs we shouldn't have to detect
  odd usage of BOOST_* defines.

Fixed two comipler warnings:
- Applied nvcc.path from https://svn.boost.org/trac/boost/ticket/9392
- Commented out unsed parameter

These warnings should be fixed in future boost versions but if not and
the changes have to be reapplied, check Gerrit for the diff. Patch Set 1
contains the unmodified version without the fixes.

Change-Id: I09179891f9979a5997171a579a6a552621a693df

Merge "Merge release-4-6 into master"

Added new amber99sb-ildn-bsc0 force field.

It extends the original amber99sb-ildn with new
parameters for nucleic acids, while the protein
parameters are the same in both force fields.

Change-Id: I665ea598015428ced103ce18f6665afd0d5ede8c

portability aspects in install guide + minor tweaks

Added information on portability aspects related to CPU instruction
sets, related to #1428.

Additionally, made several minor updates and tweaks related to
compilers, platforms, cmake, etc.

Change-Id: I621262c939c119e5bdd5e7c91dda0ae3ffc60b7b

Merge release-4-6 into master

Changed BUILD_CPU_ACCELERATION to GMX_SIMD in
CMakeLists.txt.

Reverts change to src/mdlib/pme_simd4.h. Corresponding
fix for master is under review as I5e08b0a09.

Removed reference to PD in message in shellfc.c.

Conflicts (all trivial):
CMakeLists.txt
src/contrib/fftw/CMakeLists.txt
src/gromacs/gmxlib/gmx_detect_hardware.c
src/gromacs/legacyheaders/shellfc.h
        src/gromacs/gmxlib/copyrite.cpp

Change-Id: I2d43520ae8833055eff51d431e47cd8b90d3a687

Reinstate shell code with DD

Further work on the complex/sw test case in the 5.0 regressiontests
branch reveals that the initial conditions may have been the reason
for the problems observed with DD and more than one node, rather than
the implementation.

Refs #1429

Change-Id: I26ff6d9f8c79605afa794cae4761b5643b712124

Remove particle decomposition

The paths that are eliminated are those for which MD_PARTDECOMP was
needed, those triggered by PARTDECOMP(cr), and any for which PAR(cr)
&& !cr->dd.

Reviewers, please note
* that for the purposes of this patch, OpenMP is not a form of
  parallelism,
* the definition of DOMAINDECOMP(cr) in commrec.h,
* that TPI and NM can run in parallel, but use neither PD nor DD,
* multi-simulations run in parallel but need not use DD, so
* we still need both PAR(cr) and DOMAINDECOMP(cr), and
* I have generally left the indenting alone to make for easy review,
  but we will uncrustify before merging

Summary of changes in this patch:

Removed
* code triggered by MD_PARTDECOMP
* code triggered by PARTDECOMP(cr)
* anything decomposition-related with "pd" in the name (but the name
  clash with x86 vector intrinsics is unfortunate)
* t_mdatoms field called start (DD does atom number remapping)
Note that bPDvsitecomm was never set anyway!
* manual section

Renamed two functions with "pd" in their name changed to "serial"
because they are still needed there (and moved them?)

Deleted files, functions and function paramters that only supported PD
code paths.

Cleaned up
* use of DOMAINDECOMP(cr)
* a bit of Generalized Born code that should have been static or
  was unused
* bcast_state machinery can now be much simpler; still does the right
  thing, but does it further behind the scenes and easier to
  understand
* made explicit the assumption of replica exchange that the
  integrator is dynamical, which makes for some minor
  simplifications

Refs #1292

Change-Id: If029f16e6b4b06d58d465afe072a3cde6481479e

added gmx_is_{single,double}_precision

* allows easy detection of the precision for cmake, autotools
  without parsing the output of gmx_print_version_info (no
  cross-compile support) or the output of strings command (unix
  only)
* linking against libgmx with/out -DGMX_DOUBLE will lead to
  unpredictable segfaults

Change-Id: I472f10ae374a1f42c94c55e156b53f8905bdf098

Fix aligned store to unaligned memory

Also fixes that unaligned store was used when not necessary.

Change-Id: I44bb222a07ec0af65198667787b8673b3c6cd2e7

avoid mdrun crash when rdtscp is not supported

When using rdtscp, mdrun now detects at runtime whether the CPU supports
this instruction and if this is not the case, it issues a fatal error
and instructs the user to recompile mdrun for the compute host. Note
that this will happen rarely, only when cross-compiling from a newer
host for a rather old one.

Additionally, when the user manually picks AVX, we also turn on RDTSCP
as all AVX-capable CPUs support it.

Also made CMake advanced cache option for GMX_USE_RDTSCP. This replaces
the previously hidden GMX_DISTRIBUTABLE_BUILD option.

Fixes #1428

Change-Id: I8bc884ef9ea8ea4661626b60490182ae2b302648