| 162 | 
  | 
    return pressure; | 
| 163 | 
  | 
  } | 
| 164 | 
  | 
 | 
| 165 | 
+ | 
  double Thermo::getPressure(int direction) { | 
| 166 | 
+ | 
 | 
| 167 | 
+ | 
    // Relies on the calculation of the full molecular pressure tensor | 
| 168 | 
+ | 
 | 
| 169 | 
+ | 
           | 
| 170 | 
+ | 
    Mat3x3d tensor; | 
| 171 | 
+ | 
    double pressure; | 
| 172 | 
+ | 
 | 
| 173 | 
+ | 
    tensor = getPressureTensor(); | 
| 174 | 
+ | 
 | 
| 175 | 
+ | 
    pressure = OOPSEConstant::pressureConvert * tensor(direction, direction); | 
| 176 | 
+ | 
 | 
| 177 | 
+ | 
    return pressure; | 
| 178 | 
+ | 
  } | 
| 179 | 
+ | 
 | 
| 180 | 
+ | 
 | 
| 181 | 
+ | 
 | 
| 182 | 
  | 
  Mat3x3d Thermo::getPressureTensor() { | 
| 183 | 
  | 
    // returns pressure tensor in units amu*fs^-2*Ang^-1 | 
| 184 | 
  | 
    // routine derived via viral theorem description in: |