# | 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) { |
– | Removed lines |
+ | Added lines |
< | Changed lines |
> | Changed lines |