| 57 |
|
DefineOptionalParameter(Ensemble, "ensemble"); |
| 58 |
|
DefineOptionalParameter(Dt, "dt"); |
| 59 |
|
DefineOptionalParameter(RunTime, "runTime"); |
| 60 |
< |
DefineOptionalParameter(InitialConfig, "initialConfig"); |
| 60 |
> |
//DefineOptionalParameter(InitialConfig, "initialConfig"); |
| 61 |
|
DefineOptionalParameter(FinalConfig, "finalConfig"); |
| 62 |
|
DefineOptionalParameter(SampleTime, "sampleTime"); |
| 63 |
|
DefineOptionalParameter(ResetTime, "resetTime"); |
| 99 |
|
DefineOptionalParameter(CutoffPolicy, "cutoffPolicy"); |
| 100 |
|
DefineOptionalParameter(SwitchingFunctionType, "switchingFunctionType"); |
| 101 |
|
DefineOptionalParameter(HydroPropFile, "HydroPropFile"); |
| 102 |
+ |
DefineOptionalParameter(Viscosity, "viscosity"); |
| 103 |
+ |
DefineOptionalParameter(BeadSize, "beadSize"); |
| 104 |
+ |
DefineOptionalParameter(FrozenBufferRadius, "frozenBufferRadius"); |
| 105 |
+ |
DefineOptionalParameter(LangevinBufferRadius, "langevinBufferRadius"); |
| 106 |
+ |
|
| 107 |
|
DefineOptionalParameterWithDefaultValue(UsePeriodicBoundaryConditions, "usePeriodicBoundaryConditions", true); |
| 108 |
|
DefineOptionalParameterWithDefaultValue(UseInitalTime, "useInitialTime", false); |
| 109 |
|
DefineOptionalParameterWithDefaultValue(UseIntialExtendedSystemState, "useInitialExtendedSystemState", false); |
| 118 |
|
DefineOptionalParameterWithDefaultValue(OutputForceVector, "outputForceVector", 0); |
| 119 |
|
DefineOptionalParameterWithDefaultValue(SkinThickness, "skinThickness", 1.0); |
| 120 |
|
DefineOptionalParameterWithDefaultValue(StatFileFormat, "statFileFormat", "TIME|TOTAL_ENERGY|POTENTIAL_ENERGY|KINETIC_ENERGY|TEMPERATURE|PRESSURE|VOLUME|CONSERVED_QUANTITY"); |
| 121 |
+ |
DefineOptionalParameterWithDefaultValue(UseSphericalBoundaryConditions, "useSphericalBoundaryConditions", false); |
| 122 |
+ |
DefineOptionalParameterWithDefaultValue(AccumulateBoxDipole, "accumulateBoxDipole", false); |
| 123 |
|
|
| 124 |
|
|
| 125 |
|
deprecatedKeywords_.insert("nComponents"); |
| 126 |
|
deprecatedKeywords_.insert("nZconstraints"); |
| 127 |
+ |
deprecatedKeywords_.insert("initialConfig"); |
| 128 |
|
|
| 129 |
|
} |
| 130 |
|
|
| 138 |
|
|
| 139 |
|
CheckParameter(ForceField, isNotEmpty()); |
| 140 |
|
CheckParameter(TargetTemp, isPositive()); |
| 141 |
< |
CheckParameter(Ensemble, isEqualIgnoreCase("NVE") || isEqualIgnoreCase("NVT") || isEqualIgnoreCase("NPTi") || isEqualIgnoreCase("NPTf") || isEqualIgnoreCase("NPTxyz") || isEqualIgnoreCase("NPAT") || isEqualIgnoreCase("LANGEVINDYNAMICS") ); |
| 141 |
> |
CheckParameter(Ensemble, isEqualIgnoreCase("NVE") || isEqualIgnoreCase("NVT") || isEqualIgnoreCase("NPTi") || isEqualIgnoreCase("NPTf") || isEqualIgnoreCase("NPTxyz") || isEqualIgnoreCase("NPAT") || isEqualIgnoreCase("LANGEVINDYNAMICS") || isEqualIgnoreCase("LD") || isEqualIgnoreCase("NPRT") || isEqualIgnoreCase("NPGT")); |
| 142 |
|
CheckParameter(Dt, isPositive()); |
| 143 |
|
CheckParameter(RunTime, isPositive()); |
| 144 |
< |
CheckParameter(InitialConfig, isNotEmpty()); |
| 144 |
> |
//CheckParameter(InitialConfig, isNotEmpty()); |
| 145 |
|
CheckParameter(FinalConfig, isNotEmpty()); |
| 146 |
|
CheckParameter(SampleTime, isNonNegative()); |
| 147 |
|
CheckParameter(ResetTime, isNonNegative()); |
| 150 |
|
CheckParameter(SwitchingRadius, isNonNegative()); |
| 151 |
|
CheckParameter(Dielectric, isPositive()); |
| 152 |
|
CheckParameter(ThermalTime, isNonNegative()); |
| 153 |
< |
CheckParameter(TargetPressure, isPositive()); |
| 153 |
> |
// CheckParameter(TargetPressure, isPositive()); |
| 154 |
|
CheckParameter(TauThermostat, isPositive()); |
| 155 |
|
CheckParameter(TauBarostat, isPositive()); |
| 156 |
|
CheckParameter(ZconsTime, isPositive()); |
| 166 |
|
CheckParameter(MinimizerLSMaxIter, isPositive()); |
| 167 |
|
CheckParameter(ZconsGap, isPositive()); |
| 168 |
|
CheckParameter(ZconsFixtime, isPositive()); |
| 169 |
< |
CheckParameter(ThermodynamicIntegrationLambda, isPositive()); |
| 169 |
> |
CheckParameter(ThermodynamicIntegrationLambda, isNonNegative()); |
| 170 |
|
CheckParameter(ThermodynamicIntegrationK, isPositive()); |
| 171 |
|
CheckParameter(ForceFieldVariant, isNotEmpty()); |
| 172 |
|
CheckParameter(ForceFieldFileName, isNotEmpty()); |
| 173 |
|
CheckParameter(ThermIntDistSpringConst, isPositive()); |
| 174 |
|
CheckParameter(ThermIntThetaSpringConst, isPositive()); |
| 175 |
|
CheckParameter(ThermIntOmegaSpringConst, isPositive()); |
| 176 |
< |
CheckParameter(SurfaceTension, isPositive()); |
| 176 |
> |
// CheckParameter(SurfaceTension, isNonNegative()); |
| 177 |
|
CheckParameter(ElectrostaticSummationMethod, isEqualIgnoreCase("NONE") || isEqualIgnoreCase("SHIFTED_POTENTIAL") || isEqualIgnoreCase("SHIFTED_FORCE") || isEqualIgnoreCase("REACTION_FIELD")); |
| 178 |
|
CheckParameter(ElectrostaticScreeningMethod, isEqualIgnoreCase("UNDAMPED") || isEqualIgnoreCase("DAMPED")); |
| 179 |
|
CheckParameter(CutoffPolicy, isEqualIgnoreCase("MIX") || isEqualIgnoreCase("MAX") || isEqualIgnoreCase("TRADITIONAL")); |
| 185 |
|
CheckParameter(ThermIntOmegaSpringConst, isPositive()); |
| 186 |
|
CheckParameter(DampingAlpha,isNonNegative()); |
| 187 |
|
CheckParameter(SkinThickness, isPositive()); |
| 188 |
< |
|
| 188 |
> |
CheckParameter(Viscosity, isNonNegative()); |
| 189 |
> |
CheckParameter(BeadSize, isPositive()); |
| 190 |
> |
CheckParameter(FrozenBufferRadius, isPositive()); |
| 191 |
> |
CheckParameter(LangevinBufferRadius, isPositive()); |
| 192 |
|
for(std::vector<Component*>::iterator i = components_.begin(); i != components_.end(); ++i) { |
| 193 |
|
if (!(*i)->findMoleculeStamp(moleculeStamps_)) { |
| 194 |
|
std::ostringstream oss; |