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