ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/trunk/src/integrators/NgammaT.cpp
(Generate patch)

Comparing trunk/src/integrators/NgammaT.cpp (file contents):
Revision 1782 by gezelter, Wed Aug 22 02:28:28 2012 UTC vs.
Revision 2069 by gezelter, Thu Mar 5 16:30:23 2015 UTC

# Line 35 | Line 35
35   *                                                                      
36   * [1]  Meineke, et al., J. Comp. Chem. 26, 252-271 (2005).            
37   * [2]  Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006).          
38 < * [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008).          
38 > * [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008).          
39   * [4]  Kuang & Gezelter,  J. Chem. Phys. 133, 164101 (2010).
40   * [5]  Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011).
41   */
# Line 171 | Line 171 | namespace OpenMD {
171      RealType barostat_potential = (targetPressure * thermo.getVolume() / PhysicalConstants::pressureConvert) /PhysicalConstants::energyConvert;
172  
173      Mat3x3d hmat = snap->getHmat();
174 –    RealType hz = hmat(2, 2);
174      RealType area = hmat(0,0) * hmat(1, 1);
175  
176 <    RealType conservedQuantity = totalEnergy + thermostat_kinetic + thermostat_potential +
177 <      barostat_kinetic + barostat_potential - surfaceTension * area/ PhysicalConstants::energyConvert;
176 >    RealType conservedQuantity = totalEnergy + thermostat_kinetic
177 >      + thermostat_potential + barostat_kinetic + barostat_potential
178 >      - surfaceTension * area/ PhysicalConstants::energyConvert;
179  
180      return conservedQuantity;
181  

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines