| 52 | 
  | 
 | 
| 53 | 
  | 
namespace OpenMD { | 
| 54 | 
  | 
Globals::Globals() { | 
| 55 | 
+ | 
   | 
| 56 | 
+ | 
  flucQpars_ = new FluctuatingChargeParameters(); | 
| 57 | 
+ | 
  rnemdPars_ = new RNEMDParameters(); | 
| 58 | 
+ | 
 | 
| 59 | 
  | 
  DefineParameter(ForceField, "forceField") | 
| 60 | 
  | 
  | 
| 61 | 
  | 
  DefineOptionalParameter(TargetTemp, "targetTemp"); | 
| 95 | 
  | 
  DefineOptionalParameter(DampingAlpha, "dampingAlpha"); | 
| 96 | 
  | 
  DefineOptionalParameter(SurfaceTension, "surfaceTension"); | 
| 97 | 
  | 
  DefineOptionalParameter(PrintPressureTensor, "printPressureTensor"); | 
| 98 | 
+ | 
  DefineOptionalParameter(PrintHeatFlux, "printHeatFlux"); | 
| 99 | 
  | 
  DefineOptionalParameter(TaggedAtomPair, "taggedAtomPair"); | 
| 100 | 
  | 
  DefineOptionalParameter(PrintTaggedPairDistance, "printTaggedPairDistance"); | 
| 101 | 
  | 
  DefineOptionalParameter(SwitchingFunctionType, "switchingFunctionType"); | 
| 133 | 
  | 
  DefineOptionalParameterWithDefaultValue(UseSphericalBoundaryConditions, "useSphericalBoundaryConditions", false); | 
| 134 | 
  | 
  DefineOptionalParameterWithDefaultValue(AccumulateBoxDipole, "accumulateBoxDipole", false); | 
| 135 | 
  | 
 | 
| 131 | 
– | 
  DefineOptionalParameterWithDefaultValue(UseRNEMD, "useRNEMD", false); | 
| 132 | 
– | 
  DefineOptionalParameterWithDefaultValue(RNEMD_exchangeTime, "RNEMD_exchangeTime", 100.0); | 
| 133 | 
– | 
  DefineOptionalParameterWithDefaultValue(RNEMD_nBins, "RNEMD_nBins", 16); | 
| 134 | 
– | 
  DefineOptionalParameterWithDefaultValue(RNEMD_logWidth, "RNEMD_logWidth", 16); | 
| 135 | 
– | 
  DefineOptionalParameterWithDefaultValue(RNEMD_exchangeType, "RNEMD_exchangeType", "KineticScale"); | 
| 136 | 
– | 
  DefineOptionalParameterWithDefaultValue(RNEMD_targetFlux, "RNEMD_targetFlux", 0.0); | 
| 137 | 
– | 
  DefineOptionalParameterWithDefaultValue(RNEMD_objectSelection, "RNEMD_objectSelection", "select all"); | 
| 138 | 
– | 
  DefineOptionalParameterWithDefaultValue(RNEMD_binShift, "RNEMD_binShift", false); | 
| 139 | 
– | 
  DefineOptionalParameterWithDefaultValue(RNEMD_outputDimensionalTemperature, "RNEMD_outputDimensionalTemperature", false); | 
| 136 | 
  | 
  DefineOptionalParameterWithDefaultValue(UseRestraints, "useRestraints", false); | 
| 137 | 
  | 
  DefineOptionalParameterWithDefaultValue(Restraint_file, "Restraint_file", "idealCrystal.in"); | 
| 138 | 
  | 
  DefineOptionalParameterWithDefaultValue(UseThermodynamicIntegration, "useThermodynamicIntegration", false); | 
| 139 | 
  | 
  DefineOptionalParameterWithDefaultValue(HULL_Method,"HULL_Method","Convex"); | 
| 140 | 
  | 
 | 
| 145 | 
– | 
 | 
| 146 | 
– | 
 | 
| 141 | 
  | 
  deprecatedKeywords_.insert("nComponents"); | 
| 142 | 
  | 
  deprecatedKeywords_.insert("nZconstraints"); | 
| 143 | 
  | 
  deprecatedKeywords_.insert("initialConfig"); | 
| 203 | 
  | 
  CheckParameter(FrozenBufferRadius, isPositive()); | 
| 204 | 
  | 
  CheckParameter(LangevinBufferRadius, isPositive()); | 
| 205 | 
  | 
  CheckParameter(NeighborListNeighbors, isPositive()); | 
| 212 | 
– | 
  CheckParameter(RNEMD_exchangeTime, isPositive()); | 
| 213 | 
– | 
  CheckParameter(RNEMD_nBins, isPositive() && isEven()); | 
| 214 | 
– | 
  CheckParameter(RNEMD_exchangeType, isEqualIgnoreCase("KineticSwap") || isEqualIgnoreCase("KineticScale") || isEqualIgnoreCase("Px") || isEqualIgnoreCase("Py") || isEqualIgnoreCase("Pz") || isEqualIgnoreCase("PxScale") || isEqualIgnoreCase("PyScale") || isEqualIgnoreCase("PzScale")); | 
| 206 | 
  | 
  CheckParameter(HULL_Method, isEqualIgnoreCase("Convex") || isEqualIgnoreCase("AlphaShape"));  | 
| 207 | 
  | 
  CheckParameter(Alpha, isPositive());  | 
| 217 | 
– | 
 | 
| 208 | 
  | 
   | 
| 209 | 
  | 
  for(std::vector<Component*>::iterator i = components_.begin(); i != components_.end(); ++i) { | 
| 210 | 
  | 
    if (!(*i)->findMoleculeStamp(moleculeStamps_)) { | 
| 230 | 
  | 
    return true; | 
| 231 | 
  | 
} | 
| 232 | 
  | 
 | 
| 233 | 
+ | 
bool Globals::addFluctuatingChargeParameters(FluctuatingChargeParameters* fqp) { | 
| 234 | 
+ | 
  if (flucQpars_ != NULL) | 
| 235 | 
+ | 
    delete flucQpars_; | 
| 236 | 
+ | 
     | 
| 237 | 
+ | 
  flucQpars_ = fqp; | 
| 238 | 
+ | 
  return true; | 
| 239 | 
+ | 
} | 
| 240 | 
+ | 
 | 
| 241 | 
+ | 
bool Globals::addRNEMDParameters(RNEMDParameters* rnemdPars) { | 
| 242 | 
+ | 
  if (rnemdPars_ != NULL) | 
| 243 | 
+ | 
    delete rnemdPars_; | 
| 244 | 
+ | 
     | 
| 245 | 
+ | 
  rnemdPars_ = rnemdPars; | 
| 246 | 
+ | 
  return true; | 
| 247 | 
+ | 
} | 
| 248 | 
+ | 
 | 
| 249 | 
  | 
bool Globals::addMoleculeStamp(MoleculeStamp* molStamp) { | 
| 250 | 
  | 
    std::string molStampName = molStamp->getName(); | 
| 251 | 
  | 
    std::map<std::string, MoleculeStamp*>::iterator i; |