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 983 by gezelter, Tue Jun 6 17:43:28 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 +    std::cout << "GBM = " <<  fortranForceOptions.GayBerneMu << " GBN = " << fortranForceOptions.GayBerneNu  << std::endl;
73 +
74 +    std::string DistanceMix = this->getDistanceMixingRule();
75 +    toUpper(DistanceMix);
76 +    if(DistanceMix == "ARITHMETIC"){
77 +      fortranForceOptions.DistanceMixingRule = ARITHMETIC_MIXING_RULE;
78 +    }else if(DistanceMix == "GEOMETRIC"){
79 +      fortranForceOptions.DistanceMixingRule = GEOMETRIC_MIXING_RULE;
80 +    }else if(DistanceMix == "CUBIC"){
81 +      fortranForceOptions.DistanceMixingRule = CUBIC_MIXING_RULE;
82 +    }else{
83 +      fortranForceOptions.DistanceMixingRule = ARITHMETIC_MIXING_RULE;
84 +    }
85 +    std::string EnergyMix = this->getEnergyMixingRule();  
86 +    toUpper(EnergyMix);
87 +    if(EnergyMix == "ARITHMETIC"){
88 +      fortranForceOptions.EnergyMixingRule = ARITHMETIC_MIXING_RULE;
89 +    }else if(EnergyMix == "GEOMETRIC"){
90 +      fortranForceOptions.EnergyMixingRule = GEOMETRIC_MIXING_RULE;
91 +    }else if(EnergyMix == "HHG"){
92 +      fortranForceOptions.EnergyMixingRule = HHG_MIXING_RULE;
93 +    }else{
94 +      fortranForceOptions.EnergyMixingRule = GEOMETRIC_MIXING_RULE;
95 +    }
96 +
97 +
98 +  }
99 +
100   }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines