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

Comparing trunk/OOPSE-4/src/applications/staticProps/StaticProps.cpp (file contents):
Revision 3009 by gezelter, Wed Sep 20 20:13:40 2006 UTC vs.
Revision 3076 by chuckv, Fri Nov 10 18:45:29 2006 UTC

# Line 64 | Line 64
64   #if defined(HAVE_FFTW_H) || defined(HAVE_DFFTW_H) || defined(HAVE_FFTW3_H)
65   #include "applications/staticProps/Hxy.hpp"
66   #endif
67 + #include "applications/staticProps/RhoR.hpp"
68  
69   using namespace oopse;
70  
# Line 163 | Line 164 | int main(int argc, char* argv[]){
164  
165    //parse md file and set up the system
166    SimCreator creator;
167 +  std::cout << "dumpFile = " << dumpFileName << "\n";
168    SimInfo* info = creator.createSim(dumpFileName);
169  
170    RealType maxLen;
# Line 202 | Line 204 | int main(int argc, char* argv[]){
204    } else if (args_info.rp2_given){
205      analyser = new RippleOP(info, dumpFileName, sele1, sele2);
206    } else if (args_info.bo_given){
207 <    if (args_info.rcut_given && args_info.LegendreL_given) {
207 >    if (args_info.rcut_given) {
208        analyser = new BondOrderParameter(info, dumpFileName, sele1,
209                                          args_info.rcut_arg,
208 –                                        args_info.LegendreL_arg,
210                                          args_info.nbins_arg);
211      } else {
212        sprintf( painCave.errMsg,
213 <               "Both the cutoff radius (rcut) and LegendreL must be specified when calculating Bond Order Parameters");
213 >               "A cutoff radius (rcut) must be specified when calculating Bond Order Parameters");
214        painCave.severity = OOPSE_ERROR;
215        painCave.isFatal = 1;
216        simError();
# Line 233 | Line 234 | int main(int argc, char* argv[]){
234      analyser = new Hxy(info, dumpFileName, sele1, args_info.nbins_x_arg,
235                         args_info.nbins_y_arg, args_info.nbins_arg);
236   #endif
237 +  }else if (args_info.rho_r_given) {
238 +    if (args_info.radius_given){
239 +      analyser = new RhoR(info, dumpFileName, sele1, maxLen,args_info.nbins_arg,args_info.radius_arg);
240 +    }else{
241 +      sprintf( painCave.errMsg,
242 +               "A particle radius (radius) must be specified when calculating Rho(r)");
243 +      painCave.severity = OOPSE_ERROR;
244 +      painCave.isFatal = 1;
245 +      simError();
246 +    }
247    }
248    
249    if (args_info.output_given) {

Diff Legend

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