70 |
|
(p_convert*NkBT*tb2)); |
71 |
|
} |
72 |
|
|
73 |
+ |
template<typename T> void NPTi<T>::calcVelScale(void) { |
74 |
+ |
vScale = chi + eta; |
75 |
+ |
} |
76 |
+ |
|
77 |
|
template<typename T> void NPTi<T>::getVelScaleA(double sc[3], double vel[3]) { |
78 |
|
int i; |
79 |
|
|
80 |
< |
for(i=0; i<3; i++) sc[i] = vel[i] * ( chi + eta ); |
80 |
> |
for(i=0; i<3; i++) sc[i] = vel[i] * vScale; |
81 |
|
} |
82 |
|
|
83 |
|
template<typename T> void NPTi<T>::getVelScaleB(double sc[3], int index ){ |
84 |
|
int i; |
85 |
|
|
86 |
< |
for(i=0; i<3; i++) sc[i] = oldVel[index*3 + i] * ( chi + eta ); |
86 |
> |
for(i=0; i<3; i++) sc[i] = oldVel[index*3 + i] * vScale; |
87 |
|
} |
88 |
|
|
89 |
|
|