| 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 |
|
|