ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/trunk/src/io/ForceFieldOptions.cpp
(Generate patch)

Comparing trunk/src/io/ForceFieldOptions.cpp (file contents):
Revision 805 by gezelter, Thu Dec 8 22:05:16 2005 UTC vs.
Revision 1005 by gezelter, Sun Jul 16 02:48:00 2006 UTC

# Line 38 | Line 38
38   * University of Notre Dame has been advised of the possibility of
39   * such damages.
40   */
41 <
41 > #define __C
42   #include "io/ForceFieldOptions.hpp"
43 +
44   namespace oopse {
45  
46      ForceFieldOptions::ForceFieldOptions() {
47        DefineOptionalParameter(Name, "Name");
48 +      DefineOptionalParameter(CutoffPolicy, "cutoffPolicy");
49        DefineOptionalParameterWithDefaultValue(vdWtype, "vdWtype", "Lennard-Jones");
50        DefineOptionalParameterWithDefaultValue(DistanceMixingRule, "DistanceMixingRule", "arithmetic");
51        DefineOptionalParameterWithDefaultValue(DistanceType, "DistanceType", "sigma");
# Line 51 | Line 53 | namespace oopse {
53        DefineOptionalParameterWithDefaultValue(EnergyUnitScaling, "EnergyUnitScaling", 1.0);
54        DefineOptionalParameterWithDefaultValue(DistanceUnitScaling, "DistanceUnitScaling", 1.0);
55        DefineOptionalParameterWithDefaultValue(AngleUnitScaling, "AngleUnitScaling", 1.0);
56 <      DefineOptionalParameterWithDefaultValue(TorsionAngleConvention, "TorsionAngleConvention", "180 is trans");
56 >      DefineOptionalParameterWithDefaultValue(TorsionAngleConvention, "TorsionAngleConvention", "180_is_trans");
57        DefineOptionalParameterWithDefaultValue(vdw14scale, "vdW-14-scale", 0.0);
58 <      DefineOptionalParameterWithDefaultValue(electrostatic14scale, "electrostatic-14-scale", 0.0);
58 >      DefineOptionalParameterWithDefaultValue(electrostatic14scale, "electrostatic-14-scale", 0.0)
59        DefineOptionalParameterWithDefaultValue(dielectric, "dielectric", 1.0);
60 +      DefineOptionalParameterWithDefaultValue(GayBerneMu, "GayBerneMu", 2.0);
61 +      DefineOptionalParameterWithDefaultValue(GayBerneNu, "GayBerneNu", 1.0);
62      }
63  
64 +
65 +  void ForceFieldOptions::makeFortranOptions(ForceOptions& fortranForceOptions){
66 +    
67 +    fortranForceOptions.vdw14scale = this->getvdw14scale();
68 +    fortranForceOptions.electrostatic14scale = this->getelectrostatic14scale();
69 +    fortranForceOptions.GayBerneMu = this->getGayBerneMu();
70 +    fortranForceOptions.GayBerneNu = this->getGayBerneNu();
71 +
72 +
73 +    std::string DistanceMix = this->getDistanceMixingRule();
74 +    toUpper(DistanceMix);
75 +    if(DistanceMix == "ARITHMETIC"){
76 +      fortranForceOptions.DistanceMixingRule = ARITHMETIC_MIXING_RULE;
77 +    }else if(DistanceMix == "GEOMETRIC"){
78 +      fortranForceOptions.DistanceMixingRule = GEOMETRIC_MIXING_RULE;
79 +    }else if(DistanceMix == "CUBIC"){
80 +      fortranForceOptions.DistanceMixingRule = CUBIC_MIXING_RULE;
81 +    }else{
82 +      fortranForceOptions.DistanceMixingRule = ARITHMETIC_MIXING_RULE;
83 +    }
84 +    std::string EnergyMix = this->getEnergyMixingRule();  
85 +    toUpper(EnergyMix);
86 +    if(EnergyMix == "ARITHMETIC"){
87 +      fortranForceOptions.EnergyMixingRule = ARITHMETIC_MIXING_RULE;
88 +    }else if(EnergyMix == "GEOMETRIC"){
89 +      fortranForceOptions.EnergyMixingRule = GEOMETRIC_MIXING_RULE;
90 +    }else if(EnergyMix == "HHG"){
91 +      fortranForceOptions.EnergyMixingRule = HHG_MIXING_RULE;
92 +    }else{
93 +      fortranForceOptions.EnergyMixingRule = GEOMETRIC_MIXING_RULE;
94 +    }
95 +
96 +
97 +  }
98 +
99   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines