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

Comparing trunk/src/io/ForceFieldOptions.hpp (file contents):
Revision 820 by chuckv, Fri Dec 23 23:04:25 2005 UTC vs.
Revision 1150 by chuckv, Fri Jul 6 18:14:35 2007 UTC

# Line 45 | Line 45
45   #include "utils/ParameterManager.hpp"
46   #include "utils/StringUtils.hpp"
47   #include "io/ParamConstraint.hpp"
48 + #define __C
49 + #include "UseTheForce/fForceOptions.h"
50  
51   namespace oopse {
52    
# Line 54 | Line 56 | namespace oopse {
56      DeclareParameter(DistanceMixingRule, std::string);
57      DeclareParameter(DistanceType, std::string);
58      DeclareParameter(EnergyMixingRule, std::string);
59 <    DeclareParameter(EnergyUnitScaling, double);
60 <    DeclareParameter(DistanceUnitScaling, double);
61 <    DeclareParameter(AngleUnitScaling, double);
59 >    DeclareParameter(CutoffPolicy, std::string);
60 >    DeclareParameter(EnergyUnitScaling, RealType);
61 >    DeclareParameter(MetallicEnergyUnitScaling, RealType);
62 >    DeclareParameter(DistanceUnitScaling, RealType);
63 >    DeclareParameter(AngleUnitScaling, RealType);
64      DeclareParameter(TorsionAngleConvention, std::string);
65 <    DeclareParameter(vdw14scale, double);
66 <    DeclareParameter(electrostatic14scale, double);
67 <    DeclareParameter(dielectric, double);
65 >    DeclareParameter(vdw14scale, RealType);
66 >    DeclareParameter(electrostatic14scale, RealType);
67 >    DeclareParameter(GayBerneMu, RealType);
68 >    DeclareParameter(GayBerneNu, RealType);
69      
70    public:
71      ForceFieldOptions();
# Line 69 | Line 74 | namespace oopse {
74      
75      void validateOptions() {
76        CheckParameter(vdWtype, isEqualIgnoreCase(std::string("Lennard-Jones")));
77 <      CheckParameter(DistanceMixingRule, isEqualIgnoreCase(std::string("arithmetic")) || isEqualIgnoreCase(std::string("geometric")));
77 >      CheckParameter(DistanceMixingRule, isEqualIgnoreCase(std::string("arithmetic")) || isEqualIgnoreCase(std::string("geometric")) || isEqualIgnoreCase(std::string("cubic")));
78        CheckParameter(DistanceType, isEqualIgnoreCase(std::string("sigma")) || isEqualIgnoreCase(std::string("Rmin")));
79 <      CheckParameter(EnergyMixingRule, isEqualIgnoreCase(std::string("arithmetic")) || isEqualIgnoreCase(std::string("geometric")));
79 >      CheckParameter(EnergyMixingRule, isEqualIgnoreCase(std::string("arithmetic")) || isEqualIgnoreCase(std::string("geometric")) || isEqualIgnoreCase(std::string("hhg")));
80        CheckParameter(TorsionAngleConvention, isEqualIgnoreCase(std::string("180 is trans")) || isEqualIgnoreCase(std::string("0 is trans")));
81 <    }
81 >      CheckParameter(CutoffPolicy, isEqualIgnoreCase(std::string("MIX")) || isEqualIgnoreCase(std::string("MAX")) || isEqualIgnoreCase(std::string("TRADITIONAL")));
82 >   }
83      
84      bool setData(const std::string& keyword, const std::string& value) {
85        bool result;
86        ParamMap::iterator i =parameters_.find(keyword);
87        if (i != parameters_.end()) {
88 <        if(isType<int>(value)){
88 >        if(isInteger(value)){
89            int ival = lexi_cast<int>(value);
90            result = i->second->setData(ival);
91          }      
92 <        else if (isType<double>(value)){
93 <          double dval = lexi_cast<double>(value);
92 >        else if (isType<RealType>(value)){
93 >          RealType dval = lexi_cast<RealType>(value);
94            result = i->second->setData(dval);
95          } else{
96            result = i->second->setData(value);
# Line 97 | Line 103 | namespace oopse {
103        
104        return result;
105      }
106 <
107 <    void tellFortranOptions();
106 >
107 >    void makeFortranOptions(ForceOptions & fortranForceOptions);
108    private:
109      typedef std::map<std::string, ParameterBase*> ParamMap;
110      ParamMap parameters_;                  

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines