Roland Schulz [Wed, 26 Nov 2014 01:40:30 +0000 (20:40 -0500)]
Use bitmask for bonded
Also move bitmask to utility and f_thread_t to forcerec-threading.h.
This introduces a cycle because listed-forces now depends on
mdlib/forcerec-threading.h. There is not yet a module where
forcerec-threading.h can go.
Change-Id: Ic220be2880f7263aa0e50b13f07ed509eb86e1b3
David van der Spoel [Thu, 11 Dec 2014 19:54:05 +0000 (20:54 +0100)]
Four new quotes from around the world.
Change-Id: Ib2f597c1f0e55160365d9e074b02b7eacbb84298
Mark Abraham [Thu, 3 Oct 2013 19:05:00 +0000 (21:05 +0200)]
Remove heuristic group-scheme neighbour-list updates
-1 == nstlist used to trigger the use of an algorithm that did
neighbour searching only when particles had moved far enough that
there might be a need to update the list. This supported obtaining
better energy conservation with the group cut-off scheme. It has been
superseded by the Verlet cut-off scheme.
Part of #1292
Change-Id: I074790c1c5670e874b8b587a1d46e62e2edc8961
Mark Abraham [Fri, 23 May 2014 19:07:54 +0000 (21:07 +0200)]
Fix latent bug in OpenMP initialization
mdrun worked fine because it only tries to call
gmx_omp_nthreads_init() once, but anybody trying to run thread-MPI
mdrun (e.g.) as a library or in tests would see a hang, because the
master rank never gets to the barrier because of the premature return.
Extracted a function, so that the return works as intended, which will
help us reason better about this code. Used quick return statements
for control flow, rather than nested conditionals, but there's no
other change to the functionality.
mdrun with OpenMP can now be called more than once. As a side effect,
even without OpenMP, the subsequent calls to mdrun now issue the
over-subscription warning - and all OpenMP runs report the OpenMP
settings.
Change-Id: Ie7019565ef9a513cf4dfb530e90452d4ef48289e
Mark Abraham [Sun, 18 May 2014 10:48:46 +0000 (12:48 +0200)]
Add StringFormatter and formatAndJoin to stringutil
The StringFormatter functor wraps formatString so that we can write
one-liners that take a container, format its objects to strings, and
join them with a separator into one string for e.g. logging output.
Re-implemented joinString in terms of formatAndJoin, by using an
IdentitfyFormatter, thus eliminating duplication of implementation.
Added some tests for joinString also.
Applied new functionality to reporting GPU status information.
Change-Id: I424543a1ca4e214c438cc39f104f087e7e70030a
Teemu Murtola [Thu, 26 Jun 2014 11:59:07 +0000 (14:59 +0300)]
Convert gmx rdf to the C++ framework
Convert the RDF calculation tool to the new framework. Main changed
beyond mechanical conversion:
- Remove -com and -rdf options in favor of using selections directly.
Both functionality, and a lot more, is now directly supported with
proper selections.
- Remove -hq and -fade options, as both are just post-processing that
can be trivially done based on the output xvg file alone.
- Normalization is reworked to work with dynamic selections.
- Add an -rmax parameter that allows specifying an upper cutoff for the
RDF. Combined with the use of the neighborhood search routines, this
potentially makes the tool a lot faster if the RDF is not required up
to half the box length.
Add basic unit/regression tests for a simple water box. The reference
values are not verified in any way, so they mainly catch regression
bugs.
Change-Id: Ieed19fc39a84e299291232190333652ced26abee
Teemu Murtola [Sat, 6 Dec 2014 05:14:15 +0000 (07:14 +0200)]
Fix misc. analysis neighborhood search issues
- Avoid an overflow fp-exception in case of zero cutoff (original fix by
Roland).
- If a bounding box of the reference positions has a zero size in
some dimension, fall back to a single grid cell and simple searching
in that dimension (this could possibly be improved after the search
uses a rectangular grid).
- If the bounding box (or the PBC box) is very flat or has a forced
single grid cell in some dimension, adapt the grid dimensions such
that the cells really have ~10 particles on average.
Change-Id: I4f8e5df9c4c479e480052e575354c54319bac22f
Teemu Murtola [Wed, 5 Nov 2014 19:32:46 +0000 (21:32 +0200)]
Track locations for selection parsing products
SelectionTreeElement and all intermediate structures used for selection
method parameter parsing now track the location in the input string.
By itself, this change does nothing, but it enables the subsequent
change that switches to exceptions to use this information for much
better context information in the error messages.
The location is kept as a SelectionLocation instead of as a string to
avoid spending a lot of memory for these strings with complex
selections: in general, each node in the tree contains the concatenation
of the strings from its child nodes, and there can be quite a few nodes
at the top that all contain have the full selection as their location.
Related to #655.
Change-Id: I7b1a317e043c6ff791b25f5fa357d8166c577fae
Mark Abraham [Mon, 8 Dec 2014 23:14:56 +0000 (00:14 +0100)]
Merge "Merge release-5-0 into master"
Teemu Murtola [Fri, 24 Oct 2014 17:45:21 +0000 (20:45 +0300)]
Rectangular grid for analysis nbsearch
Previously, the grid-based analysis neighborhood search used triclinic
grid cells covering the triclinic unit cell. This made PBC handling
easy, but all other math is more complicated, and the volume where pairs
are searched is sub-optimal for triclinic boxes. Now, the grid is
composed of rectangular cells that cover the rectangular unit cell.
The most significant change to support this is that most computation for
what cells to loop over is now done in fractional cell coordinates to
support for non-integral shift of the cell grid at the edges of
triclinic boxes. The fixed list of cell neighbors is replaced with
on-the-fly computation, since the non-integral shifts mean that the list
of neighbors is not always the same. This approach also allows the list
of cells in the loop to easily depend on the test position. Any cells
that are completely outside the cutoff from the current test position
are skipped completely, improving the efficiency also for rectangular
boxes.
Currently, the list of cells looped over is still a rectangular
box, but with rectangular cells, the math to exclude corners should be
straightforward to add in subsequent changes.
Change-Id: I8bd4fa67e26e282cd5bd82dee5edfaf651e59a19
Mark Abraham [Mon, 31 Mar 2014 13:54:12 +0000 (15:54 +0200)]
Add basic GPU support for gmx tune-pme
New functionality
* permits the user to specify a per-node list of eligible GPU ids (as
distinct from mdrun, which permits the user to specify the array of
actual GPU ids to use), which it will then map to the PP ranks
in round-robin style during the optimization
* detects whether the mdrun is GPU-aware if the user specifies GPU ids
* handles error from mdrun when number of PP ranks does not match the
GPU id string provided
* handles error from mdrun when GPU with given id is invalid
* runs check that a .tpr will work mdrun with -npme 0, which is less
likely to run into issues at small scale with eligible partitioning
issues than -npme -1; noted TODO to fix this for large scale also
This code does not do anything active about optimizing the use of
OpenMP, as current-generation tune-pme is basically designed to have
no ability to change the number of MPI ranks at all. However, if the
MPI system is set up such that mdrun's automatic allocation of OpenMP
threads to ranks will work, then the optimization can take place over
runs that actually trigger the use of OpenMP. (This was true before
this patch.)
Some TODOs noted inline for future clean up / feature addition.
Refs #1119
Change-Id: Ic1ea4e7524f62e7e7480b2e61a0c0afad41ca710
Roland Schulz [Fri, 5 Dec 2014 23:31:17 +0000 (18:31 -0500)]
Merge release-5-0 into master
Conflicts (trivial):
src/gromacs/gmxlib/gmx_thread_affinity.c
src/gromacs/mdlib/clincs.cpp
Manual changes:
Removed config.h because not necessary anymore:
src/gromacs/mdlib/nbnxn_search_simd_4xn.h
src/gromacs/mdlib/nbnxn_search.c
Removed config.h supression:
docs/doxygen/suppressions.txt
uncrustify:
src/gromacs/mdlib/clincs.cpp
Change-Id: I00064120f12f609fabecac9b3c823d33c015e59c
Teemu Murtola [Mon, 27 Oct 2014 18:17:32 +0000 (20:17 +0200)]
Simplify analysisdata parallelism model
Instead of trying to hide all parallelism within the AnalysisData and
AnalysisDataStorage classes, add an AnalysisData::finishFrameSerial()
call and make higher-level code responsible of calling that in suitable
context to do serial processing for data frames. This then needs to be
propagated to parallel-aware data modules.
This substantially simplifies the implementation and makes it easier to
argue about correctness, since now nearly all potentially parallel code
will operate on distinct data, reducing the need to think about
synchronization. The drawback is that now the calling code needs to
think about this, but there will be very little code that this actually
affects, and that is in the framework, not in every analysis tool.
Further, this actually makes the parallelism easy to express in terms
of, e.g., tbb::parallel_pipeline(), so the clarity can actually improve
without too much complexity.
Update and improve related high-level documentation to make these
aspects easier to explain.
Change-Id: Ibbfb4b7aa3acf7fcda0bfdb4da7212f5b11c1fcb
Jiri Kraus [Wed, 22 Oct 2014 15:08:51 +0000 (08:08 -0700)]
Add optional NVIDIA Management Library (NVML) Integration
NVML integration allows to automatically control GPUBoost from GROMACS
directly on supported GPUs. With this GROMACS either changes
application clocks automatically to the best performing setting or
informs the user that he should do it. The latter is done if permissions
do not allow to change application clocks.
Change-Id: I998e45596e2f556a5be7fa146fcd216ecae11452
Mark Abraham [Sun, 12 Oct 2014 16:35:16 +0000 (18:35 +0200)]
Clean up non-PME part of ewald module
ewald/ewald and ewald/ewald-util code had several different kinds of
stuff in it. Separated code that does Ewald long-range and
Ewald-family charge correction from the code specific to the group
scheme.
Moved general-purpose routines that are used in a few other places to
calculate Ewald splitting parameters to the math module.
Minimized header dependencies.
Removed unused things: FLBS FLBSZ
Removed use of typedef for existing opaque struct ewald_tab, per
policy in #1490. Renamed to gmx_ewald_tab_t.
Change-Id: I1394bbd02aa92e6581d011e52c5bee12406a0144
Mark Abraham [Sat, 17 May 2014 13:35:56 +0000 (15:35 +0200)]
Make mdrun print the list of compatible GPUs
This extends functionality so that the available, compatible, and
actually-used GPU IDs are reported when that makes sense. This will
help (e.g.) gmx tune-pme optimize resource usage by running mdrun
to probe for the available hardware.
The old implementation and documentation had a few issues that have
been corrected or updated. The major issue was that the list of
compatible device IDs was formed and then replaced by the list of
device IDs to use, so the former was no longer available during
subsequent output.
Change-Id: I3ca3e88054bd7d288f556af7e5d7704145ab6fca
Mark Abraham [Sun, 14 Sep 2014 16:30:55 +0000 (18:30 +0200)]
Move threading support into listed-forces module
Added Doxygen, converted a macro to a const int.
Change-Id: I1e5579b1f80aecaf46827704cf5fd3b6d6285d80
Mark Abraham [Mon, 15 Sep 2014 07:07:18 +0000 (09:07 +0200)]
Created do_force_listed and moved code to the module
This commit is almost pure code motion and renaming, modulo things
like changes to predicates because they've been handled earlier in the
call tree.
Converted the logic to return immediately if GMX_FORCE_LISTED is
false, to make understanding and indenting easier.
Also migrated position-restraints code to bonded module, and called it
from a sensible place together with the other code inside
GMX_FORCE_LISTED. Elapsed time for position restraints now accrues to
the listed sub-timer, instead of a catch-all.
Reorganized the foreign-lambda part of position restraints to be
called together with the other foreign-lambda bonded
interactions. Fixed benign bug introduced by the copy-pasting in the
Verlet-scheme monster commit, where dvdl was re-used for
foreign-lambda position restraints in the then-newly-created
posres_wrapper(), instead of the dummy dvdl_dum.
Made some more const correctness of x and pbc.
Cleaned up some include-dependency logic in files in the module.
Change-Id: I118af1d66da210e08d9ea065d1f431b7470742cc
Berk Hess [Mon, 1 Dec 2014 21:40:40 +0000 (22:40 +0100)]
Make domdec a proper module
Now has
* Doxygen group definition
* initial Doxygen comments
* correct #include guards
* no export "C" for header used only in C++ code
* all relevant current files are declared as such in the module
* new domdec_constraints.h and domdec_internal.h
* split off domdec_specatomcomm.cpp and domdec_vsite.c from
domdec_constraints.cpp
Change-Id: I7ed79c1654850df31e242aa1e907afd124e3ae1b
Berk Hess [Tue, 2 Dec 2014 08:08:45 +0000 (09:08 +0100)]
Clarified grompp pull error messages
Change-Id: If0df3fcd0d034c4b4ced04eb32ed7bb74e596776
Teemu Murtola [Sun, 2 Nov 2014 13:45:21 +0000 (15:45 +0200)]
Improve selection parsing error reporting
Now the selection parser tracks the part of selection text that
contributes to the current parsing context, and uses that for more
useful error messages. At least in some cases, there is now also
context information for syntax errors, although this could likely be
improved further.
Exception handling during selection parsing is also improved:
- The above context is now added to the exception, similar to what the
old reporting mechanism does.
- The full selection is added to the exception as context, also for
non-interactive parsing.
- More consistent exception handling, e.g., in cases where an exception
is thrown during error handling.
Improve documentation for related functions.
Part of #655.
Change-Id: Ib276ce4e219692c9819e45a1637e3660d958803f
Berk Hess [Wed, 26 Nov 2014 20:17:12 +0000 (21:17 +0100)]
Avoid race on dvdl with Verlet+OpenMP+LINCS+FE+VV
Also restructured the dH/dlambda reduction in do_lincs (used for
coordinates and not affected by the race issue) to work similar
do the do_lincsp code and properly use thread parallelization.
Fixes #1647.
Change-Id: I4eeb131018abca88b3635932491d99a779e16037
Teemu Murtola [Tue, 26 Aug 2014 04:19:40 +0000 (07:19 +0300)]
Non-periodic nbsearch with grid based on bounding box
For cases where there is no full PBC, the analysis neighborhood search
routines now allow constructing the grid based on the bounding box of
the reference positions. This frees the caller from the need to provide
a reasonably sized box even without PBC to benefit from the grid
searching.
Change-Id: Ia14b7c584a6ca84a5f698c12dc1491dde78bfe6d
Teemu Murtola [Sun, 23 Nov 2014 11:27:30 +0000 (13:27 +0200)]
Move/add uncrustify documentation to dev manual
Add some basic description of code indentation/formatting guidelines,
and move documentation of uncrustify and related scripts to the
developer manual. Improve the documentation at the same time, and make
all script reference the developer manual for additional documentation.
Make uncrustify.sh work without explicitly setting the UNCRUSTIFY
environment variable if the pre-commit hook is already configured.
Change-Id: I81ae71fbcdae0e177cc7b05fe278f8c20a19fe7a
Teemu Murtola [Tue, 23 Sep 2014 14:24:36 +0000 (17:24 +0300)]
Add 'gmx pairdist' tool as partial g_mindist replacement
Add a new tool that can compute pairwise distances from one reference
selections to one or more calculation selections. The tool can do
minimum or maximum distance over the whole selection, as well as between
each residue/molecule pair in the selections. This should provide
(nearly) all the distance outputs that g_mindist has, and with a bit of
effort also be able to replace g_mdmat.
The tool intentionally does not compute contact counts; those should be
added into a separate tool to keep things simpler.
Update some documentation references.
Change-Id: I553c3477ae808a133a3cec2b46feaccc3420022f
Teemu Murtola [Sun, 24 Aug 2014 04:16:41 +0000 (07:16 +0300)]
Analysis grid nbsearch for more cases
Add support for grid neighborhood searching also for epbcXY and
epbcNONE, as well as the XY mode in the analysis neighborhood search.
The grid is done based on the input box size also for the non-periodic
dimensions, or a single grid cell is used if the box size is unknown in
that dimension. In the future, an option could be added to determine
the grid automatically from the bounding box, but that is outside the
scope of this change. For the XY mode, a single grid cell is always
used in Z.
Improved the logic of determining when grid search is possible and/or
reasonable and implemented the option to force the grid
searching (since the test data currently would not trigger it
otherwise).
Change-Id: I0c43bf95d7d5dd46e57923cdf0d6864fbb3dfe59
Berk Hess [Thu, 13 Nov 2014 11:30:33 +0000 (12:30 +0100)]
Added warning for unnecessary soft-core
The sc-coul mdp option is, presently, only active with lambda states.
grompp now issues a warning when using soft-core without Van der Waals
decoupling without lambda states.
Also fixed an incorrect twin-range grompp check with PME.
Refs #1421.
Change-Id: I0605fe0f735d69f96f478612a00434eccef6232f
Teemu Murtola [Wed, 29 Oct 2014 18:44:07 +0000 (20:44 +0200)]
Basic support for 'z of ...' selections
Add basic support for selections of type 'z of ...', where "z" can in
principle be any keyword and "..." any expression. However, currently
the only syntax that is actually supported is '[xyz] of co[gm] of ...',
because
1) this is the original case that was discussed in gmx-users some time
ago,
2) only the x, y, and z keywords satisfy the constraints that the
simplest possible implementation requires,
3) there is no context that would accept multiple values as produced by
an expression like 'z of resnr 1 to 3', and
4) there is no mechanism to detect that 'z of atomnr 1' actually
evaluates to a single value.
Some of these limitations can be lifted in the future with reasonable
effort, though.
Change-Id: I6b87ce8d5c93e1ad05722dd59a66375b52b8e363
Szilárd Páll [Sat, 8 Nov 2014 02:05:31 +0000 (03:05 +0100)]
Add kernel and compiler suport for CC 3.7 devices
On compute capability 3.7 NVIDIA GPUs we can make use of the increased
register size by running 128 threads/block with keeping the minimum
number of blocks per multiprocessor at 16.
Change-Id: I84ec179a409668fe44fb9183cf3485c21bd53254
Christian Wennberg [Mon, 24 Nov 2014 13:55:47 +0000 (14:55 +0100)]
Fixed typo in amber99sb-ildn force field
Torsion calculations for N- and C-terminus ILE tries to find a non-existing atom type
Change-Id: Ibdc6ecc88e0c4da304b97f092d0a908451e12ac4
Teemu Murtola [Sat, 1 Nov 2014 10:57:00 +0000 (12:57 +0200)]
Clarify memory management for selection values
Add asserts for cases where memory would be leaked when assigning the
storage pointer for gmx_ana_selvalue_t. Add functions for cases where
the assignment is combined with some other operation that takes care of
the memory deallocation or ownership transfer to make the asserts not
give false positives.
The above clarifies the memory management and makes it less prone to
errors, but does not fix the memory leak that was revealed by "z of ..."
implementation. To actually fix that, additionally split the logic for
freeing memory for gmx_ana_selparam_t, and call that also from the place
that earlier set the value pointer to NULL.
Move Doxygen documentation for some of the affected methods and their
neighbors to the header per newer guideline.
Add some extra debug output.
Change-Id: I4e5dfa1248a20cab8a242be7209b5b0779204e64
Teemu Murtola [Tue, 2 Sep 2014 17:20:21 +0000 (20:20 +0300)]
Helper function for grouping selection positions
Add a helper function that can be used to divide selection positions
into groups based on residues and/or molecules.
Currently, 'gmx sasa' uses a construct like this, and 'gmx rdf -surf'
will require it as well, as will a replacement for 'gmx mindist', so it
makes sense to have a common helper function. The helper function can
then do extensive checking on the input just in one place, and it keeps
the tool code more readable by removing some complicated loops from
there. Also, a single set of tests is necessary to cover all the paths,
and the tool tests can focus on the tool-specific functionality, without
needing to test also the different groupings.
Improve related documentation.
Change-Id: Iaa6904c4b4b7ff5b824f37d264c327bdf52da273
David van der Spoel [Fri, 28 Nov 2014 06:40:42 +0000 (08:40 +0200)]
Fixes recently introduced bug in gmx analyze.
Since the introduction of a separate module for autocorrelation
and related algorithms the program gmx analyze did not work anymore
due to a typo introduced, such that the acf related command line
arguments were not properly added to the others.
Change-Id: If79adc761de397e42c387d1aa6bd53d565b02820
Mark Abraham [Thu, 4 Sep 2014 21:39:15 +0000 (23:39 +0200)]
Convert remaining callers of do_force to C++
Major callers are converted in other patches.
Removed unused variables, used std::min/max, named constants, put
maximum length on sscanf floats.
Introduced an instance of bMaster. Initialized gstat to and later
checked for NULL before use. These help static analysis reason
correctly. Previous code functioned correctly in practice, but correct
analysis would require cross-function analysis that cr, when passed as
non-const, is not actually modified. Fixed a warning about stepsize
not being used in L-BFGS, and added more comments to this algorithm
while editing the code to help us understand what it does.
The predict_shells() function has been modified to use the
gmx_mtop_atomnr_to_atom() call for case 2; the previous version would
use an uninitialized pointer.
Fixes #1593.
Change-Id: I54f238cedc78aadb0dad080ea3b28f001dce8d94
Mark Abraham [Wed, 12 Nov 2014 02:41:15 +0000 (02:41 +0000)]
Remove use of interaction_mask_indices on BG/Q
This field was degrading cache performance ~1% on x86. It probably
made little difference on BG/Q, because the extra integer operations
can use the second instruction-issue port, assuming the use of OpenMP
to use more than one hardware thread per core. Overall, this code is
about 1% faster on BG/Q.
Minor fix to the gmx_load_simd_4xn_interactions() function that looks
up the exclusion masks, so that new non-x86 platforms won't silently
fail for want of an implementation of this function.
Minor simplication to always pass simd_interaction_indices to
gmx_load_simd_4xn_interactions(), since it is only used on BG/Q and
then it is non-null.
Change-Id: I140a11607810e9cf08b702cae0b48426c3592fec
Mark Abraham [Thu, 27 Nov 2014 17:04:46 +0000 (18:04 +0100)]
Merge branch release-4-6 into release-5-0
Conflicts:
src/gromacs/gmxpreprocess/readir.c
(added content at same point in both branches)
Change-Id: I9df7cef8c2e4920461fb3d84806c58cf07c645e2
Teemu Murtola [Fri, 10 Oct 2014 04:06:47 +0000 (07:06 +0300)]
Clean up gmx solvate
- Remove unused code.
- Remove help text that did not have corresponding implementation.
- Remove paths in addconf.c that were no longer triggered now that
gmx insert-molecules doesn't use it (will remove the whole file next,
though, so careful review may not be warranted).
- Fix most memory leaks outside add_conf (add_conf still leaks several
MB during the unit tests...).
- Consolidated memory management and reduced the scope of some
variables.
Change-Id: Ibf570247a7e91bc028401165000142b03baebe4d
Teemu Murtola [Thu, 9 Oct 2014 18:59:57 +0000 (21:59 +0300)]
Port gmx-solvate tests to use CommandLineTestBase
This makes the tests relatively easily extendable to support checking
the output.
Change-Id: Icf835aa143d3832a340d156327c0aee6c80ef4ea
Mark Abraham [Thu, 27 Nov 2014 16:55:56 +0000 (17:55 +0100)]
Merge "Merge release-5-0 into master"
Teemu Murtola [Fri, 10 Oct 2014 04:03:54 +0000 (07:03 +0300)]
More control over file backups
Add a gmx_set_max_backup_count() to specify how many backups should be
made, instead of just relying on the environment variable. Make the
default to be no backups: now tests and FAHCORE do not need special
treatment. Instead, make the command-line initialization specify the
maximum backup count, and add a command-line option that also allows
disabling backups as an alternative to the environment variable.
This also allows removing a valgrind suppression for memory leaks from
the test code.
Change-Id: Ib512eb2a491691c76404aacbf93b67087a1afe41
Teemu Murtola [Tue, 30 Sep 2014 03:55:37 +0000 (06:55 +0300)]
Further insert-molecules cleanup
- Run insert-molecules tests under valgrind, and fix (most) memory
leaks.
- Clarify memory management by moving memory allocation of some arrays
into functions where they are actually used.
- More sensible status output.
- Additional tests.
Change-Id: I7bc2f7b4b155ac650235da00b414561b963d347c
Christoph Junghans [Wed, 26 Nov 2014 22:18:01 +0000 (15:18 -0700)]
Merge release-5-0 into master
Change-Id: If41f9c458e5d3aa5f841afcd71bc7f72a1f312c7
Roland Schulz [Wed, 26 Nov 2014 21:19:18 +0000 (16:19 -0500)]
Fix clang with boost 1.46.1 and 1.47
Change-Id: I53b9fcb1bcf553ea128f7edb2d15bec87fd4f9e9
Teemu Murtola [Mon, 17 Nov 2014 19:07:19 +0000 (21:07 +0200)]
Convert surface area calculation to C++
Minimal conversion to make min/max work, and remove a few unused
variables. This enables writing a C++ wrapper over the functionality,
as well as using the analysis neighborhood search internally to fix
several issues with the original algorithm.
Change-Id: If5089a749c8db9630e33edab3d2b79968a28a921
Teemu Murtola [Fri, 28 Feb 2014 04:37:13 +0000 (06:37 +0200)]
Add some unit tests for surface area calculation
There are a few real unit tests that compare the results against
analytical calculation for trivial cases, but the rest of the tests are
more of regression tests, designed to catch changes in the output.
I have not actually checked that the reference values for the more
complex values are correct.
Put the tests in surfacearea.cpp instead of nsc.cpp and renamed the
implementation to match in preparation for future cleanup.
The volume computation with PBC is not tested, since it is broken.
Change-Id: I18e3baa899b1c10fec0c8b43d8dd76a357297005
Mark Abraham [Sat, 14 Jun 2014 01:07:43 +0000 (21:07 -0400)]
Fix trjconv modifying fr.time unnecessarily
The was a bug where trjconv modified the input frame time, which
caused issues with xtc seeking. This has recently been fixed
by commit
81d10b25, but it's better not to mess with the input
frame and only modify the output frame.
Refs #1405, #1406
Change-Id: Ibf7c8655e916c12f8cc51523a67556cc7988e6e6
Ake Sandgren [Thu, 2 Oct 2014 15:29:46 +0000 (17:29 +0200)]
Fix problem with mixed affinity mask on different nodes.
If task distribution (with slurm for instance) causes both fully
allocated and not-fully allocated nodes to be assigned to the job then
there may be tasks with a all-cores affinity mask and tasks with a
not-all-cores affinity masks.
Fixes #1613
Change-Id: I71c0daa43a5dd42da57bfd09037806ce1d9334b5
Berk Hess [Fri, 24 Oct 2014 13:42:38 +0000 (15:42 +0200)]
Fixed triclinic 1xNx1 domain decomposition
With the Verlet scheme, 1D triclinic domain decomposition along
dimension y produces incorrect bounding boxes for the non-bonded grid.
This led to a lot of missing non-bonded interactions, which quickly
crashes any simulation affected by this.
Fixes #1631.
Change-Id: I9bd1fc9d983be839e0c9a8e62d47f6cf17684a03
Roland Schulz [Wed, 26 Nov 2014 01:10:03 +0000 (20:10 -0500)]
Merge release-5-0 into master
Conflicts:
src/gromacs/mdlib/domdec.cpp (std::max)
Change-Id: Iaa77c561f52cf5815c35cb08008cf578d1be2314
Mark Abraham [Tue, 18 Nov 2014 19:46:26 +0000 (20:46 +0100)]
Move code to domdec module folder
It's not a module yet, but here's where it will be. Updated
suppressions for module cyclic dependency warnings for domdec and
mdlib.
Change-Id: Ibe449ec8078fb7f25c4e95ebae28bac74db4f577
Berk Hess [Mon, 24 Nov 2014 14:36:13 +0000 (15:36 +0100)]
Fixed twin-range + freeze + constraints
With twin-range cut-offs, atoms which are both frozen and constrained
would experience very large or NaN forces, leading to constraint
warnings and errors.
Fixes #1639.
Change-Id: I1871a87054fec1149c9ed75872451df79a8ac2d2
Berk Hess [Thu, 16 Oct 2014 12:03:49 +0000 (14:03 +0200)]
Avoid NaN in rlist buffer calculation
With constraints on large masses and very low tolerance, the Verlet
buffer calculation could produce 1/0. This could lead to a somewhat
too smaller buffer.
Also added a missing scaling factor to the contribution of the third
derivative of the potential. This issue only caused a minor
overestimate for systems with constraints and without electrostatics.
Change-Id: I97e9d428a83f1b4954012ebd39bc49d397574f8c
Magnus Lundborg [Fri, 24 Oct 2014 08:24:22 +0000 (10:24 +0200)]
When writing TNG the input must contain molecule data.
When using trjcat to write TNG the input must also be TNG since
no tpr file is used.
Fixes #1626
Change-Id: I060e0174f9a471e134a4a899f83afaf02f4fce00
David van der Spoel [Sun, 23 Nov 2014 11:01:02 +0000 (12:01 +0100)]
Added warnings for ewald-geometry and surface-epsilon
ewald-geometry and surface-epsilon require the system dipole,
which will be incorrect when charge groups with net charge cross pbc.
grompp now checks and warns for this.
Refs #1645.
Change-Id: I02e317cbddb47256f942312ec53c5bab2b13be2a
Teemu Murtola [Fri, 21 Nov 2014 18:26:42 +0000 (20:26 +0200)]
Rename doc-* targets to be less generic
Now that there are a lot of other documentation targets as well, rename
the targets that generate Doxygen documentation to doxygen-*.
Also, rename the doc-check target to source-check, since it has evolved
to more than just a Doxygen documentation checker.
There is some backwards-compatibility glue for Jenkins that expects to
execute the doc-check target with the old name, but otherwise this
change should only be visible to those few developers that may generate
the Doxygen documentation themselves or run the check-source target.
Change-Id: If69b4becc2bdc17899fa17ffd12ea1516133457b
Roland Schulz [Fri, 10 Oct 2014 07:10:15 +0000 (03:10 -0400)]
Set GMX_SIMD to MIC for Xeon Phi
The reason GMX_SIMD wasn't set the same wasy as it is for other architectures,
was that we wanted it to also work for offload. But for offload it is easy to
overwrite it for a compilation unit being offloaded. Setting it normally for
native mode makes it more consistent and lets the user disable SIMD.
Change-Id: Id660a6560278039673ce0a4feaea3261647f40e8
Christoph Junghans [Tue, 9 Sep 2014 17:57:40 +0000 (11:57 -0600)]
Added .travis.yml
* make use of travis-ci.org once a commit is pushed to github
Change-Id: I96673fb392a49d842ae6d959ce844e727def0d26
Teemu Murtola [Tue, 16 Sep 2014 04:13:06 +0000 (07:13 +0300)]
Describe doc-check in developer manual
Move existing documentation for the doc-check target into the developer
manual, and extend it to cover also the include sorter. Add a separate
page that explains the include order enforced by the checker.
Various updates to make the description match better the current state
of affairs, with the module dependency graph embedded into the Doxygen
documentation.
Change-Id: I570e4f6af78453cffea89855f294a599f28dedd9
Teemu Murtola [Mon, 15 Sep 2014 18:49:54 +0000 (21:49 +0300)]
Skeleton for a developer manual
Add a skeleton for a developer manual, mostly focusing on the used tools
for now. Some existing content from Doxygen is linked from the new
pages, but subsequent changes will add more content (with initial focus
on explaining the formatting and documentation checking that Jenkins
does) and reformat the existing content to better fit in the new guide.
Additional sections can be added when a need is identified.
Change-Id: Id28daa9effe93fb66d03ef7b696b14ff2554e602
Szilárd Páll [Fri, 24 Oct 2014 00:37:34 +0000 (02:37 +0200)]
Allow increasing CUDA thread block size
This change parametrizes the CUDA kernel to allow increasing the number
of threads per block by processing multiple j-clusters concurrently on
additional pairs of warps. The change supports 1-, 2-, and 4-way
concurrent j-cluster processing resulting in 64, 128, and 256 threads
per block, respectively.
Due to register limitations, on current CUDA architectures the version
with 64 threads/block (equivalent with the original kernels) is fastest.
The latter configurations using 128 and 256 threads are 3-4% and 10-13%
slower on CC 3.5/5.2, respectively.
Change-Id: I7c12a826f9347f724827320184628d6b310c1424
Szilárd Páll [Sat, 8 Nov 2014 03:47:03 +0000 (04:47 +0100)]
Add CUDA compiler support for CC 5.0
With CUDA 6.5 and later compute capability 5.0 devices are supported, so
we generate cubin and PTX for these too and remove PTX 3.5.
This change also removes explicit optimization for CC 2.1 where
sm_20 binary code runs equally fast as sm_21.
Change-Id: I5a277c235b873afb2d1b2b12b5db64b370f1bade
Teemu Murtola [Fri, 31 Oct 2014 18:34:04 +0000 (20:34 +0200)]
Refactor position-evaluated selection keyword handling
Now the keywords themselves only return one value per position, not one
value per atom. The expansion to atoms is instead done by the caller in
evaluate.cpp. This makes it easier to reuse the methods for
implementing syntax like "z of ...", where the expansion is not
appropriate.
Also, add an additional unit test for code around this path.
Change-Id: Ida74244fe1c76eaf89c3d2e1ca9b6f6784b0476d
Szilárd Páll [Thu, 23 Oct 2014 23:49:46 +0000 (01:49 +0200)]
Start nvprof profiling at counter reset
When running in the NVIDIA profiler, to eliminate initial kernel
performance fluctuations, load balancing effects, as well as
initialization API calls from the traces, we now start the NVIDIA
profiler through the CUDA runtime API at the performance counter
resetting. This has an effect only if mdrun is started in nvprof with
profiling off.
Change-Id: Idfb3c86a96cb8b55cd874f641f4922b5517de6e3
Roland Schulz [Sat, 8 Nov 2014 16:16:08 +0000 (11:16 -0500)]
Use stat to check that paths are equivalent
Replaces the previous approach of comparing path as string.
The previous approach didn't work reliable on case insenstive filesystems
or with symlinks.
Change-Id: Iee11c172db89b98a26b3592be2f285110a85a632
Berk Hess [Fri, 31 Oct 2014 18:10:13 +0000 (19:10 +0100)]
Corrected parameter for sorting on search grid
The nbnxn grid search sorting range parameter used the whole box
height instead of the local DD cell height with domain decomposition
along z. This could not cause errors, but the search used O(N^2) time
with more than 4 domains along z, instead of O(N).
Change-Id: I324df6b15635a5ecff7c8bfeb124fdb933e5f845
Berk Hess [Fri, 17 Oct 2014 13:33:46 +0000 (15:33 +0200)]
Fix DD bonded interaction range print
A recent fix caused the DD setup printing of the maximum bonded
distance to instead print the max of the bonded and the list buffer
to the log file. This was a printing issue only.
Refs #1607.
Change-Id: I685e2e5e07f2f1a0a39c5eef4264a77ddfcecb31
Roland Schulz [Mon, 25 Nov 2013 21:06:26 +0000 (16:06 -0500)]
Add bitmask type
Adds support for arbirtrary number of OpenMP threads
Fixes #1386
Change-Id: I3747b4c7cd4b45d5901f710fcab47984c8913bc6
Carsten Kutzner [Wed, 5 Nov 2014 14:53:34 +0000 (15:53 +0100)]
Updated some .mdp file variables to use "-" instead of deprecated "_"
In the .mdp file, now a minus sign instead of an underline is
generally used. Some sections of the manual still listed .mdp
variable names in the old notation.
- three manual sections (pulling, enforced rotation and computational
electrophysiology) now use "-" instead of "_" in .mdp file variable
names
- updated output of rotation code to reflect this up-to-date notation
- no functional changes are in this patch
Change-Id: I7b00193d3e3d549ce8c9a728b307280c0bdd35ce
Justin Lemkul [Wed, 29 Oct 2014 13:10:36 +0000 (09:10 -0400)]
Removed spurious GLY dihedral in AMBER03.
Fixes #1632
Change-Id: I7302e42558454d546e14d7784077fb72c0dbfad7
Viveca Lindahl [Fri, 31 Oct 2014 17:09:10 +0000 (18:09 +0100)]
Fixed typo in cmake warning.
Change-Id: Id096ed65567fa671dce40970948b6435d9c357a0
Roland Schulz [Sun, 9 Nov 2014 07:42:58 +0000 (02:42 -0500)]
Backport valgrind suppression
Was added to master by
106ca9e6 and is now also required by
Jenkins for 5.0. Unclear what triggered that it is required.
Change-Id: I272600716fc986fe77a03fb86da827ba40d17adf
Roland Schulz [Tue, 28 Oct 2014 21:49:08 +0000 (17:49 -0400)]
Fix CMP0054 cmake 3.1.0 warnings
Change-Id: I5d91289a543240fa6f493b7f7738522f3c268c07
Roland Schulz [Wed, 5 Nov 2014 00:01:29 +0000 (19:01 -0500)]
Minor update of boost README
Change-Id: Ib6017efd89225771cea4f56d9e18994c25b8361e
Teemu Murtola [Thu, 30 Oct 2014 18:45:52 +0000 (20:45 +0200)]
Add more selection unit tests
Add tests for syntax like
res_com x < 5
res_com distance from ... < 5
as these were not covered well with earlier tests. Subsequent changes
will refactor the code responsible for evaluating these to allow syntax
like "z of com of ...".
Change-Id: I63b60cc2958fcf919a82b809f60465b8444904d2
Teemu Murtola [Thu, 30 Oct 2014 18:41:14 +0000 (20:41 +0200)]
Fix more memory leaks in default group creation
Free temporary arrays and strings allocated during the processing.
Don't know why these were not spotted in my earlier valgrind run (all of
these fixes are from working with the interactive selection tester).
Remove one unused array that was allocated and freed, but never used.
Change-Id: Ib78b6366bb5230499a77717f36afed652a1ec7fd
Teemu Murtola [Wed, 29 Oct 2014 18:40:34 +0000 (20:40 +0200)]
Fix memory leak in default index group creation
Temporary memory should be deallocated in the scope where it is
allocated, not only within some conditional paths.
Change-Id: I102d2b82ba531c25dd26728729b46c0b99c124e4
Roland Schulz [Thu, 30 Oct 2014 18:19:22 +0000 (19:19 +0100)]
Merge "Merge release-5-0 into master"
Teemu Murtola [Wed, 29 Oct 2014 18:52:16 +0000 (20:52 +0200)]
Merge release-5-0 into master
Conflicts:
src/external/gmock-1.7.0/CMakeLists.txt
- took both changes on the same line
src/gromacs/gmxpreprocess/calc_verletbuf.c
- took both changes on the same line
src/gromacs/mdlib/nbnxn_cuda/nbnxn_cuda_data_mgmt.cu
- discarded temporary common.h hack from release-5-0 in favor of
better solution in master
src/gromacs/selection/nbsearch.cpp
- resolved all conflicts in favor of master (were from a
back-ported change)
Change-Id: I34190b704bcb1c020dcc2a3586a25620fd3c3ccf
Teemu Murtola [Sat, 23 Aug 2014 03:27:24 +0000 (06:27 +0300)]
Improve analysis nbsearch grid mapping
Now the analysis neighborhood search implements its own version of
put_atoms_in_triclinic_unitcell(). While computing the index of the
correct grid cell, it is relatively easy to produce also the coordinates
that lay within that cell instead of using a separate call. This
provides two benefits:
- It avoids rare rounding problems if put_atoms_in_triclinic_unitcell()
would put the atom right at the edge of the box, but the mapping code
would consider it outside the box, causing out-of-range grid cell
index to be generated.
- It allows to customize the grid mapping more freely (e.g., to create
grids that are not periodic).
Backported from master with minor changes, fixes #1611. Kept commit
message the same; the second point will be only relevant for master.
Change-Id: Ib7602fa49a1b8f7882a63843322786b3e51e8e32
(cherry-picked from
b3e2e82 in master)
David van der Spoel [Tue, 28 Oct 2014 17:43:31 +0000 (18:43 +0100)]
Even more refreshing quotes.
From a variety of source.
Change-Id: If4f9ecc77c0941d460628443491d6b2ce4c6c1a0
Teemu Murtola [Tue, 28 Oct 2014 04:21:59 +0000 (06:21 +0200)]
Avoid common.h inclusion in CUDA code
This is solved nicer in master with I943f90f, but there is little value
in backporting that compared to the effort, since more or less all
changes from there would conflict. This may fix compilation problems
with some nvcc/boost combinations.
Change-Id: I82a26c60d366bfe5fdbcfd8a17bffdb1139ce017
Roland Schulz [Sun, 19 Oct 2014 17:37:22 +0000 (13:37 -0400)]
Add Windows support for host-/username
Change-Id: I971ffc811bddea8c8663d3d1ce95678eddfd5860
Teemu Murtola [Sun, 19 Oct 2014 09:36:02 +0000 (12:36 +0300)]
Move truncate() wrapper to futil.h
- Isolate platform-dependent file handling code better, removing
some platform-dependent #ifdefs and #includes from higher-level files.
- Make 'gmx trjconv -trunc' independent of the platform (should now also
work on Windows).
- Don't leak file handles when truncating on Windows.
Change-Id: I66162978c7b29499c8fc796d15b5f2649abf59dd
Teemu Murtola [Sun, 19 Oct 2014 06:55:37 +0000 (09:55 +0300)]
Encapsulate use of nice()
- Instead of platform-specific #ifdefs in otherwise platform-agnostic
code, encapsulate the call in sysinfo.h/.cpp.
- Make it fully CMake's responsibility to determine whether nice() is
going to be used. This could probably be simplified further (even to
the point of using HAVE_NICE directly), but since I'm not sure how all
the special cases work, left them as they were, but with TODO
comments.
- Remove some unused platform-specific #includes.
Change-Id: I9f1172ba2c02bb042574ad5d2b9ccbf5d31a3d6f
Teemu Murtola [Sun, 19 Oct 2014 04:24:36 +0000 (07:24 +0300)]
Collect more platform-specific code to sysinfo.h
Move a few existing functions of this type to sysinfo.* from
basenetwork.* and cstringutil.*. Mostly this is just code movement, but
there are two changes:
* getuid() usage is now controlled by HAVE_UNISTD_H instead of
HAVE_PWD_H, and excluded for MinGW (the latter is just a guess, but
many other similar calls are excluded for MinGW).
* Getting the current username uses getlogin_r() instead of getpwuid(),
following the implementation in tngio.cpp instead of the original
implementation in cstringutil.c. If this is OK, then we can get rid
of HAVE_PWD_H completely. Alternative would be to make the code use
getpwuid_r() if the old functionality is more desirable, but that is
a bit more effort, and requires some testing for platform support, as
that has not been used before in the code.
This allows removing platform-specific #ifdefs and #includes from
multiple files completely, making these aspects better encapsulated.
While checking the presence of #ifdefs, removed one useless #ifndef
GMX_FAHCORE from tngio.cpp; make_backup() already handles FAHCORE
internally.
Change-Id: I46876635fb254eb9558cce80b34ebe3078425375
Teemu Murtola [Sat, 18 Oct 2014 04:09:34 +0000 (07:09 +0300)]
Rename common.h to classhelpers.h
This reflects the remaining contents better. Since basedefinitions.h
was created, it is a more natural home for all truly common definitions,
so the old name no longer makes sense, either.
Change-Id: Ia2aa84e0b9696541f61490312aae9c18f719076f
Teemu Murtola [Sat, 18 Oct 2014 03:59:19 +0000 (06:59 +0300)]
Move basic macros from common.h to basedefinitions.h
Move GMX_UNUSED_VALUE and GMX_IGNORE_RETURN_VALUE from common.h to
basedefinitions.h to not force a boost dependency on code that uses only
these macros. In particular, CUDA code no longer includes common.h and
thus boost, which yet again works around broken boost/nvcc combinations.
basedefinitions.h is also a more natural place for these, as gmx_unused
is already declared there.
Fix one header that included others within an extern "C" block.
Change-Id: I943f90f068987456219ac7aad388258907018692
Roland Schulz [Sun, 19 Oct 2014 16:23:53 +0000 (12:23 -0400)]
Remove Catamount cmake
Was superseeded in 2006. No Cray XT3 should be around anymore.
Also gets rid of GMX_NO_SYSTEM which was only set by Catamount
Change-Id: Ib3c0450995e395ccdd2935b629da7930d9da5a4a
Roland Schulz [Sun, 19 Oct 2014 16:04:16 +0000 (12:04 -0400)]
Fix C++11 flags for MingW
Also replace all 'WIN32 AND NOT CYGWIN' with 'WIN32'. Since cmake
2.8.4 WIN32 doesn't include CYGWIN.
Change-Id: Ie93fa8a79daabcc3515e00cd1f95b76041417eb6
Roland Schulz [Fri, 17 Oct 2014 06:29:15 +0000 (02:29 -0400)]
Create scoped_cptr to replace scoped_ptr_sfree
Compareed to the previous scoped_ptr_sfree, scoped_cptr has a get
method, a bool-conversion, and allows to specify a custom deleter.
If just a guard is required then the scoped_guard_sfree typedef
can be used which doesn't require to specify a type.
Change-Id: I81cb50c4cf34e2d5a93e9c4b0043afcbc5feb535
Roland Schulz [Fri, 30 May 2014 05:02:23 +0000 (01:02 -0400)]
Used IWYU to partially clean up some includes
- Used trunk version
- Works well for C and simple C++ files
- Does not yet work without manual reverting some changes for more complicated
files e.g. options.cpp (see issues 105, 111, 135, and 138)
- IWYU doesn't check for conditional compiled code. I ran for
pme.c and wallcycle.c with -DDEBUG_PME/_WCYCLE
Given that one needs to fix a few issues and have to check for #ifdefs it is
still manual work, but I think quite a bit faster and more complete than doing
it manually.
Change-Id: I7c5569693a4b84f481f0f7afd85f0f01c33295cf
Berk Hess [Wed, 15 Oct 2014 11:24:14 +0000 (13:24 +0200)]
Removed f_novirsum reduction in Verlet scheme
There is a reduction over MPI ranks (dd_move_f) for fr->f_novirsum
for the exclusion correction of PME, which should not contribute
to the virial. But with the Verlet scheme this is unnecessary,
since the exclusions fully are handled in the non-bonded kernels.
Change-Id: I65b3dc90cf01de7b33e955074dfa32bd1940f781
Carsten Kutzner [Tue, 14 Oct 2014 08:56:14 +0000 (10:56 +0200)]
Fix order of ED/swap blocks during checkpoint reading.
For some reason the order of reading the swap state and the essential
dynamics (ED) state were mixed up in read_checkpoint(), which lead to
unreadable state files if both ED and swap functionality were used.
This fix changes the order of those blocks to be consistent with
write_checkpoint().
Change-Id: Ie234486b990861decef9e00fc5125f867001f814
Mark Abraham [Fri, 17 Oct 2014 15:50:12 +0000 (17:50 +0200)]
Extend gmock to turn off socket listener
On Linux, you'd generally expect this feature to work, but on Cray the
linker warns that the use of getaddrinfo() would require linking a
shared library at runtime. That never happens because we don't use the
socket listener feature, but we don't want people stressing over that
at build time.
Change-Id: I9c358a2923d7563809d471f2ea1767cf2e8bbe25
Mark Abraham [Thu, 2 Oct 2014 13:30:38 +0000 (15:30 +0200)]
Move code to what will become the ewald module
Moved PME implementation, PME load-balancing code and basic Ewald code
into new module. PME is the major component, but they share a little
bit of infrastructure and might one day pass the same tests (or test
each other) so there's some value to being together.
Moved legacyheaders/pme.h to module, split it into pme.h and
pme-internal.h, and moved the definition of splinevec into the only
source file that used it. pme-internal.h is intended as a temporary
measure - a place to put details that are internal to the module -
while developing some internal structure with more normal
module/file.{cpp,h} layout.
Moved legacyheaders/coulomb.h to the new module and named it
ewald/ewald-util.h, since it pretty much declares functions defined in
the former ewald-util.c (now ewald/ewald-util.c). Both files will
disappear in future clean up, though.
The internal structure of this module is still messy, but some
clean-up will follow. Then C++ conversion and Doxygen.
Renamed files to use '-' rather than '_', because I like using fewer
keystrokes.
Removed a useless extern on calc_ewaldcoeff_lj() declaration.
Introduced new cyclic dependency (and doc-check suppression) on mdlib,
because the PME-PP load-balancing needs to let mdlib/nbnxn_cuda know
about the changed parameters. It also uses some functions that are
currently declared in pme-internal.h, so there's no single-step
solution, even if nxn was already a module. In the long term, a PME
load-balancing module that talks to the ewald and nxn modules through
their library interfaces makes good sense, but we can't do that right
now.
Change-Id: I6ddd5579b6e52812d8fb110b4950a207b159ea4a
Mark Abraham [Sat, 20 Sep 2014 10:59:41 +0000 (12:59 +0200)]
Converted many integrator-loop files to C++
Moved code from mdlib/init.c to other locations, to avoid creating an
init.cpp file that was freaking out the ASAN build, perhaps via linker
filename clash with utility/init.cpp. Renamed gmxlib/mvdata.c to
mdlib/broadcaststructs.cpp since that makes sense now.
Converted min and max to std:: varieties, and provided explicit
template specifications where required. Used std::sqrt to help the
compiler do automatic template type deduction.
Used static_cast<int> to suppress warnings about comparison between
signed and unsigned integers, and conversion of double to int where
that is the intention of the code.
Converted a few variables between int, real and double to correctly
reflect the intention of the code.
Fixed use of double for number of steps completed, which should be
gmx_int64_t. Formally, this is a bug in release-5-0, but you'd have to
do 2^53 steps to notice it.
Used const real to define constants so that the std::pow overloads
on Windows can work smoothly.
Added some width specifiers on environment-variable scanf calls
to avoid problems.
Removed many unused variables.
Added inline cppcheck suppressions for writes to variables that are
only read when OpenMP is being used. A recent change introduced such
variables to work around compilers that have bugs if you call a
function inside the #pragma omp, but those variables are unused if the
compilation does not use OpenMP.
Added C++ guard to free-energy kernel header.
Change-Id: Id190e36758fa6aef68b14e5c9b78eacfb0a86949
David van der Spoel [Wed, 21 May 2014 14:38:25 +0000 (16:38 +0200)]
Reorganizing analysis of correlation functions.
Moved routines computing correlation functions and
fitting to those in a special subdirectory not
dependent on gmxana. The autocorrelation now functions
as it should (previously there were issues with long
time correlations).
Replaced Levenberg-Marquardt algorithm by another implementation
from http://apps.jcns.fz-juelich.de/doku/sc/lmfit,
under a FreeBSD license. This code is in src/external/lmfit.
Added Doxygen comments to all headers.
Fixed wordcount algorithm in xvgr.c
The routines many_auto_correl and many_cross_correl
are optimized using OpenMP.
Change-Id: Ifbe003437db36e35d606cd942bb38deeef86ea64
Teemu Murtola [Sun, 5 Oct 2014 18:56:07 +0000 (21:56 +0300)]
Check config.h usage in doc-check
Make the doc-check target check that config.h is included if and only if
it is needed by the source file. Add suppressions for locations where
this is currently not the case.
Change-Id: I6110e9887d4df3c4644f20003faa3120ea9048ea
Teemu Murtola [Sun, 28 Sep 2014 03:51:52 +0000 (06:51 +0300)]
Use analysis nbsearch in insert-molecules
Advantages:
- This reduces the amount of code by ~90% compared to what addconf.c
has, making it significantly easier to understand.
- Now the tool is independent of potential changes in the
mdrun-specific neighborhood search.
- Memory leaks related to addconf.c are gone.
- The neighborhood search is terminated as soon as one pair within the
cutoff is found, potentially making it faster. This likely offsets
any performance differences between the nbsearch implementations.
The unit tests are ~35% faster.
- Confusing mdrun-specific output related to the neighborhood
searching is gone. This includes notes that "This file uses the
deprecated 'group' cutoff_scheme" and references to Coulomb or VdW
tables and cutoffs.
Change-Id: Iba82858b9a2b43b6e10a49cd3964b99b22996166
Mark Abraham [Mon, 13 Oct 2014 10:02:59 +0000 (12:02 +0200)]
Merge release-4-6 into release-5-0
Change-Id: I08dc9d5e776480e7686a159cfcd981ba46c510ad