Fix for the pressure in MTTK with constraints + dispersion
commitb18b78331458910f0186f299af6c929549c79f61
authorMichael Shirts <michael.shirts@virginia.edu>
Sat, 1 Dec 2012 19:55:47 +0000 (1 14:55 -0500)
committerMark Abraham <mark.j.abraham@gmail.com>
Thu, 20 Dec 2012 13:14:43 +0000 (20 14:14 +0100)
tree98c418a24865f5348544c99ab5747ce6ae095116
parent4d3a2d48d7bdde9d1297230221ea2d088af0e8cc
Fix for the pressure in MTTK with constraints + dispersion

Fix for pressure in MTTK - broken at some point in the
merge process.  The MTTK+constraints code definitely needs to
be simplified for 5.0, it's too easy to break on code rearrangements
right now!

Testing with simple systems over long time shows that the pressure is
unbiased, while the PR pressure is off by a factor of 1/N (which is totally
fine for most realistic systems).

Fixes redmine bug #1061
Change-Id: Ie82cf8d6bdeb3e47db044ae34fc42c077ad73c75
src/kernel/md.c
src/mdlib/coupling.c