Remove overeager selection index group check

The check for sorted groups in selections was a bit too eager, and
prevented use of many index groups that would previously have been
perfectly valid for, e.g., g_bond, with gmx distance.  Remove the
check.  Instead, fix the case where unsorted group is passed as the
value of a selection, and add tests for that.

Will reinstantiate the check later in another form, as currently
unsorted index groups get silently accepted and treated in a more or
less random manner in all other cases.  However, this change provides
a usability improvement that could be nice to have, even if I don't
immediately have time to implement the check.

Change-Id: I86e935d8bd08e63a9e3d01fb5a4084d96c5b7006

Merge release-4-6 into release-5-0

Conflicts:
CMakeLists.txt
share/top/oplsaa.ff/watermodels.dat
src/gromacs/fileio/trxio.c
src/gromacs/gmxana/anadih.c
src/gromacs/gmxana/gmx_chi.c
src/gromacs/gmxana/gmx_tcaf.c
src/gromacs/gmxlib/checkpoint.c
src/gromacs/legacyheaders/tables.h
src/gromacs/mdlib/shellfc.c
src/gromacs/mdlib/tables.c
src/gromacs/pulling/pull.c
src/kernel/readpull.c -> applied to
             src/gromacs/gmxpreprocess/readpull.c
src/programs/mdrun/repl_ex.c

Change ignored because was already merged:
       src/gromacs/mdlib/sim_util.c

Deleted. No merge:
src/gmxlib/gpu_utils/dummy.cpp

Removal of INSTALL_NAME_DIR not necessary because
was already not there anymore:
src/gmxlib/CMakeLists.txt
src/tools/CMakeLists.txt
src/kernel/CMakeLists.txt
src/mdlib/CMakeLists.txt

Tools deleted - no merge:
src/tools/gmx_bond.c
src/tools/gmx_kinetics.c

Change-Id: Ibeb5a23cce94fb1765ebee271c9c92e001b83630

Always set b-state posres, even if identical to a-state

Previously we did not read the b-state position restraint
file (option -rb) if the name was identical to the a-state,
which caused the position restraints to be zero. Rather
than trying to be smart, we now always read it.

Fixes #1399.

Change-Id: I13d93ab667734f022e48143e4cd0672b7f303e1c

Update g_density to handle bilayers better

As suggested by Chris Neale, we now properly
perform the binning relative to the center rather than
merely shifting the center of the output. While at it,
I have also added a selection for the group to center on
(useful when the membrane is not in the center of
the system, say with a membrane protein), and an option
that allows us to perform the binning in relative
coordinates and output the average dimensions. The last
is useful when there are large box fluctuations. The
tool now also has an expanded help section with a few
recommendations for bilayers.

Fixes #1168.

Change-Id: Id132311f1d6cd4858f99e8f44a0266667100a406

Fix LJPME Lorentz-Berthelot combination rules

Christian found a bug in the way we evaluate
arithmetic combination rules for LJ-PME reciprocal
space. With this change, both group and verlet kernels
produce virtually identical total energies for both
types of reciprocal-space combination rules. It does
not affect the (default) geometric LJPME case.

Change-Id: I4bb0de1df9588dce3ca94e39a370a0b994b25b53

Changed linking order to make Extrae work.

In addition removed redundant GMX_EXTRA_LIBRARIES linkage.
The list is already linked in libgromacs and there's no need
to include it again for mdrun and the tools.

Change-Id: I12e8d32f1f0fb2abceea2a5bd929fad563c1b609

Improve PDB I/O (keep occupancy, b-factor, element)

Retain the extra PDB info such as occupancy, b-factors,
and element name columns with trjconv when providing a
PDB structure file. If you use a TPR file, we still
write the (given atomnumbers being present).
However, we deliberately don't store the pdbinfo in the
TPR file since that is a specification for a particular
experiment, rather than something that's valid in the
simulation.

Fixes #917, #1307.

Change-Id: I80cb8f5a250ba094fe81f32c58b4eb0298164053

Remove obsolete pull-geometry option 'position' from documentation.

Change-Id: Ie7e3e962e1b814c92ebba5de499167b3310d713f

Fix Phi Platform file for MPI

The Platform file was only working without MPI. One could still
compile with MPI without using the toolchain. But now this method
described in the installguide works for MPI too.

Change-Id: I188e3e970a2d4e777ba98cb9e0ea83f6c1276c62

Avoid warnings about coupling when free_energy=no

At least during free energy testing it is irritating to
get a new warning about the molecule listed in couple-mol not
being decoupled when we disable free energy - that is why we
are disabling it. Changing the warning to a note, since it is
really a note to to help, not something that will result in a
bad simulation.

Change-Id: If56061be0ff784aeedfc31b92532927ab562392c

Added information about dihedral restraints.

Fixes #1149.

Change-Id: Ifc6ab09d91116f3059addded11f8b0ed4f200545

Remove install_manifest.txt during make clean.

install_manifest.txt is generated in the build tree during install.
If installation is done as root, subsequent 'make install' as
normal user will fail due to permissions issue to overwrite the file.

Change-Id: Ie29b05db3a6b4f7e07ecd88f02e61c92a99ff003

Fix Visual Studio Project build

Fix two issues with project builds (the default cmake Generator on
Windows, jenkins uses jom):
- build configuration linker flags have to be defined
- only c++ flags are used

Change-Id: I773f45d0937f6bad9a2008ef1cfcbfc6c58a13a0

Merge "Merge release-4-6 into release-5-0" into release-5-0

Merge release-4-6 into release-5-0

This merges commit dced970, which changes many free-energy,
modifier and table-generation code paths, and its fix 349d8056.
That patch 349d8056 contains fixes to potential-shift and
potential-switch, as well as the shift/switch interactions
in combination with free energy. Since 5.0 has undergone
changes in the same areas (both nbnxn free energy, LJ-PME
and force-switch), this commit is a likely place for
bugs to have been introduced, so we keep it as a separate
commit.

Uncrustified the result of the merge.

Conflicts:
src/gmxlib/nonbonded/nb_free_energy.c
Resolved in favour of whichever branch seemed most right; changes from
dced970, 99aa704d and 5f59569a8 were all relevant here. We have
introduced some new LJ-PME-related variables so that code path is
reasonably similar to the coulomb path. We have also fixed a small
bug where the LJPME self-energy (i==j for verlet kernels) was not
multiplied by 0.5.

src/gromacs/gmxlib/nonbonded/nonbonded.c
Resolved as for dced970

src/gromacs/gmxpreprocess/readir.c
Resolved as for dced970

src/gromacs/mdlib/forcerec.c
Resolved from both branches

src/gromacs/mdlib/sim_util.c
Resolved from both branches, and from 349d8056

src/gromacs/mdlib/tables.c
Resolved from both branches, and added a few lines of code
to make LJ-PME work with shift modifiers.

As noted above, to avoid breaking the 5.0 branch, we have
manually added the changes corresponding to 349d8056 to make
sure force-switch (same as vdwtype=shift) results in correct
dispersion correction, and we have added a fix for the sign
of the LJPME grid c6 term in the generic nonbonded kernels.
This means 349d8056 should not be merged in again to 5.0 later.

Change-Id: Ida29b143a1bcb727ff38f9c63bf133bf749477b1

Fix MacOSX rpath

d9d2b0aa74 added rpath for MacOSX but the install rpath
was wrong for MacOSX.

Change-Id: I33984696a0800baf0c58bc13ffbdddc514008750

Updated Amber GBSA parameters

Updated parameters according to Per Larsson's post on
gmx-users in November 2012 (076230.html) that found energies
to agree to at least within 0.2kJ/mol (OK, not great).

Fixes #1161.

Change-Id: I1666238c332fb9d8cc82ab6faea672e016d78757

Make source compile under Linux X32 ABI

The X32 ABI uses 32-bit pointers in combination with all
other features of x86-64 to save memory. This requires two
minor changes to x86 inline assembly macro checks, where we
need check for __LP64__ (64-bit pointers) rather than
__x86_64__ (which is set for X32). The resulting binaries
run approximately 2% faster, but since it requires all
dependent libraries to be available in X32 format it is
not something we will enable with a simple setting for now.
To use it manually, set CFLAGS and CXXFLAGS to "-O3 -mx32".

Fixes #1263.

Change-Id: I2555b2c68d4bba24f4e581268c472a47fe7d1531

Make it harder for compiler to optimize away SIMD CMake tests

On at least one old version of Linux, a new compiler in
combination with an old assembler led to the compiler
understanding the SIMD code but optimizing it away, which made
the test pass even though the assembler could not handle it.
This changes the return value of the CMake tests to be based on
the SIMD operations, which should make them a lot more
difficult to optimize away.

Fixes #1493.

Change-Id: I3e021c3c718cf54afaadf131c5fa911b3933f61e

Added check for domain decomposition and shells.

Fixes #1376 in as much as that it prevents incorrect calculations.
Can be left out when merging to 5.0 that does not support
particle decomposition. In 5.0 we can write a message about OpenMP
if that works with shells.
Change-Id: Ic07ae70cd3cd5298b6b2d579719da4d267b2d8fb

Manual fixes

* updates for cut-off schemes
* used more Roman fonts for abbreviations that are not algebraic
* fixed layout and bounding boxes of figures
* removed references to removed GROMACS-87 force fields
* changed formatting of column of [moleculetype] parameter-listing
  table because we now have some interactions with really long names
  (and some consequent changes to how footnotes and labels
  have to work)
* removed out-of-date appendix B section on Modifying GROMACS
* fixed aspects of CHARMM description, including adding the second
  reference for nucleic acids, and adding a link to the CHARMM36 port

Change-Id: Icb8c0a6087af282a653bb528d71f77e1bb04fae1

Improved check to decide if to add chain/residue.

Together with the latest TNG library updates this
fixes #1507

Change-Id: Iea7160cc1194fc398f9c4db0cc66f4bd43ce724a

Added check in grompp for shells and inputrec issues.

Some combinations of inputrec settings do not work with shells,
in particular nstcalcenergy > 1 or use of a twin range cutoff.
This is now checked for in grompp. Fixes #1376.

Change-Id: I4382bcf5231920b22c725a7cf8e7b4e17c2526d9

Fixed the sign of mass perturbed contribution

Corrected the sign of the kinetic energy contribution to dHdl and dH.

Fixes #1527.

Change-Id: Ia1a2f708cdbba9a5a120a2ba8510cc2498f4d054

Fixed GROMOS assignment of PRO-NH terminus torsions

This patch also adjusts the parameters for torsions
over N-CA bonds in the termini, according to
suggestions by Justin Lemkul & Drazen Petrov.

Fixes #1022.

Change-Id: Ib00d69eeddf924b1e7e2de33475d774d7da54899

Make sure duplicates bondeds are set from hackblocks

Previously the rtp file could end up containing two
entries when the N/C termini databases added interactions.
These were then quicksorted, which could make the
assignment random. We now check for duplicates and
discard the rtp entry in favor of the hackblock.

Fixes #1395.

Change-Id: I505cfbb7e6c8e8050e6e2935cc417e49a515d8d8

Clean up CMAP placement in parameter list

Some of the CMAP variables were always placed into and
read from the first entry of the parameter list. While this
did not result in any errors, this patch now places them
correctly in the F_CMAP position.

Fixes #1345.

Change-Id: Ic6e09e46c352976cfc1f57ff903c082b4ec43df8

Added multisimulation info to log/output

mdrun now reports that it is part of a multisimulation
setup, how many simulations there are, and that the
detailed MPI/OpenMP data applies to the individual
simulation.

Fixes #1204.

Change-Id: I04f1edf9a7fd1eda202b83713bdd877807242dca

Fixed recent bug with disp.corr. and vdwtype=shift

The very recent commit dced970a introduced a small error
in the dispersion correction with vdwtype=shift.

Change-Id: I23c082ed190a5eda096504a23b277fb879c5feb4

Extend version checking/warnings for checkpoint continuation

Detect if the major/minor versions differ, and write a more
explicit warning that both explains this type of continuation
is unsupported, and recommend to use the -noappend option.
In contrast, for cases where only the patch level differs,
we now only issue a note that the binary has changed, since
this is fully supported.

Fixes #1230.

Change-Id: Ifec509de0e5a05d170ada89fd000c961e311c9ca

Enable terminus-specific atom translation

By adding protein-nterm and protein-cterm classes for atom
name translation in xlateat.dat we now avoid replacing
names such as O2 in non-standard amino acids. This patch
also corrects a cosmetic issue in the number of residues
claimed to be found in each chain by pdb2gmx.

Fixes #567.

Change-Id: I5956991e3ff3e5e1202330a857fd773f9e9b6945

Added tip5p to the list of watermodels in some FF.

Added support for flexible tip5p to OPLSAA and parameters
were taken from the Amber/Charmm tip5p.itp files.

There are, however, various problems with that water model.
This commit doesn't resolve them so the bug report below
will need to be followed up at later stage.

watermodels.dat now includes a note about these issues
and refers users to the redmine report.

Refs #1348.

Change-Id: I3adaae21f35efbb4240716e91f1bb5d256c8853e

Enable 4-letter resname in PDB output, keeps more pdbinfo.

This still fully adheres to the PDB standard since column 21
is not used by the standard. All common programs (PyMol, VMD, etc)
understand the 4-letter format, and programs that only read three
letters will still read the same filename as they used to. In
particular, this conserves most residue names during pdb<->gro
format conversions. We have also killed the non-standard
wide pdb format to avoid writing broken PDB files.

Fixes #725. Refs #917.

Change-Id: I9b6b8f2e191acdfb65ca2b5d96f39249cd71ea98

Detect non-fatal warnings when trying CFlags for source

We now first check the flag itself (and catch warnings),
after which we check if our SIMD source compiles with it.

Fixes #1461.

Change-Id: Ic6ae99f62c5657195a59dc83ab08399d3491249a

Enable xvgr output and labels for g_principal

g_principal now writes proper xgr files instead of raw
data files, and the output has more explanatory axis
labels as well as legends for the data sets.

Fixes #609.

Change-Id: I1eb7914a25661734e6a95e70413194ef723f9635

Disallow extracting V/F from XTC files in g_traj

We now issue a fatal error if the user asks for velocities or
forces to be written to as output from g_traj when the input
XTC file does not contain these quantities.

Fixes #1298.

Change-Id: I45a12fb1e5db5b2b2e1b072444e5430373fbaaa3

Warn if using nsteps=0 without continuation=yes

To achieve true single-point energy evaluations for a structure
we should not constrain the initial conformation. This is
typically done by settings nsteps=0, so grompp now warns if this
is used while keeping continuation = no (the default).

Fixes #1231.

Change-Id: I4e65b66f6535da4206e7dbd5f7769a33c94ec12d

Implemented fatal error for nstcalcenergy!=1 with shells.

The combination of using shell particles with a value for
nstcalcenergy that is different from 1 leads to energies
not being communicated around on a parallel machine and hence
a meaningless shell relaxation procedure. This patch will just
bail out with a fatal error if the combination is detected.
Fixes #1376.

Change-Id: Ia54f1480b2358355438acb5de31dfa56b7f8d603

More clear message for undefined pullgroups.

When a pull group was not defined in the index file,
grompp didn't report in a clear way where the problem was.

Fixes #1446.

Change-Id: If571b46b91f6ad28213808b6e18d96e1da0be79d

Make log file option listing similar to mdp names

Modify names to be consistent with the mdp file spelling
and abbreviations, to the extent possible. This patch also
removes the unused ndelta option.

Fixes #1514.

Change-Id: I3e56ee3dc773d027b2c0c5ca060fa34e4930a4d6

Use ${CMAKE_INSTALL_LIBDIR} as lib directory for tng_io library

Directory lib was used as library directory for tng_io library. This
causes problems on Debian wheezy (and future versions), which expect
libraries to be installed under lib/<arch>. GROMACS already uses
${CMAKE_INSTALL_LIBDIR} for libgromacs.so, and this patch makes
libtng_io and libtng_compress install in the same location.

This was also addressed in the TNG repository and this commit
corresponds to commit e8fa7bf3abfc74acf60a048cdfe1f65fdbc0dc2d in
the TNG repository.

This patch also bumps minimum CMake version required for building
tng_io from 2.8 to 2.8.8 due to inclusion of GNUInstallDirs module.

Change-Id: I68b2196ab77e3e18c3a50717327e59b3482f2a29

Make sure masses are copied from termini databases

Masses should always be copied from termini databases,
where they are mandatory since they are specified before
the charge, rather than picking them from atomtypes.atp.
This patch also adjusts a few of the default CHARMM27
charges in atomtypes.atp that deviated in the 3rd decimal.

Fixes #1207.

Change-Id: I6f0c87989ae285ff435cd56c64432988c6634824

Fix g_mindist with pbc=XY

Fixes #1189.

Change-Id: I7c8242df6f2644bebc49c13c86412f8bc1d37428

Fixed typo in oplsaa ASPH charges.

Part of #1395.

Change-Id: I03c88c8377affadc629c7b9a04d5e8ddd1cb946a

Fix that seeking was wrong if frame time was modified

The trjconv -t0 option causes fr->time to be modified and
this then causes the skipping to be incorrect. Alternative
we could document that fr->time shouldn't be modified
and change that trjconv makes a copy of fr before changing
time. This change has the advantage it is less likely
to reintroduce the problem.

Fixes #1405, #1406

Change-Id: Ibd692dfecdf175450fc43f7f35b774bf4446b316

Removed unused parameters and cleanup of inputrec.

Removed optimize_fft, dihre_fc, nstcheckpoint.

Fixes #820, related to #552.

Change-Id: I8d1c3700834e5575972165410b3b29cbb56cba53

Fixed shells with particle decomposition

In the spirit of version 5.0, the recent Verlet scheme fix for shells,
commit debd482d, removed the support for PD. Now it's back.

Refs #1429.

Change-Id: I686d4287ec8946e418aa98e739a1a81a0b7f7055

Fix ref error in constraint pull code

With multiple pull constraints, the convergence check would use
the reference distance of the last constraint for all constraints.

Change-Id: Ic15d5dd618acc6d3aa5ca0f1e009c6f8096ed0d1

Small typos in surface tension description.

Fixes #816.

Change-Id: I9273f186ad17b5c373d4a3e87fcf0a981dba40e3

Bumped the TNG library to the latest version.

This is commit f7083045c6f7d23bac47b8302e413becd53fc030
in the TNG repository.

Change-Id: I1f5a9edbaf1588b6d1cee6cf86a7ea254aad23a1

Disable maybe-uninitialized warning

Causes false positives with old gmx_fatal. Don't merge
into master branch.

Change-Id: I48e9d7798bbf1990cfeb79efaeabc354ab469746

Add ASAN build configuration

Also renamed ThreadSanitizer to TSAN

Change-Id: I4373a626acc3962bfb8f163297ab146684bd5ea6

Fixed CMake 2.8.12+ CUDA dylib bugs and warnings on OS X

This patch fixes up a number of minor compilation issues
with CUDA on OS X. The __STRICT_ANSI__ define has been
added to make CUDA work with gcc version 4.8 and 4.9,
a dummy variable has been added to the dummy c++ gpu utils
file to avoid warnings about empty objects, and we now
properly handle both RPATH (newer CMake) and INSTALL_NAME_DIR
(older CMake) without warnings. CMake-2.8.12 introduced
a bug that caused Gromacs executables to be malformed since
the rpath to the CUDA libraries was added twice. This is
caused by fairly deep internal changes in CMake, but it has
been worked around by removing the rpath (which FindCUDA.cmake
only adds on OS X) manually for CMake>2.8.11.

Fixes #1471.

Change-Id: I8b01204dc07678bc305821162f20db65e0e7b88d

Added reference to cmap in pdb2gmx.

Fixes #886.

Change-Id: Ib4c84a50962ed90ce145ad001322a9a46c8a2ab0

Fixed sign error in posres dvdl

For position restraints the sign of the contribution of the change
in reference location was incorrect. Note that the contribution of
the change in force constant was and is correct.

Fixes #1408.

Change-Id: I75210c3855d13fdd58e0a698d6c5d7c7116f456b

Made mdrun warn when generic kernels are used

mdrun now issues a warning if the generic group-scheme, or plain-C
Verlet-scheme kernels actually did some work. Such kernels should not
be called during routine use of mdrun.

It would be possible to do the check for the use of the generic
group-scheme kernel earlier in mdrun, but it's not easy to ensure we
only write the output at the first such step. Given that the group
scheme is deprecated, warning only at the end is sufficiently nice
behaviour.

Change-Id: Ib9b15165e53510f537687a4b448aaa01f1418522

Disabled replica exchange when T not in order

There are issues when the replica property is not increasing
with the replica index. This patch has a partial fix. But since
there are still issues, unordered replicas now lead to a fatal error.

Fixes #1377.

Change-Id: I272d9831df85fdc18c1b2fe5724cbfc06e73444e

Fix a Replica Exchange output problem

when running NPT simulations a "b" was added at the beginning
of the Repl Ex line thus breaking demux.pl

Change-Id: I72781d40187d94216e844a8ea522c5d1832bba09

Updated mdrun -npme documentation

The number of nodes at which PME nodes are used has increased.
Added note on PME nodes not being selected automatically with GPUs.

Fixes #1374.

Change-Id: Ie1de87abd3d1204d99af8b4f8e6809e7806f5c08

Fix that skipping fails if target is 2nd frame

With bSeekForwardOnly the current file position should be an inclusive boundary,
but because the binary search stops if the search region (high-low) becomes
header_size, the current frame wasn't found. Because trjconv always reads the
first frame before calling seek, it failed if the target was the 2nd frame.

Fixes #1154

Change-Id: Id4151f0106abdb68d6067f02fe2927a4ec6d77a0

Removed rounding issue in nbnxn Ewald table

The nbnxn Ewald table could got the spacing dx passed as a real.
Changing this to double improves the accuracy of the single precision
table by two orders of magnitude, because rounding errors are avoided.

Change-Id: I11b6c27bd3f4e306c412da8f01761b03dc337ce7

Small changes in appendix A of the manual.

Final fixes #552.

Change-Id: I8e7939ae70eb14d2e1fc0a25b50774d472889bd1

Fixed PME bug with #OpenMP-threads a large prime

With hybrid MPI+OpenMP parallelization and the rank local FFT grid
size for a dimension not divisible by the number of OpenMP threads
in that dimension, a very small amount of the FFT grid overlap part
could set/added twice. This would only occur at low-medium MPI
parallelization with OpenMP thread counts with large prime factors,
which is practice means almost never. Even when it occured, no actual
differences in PME energies or forces were observed.
This issue was due to a leftover from when space was uniformly
divided over the grids iso assigning whole grid lines.

Fixes #1388.

Change-Id: I22904c7f55d2e96fc4b8cd1498af2087eaed47ac

Added tip5pe water model to oplsa.

The model introduces a small modification to the tip5p model
to reproduce better experimental data for simulations using
Ewald summation.

Reference publication J. Chem. Phys. 120 (2004) 6085-6093.

Fixes #1016.

Change-Id: I1afba614d65d1e4e8eb8594e99eea41e31f1e391

Make sure figure legends are adapted to xmgr/xmgrace

A few programs were not using the xvgr_legend call. g_mdmat
has been modified to use it, and for g_analyse and g_enemat
I at least fixed the existing direct-writing code by adding
a conditional.

Fixes #783.

Change-Id: I782af539b93737b570f266ba744e718fbd846e5f

Changed appending file name error message

The error message now mentiones potential names issues.

Fixes #1497.

Change-Id: I0c0195abd21c5012b49119024505b18a906078be

Fix memory-usage error in g_chi

Out of array bounds writing error as pointed out by valgrind
is fixed by this patch, preventing a SEGV.

Fixes #1503

Change-Id: I57959f4751f3eb645ff66bbcd4418607af56c1a4

Make nbnxn template and generated files match

Change-Id: I80b80e242d920862a64db22ae2321dcdb532faa1

Force g_bond to use structure file

Previous versions of Gromacs never wrote broken molecules, but
this is no longer true for domain decomposition. Box dimensions
are set from the trajectory, but the molecular connectivity
from the structure file.

Fixes #834.

Change-Id: I61ecb4957ef11e03bfeef10095a406cbddd23a34

Made shells work with the Verlet scheme

The issue was that atoms iso charge group should be put in the box.
Additionally water_pol only worked within a charge group, now fixed.
Note that with shells working across charge groups, which is always
the case with the Verlet scheme, no shell prediction is done.
We should implement prediction, since it improves performance.

Fixes #1429.

Change-Id: I2ebfc2d91fc161167f8f2573b61e1f519cf11fd8

Cleanup environment variables.

Renamed variables to more descriptive names and
removed unused ones.

Fixes #976.

Change-Id: I0c17443c93efe581bf728f783893181e5a7f8ed8

Added dummy mass for charmm TYR aromatic vsites

pdb2gmx was not able to create aromatic vsites
for charmm27 since this dummy mass center was missing.
Note that aromatic vsites are still largely untested,
and not recommended for general usage.

Fixes #587.

Change-Id: Ic7d2976536c6269d1dcf12b7d33bddab096c353f

Fixed perturbed wall interactions

With free-energy calculations, the B-state wall force and potential
were set to zero.
Fixes #1501

Change-Id: Ib53073c19fee19979ade96455b6505001187a2fa

Made gmxcheck -rmsd work again with forces

Change-Id: I776b282b17ee7c6cc4e56130dd804aa0ea04bce6

Make sure genrestr uses the disre_up2 parameter

Previously, a constant value of 1.0 was used instead of the
margin read from the command line.

Fixes #1357.

Change-Id: I65ef225159e7b6ef818c71d062b619e78fa7ea22

Removed undocumented -map and -mdata options from g_sham.

Both options and code were completely undocumented, and
nobody seemed interested in helping document them for the
last two years, so for now they are gone.

Fixes #812.

Change-Id: I4d60cf3d58da89d76a4461cdc2e75f031676d2c9

Always write output from genion

To facilitate scripting we now copy the input conformation
and write it to the output, even if no ions are to be added.

Refs #1458.

Change-Id: I3e467169e30bc7448758c2d506ca35f49b001b50

Decoupled repl_ex_nst from nstcalcenergy

The replica exchange frequency is automatically changed by mdrun
to a multiple of nstcalcenergy, which is annoying. It turns out that
it doesn't need to be a multiple, so this changing has been removed.
Fixes #1494

Change-Id: I00833f92fd468924f61879aff8b7c85fe79d3c2e

Mapped HIS1 to HSD in charmm27.

Fixes #1133.

Change-Id: I6c289af9b58c8e7cd66f36982bf118505aa38f7d

Made repl_ex time step more readable

A recent change made the repl_ex time print %e, now it's %.5f,
as in the log file.

Refs #1486.

Change-Id: I2e669d91a61e6a696ce6d123c701fa8b773a3c2f

Avoid writing xvgr formatting with -xvg none

Several tools were writing xvgr formatting code directly
to output files, even when the users selects -xvg none as
a command line option.

Fixes #1407, #1479.

Change-Id: I1db10c2ad7455332e2a415a922885cddaec5efd1

Print in higher precision repl_ex timestep.

Fixes #1486.

Change-Id: Ibe5bcee18f761bea96e1b4a432774a4472a40f47

Process negative sigma correctly with combrule 2 or 3

The initial combination rule code (for sigma/epsilon)
did not take negative sigma rules into account, which
caused segfaults instead of the values reaching the code
in convparm.c that handles sigma<0 during conversion
to c6/c12 (where it is used to signal c6=0).

Fixes #1391.

Change-Id: I437f06d67c5ecfb58d236590288ad122bcdf2df0

Replaced gmx_hostname_num by gmx_physicalnode_id_hash

gmx_hostname_num was only used to distinguish physical nodes.
Since it only worked for hostnames with numbers, we replaced it
by a hash of the hostname.

Fixes #1513

Change-Id: I8e60757707386f43269afe0bb38e8500decefcd6

Issue a warning for using gmx_rms -prev with large trajectories.

Refs #716.

Change-Id: I607d6a63259e0dd9de7c8c99e4917bbe9f37ff8e

Install libraries in correct directories on multiarch systems.

NB: this patch may change the default install dirs on some systems,
and will _not_ overwrite previous installations. Thus make sure to
link to the correct library versions of e.g. libgromacs.

Fixes #591 and #623.

Change-Id: Ibd1d7d216236a0f5841ba9f0ac622a3f8f1a8718

Set make_ndx max name length and count to 1024

Limits were 30, which was a bit on the low side.
Fixes #1512

Change-Id: I62d5146fdac06b7b5b7f56925addeca561c83173

Fixed the output format of g_principal

Gromacs-4.6 and earlier versions had the output transposed. This
patch fixes it such that paxisN.dat contains the x/y/z components
of the N:th prinicipal axis. The default file output names have
been changed to paxisN.dat to increase the probability that users
who rely on old scripts will need to read the help text and find
the changed format.

Fixes #1441, related to #609.

Change-Id: Ic1ed9370145d3389ae8f43c2a419765dabf3a66f

Fixed typo in description of conversion factors.

Change-Id: I5f9875ef94e8d5526f84c029f539fab642319a04

Modify single/double EPS to agree with ISO C

The previous Gromacs values were a factor 2 too tight. While
we're at it, we also added float modifiers to the single-precision
constants, increased the accuracy in the double max/min values,
and improved the least significant bit in one of the single-precision
values. The previous relative accuracies for the PME correction
analytical approximations where set fairly arbitrary by dividing a
tightly chose constant with GMX_REAL_EPS, and since this is now a
larger number the unit test tolerance for those corrections has been
increased slightly to avoid triggering false failures.

Fixes #1356.

Change-Id: If29d2ff4141d86fd019e0dbdbfb7737ead28c3d4

Added package version to tng_io CMakelists.txt

Without a package version (for tng_io), the CMake 3 beta
versions will issue warnings.

Fixes #1477.

Change-Id: I4a8f64f70e7372118eecfebb69b4a24cd18835cb

Merge release-4-6 into release-5-0

Conflicts:
src/gromacs/gmxana/gmx_energy.c
Eliminated time parameter (removed in release-5-0),
used gmx_ffclose (fixed in release-5-0).

src/gromacs/mdlib/update.c
Eliminated fplog parameter (removed in release-5-0), otherwise added
new MTS-fixing functionality from release-4-6.

src/programs/mdrun/md.c
Added new code path for MTS (from release-4-6), eliminated
wcycle and bInitStep parameters (from release-5-0).

Change-Id: I303bb649ed2373656aa4420d0afb27f738a9ea30

Fixed shift and switch modifiers, particularly for free-energy

When using tabulated interactions (historically with PME-Switch), the
previous free-energy kernels used tabulated interactions which gave
correct results. However, as we have moved to using the new
interaction modifiers, Ewald short-ranged interactions are computed
analytically. To extend the range over which we apply the soft-core
interaction, the free-energy kernels evaluated interactions by
subtracting the reciprocal-space component, and then applying the
free-energy evaluation to the Coulomb (1/r) short-range
interaction. This works fine for vanilla PME, but led to problems when
combined with a switch modifier, since we are switching a different
function compared to the non-free-energy kernels. This could lead to
large artefacts where the free energy was 100x off if we were applying
the cutoff to r while the switch was applied to the scaled soft-core
radius.

This patch modifies the free-energy kernel so that the vanilla, shift,
and exact-cutoff versions still use the compensation trick, while the
switch modifier always operates on the traditional short-range Ewald
functional form.

The (very small) Ewald shift has also been added when computing free
energy in combination with Ewald summation and potential-shift
modifiers. As the perturbation goes to zero, the interaction will also
approach the non-free-energy interactions. Tested to match the
non-free-energy kernel to with 1e-8 in the fully coupled state, it
conserves energy, and produces reasonable free energies for ethanol in
water.

This also modifies table-generation, table-usage, and
dispersion-correction code to use shift/switch forms (and correctly),
when that has been selected in the interaction modifiers. This
provides much more accurate results for our new shifted interactions.

Correct (unmodified) tables are now generated for 1-4 interactions
in a few corner cases in the presence of modifiers for non-bonded
interactions.

Code paths for using exact cutoffs now work correctly when
rcoulomb-switch != rvdw-switch, or if only one kind of switch is
active.

Free-energy calculations using a plain Coulomb interaction now
incorporate a potential shift if one exists.

The GMX_NB_GENERIC environment variable can now be used to specify the
use of the generic kernel even with shifts or switches active.

Fixes #1463.

Change-Id: Ia63a1ed7d6c9cdf9cd9e6209b6326a49043060ec

Fix constraint virial with multiple time stepping

With multiple time stepping the additional nstcalclr-1 force
contribution was constrained to remove it from the virial.
This procedure neglected the non-linear contribution due to
the rotation of constraints. Now the contribution of this force
component to the coordinate update is constrained instead and
the corresponding virial contribution is subtracted from the
constraint virial.

Fixes #1400

Change-Id: If3217f52808bf7491998324f8dc3161bc003ec1b

New quotes

Change-Id: Ia8494c51d7bce41594152644a4a25b9188d41da8

Updated C-/N-terminal partial charges in Amber03.ff.

At the time of porting the AmberFFs were validated against
AMBER 8 and the results have matched precisely. However,
that specific AMBER version had a bug due to which CT/NT charges in ff03
were in fact using ff94 charges. The bug correspondingly propagated
to the Gromacs ports. In newer versions of AMBER this has been fixed.

The current GROMACS patch uses charges as specified in the
all_aminoct03.lib and all_aminont03.lib files as
taken from the AmberTools14 distribution.

In that distribution (14) seem to be no updates to the ff9x parameters.

Fixes #1466.

Change-Id: Ie6cfea5702500ff6cd5019edb22f224d29135425

Use RNG correctly for v-rescale thermostat

Two integers were passed in the wrong order. I suspect from the
construction of the RNG that the only effect of this is to permit a
rare re-use of a random number in a different RNG stream (i.e. no
effect in practice).

Change-Id: I1fb8eddbe7c8b029dc3686be80f3f083108fc28c

Fixed g_energy Einstein viscosity

The g_energy -vis Einstein viscosity output was obviously incorrect.
This bug has been introduced in version 4.5.
Fixes #1516

Change-Id: I5ffc1a232f0c64769cc438c977b757e8d8b55b98

Added grompp checks for zero atom masses in state B

Grompp already issued fatal errors if atoms had zero mass, or
if vsites had non-zero mass, in state A. This patch adds
state B to the same check.

Refs #1433.

Change-Id: Iec3d764ff7db6e7bdc3901d51b475cf2c35fe137

Add quote on the kT-kj/mol conversion factor

Great advice to all those struggling with remembering numbers by Anders
Gabrielsson.

Change-Id: I65fe5247b9342f6f48de709ad5bdbc69a52abc9a