39 |
|
bool haveMaxLength; |
40 |
|
double maxLength; |
41 |
|
|
42 |
+ |
int nBins; |
43 |
+ |
|
44 |
|
vector<PairCorrList> theList; |
45 |
|
pairCorrEnum pairType; |
46 |
|
|
74 |
|
haveMaxLength = false; |
75 |
|
|
76 |
|
maxLength = 1.0; |
77 |
+ |
nBins = 100; |
78 |
|
|
79 |
|
programName = argv[0]; /*save the program name in case we need it*/ |
80 |
|
|
211 |
|
if( optionError ){ |
212 |
|
sprintf( painCave.errMsg, |
213 |
|
"\n" |
214 |
< |
"The -i flag should end an option sequence.\n" |
215 |
< |
" example: -ri <int> *NOT* -ir <int>\n" ); |
214 |
> |
"The -l flag should end an option sequence.\n" |
215 |
> |
" example: -sl <double> *NOT* -ls <double>\n" ); |
216 |
|
usage(); |
217 |
|
painCave.isFatal = 1; |
218 |
|
simError(); |
222 |
|
if( i>=argC ){ |
223 |
|
sprintf( painCave.errMsg, |
224 |
|
"\n" |
225 |
< |
"not enough arguments for -i\n"); |
225 |
> |
"not enough arguments for -l\n"); |
226 |
|
usage(); |
227 |
|
painCave.isFatal = 1; |
228 |
|
simError(); |
229 |
|
} |
230 |
|
|
231 |
< |
maxLength = strtod( argV[i], &conversionCheck, 10 ); |
231 |
> |
maxLength = strtod( argV[i], &conversionCheck); |
232 |
|
if( conversionCheck == argV[i] ) conversionError = true; |
233 |
|
if( *conversionCheck != '\0' ) conversionError = true; |
234 |
|
|
245 |
|
|
246 |
|
break; |
247 |
|
|
248 |
+ |
case 'n': |
249 |
+ |
// -n <int> set <int> to the nBins |
250 |
+ |
|
251 |
+ |
j++; |
252 |
+ |
currentFlag = argV[i][j]; |
253 |
+ |
|
254 |
+ |
if( currentFlag != '\0' ) optionError = true; |
255 |
+ |
|
256 |
+ |
if( optionError ){ |
257 |
+ |
sprintf( painCave.errMsg, |
258 |
+ |
"\n" |
259 |
+ |
"The -n flag should end an option sequence.\n" |
260 |
+ |
" example: -sn <int> *NOT* -ns <int>\n" ); |
261 |
+ |
usage(); |
262 |
+ |
painCave.isFatal = 1; |
263 |
+ |
simError(); |
264 |
+ |
} |
265 |
+ |
|
266 |
+ |
i++; |
267 |
+ |
if( i>=argC ){ |
268 |
+ |
sprintf( painCave.errMsg, |
269 |
+ |
"\n" |
270 |
+ |
"not enough arguments for -n\n"); |
271 |
+ |
usage(); |
272 |
+ |
painCave.isFatal = 1; |
273 |
+ |
simError(); |
274 |
+ |
} |
275 |
+ |
|
276 |
+ |
nBins = strtol( argV[i], &conversionCheck, 10 ); |
277 |
+ |
if( conversionCheck == argV[i] ) conversionError = true; |
278 |
+ |
if( *conversionCheck != '\0' ) conversionError = true; |
279 |
+ |
|
280 |
+ |
if( conversionError ){ |
281 |
+ |
sprintf( painCave.errMsg, |
282 |
+ |
"Error converting \"%s\" to an int for nBins.\n", |
283 |
+ |
argV[i] ); |
284 |
+ |
usage(); |
285 |
+ |
painCave.isFatal = 1; |
286 |
+ |
simError(); |
287 |
+ |
} |
288 |
+ |
|
289 |
+ |
if( nBins < 1 ){ |
290 |
+ |
sprintf( painCave.errMsg, |
291 |
+ |
"nBins error: nBins = %d is less than 1.\n", |
292 |
+ |
nBins ); |
293 |
+ |
usage(); |
294 |
+ |
painCave.isFatal = 1; |
295 |
+ |
simError(); |
296 |
+ |
} |
297 |
+ |
|
298 |
+ |
done = true; |
299 |
+ |
|
300 |
+ |
break; |
301 |
+ |
|
302 |
|
default: |
303 |
|
|
304 |
|
sprintf(painCave.errMsg, |
388 |
|
printf("Initializing the pair correlations..." ); |
389 |
|
fflush(stdout); |
390 |
|
|
391 |
+ |
if(haveMaxLength) |
392 |
+ |
theCorrs.setMaxLength( maxLength ); |
393 |
+ |
|
394 |
+ |
theCorrs.setNbins( nBins ); |
395 |
|
theCorrs.setFrames( infoArray, nFrames, reader ); |
396 |
|
theCorrs.setPairCorrList( theList ); |
397 |
|
theCorrs.initCorrelations( outPrefix ); |
419 |
|
" short:\n" |
420 |
|
" ------\n" |
421 |
|
" -o <name> The output prefix\n" |
422 |
+ |
" -n <nBins> Set the number of bins in the Histogram.\n" |
423 |
+ |
" *Defaults to 100\n" |
424 |
|
" -l <maxLength> set the maximum value of r\n" |
425 |
|
" *Defaults to 1/2 smallest length of first frame.\n" |
426 |
|
" -s Turn on separate output files\n" |