92 |
|
void NPTi::getPosScale(const Vector3d& pos, const Vector3d& COM, |
93 |
|
int index, Vector3d& sc){ |
94 |
|
/**@todo*/ |
95 |
< |
sc = (oldPos[index] + pos)/2.0 -COM; |
95 |
> |
sc = (oldPos[index] + pos)/(RealType)2.0 -COM; |
96 |
|
sc *= eta; |
97 |
|
} |
98 |
|
|
99 |
|
void NPTi::scaleSimBox(){ |
100 |
|
|
101 |
< |
double scaleFactor; |
101 |
> |
RealType scaleFactor; |
102 |
|
|
103 |
|
scaleFactor = exp(dt*eta); |
104 |
|
|
123 |
|
return ( fabs(prevEta - eta) <= etaTolerance ); |
124 |
|
} |
125 |
|
|
126 |
< |
double NPTi::calcConservedQuantity(){ |
126 |
> |
RealType NPTi::calcConservedQuantity(){ |
127 |
|
|
128 |
|
chi= currentSnapshot_->getChi(); |
129 |
|
integralOfChidt = currentSnapshot_->getIntegralOfChiDt(); |
138 |
|
// of freedom). |
139 |
|
fkBT = info_->getNdf()*OOPSEConstant::kB *targetTemp; |
140 |
|
|
141 |
< |
double conservedQuantity; |
142 |
< |
double Energy; |
143 |
< |
double thermostat_kinetic; |
144 |
< |
double thermostat_potential; |
145 |
< |
double barostat_kinetic; |
146 |
< |
double barostat_potential; |
141 |
> |
RealType conservedQuantity; |
142 |
> |
RealType Energy; |
143 |
> |
RealType thermostat_kinetic; |
144 |
> |
RealType thermostat_potential; |
145 |
> |
RealType barostat_kinetic; |
146 |
> |
RealType barostat_potential; |
147 |
|
|
148 |
|
Energy =thermo.getTotalE(); |
149 |
|
|