ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE/libmdtools/Thermo.cpp
(Generate patch)

Comparing trunk/OOPSE/libmdtools/Thermo.cpp (file contents):
Revision 708 by tim, Wed Aug 20 22:23:34 2003 UTC vs.
Revision 755 by mmeineke, Tue Sep 9 20:35:25 2003 UTC

# Line 166 | Line 166 | double Thermo::getPressure() {
166    return pressure;
167   }
168  
169 + double Thermo::getPressureX() {
170  
171 +  // Relies on the calculation of the full molecular pressure tensor
172 +  
173 +  const double p_convert = 1.63882576e8;
174 +  double press[3][3];
175 +  double pressureX;
176 +
177 +  this->getPressureTensor(press);
178 +
179 +  pressureX = p_convert * press[0][0];
180 +
181 +  return pressureX;
182 + }
183 +
184 + double Thermo::getPressureY() {
185 +
186 +  // Relies on the calculation of the full molecular pressure tensor
187 +  
188 +  const double p_convert = 1.63882576e8;
189 +  double press[3][3];
190 +  double pressureY;
191 +
192 +  this->getPressureTensor(press);
193 +
194 +  pressureY = p_convert * press[1][1];
195 +
196 +  return pressureY;
197 + }
198 +
199 + double Thermo::getPressureZ() {
200 +
201 +  // Relies on the calculation of the full molecular pressure tensor
202 +  
203 +  const double p_convert = 1.63882576e8;
204 +  double press[3][3];
205 +  double pressureZ;
206 +
207 +  this->getPressureTensor(press);
208 +
209 +  pressureZ = p_convert * press[2][2];
210 +
211 +  return pressureZ;
212 + }
213 +
214 +
215   void Thermo::getPressureTensor(double press[3][3]){
216    // returns pressure tensor in units amu*fs^-2*Ang^-1
217    // routine derived via viral theorem description in:

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines