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