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

Comparing trunk/src/io/Globals.hpp (file contents):
Revision 609 by chuckv, Fri Sep 16 23:04:47 2005 UTC vs.
Revision 710 by chrisfen, Wed Nov 2 20:36:25 2005 UTC

# Line 53 | Line 53
53   #include "types/Component.hpp"
54   #include "types/MakeStamps.hpp"
55   #include "types/ZconStamp.hpp"
56 + #include "utils/CaseConversion.hpp"
57  
58 + template<typename T>
59 + struct ParameterTraits;
60  
61 + //string
62 + template<>                    
63 + struct ParameterTraits<std::string>{
64 +  typedef std::string RepType;       // Representation type of the value
65  
66 < /**
67 < * @class Globals Globals.hpp "io/Globals.hpp"
68 < * @brief parsing and storing global parameters for simulation
69 < * @todo need refactorying
70 < */
71 < class Globals{
72 <  
73 < public:
74 <  
75 <  Globals();
76 <  ~Globals();
66 > template<typename T> static bool    convert(T v, RepType& r){return false;} // !NB everything is ok
67 > template<typename T> static RepType convert(T v)            {RepType tmp; convert(v,tmp);return tmp;}
68 >  static bool convert(RepType v, RepType& r) { r = v; return true;}
69 >  static std::string getParamType() { return "string";}  
70 > };
71 > //bool
72 > template<>                    
73 > struct ParameterTraits<bool>{
74 >  typedef bool RepType;
75 >  template<typename T> static bool    convert(T, RepType&){return false;}
76 >  template<typename T> static RepType convert(T v)        {RepType tmp; convert(v,tmp);return tmp;}
77 >  static bool convert(std::string v, RepType& r) {
78 >    oopse::toLower(v);
79 >    bool result = false;
80 >    if (v == "true") {
81 >        r = true;
82 >        result = true;
83 >    } else if (v == "false") {
84 >        r = false;
85 >        result = true;
86 >    }
87 >    
88 >     return result;
89 >  }
90 >  static std::string getParamType() { return "bool";}  
91 > };
92  
93 <  void initalize();
94 <  
95 <  int newComponent( event* the_event );
96 <  int componentAssign( event* the_event );
97 <  int componentEnd( event* the_event );
93 > //int  
94 > template<>
95 > struct ParameterTraits<int>{
96 >    typedef int RepType;
97 >    template<typename T> static bool    convert(T, RepType&){return false;}
98 >    template<typename T> static RepType convert(T v)        {RepType tmp; convert(v,tmp);return tmp;}
99 >    static bool convert(RepType v, RepType& r)            { r=v; return true;}
100 >    static std::string getParamType() { return "int";}  
101 > };
102  
103 <  int newZconstraint( event* the_event );
104 <  int zConstraintAssign( event* the_event );
105 <  int zConstraintEnd( event* the_event );
106 <  
107 <  int globalAssign( event* the_event );
108 <  int globalEnd( event* the_event );
109 <
110 <  char*  getForceField( void )      { return force_field; }
111 <  int    getNComponents( void )     { return n_components; }
112 <  double getTargetTemp( void )      { return target_temp; }
87 <  double getTargetPressure( void )  { return target_pressure; }
88 <  double getQmass( void )           { return q_mass; }
89 <  double getTauThermostat( void )   { return tau_thermostat; }
90 <  double getTauBarostat( void )     { return tau_barostat; }
91 <  char*  getEnsemble( void )        { return ensemble; }
92 <  double getDt( void )              { return dt; }
93 <  double getRunTime( void )         { return run_time; }
94 <
95 <  int    getNzConstraints( void )   { return n_zConstraints; }
96 <  char*  getInitialConfig( void )   { return initial_config; }
97 <  char*  getFinalConfig( void )     { return final_config; }
98 <  int    getNMol( void )            { return n_mol; }
99 <  double getDensity( void )         { return density; }
100 <  double getBox( void )             { return box; }
101 <  double getBoxX( void )            { return box_x; }
102 <  double getBoxY( void )            { return box_y; }
103 <  double getBoxZ( void )            { return box_z; }
104 <  double getSampleTime( void )      { return sample_time; }
105 <  double getStatusTime( void )      { return status_time; }
106 <  double getResetTime( void )       { return resetTime; }
107 <  double getThermalTime( void )     { return thermal_time; }
108 <  double getDielectric( void )      { return dielectric; }
109 <  double getRcut( void )            { return rcut; }
110 <  double getRsw( void )             { return rsw; }
111 <  int    getTempSet( void )         { return tempSet; }
112 <  int    getUseInitTime( void )     { return useInitTime; }
113 <  int    getUseInitXSstate( void )  { return useInitXSstate; }
114 <  double getOrthoBoxTolerance(void) { return orthoBoxTolerance; }
115 <  int    getPBC( void )             { return usePBC; }
116 <  char*  getMixingRule( void)       { return mixingRule; }
117 <  double getZconsTime(void)         { return zcons_time; }
118 <  double getZconsTol(void)          { return zcons_tol; }
119 <  char*  getZconsForcePolicy(void)  { return zconsForcePolicy; }
120 <  double getZconsGap(void)          { return zcons_gap; }
121 <  double getZconsFixtime(void)      { return zcons_fixtime; }
122 <  int    getZconsUsingSMD(void)     { return zcons_using_smd; }
123 <  int    getSeed(void)              { return seed; }
124 <  char*  getMinimizer(void)         { return minimizer_name; }
125 <  int    getMinMaxIter(void)        { return minimizer_maxiteration; }
126 <  int    getMinWriteFrq(void)       { return minimizer_writefrq; }
127 <  double getMinStepSize(void)       { return minimizer_stepsize; }
128 <  double getMinFTol(void)           { return minimizer_ftol; }
129 <  double getMinGTol(void)           { return minimizer_gtol; }
130 <  double getMinLSTol(void)          { return minimizer_ls_tol; }
131 <  int    getMinLSMaxIter(void)      { return minimizer_ls_maxiteration; }
132 <  int    getUseSolidThermInt(void)  { return useSolidThermInt; }
133 <  int    getUseLiquidThermInt(void) { return useLiquidThermInt; }
134 <  double getThermIntLambda(void)    { return thermodynamic_integration_lambda; }
135 <  double getThermIntK(void)         { return thermodynamic_integration_k; }
136 <  char*  getForceFieldVariant( void ) { return forcefield_variant; }
137 <  char*  getForceFieldFileName()    { return forcefield_filename; }
138 <  double getDistSpringConst(void)   { return therm_int_dist_spring; }
139 <  double getThetaSpringConst(void)  { return therm_int_theta_spring; }
140 <  double getOmegaSpringConst(void)  { return therm_int_omega_spring; }
141 <  double getSurfaceTension(void)    { return surface_tension; }
142 <  bool   getPrintPressureTensor(void) { return print_pressure_tensor;}
143 <  char*  getElectrostaticSummationMethod(void) { return electrostaticSummationMethod; }
144 <  double getDampingAlpha(void)      { return dampingAlpha; }
145 <  char*  getCutoffPolicy(void)      { return cutoffPolicy; }
146 <  bool   getCompressDumpFile(void)  { return compressDumpFile;}
147 <  
148 <  short int haveDt( void )            { return have_dt; }
149 <  short int haveRunTime( void )       { return have_run_time; }
150 <  short int haveEnsemble( void )      { return have_ensemble; }
151 <  short int haveTargetTemp( void )    { return have_target_temp; }
152 <  short int haveInitialConfig( void ) { return have_initial_config; }
153 <  short int haveFinalConfig( void )   { return have_final_config; }
154 <  short int haveNMol( void )          { return have_n_mol; }
155 <  short int haveDensity( void )       { return have_density; }
156 <  short int haveBox( void )           { return have_box; }
157 <  short int haveBoxX( void )          { return have_box_x; }
158 <  short int haveBoxY( void )          { return have_box_y; }
159 <  short int haveBoxZ( void )          { return have_box_z; }
160 <  short int haveSampleTime( void )    { return have_sample_time; }
161 <  short int haveResetTime( void )     { return have_reset_time; }
162 <  short int haveStatusTime( void )    { return have_status_time; }
163 <  short int haveThermalTime( void )   { return have_thermal_time; }
164 <  short int haveRcut( void )          { return have_rcut; }
165 <  short int haveRsw( void )           { return have_rsw; }
166 <  short int haveDielectric( void )    { return have_dielectric; }
167 <  short int haveTempSet( void )       { return have_tempSet; }
168 <  short int haveTargetPressure( void ){ return have_target_pressure; }
169 <  short int haveQmass( void )         { return have_q_mass; }
170 <  short int haveTauThermostat( void ) { return have_tau_thermostat; }
171 <  short int haveTauBarostat( void )   { return have_tau_barostat; }
172 <  short int haveZconstraintTime(void) { return have_zcons_time; }
173 <  short int haveZconstraints( void )  { return have_zConstraints; }
174 <  short int haveZconsTol(void)        { return have_zcons_tol; }
175 <  short int haveZconsForcePolicy(void){ return have_zcons_force_policy; }
176 <  short int haveZConsGap(void)        { return have_zcons_gap; }
177 <  short int haveZConsFixTime(void)    { return have_zcons_fixtime; }
178 <  short int haveZConsUsingSMD(void)   { return have_zcons_using_smd; }  
179 <  short int haveSeed(void)            { return have_seed; }
180 <  short int haveMinimizer(void)       { return have_minimizer; }
181 <  short int haveMinMaxIter(void)      { return have_minimizer_maxiteration; }
182 <  short int haveMinWriteFrq(void)     { return have_minimizer_writefrq; }
183 <  short int haveMinStepSize(void)     { return have_minimizer_stepsize; }
184 <  short int haveMinFTol(void)         { return have_minimizer_ftol; }
185 <  short int haveMinGTol(void)         { return have_minimizer_gtol; }
186 <  short int haveMinLSTol(void)        { return have_minimizer_ls_tol; }
187 <  short int haveMinLSMaxIter(void)    { return have_minimizer_ls_maxiteration;}
188 <  short int haveThermIntLambda(void)  { return have_thermodynamic_integration_lambda; }
189 <  short int haveThermIntK(void)       { return have_thermodynamic_integration_k; }
190 <  short int haveForceFieldVariant(void)  { return have_forcefield_variant; }
191 <  short int haveForceFieldFileName(void) { return have_forcefield_filename; }
192 <  short int haveDistSpringConst(void)    { return have_dist_spring_constant; }
193 <  short int haveThetaSpringConst(void)   { return have_theta_spring_constant; }
194 <  short int haveOmegaSpringConst(void)   { return have_omega_spring_constant; }
195 <  short int haveSurfaceTension(void)     { return have_surface_tension; }
196 <  short int havePrintPressureTensor(void) {return have_print_pressure_tensor;}
197 <  short int haveElectrostaticSummationMethod(void) {return have_electro_sum_method;}
198 <  short int haveDampingAlpha(void)       { return have_damping_alpha; }
199 <  short int haveCutoffPolicy(void)       { return have_cutoff_policy; }
200 <  short int haveCompressDumpfile(void)   { return have_compress_dumpfile; }
103 > //double
104 > template<>                    
105 > struct ParameterTraits<double>{
106 >    typedef double RepType;
107 >    template<typename T> static bool    convert(T, RepType&){return false;}
108 >    template<typename T> static RepType convert(T v)        {RepType tmp; convert(v,tmp);return tmp;}
109 >    static bool convert(RepType v, RepType& r)            {r=v; return true;}
110 >    static bool convert(int v, RepType& r)                {r = static_cast<double>(v); return true;}
111 >    static std::string getParamType() { return "double";}    
112 > };
113  
202 –  /* other accessors */
203 –  Component** getComponents( void )   { return components; }
204 –  ZconStamp** getZconStamp( void )    { return zConstraints; }
205 –  
206 – private:
207 –  
114  
115 <  typedef std::map<std::string, int> CommandMapType;
116 <  CommandMapType command_table;
115 > class ParameterBase {
116 >  public:    
117 >    ParameterBase() : keyword_(), optional_(false), defaultValue_(false), empty_(true) {}
118 >    virtual ~ParameterBase() {}
119 >    bool isOptional() {return optional_;}
120 >    void setOptional(bool optional) {optional_ = optional;}
121 >    bool hasDefaultValue() {return defaultValue_;}
122 >    virtual bool isValid() { return true;}
123 >    const std::string& getKeyword() {return keyword_;}
124 >    void setKeyword(const std::string& keyword) { keyword_ = keyword;}
125 >    bool empty() {return empty_;}
126 >    virtual bool setData(std::string) = 0;
127 >    virtual bool setData(int) = 0;
128 >    virtual bool setData(double) = 0;
129 >    virtual std::string getParamType() = 0;
130 >  protected:
131 >    std::string keyword_;
132 >    bool optional_;
133 >    bool defaultValue_;
134 >    bool empty_;
135 > };
136  
137 <  
138 <  char* checkMe( void );
139 <  
140 <  Component* current_component;
141 <  Component** components; // the array of components
142 <
218 <  ZconStamp* current_zConstraint;
219 <  ZconStamp** zConstraints; // the array of zConstraints
137 > template<class ParamType>
138 > class Parameter : public ParameterBase{
139 >   public:    
140 >    typedef ParameterTraits<ParamType> ValueType;
141 >     void setDefaultValue(const ParamType& value) {data_ = value; defaultValue_ = true;}
142 >     ParamType getData() { return data_;}
143  
144 <  char force_field[100];
145 <  int n_components;
146 <  int n_zConstraints;
147 <  double target_temp;
148 <  double target_pressure;
149 <  char ensemble[100];
150 <  char mixingRule[100];
151 <  double dt;
152 <  double run_time;
153 <  char initial_config[120];
154 <  char final_config[120];
155 <  int n_mol;
156 <  double density;
157 <  double box;
158 <  double box_x, box_y, box_z;
159 <  double sample_time;
160 <  double status_time;
161 <  double resetTime;
162 <  double orthoBoxTolerance;
163 <  double thermal_time;
164 <  double rcut;
165 <  double rsw;
166 <  double dielectric;
167 <  int tempSet;
168 <  int useInitTime;
169 <  int useInitXSstate;
247 <  int usePBC;
248 <  double q_mass;
249 <  double tau_thermostat;
250 <  double tau_barostat;
251 <  double zcons_time;    
252 <  double zcons_tol;
253 <  char zconsForcePolicy[100];
254 <  double zcons_gap;
255 <  double zcons_fixtime;
256 <  int zcons_using_smd;
257 <  
258 <  int seed;
259 <  char minimizer_name[100];
260 <  int minimizer_maxiteration;
261 <  int minimizer_writefrq;
262 <  double minimizer_stepsize;
263 <  double minimizer_ftol;
264 <  double minimizer_gtol;
265 <  double minimizer_ls_tol;
266 <  int minimizer_ls_maxiteration;
267 <  int useSolidThermInt;
268 <  int useLiquidThermInt;
269 <  double thermodynamic_integration_lambda;
270 <  double thermodynamic_integration_k;
271 <  char forcefield_variant[100];
272 <  char forcefield_filename[100];
273 <  double therm_int_dist_spring;
274 <  double therm_int_theta_spring;
275 <  double therm_int_omega_spring;
276 <  double surface_tension;
277 <  bool print_pressure_tensor;
278 <  char electrostaticSummationMethod[100];
279 <  double dampingAlpha;
280 <  char cutoffPolicy[100];
281 <  bool compressDumpFile;
282 <  
283 <  //required arguments
284 <  short int have_force_field, have_n_components, have_target_temp;
285 <  short int have_target_pressure, have_ensemble, have_dt, have_run_time;
286 <  
287 <  // optional arguments
288 <  short int have_initial_config, have_final_config, have_n_mol;
289 <  short int have_density, have_box, have_box_x, have_box_y, have_box_z;
290 <  short int have_sample_time, have_status_time, have_rcut, have_dielectric;
291 <  short int have_tempSet, have_thermal_time, have_rsw, have_q_mass;
292 <  short int have_tau_thermostat, have_tau_barostat;
293 <  short int have_zcons_time, have_zConstraints, have_n_zConstraints;
294 <  short int have_zcons_tol, have_seed;
295 <  short int have_zcons_force_policy, have_reset_time;
296 <  short int have_zcons_gap, have_zcons_fixtime;
297 <  short int have_zcons_using_smd;
298 <  short int have_minimizer, have_minimizer_maxiteration;
299 <  short int have_minimizer_writefrq, have_minimizer_stepsize;
300 <  short int have_minimizer_ftol, have_minimizer_gtol;
301 <  short int have_minimizer_ls_tol, have_minimizer_ls_maxiteration;
302 <  short int have_thermodynamic_integration_lambda;
303 <  short int have_thermodynamic_integration_k;
304 <  short int have_forcefield_variant;
305 <  short int have_forcefield_filename;  
306 <  short int have_dist_spring_constant;
307 <  short int have_theta_spring_constant;
308 <  short int have_omega_spring_constant;
309 <  short int have_surface_tension;
310 <  short int have_print_pressure_tensor;
311 <  short int have_electro_sum_method;
312 <  short int have_damping_alpha;
313 <  short int have_cutoff_policy;
314 <  short int have_compress_dumpfile;
144 >    virtual bool setData(std::string sval) {
145 >        return internalSetData<std::string>(sval);
146 >    }
147 >    virtual bool setData(int ival) {
148 >        return internalSetData<int>(ival);
149 >    }
150 >    
151 >    virtual bool setData(double dval) {
152 >        return internalSetData<double>(dval);
153 >    }
154 >
155 >    virtual std::string getParamType() { return ParameterTraits<ParamType>::getParamType();}
156 >   private:
157 >     template<class T> bool internalSetData(T data) {
158 >        ParamType tmp;
159 >        bool result = ValueType::convert(data, tmp);
160 >        if (result) {
161 >            empty_ = false;
162 >            data_ = tmp;
163 >        }
164 >        return result;
165 >     }
166 >    
167 >   private:
168 >     ParamType data_;
169 >      
170   };
171  
172 + #define DeclareParameter(NAME, TYPE)         \
173 +    private:                                                   \
174 +      Parameter<TYPE> NAME;                                     \
175 +    public:                                                      \
176 +      bool have##NAME() { return !NAME.empty();}  \
177 +      TYPE get##NAME() { return NAME.getData();}
178 +
179 +
180 + class Globals {
181 +  public:
182 +    Globals();
183 +    
184 +  DeclareParameter(ForceField, std::string);
185 +  DeclareParameter(NComponents, int);  
186 +  DeclareParameter(TargetTemp, double);
187 +  DeclareParameter(Ensemble, std::string);
188 +  DeclareParameter(Dt, double);
189 +  DeclareParameter(RunTime, double);
190 +  DeclareParameter(InitialConfig, std::string);
191 +  DeclareParameter(FinalConfig, std::string);
192 +  DeclareParameter(NMol, int);
193 +  DeclareParameter(Density, double);
194 +  DeclareParameter(Box, double);
195 +  DeclareParameter(BoxX, double);
196 +  DeclareParameter(BoxY, double);
197 +  DeclareParameter(BoxZ, double);
198 +  DeclareParameter(SampleTime, double);
199 +  DeclareParameter(ResetTime, double);
200 +  DeclareParameter(StatusTime, double);
201 +  DeclareParameter(CutoffRadius, double);
202 +  DeclareParameter(SwitchingRadius, double);
203 +  DeclareParameter(Dielectric, double);
204 +  DeclareParameter(TempSet, bool);
205 +  DeclareParameter(ThermalTime, double);
206 +  DeclareParameter(MixingRule, std::string);
207 +  DeclareParameter(UsePeriodicBoundaryConditions, bool);
208 +  DeclareParameter(TargetPressure, double);
209 +  DeclareParameter(TauThermostat, double);
210 +  DeclareParameter(TauBarostat, double);
211 +  DeclareParameter(ZconsTime, double);
212 +  DeclareParameter(NZconstraints, int);
213 +  DeclareParameter(ZconsTol, double);
214 +  DeclareParameter(ZconsForcePolicy, std::string);
215 +  DeclareParameter(Seed, int);
216 +  DeclareParameter(UseInitalTime, bool);
217 +  DeclareParameter(UseIntialExtendedSystemState, bool);
218 +  DeclareParameter(OrthoBoxTolerance, double);
219 +  DeclareParameter(Minimizer, std::string);
220 +  DeclareParameter(MinimizerMaxIter, double);
221 +  DeclareParameter(MinimizerWriteFrq, int);
222 +  DeclareParameter(MinimizerStepSize, double);
223 +  DeclareParameter(MinimizerFTol, double);
224 +  DeclareParameter(MinimizerGTol, double);
225 +  DeclareParameter(MinimizerLSTol, double);
226 +  DeclareParameter(MinimizerLSMaxIter, int);
227 +  DeclareParameter(ZconsGap, double);
228 +  DeclareParameter(ZconsFixtime, double);
229 +  DeclareParameter(ZconsUsingSMD, bool);
230 +  DeclareParameter(UseSolidThermInt, bool);
231 +  DeclareParameter(UseLiquidThermInt, bool);
232 +  DeclareParameter(ThermodynamicIntegrationLambda, double);
233 +  DeclareParameter(ThermodynamicIntegrationK, double);
234 +  DeclareParameter(ForceFieldVariant, std::string);
235 +  DeclareParameter(ForceFieldFileName, std::string);
236 +  DeclareParameter(ThermIntDistSpringConst, double);
237 +  DeclareParameter(ThermIntThetaSpringConst, double);
238 +  DeclareParameter(ThermIntOmegaSpringConst, double);
239 +  DeclareParameter(SurfaceTension, double);
240 +  DeclareParameter(PrintPressureTensor, bool);
241 +  DeclareParameter(ElectrostaticSummationMethod, std::string);
242 +  DeclareParameter(ScreeningMethod, std::string);
243 +  DeclareParameter(DampingAlpha, double);
244 +  DeclareParameter(CutoffPolicy, std::string);
245 +  DeclareParameter(CompressDumpFile, bool);
246 +  DeclareParameter(SkinThickness, double);
247 +  DeclareParameter(StatFileFormat, std::string);    
248 +
249 +  private:
250 +    typedef std::map<std::string, ParameterBase*> ParamMap;
251 +    ParamMap parameters_;
252 +    
253 +    Component* current_component;
254 +    Component** components; // the array of components
255 +
256 +    ZconStamp* current_zConstraint;
257 +    ZconStamp** zConstraints; // the array of zConstraints
258 +
259 +    char* checkMe();
260 +
261 +  public:
262 +    int newComponent( event* the_event );
263 +    int componentAssign( event* the_event );
264 +    int componentEnd( event* the_event );
265 +
266 +    int newZconstraint( event* the_event );
267 +    int zConstraintAssign( event* the_event );
268 +    int zConstraintEnd( event* the_event );
269 +  
270 +    int globalAssign( event* the_event );
271 +    int globalEnd( event* the_event );    
272 +
273 +    ZconStamp** getZconStamp() {return zConstraints;}
274 +    Component** getComponents() {return components;}
275 + };
276   #endif
277 +

Diff Legend

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