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