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

Comparing trunk/src/applications/staticProps/StaticProps.cpp (file contents):
Revision 1513 by gezelter, Tue Oct 19 18:40:54 2010 UTC vs.
Revision 1585 by gezelter, Sun Jun 19 17:51:15 2011 UTC

# Line 69 | Line 69
69   #include "applications/staticProps/pAngle.hpp"
70   #include "applications/staticProps/BondAngleDistribution.hpp"
71   #include "applications/staticProps/NanoVolume.hpp"
72 + #include "applications/staticProps/NanoLength.hpp"
73   #if defined(HAVE_FFTW_H) || defined(HAVE_DFFTW_H) || defined(HAVE_FFTW3_H)
74   #include "applications/staticProps/Hxy.hpp"
75   #endif
76   #include "applications/staticProps/RhoR.hpp"
77 + #include "applications/staticProps/AngleR.hpp"
78 + #include "applications/staticProps/RhoAngleR.hpp"
79 + #include "applications/staticProps/TetrahedralityParam.hpp"
80  
81   using namespace OpenMD;
82  
# Line 187 | Line 191 | int main(int argc, char* argv[]){
191      maxLen = std::min(std::min(hmat(0, 0), hmat(1, 1)), hmat(2, 2)) /2.0;
192      zmaxLen = hmat(2,2);    
193    }    
194 <
194 >  
195    StaticAnalyser* analyser;
196    if (args_info.gofr_given){
197      analyser= new GofR(info, dumpFileName, sele1, sele2, maxLen,
198 <                       args_info.nbins_arg);        
198 >                       args_info.nbins_arg);        
199    } else if (args_info.gofz_given) {
200      analyser= new GofZ(info, dumpFileName, sele1, sele2, maxLen,
201 <                       args_info.nbins_arg);
201 >                       args_info.nbins_arg);
202    } else if (args_info.r_z_given) {
203      analyser  = new GofRZ(info, dumpFileName, sele1, sele2, maxLen, zmaxLen,
204 <                          args_info.nbins_arg, args_info.nbins_z_arg);
204 >                          args_info.nbins_arg, args_info.nbins_z_arg);
205    } else if (args_info.r_theta_given) {
206      analyser  = new GofRTheta(info, dumpFileName, sele1, sele2, maxLen,
207 <                              args_info.nbins_arg, args_info.nanglebins_arg);
207 >                              args_info.nbins_arg, args_info.nanglebins_arg);
208    } else if (args_info.r_omega_given) {
209      analyser  = new GofROmega(info, dumpFileName, sele1, sele2, maxLen,
210 <                              args_info.nbins_arg, args_info.nanglebins_arg);
210 >                              args_info.nbins_arg, args_info.nanglebins_arg);
211    } else if (args_info.theta_omega_given) {
212      analyser  = new GofAngle2(info, dumpFileName, sele1, sele2,
213 <                              args_info.nanglebins_arg);
213 >                              args_info.nanglebins_arg);
214    } else if (args_info.gxyz_given) {
215      if (args_info.refsele_given) {
216        analyser= new GofXyz(info, dumpFileName, sele1, sele2,args_info.refsele_arg,
217 <                           maxLen, args_info.nbins_arg);        
217 >                           maxLen, args_info.nbins_arg);        
218      } else {
219        sprintf( painCave.errMsg,
220 <               "--refsele must set when --gxyz is used");
220 >               "--refsele must set when --gxyz is used");
221        painCave.severity = OPENMD_ERROR;
222        painCave.isFatal = 1;
223        simError();  
# Line 221 | Line 225 | int main(int argc, char* argv[]){
225    } else if (args_info.twodgofr_given){
226      if (args_info.dz_given) {
227        analyser= new TwoDGofR(info, dumpFileName, sele1, sele2, maxLen,
228 <                             args_info.dz_arg, args_info.nbins_arg);        
228 >                             args_info.dz_arg, args_info.nbins_arg);        
229      } else {
230        sprintf( painCave.errMsg,
231 <               "A slab width (dz) must be specified when calculating TwoDGofR");
231 >               "A slab width (dz) must be specified when calculating TwoDGofR");
232        painCave.severity = OPENMD_ERROR;
233        painCave.isFatal = 1;
234        simError();
235 <    }
235 >    }    
236    } else if (args_info.p2_given) {
237 <    analyser  = new P2OrderParameter(info, dumpFileName, sele1, sele2);
237 >    if (args_info.sele1_given) {    
238 >      if (args_info.sele2_given)
239 >        analyser  = new P2OrderParameter(info, dumpFileName, sele1, sele2);
240 >      else
241 >        analyser  = new P2OrderParameter(info, dumpFileName, sele1);
242 >    } else {
243 >      sprintf( painCave.errMsg,
244 >               "At least one selection script (--sele1) must be specified when calculating P2 order parameters");
245 >      painCave.severity = OPENMD_ERROR;
246 >      painCave.isFatal = 1;
247 >      simError();
248 >    }
249    } else if (args_info.rp2_given){
250      analyser = new RippleOP(info, dumpFileName, sele1, sele2);
251    } else if (args_info.bo_given){
252      if (args_info.rcut_given) {
253        analyser = new BondOrderParameter(info, dumpFileName, sele1,
254 <                                        args_info.rcut_arg,
255 <                                        args_info.nbins_arg);
254 >                                        args_info.rcut_arg,
255 >                                        args_info.nbins_arg);
256      } else {
257        sprintf( painCave.errMsg,
258 <               "A cutoff radius (rcut) must be specified when calculating Bond Order Parameters");
258 >               "A cutoff radius (rcut) must be specified when calculating Bond Order Parameters");
259        painCave.severity = OPENMD_ERROR;
260        painCave.isFatal = 1;
261        simError();
262      }
263 +    
264 +  } else if (args_info.tet_param_given) {
265 +    if (args_info.rcut_given) {  
266 +      analyser = new TetrahedralityParam(info, dumpFileName, sele1,
267 +                                         args_info.rcut_arg,
268 +                                         args_info.nbins_arg);
269 +    } else {
270 +      sprintf( painCave.errMsg,
271 +               "A cutoff radius (rcut) must be specified when calculating Tetrahedrality Parameters");
272 +      painCave.severity = OPENMD_ERROR;
273 +      painCave.isFatal = 1;
274 +      simError();
275 +    }
276    } else if (args_info.bor_given){
277      if (args_info.rcut_given) {
278        analyser = new BOPofR(info, dumpFileName, sele1, args_info.rcut_arg,
279 <                            args_info.nbins_arg, maxLen);
279 >                            args_info.nbins_arg, maxLen);
280      } else {
281        sprintf( painCave.errMsg,
282 <               "A cutoff radius (rcut) must be specified when calculating Bond Order Parameters");
282 >               "A cutoff radius (rcut) must be specified when calculating Bond Order Parameters");
283        painCave.severity = OPENMD_ERROR;
284        painCave.isFatal = 1;
285        simError();
# Line 259 | Line 287 | int main(int argc, char* argv[]){
287    } else if (args_info.bad_given){
288      if (args_info.rcut_given) {
289        analyser = new BondAngleDistribution(info, dumpFileName, sele1, args_info.rcut_arg,
290 <                                           args_info.nbins_arg);
290 >                                           args_info.nbins_arg);
291      } else {
292        sprintf( painCave.errMsg,
293 <               "A cutoff radius (rcut) must be specified when calculating Bond Angle Distributions");
293 >               "A cutoff radius (rcut) must be specified when calculating Bond Angle Distributions");
294        painCave.severity = OPENMD_ERROR;
295        painCave.isFatal = 1;
296        simError();
297 <    }
297 >      }
298    } else if (args_info.scd_given) {
299      if (batchMode) {
300        analyser  = new SCDOrderParameter(info, dumpFileName, args_info.molname_arg,
301 <                                        args_info.begin_arg, args_info.end_arg);
301 >                                        args_info.begin_arg, args_info.end_arg);
302      } else{
303        std::string sele3 = args_info.sele3_arg;
304        analyser  = new SCDOrderParameter(info, dumpFileName, sele1, sele2, sele3);
305      }
306    }else if (args_info.density_given) {
307      analyser= new DensityPlot(info, dumpFileName, sele1, sele2, maxLen,
308 <                              args_info.nbins_arg);  
308 >                              args_info.nbins_arg);  
309    } else if (args_info.count_given) {
310      analyser = new ObjectCount(info, dumpFileName, sele1 );
311    } else if (args_info.slab_density_given) {
# Line 287 | Line 315 | int main(int argc, char* argv[]){
315   #if defined(HAVE_FFTW_H) || defined(HAVE_DFFTW_H) || defined(HAVE_FFTW3_H)
316    }else if (args_info.hxy_given) {
317      analyser = new Hxy(info, dumpFileName, sele1, args_info.nbins_x_arg,
318 <                       args_info.nbins_y_arg, args_info.nbins_arg);
318 >                       args_info.nbins_y_arg, args_info.nbins_arg);
319   #endif
320    }else if (args_info.rho_r_given) {
321      if (args_info.radius_given){
322        analyser = new RhoR(info, dumpFileName, sele1, maxLen,args_info.nbins_arg,args_info.radius_arg);
323      }else{
324        sprintf( painCave.errMsg,
325 <               "A particle radius (radius) must be specified when calculating Rho(r)");
325 >               "A particle radius (radius) must be specified when calculating Rho(r)");
326        painCave.severity = OPENMD_ERROR;
327        painCave.isFatal = 1;
328        simError();
329      }
330 <  }else if (args_info.hullvol_given) {
330 >  } else if (args_info.hullvol_given) {
331      analyser = new NanoVolume(info, dumpFileName, sele1);
332 +  } else if (args_info.rodlength_given) {
333 +    analyser = new NanoLength(info, dumpFileName, sele1);
334 +  } else if (args_info.angle_r_given) {
335 +    analyser = new AngleR(info, dumpFileName, sele1, maxLen,args_info.nbins_arg);
336    }
337 <  
337 >    
338    if (args_info.output_given) {
339      analyser->setOutputName(args_info.output_arg);
340    }
341    if (args_info.step_given) {
342      analyser->setStep(args_info.step_arg);
343    }
344 <
344 >  
345    analyser->process();
346 <
346 >  
347    delete analyser;    
348    delete info;
349  

Diff Legend

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