--- branches/development/src/nonbonded/SHAPES.cpp 2011/09/15 16:24:03 1635 +++ branches/development/src/nonbonded/SHAPES.cpp 2013/05/15 15:09:35 1874 @@ -35,8 +35,9 @@ * * [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). * [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). - * [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008). - * [4] Vardeman & Gezelter, in progress (2009). + * [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008). + * [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). + * [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). */ #include @@ -266,6 +267,8 @@ namespace OpenMD { RealType proji = sqrt(r * 1.0e-12); Vector3d dcpidx(1.0 / proji, 0.0, + + // pickup the ball here! dcpidx = 1.0_dp / proji dcpidy = 0.0_dp