Refactored 'solvate' and 'insert-molecules'
commitc7a1e248090ccd3f7e1501bec4a687bac30232fb
authorMark Abraham <mark.j.abraham@gmail.com>
Mon, 3 Feb 2014 00:13:27 +0000 (3 01:13 +0100)
committerRoland Schulz <roland@rschulz.eu>
Wed, 12 Feb 2014 02:42:21 +0000 (12 03:42 +0100)
tree89b643ea11a70780d6f35f8c606e6c2a6cb3b533
parent8fa45f0b685bfe8b2faeaf40c3f70293a1d4673d
Refactored 'solvate' and 'insert-molecules'

The use of "VDW radii" in these tools was misleading in both
documentation and code. VDW radii are used if they are found, but they
used to be scaled down via a hidden option, and the user had to
magically know what kind of number to give to -vdwd. Documented what
is going on here.

Made that hidden option public, and renamed -vdwd to -radius. Renamed
variables named r to be more reflective of their use as additive
inter-atomic exclusion half-distances. This clarified what "mk_vdw"
was really doing, and that a call to it should not be living at the
end of a function that reads in a conformation from a file. Renamed
lots of things and documented more behaviour. Freed some memory
that should have been freed all along.

This exposed some horrible use of t_atoms. Fixed the code to always
allocate on the heap, call the initializer and call the destructor.

Change-Id: I3f8bbbed07442ccc6cd8e563168f5a68736c6b40
src/gromacs/gmxpreprocess/insert-molecules.cpp
src/gromacs/gmxpreprocess/read-conformation.cpp
src/gromacs/gmxpreprocess/read-conformation.h
src/gromacs/gmxpreprocess/solvate.cpp