queso-0.51.1
SequenceStatisticalOptions.C
Go to the documentation of this file.
1 //-----------------------------------------------------------------------bl-
2 //--------------------------------------------------------------------------
3 //
4 // QUESO - a library to support the Quantification of Uncertainty
5 // for Estimation, Simulation and Optimization
6 //
7 // Copyright (C) 2008,2009,2010,2011,2012,2013 The PECOS Development Team
8 //
9 // This library is free software; you can redistribute it and/or
10 // modify it under the terms of the Version 2.1 GNU Lesser General
11 // Public License as published by the Free Software Foundation.
12 //
13 // This library is distributed in the hope that it will be useful,
14 // but WITHOUT ANY WARRANTY; without even the implied warranty of
15 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 // Lesser General Public License for more details.
17 //
18 // You should have received a copy of the GNU Lesser General Public
19 // License along with this library; if not, write to the Free Software
20 // Foundation, Inc. 51 Franklin Street, Fifth Floor,
21 // Boston, MA 02110-1301 USA
22 //
23 //-----------------------------------------------------------------------el-
24 
25 #include <queso/Defines.h>
26 #ifdef QUESO_USES_SEQUENCE_STATISTICAL_OPTIONS
27 
28 #include <queso/SequenceStatisticalOptions.h>
29 #include <queso/Miscellaneous.h>
30 
31 namespace QUESO {
32 
33 SsOptionsValues::SsOptionsValues()
34  :
35  m_initialDiscardedPortions(0 ),//,0.),
36 #ifdef QUESO_COMPUTES_EXTRA_POST_PROCESSING_STATISTICS
37  m_meanMonitorPeriod (UQ_SEQUENCE_MEAN_MONITOR_PERIOD_ODV ),
38  m_bmmRun (UQ_SEQUENCE_BMM_RUN_ODV ),
39  m_bmmLengths (0 ),//,0),
40  m_fftCompute (UQ_SEQUENCE_FFT_COMPUTE_ODV ),
41  m_fftParamId (UQ_SEQUENCE_FFT_PARAM_ID_ODV ),
42  m_fftSize (UQ_SEQUENCE_FFT_SIZE_ODV ),
43  m_fftTestInversion (UQ_SEQUENCE_FFT_TEST_INVERSION_ODV ),
44  m_fftWrite (UQ_SEQUENCE_FFT_WRITE_ODV ),
45  m_psdCompute (UQ_SEQUENCE_PSD_COMPUTE_ODV ),
46  m_psdNumBlocks (UQ_SEQUENCE_PSD_NUM_BLOCKS_ODV ),
47  m_psdHopSizeRatio (UQ_SEQUENCE_PSD_HOP_SIZE_RATIO_ODV ),
48  m_psdParamId (UQ_SEQUENCE_PSD_PARAM_ID_ODV ),
49  m_psdWrite (UQ_SEQUENCE_PSD_WRITE_ODV ),
50  m_psdAtZeroCompute (UQ_SEQUENCE_PSD_AT_ZERO_COMPUTE_ODV ),
51  m_psdAtZeroNumBlocks (0 ),//,0),
52  m_psdAtZeroHopSizeRatio (UQ_SEQUENCE_PSD_AT_ZERO_HOP_SIZE_RATIO_ODV ),
53  m_psdAtZeroDisplay (UQ_SEQUENCE_PSD_AT_ZERO_DISPLAY_ODV ),
54  m_psdAtZeroWrite (UQ_SEQUENCE_PSD_AT_ZERO_WRITE_ODV ),
55  m_gewekeCompute (UQ_SEQUENCE_GEWEKE_COMPUTE_ODV ),
56  m_gewekeNaRatio (UQ_SEQUENCE_GEWEKE_NA_RATIO_ODV ),
57  m_gewekeNbRatio (UQ_SEQUENCE_GEWEKE_NB_RATIO_ODV ),
58  m_gewekeDisplay (UQ_SEQUENCE_GEWEKE_DISPLAY_ODV ),
59  m_gewekeWrite (UQ_SEQUENCE_GEWEKE_WRITE_ODV ),
60  m_meanStaccCompute (UQ_SEQUENCE_MEAN_STACC_COMPUTE_ODV ),
61  m_histCompute (UQ_SEQUENCE_HIST_COMPUTE_ODV ),
62  m_histNumInternalBins (UQ_SEQUENCE_HIST_NUM_INTERNAL_BINS_ODV ),
63  m_cdfStaccCompute (UQ_SEQUENCE_CDF_STACC_COMPUTE_ODV ),
64  m_cdfStaccNumEvalPositions(UQ_SEQUENCE_CDF_STACC_NUM_EVAL_POSITIONS_ODV),
65 #endif
66  m_autoCorrComputeViaDef (UQ_SEQUENCE_AUTO_CORR_COMPUTE_VIA_DEF_ODV ),
67  m_autoCorrComputeViaFft (UQ_SEQUENCE_AUTO_CORR_COMPUTE_VIA_FFT_ODV ),
68  m_autoCorrSecondLag (UQ_SEQUENCE_AUTO_CORR_SECOND_LAG_ODV ),
69  m_autoCorrLagSpacing (UQ_SEQUENCE_AUTO_CORR_LAG_SPACING_ODV ),
70  m_autoCorrNumLags (UQ_SEQUENCE_AUTO_CORR_NUM_LAGS_ODV ),
71  m_autoCorrDisplay (UQ_SEQUENCE_AUTO_CORR_DISPLAY_ODV ),
72  m_autoCorrWrite (UQ_SEQUENCE_AUTO_CORR_WRITE_ODV ),
73  m_kdeCompute (UQ_SEQUENCE_KDE_COMPUTE_ODV ),
74  m_kdeNumEvalPositions (UQ_SEQUENCE_KDE_NUM_EVAL_POSITIONS_ODV ),
75  m_covMatrixCompute (UQ_SEQUENCE_COV_MATRIX_COMPUTE_ODV ),
76  m_corrMatrixCompute (UQ_SEQUENCE_CORR_MATRIX_COMPUTE_ODV )
77 {
78 }
79 
80 SsOptionsValues::~SsOptionsValues()
81 {
82 }
83 
84 SsOptionsValues::SsOptionsValues(const SsOptionsValues& src)
85 {
86  this->copy(src);
87 }
88 
89 SsOptionsValues&
90 SsOptionsValues::operator=(const SsOptionsValues& rhs)
91 {
92  this->copy(rhs);
93  return *this;
94 }
95 
96 void
97 SsOptionsValues::copy(const SsOptionsValues& src)
98 {
99  m_initialDiscardedPortions = src.m_initialDiscardedPortions;
100 #ifdef QUESO_COMPUTES_EXTRA_POST_PROCESSING_STATISTICS
101  m_meanMonitorPeriod = src.m_meanMonitorPeriod;
102  m_bmmRun = src.m_bmmRun;
103  m_bmmLengths = src.m_bmmLengths;
104  m_fftCompute = src.m_fftCompute;
105  m_fftParamId = src.m_fftParamId;
106  m_fftSize = src.m_fftSize;
107  m_fftTestInversion = src.m_fftTestInversion;
108  m_fftWrite = src.m_fftWrite;
109  m_psdCompute = src.m_psdCompute;
110  m_psdNumBlocks = src.m_psdNumBlocks;
111  m_psdHopSizeRatio = src.m_psdHopSizeRatio;
112  m_psdParamId = src.m_psdParamId;
113  m_psdWrite = src.m_psdWrite;
114  m_psdAtZeroCompute = src.m_psdAtZeroCompute;
115  m_psdAtZeroNumBlocks = src.m_psdAtZeroNumBlocks;
116  m_psdAtZeroHopSizeRatio = src.m_psdAtZeroHopSizeRatio;
117  m_psdAtZeroDisplay = src.m_psdAtZeroDisplay;
118  m_psdAtZeroWrite = src.m_psdAtZeroWrite;
119  m_gewekeCompute = src.m_gewekeCompute;
120  m_gewekeNaRatio = src.m_gewekeNaRatio;
121  m_gewekeNbRatio = src.m_gewekeNbRatio;
122  m_gewekeDisplay = src.m_gewekeDisplay;
123  m_gewekeWrite = src.m_gewekeWrite;
124  m_meanStaccCompute = src.m_meanStaccCompute;
125  m_histCompute = src.m_histCompute;
126  m_histNumInternalBins = src.m_histNumInternalBins;
127  m_cdfStaccCompute = src.m_cdfStaccCompute;
128  m_cdfStaccNumEvalPositions = src.m_cdfStaccNumEvalPositions;
129 #endif
130  m_autoCorrComputeViaDef = src.m_autoCorrComputeViaDef;
131  m_autoCorrComputeViaFft = src.m_autoCorrComputeViaFft;
132  m_autoCorrSecondLag = src.m_autoCorrSecondLag;
133  m_autoCorrLagSpacing = src.m_autoCorrLagSpacing;
134  m_autoCorrNumLags = src.m_autoCorrNumLags;
135  m_autoCorrDisplay = src.m_autoCorrDisplay;
136  m_autoCorrWrite = src.m_autoCorrWrite;
137  m_kdeCompute = src.m_kdeCompute;
138  m_kdeNumEvalPositions = src.m_kdeNumEvalPositions;
139  m_covMatrixCompute = src.m_covMatrixCompute;
140  m_corrMatrixCompute = src.m_corrMatrixCompute;
141 
142  return;
143 }
144 
145 SequenceStatisticalOptions::SequenceStatisticalOptions( const BaseEnvironment& env, const std::string& prefix)
148  :
149  m_ov (),
150  m_prefix ((std::string)(prefix) + "stats_"),
151  m_env (env),
152  m_optionsDesc (new po::options_description("Chain statistical options")),
153  m_option_help (m_prefix + "help" ),
154  m_option_initialDiscardedPortions (m_prefix + "initialDiscardedPortions" ),
155 #ifdef QUESO_COMPUTES_EXTRA_POST_PROCESSING_STATISTICS
156  m_option_mean_monitorPeriod (m_prefix + "mean_monitorPeriod" ),
157  m_option_bmm_run (m_prefix + "bmm_run" ),
158  m_option_bmm_lengths (m_prefix + "bmm_lengths" ),
159  m_option_bmm_display (m_prefix + "bmm_display" ),
160  m_option_bmm_write (m_prefix + "bmm_write" ),
161  m_option_fft_compute (m_prefix + "fft_compute" ),
162  m_option_fft_paramId (m_prefix + "fft_paramId" ),
163  m_option_fft_size (m_prefix + "fft_size" ),
164  m_option_fft_testInversion (m_prefix + "fft_testInversion" ),
165  m_option_fft_write (m_prefix + "fft_write" ),
166  m_option_psd_compute (m_prefix + "psd_compute" ),
167  m_option_psd_numBlocks (m_prefix + "psd_numBlocks" ),
168  m_option_psd_hopSizeRatio (m_prefix + "psd_hopSizeRatio" ),
169  m_option_psd_paramId (m_prefix + "psd_paramId" ),
170  m_option_psd_write (m_prefix + "psd_write" ),
171  m_option_psdAtZero_compute (m_prefix + "psdAtZero_compute" ),
172  m_option_psdAtZero_numBlocks (m_prefix + "psdAtZero_numBlocks" ),
173  m_option_psdAtZero_hopSizeRatio (m_prefix + "psdAtZero_hopSizeRatio" ),
174  m_option_psdAtZero_display (m_prefix + "psdAtZero_display" ),
175  m_option_psdAtZero_write (m_prefix + "psdAtZero_write" ),
176  m_option_geweke_compute (m_prefix + "geweke_compute" ),
177  m_option_geweke_naRatio (m_prefix + "geweke_naRatio" ),
178  m_option_geweke_nbRatio (m_prefix + "geweke_nbRatio" ),
179  m_option_geweke_display (m_prefix + "geweke_display" ),
180  m_option_geweke_write (m_prefix + "geweke_write" ),
181  m_option_meanStacc_compute (m_prefix + "meanStacc_compute" ),
182  m_option_hist_compute (m_prefix + "hist_compute" ),
183  m_option_hist_numInternalBins (m_prefix + "hist_numInternalBins" ),
184  m_option_cdfStacc_compute (m_prefix + "cdfStacc_compute" ),
185  m_option_cdfStacc_numEvalPositions(m_prefix + "cdfStacc_numEvalPositions"),
186 #endif
187  m_option_autoCorr_computeViaDef (m_prefix + "autoCorr_computeViaDef" ),
188  m_option_autoCorr_computeViaFft (m_prefix + "autoCorr_computeViaFft" ),
189  m_option_autoCorr_secondLag (m_prefix + "autoCorr_secondLag" ),
190  m_option_autoCorr_lagSpacing (m_prefix + "autoCorr_lagSpacing" ),
191  m_option_autoCorr_numLags (m_prefix + "autoCorr_numLags" ),
192  m_option_autoCorr_display (m_prefix + "autoCorr_display" ),
193  m_option_autoCorr_write (m_prefix + "autoCorr_write" ),
194  m_option_kde_compute (m_prefix + "kde_compute" ),
195  m_option_kde_numEvalPositions (m_prefix + "kde_numEvalPositions" ),
196  m_option_covMatrix_compute (m_prefix + "covMatrix_compute" ),
197  m_option_corrMatrix_compute (m_prefix + "corrMatrix_compute" )
198 {
199  if (m_env.subDisplayFile()) {
200  *m_env.subDisplayFile() << "Entering SequenceStatisticalOptions::constructor(1)"
201  << ", prefix = " << m_prefix
202  << std::endl;
203  }
204 
205  UQ_FATAL_TEST_MACRO(m_env.optionsInputFileName() == "",
206  m_env.worldRank(),
207  "SequenceStatisticalOptions::constructor(1)",
208  "this constructor is incompatible with the abscense of an options input file");
209 
210  defineMyOptions (*m_optionsDesc);
211  m_env.scanInputFileForMyOptions(*m_optionsDesc);
212  getMyOptionValues (*m_optionsDesc);
213 
214  if (m_env.subDisplayFile()) {
215  *m_env.subDisplayFile() << "After reading values of options with prefix '" << m_prefix
216  << "', state of SequenceStatisticalOptions object is:"
217  << "\n" << *this
218  << std::endl;
219  }
220 
221  if (m_env.subDisplayFile()) {
222  *m_env.subDisplayFile() << "Leaving SequenceStatisticalOptions::constructor(1)"
223  << ", prefix = " << m_prefix
224  << std::endl;
225  }
226 }
227 
228 SequenceStatisticalOptions::SequenceStatisticalOptions( const BaseEnvironment& env, const std::string& prefix,
231  const SsOptionsValues& alternativeOptionsValues)
232  :
233  m_ov (alternativeOptionsValues),
234  m_prefix ((std::string)(prefix) + "stats_"),
235  m_env (env),
236  m_optionsDesc (NULL),
237  m_option_help (m_prefix + "help" ),
238  m_option_initialDiscardedPortions (m_prefix + "initialDiscardedPortions" ),
239 #ifdef QUESO_COMPUTES_EXTRA_POST_PROCESSING_STATISTICS
240  m_option_mean_monitorPeriod (m_prefix + "mean_monitorPeriod" ),
241  m_option_bmm_run (m_prefix + "bmm_run" ),
242  m_option_bmm_lengths (m_prefix + "bmm_lengths" ),
243  m_option_bmm_display (m_prefix + "bmm_display" ),
244  m_option_bmm_write (m_prefix + "bmm_write" ),
245  m_option_fft_compute (m_prefix + "fft_compute" ),
246  m_option_fft_paramId (m_prefix + "fft_paramId" ),
247  m_option_fft_size (m_prefix + "fft_size" ),
248  m_option_fft_testInversion (m_prefix + "fft_testInversion" ),
249  m_option_fft_write (m_prefix + "fft_write" ),
250  m_option_psd_compute (m_prefix + "psd_compute" ),
251  m_option_psd_numBlocks (m_prefix + "psd_numBlocks" ),
252  m_option_psd_hopSizeRatio (m_prefix + "psd_hopSizeRatio" ),
253  m_option_psd_paramId (m_prefix + "psd_paramId" ),
254  m_option_psd_write (m_prefix + "psd_write" ),
255  m_option_psdAtZero_compute (m_prefix + "psdAtZero_compute" ),
256  m_option_psdAtZero_numBlocks (m_prefix + "psdAtZero_numBlocks" ),
257  m_option_psdAtZero_hopSizeRatio (m_prefix + "psdAtZero_hopSizeRatio" ),
258  m_option_psdAtZero_display (m_prefix + "psdAtZero_display" ),
259  m_option_psdAtZero_write (m_prefix + "psdAtZero_write" ),
260  m_option_geweke_compute (m_prefix + "geweke_compute" ),
261  m_option_geweke_naRatio (m_prefix + "geweke_naRatio" ),
262  m_option_geweke_nbRatio (m_prefix + "geweke_nbRatio" ),
263  m_option_geweke_display (m_prefix + "geweke_display" ),
264  m_option_geweke_write (m_prefix + "geweke_write" ),
265  m_option_meanStacc_compute (m_prefix + "meanStacc_compute" ),
266  m_option_hist_compute (m_prefix + "hist_compute" ),
267  m_option_hist_numInternalBins (m_prefix + "hist_numInternalBins" ),
268  m_option_cdfStacc_compute (m_prefix + "cdfStacc_compute" ),
269  m_option_cdfStacc_numEvalPositions(m_prefix + "cdfStacc_numEvalPositions"),
270 #endif
271  m_option_autoCorr_computeViaDef (m_prefix + "autoCorr_computeViaDef" ),
272  m_option_autoCorr_computeViaFft (m_prefix + "autoCorr_computeViaFft" ),
273  m_option_autoCorr_secondLag (m_prefix + "autoCorr_secondLag" ),
274  m_option_autoCorr_lagSpacing (m_prefix + "autoCorr_lagSpacing" ),
275  m_option_autoCorr_numLags (m_prefix + "autoCorr_numLags" ),
276  m_option_autoCorr_display (m_prefix + "autoCorr_display" ),
277  m_option_autoCorr_write (m_prefix + "autoCorr_write" ),
278  m_option_kde_compute (m_prefix + "kde_compute" ),
279  m_option_kde_numEvalPositions (m_prefix + "kde_numEvalPositions" ),
280  m_option_covMatrix_compute (m_prefix + "covMatrix_compute" ),
281  m_option_corrMatrix_compute (m_prefix + "corrMatrix_compute" )
282 {
283  if (m_env.subDisplayFile()) {
284  *m_env.subDisplayFile() << "Entering SequenceStatisticalOptions::constructor(2)"
285  << ", prefix = " << m_prefix
286  << std::endl;
287  }
288 
289  UQ_FATAL_TEST_MACRO(m_env.optionsInputFileName() != "",
290  m_env.worldRank(),
291  "SequenceStatisticalOptions::constructor(2)",
292  "this constructor is incompatible with the existence of an options input file");
293 
294  if (m_env.subDisplayFile()) {
295  *m_env.subDisplayFile() << "After setting values of options with prefix '" << m_prefix
296  << "', state of SequenceStatisticalOptions object is:"
297  << "\n" << *this
298  << std::endl;
299  }
300 
301  if (m_env.subDisplayFile()) {
302  *m_env.subDisplayFile() << "Leaving SequenceStatisticalOptions::constructor(2)"
303  << ", prefix = " << m_prefix
304  << std::endl;
305  }
306 }
307 
308 SequenceStatisticalOptions::~SequenceStatisticalOptions()
309 {
310  if (m_optionsDesc) delete m_optionsDesc;
311 }
312 
313 void
314 SequenceStatisticalOptions::defineMyOptions(
315  po::options_description& optionsDesc) const
316 {
317  optionsDesc.add_options()
318  (m_option_help.c_str(), "produce help message for chain statistical options" )
319  (m_option_initialDiscardedPortions.c_str(), po::value<std::string >()->default_value(UQ_SEQUENCE_INITIAL_DISCARDED_PORTIONS_ODV ), "list of initial discarded portions for chain statistics" )
320 #ifdef QUESO_COMPUTES_EXTRA_POST_PROCESSING_STATISTICS
321  (m_option_mean_monitorPeriod.c_str(), po::value<unsigned int>()->default_value(UQ_SEQUENCE_MEAN_MONITOR_PERIOD_ODV ), "period for monitoring mean" )
322  (m_option_bmm_run.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_BMM_RUN_ODV ), "compute variance of sample mean with batch means method" )
323  (m_option_bmm_lengths.c_str(), po::value<std::string >()->default_value(UQ_SEQUENCE_BMM_LENGTHS_ODV ), "list of batch lenghts for BMM" )
324  (m_option_fft_compute.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_FFT_COMPUTE_ODV ), "compute fft" )
325  (m_option_fft_paramId.c_str(), po::value<unsigned int>()->default_value(UQ_SEQUENCE_FFT_PARAM_ID_ODV ), "parameter id for fft computations" )
326  (m_option_fft_size.c_str(), po::value<unsigned int>()->default_value(UQ_SEQUENCE_FFT_SIZE_ODV ), "fft size" )
327  (m_option_fft_testInversion.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_FFT_TEST_INVERSION_ODV ), "test fft inversion" )
328  (m_option_fft_write.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_FFT_WRITE_ODV ), "write fft" )
329  (m_option_psd_compute.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_PSD_COMPUTE_ODV ), "compute psd" )
330  (m_option_psd_numBlocks.c_str(), po::value<unsigned int>()->default_value(UQ_SEQUENCE_PSD_NUM_BLOCKS_ODV ), "number of blocks for computation of psd" )
331  (m_option_psd_hopSizeRatio.c_str(), po::value<double >()->default_value(UQ_SEQUENCE_PSD_HOP_SIZE_RATIO_ODV ), "hop size ratio for psd" )
332  (m_option_psd_paramId.c_str(), po::value<unsigned int>()->default_value(UQ_SEQUENCE_PSD_PARAM_ID_ODV ), "parameter id for psd computations" )
333  (m_option_psd_write.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_PSD_WRITE_ODV ), "write psd" )
334  (m_option_psdAtZero_compute.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_PSD_AT_ZERO_COMPUTE_ODV ), "compute power spectral densities" )
335  (m_option_psdAtZero_numBlocks.c_str(), po::value<std::string >()->default_value(UQ_SEQUENCE_PSD_AT_ZERO_NUM_BLOCKS_ODV ), "list of numbers of blocks for computation of psd at zero" )
336  (m_option_psdAtZero_hopSizeRatio.c_str(), po::value<double >()->default_value(UQ_SEQUENCE_PSD_AT_ZERO_HOP_SIZE_RATIO_ODV ), "hop size ratio for psd at zero" )
337  (m_option_psdAtZero_display.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_PSD_AT_ZERO_DISPLAY_ODV ), "display computed psd at frequency zero on screen" )
338  (m_option_psdAtZero_write.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_PSD_AT_ZERO_WRITE_ODV ), "write computed psd at frequency zero to the output file" )
339  (m_option_geweke_compute.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_GEWEKE_COMPUTE_ODV ), "compute Geweke coefficients" )
340  (m_option_geweke_naRatio.c_str(), po::value<double >()->default_value(UQ_SEQUENCE_GEWEKE_NA_RATIO_ODV ), "ratio NA for Geweke" )
341  (m_option_geweke_nbRatio.c_str(), po::value<double >()->default_value(UQ_SEQUENCE_GEWEKE_NB_RATIO_ODV ), "ratio NB for Geweke" )
342  (m_option_geweke_display.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_GEWEKE_DISPLAY_ODV ), "display computed Geweke on screen" )
343  (m_option_geweke_write.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_GEWEKE_WRITE_ODV ), "write computed Geweke to the output file" )
344  (m_option_meanStacc_compute.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_MEAN_STACC_COMPUTE_ODV ), "compute statistical accuracy of mean" )
345  (m_option_hist_compute.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_HIST_COMPUTE_ODV ), "compute histograms" )
346  (m_option_hist_numInternalBins.c_str(), po::value<unsigned int>()->default_value(UQ_SEQUENCE_HIST_NUM_INTERNAL_BINS_ODV ), "number of internal bins" )
347  (m_option_cdfStacc_compute.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_CDF_STACC_COMPUTE_ODV ), "compute statisical accuracy of cdf" )
348  (m_option_cdfStacc_numEvalPositions.c_str(), po::value<unsigned int>()->default_value(UQ_SEQUENCE_CDF_STACC_NUM_EVAL_POSITIONS_ODV ), "number of evaluations points for statistical accuracy of cdf" )
349 #endif
350  (m_option_autoCorr_computeViaDef.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_AUTO_CORR_COMPUTE_VIA_DEF_ODV ), "compute correlations via definition" )
351  (m_option_autoCorr_computeViaFft.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_AUTO_CORR_COMPUTE_VIA_FFT_ODV ), "compute correlations via fft" )
352  (m_option_autoCorr_secondLag.c_str(), po::value<unsigned int>()->default_value(UQ_SEQUENCE_AUTO_CORR_SECOND_LAG_ODV ), "second lag for computation of autocorrelations" )
353  (m_option_autoCorr_lagSpacing.c_str(), po::value<unsigned int>()->default_value(UQ_SEQUENCE_AUTO_CORR_LAG_SPACING_ODV ), "lag spacing for computation of autocorrelations" )
354  (m_option_autoCorr_numLags.c_str(), po::value<unsigned int>()->default_value(UQ_SEQUENCE_AUTO_CORR_NUM_LAGS_ODV ), "number of lags for computation of autocorrelations" )
355  (m_option_autoCorr_display.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_AUTO_CORR_DISPLAY_ODV ), "display computed autocorrelations on the screen" )
356  (m_option_autoCorr_write.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_AUTO_CORR_WRITE_ODV ), "write computed autocorrelations to the output file" )
357  (m_option_kde_compute.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_KDE_COMPUTE_ODV ), "compute kernel density estimators" )
358  (m_option_kde_numEvalPositions.c_str(), po::value<unsigned int>()->default_value(UQ_SEQUENCE_KDE_NUM_EVAL_POSITIONS_ODV ), "number of evaluation positions" )
359  (m_option_covMatrix_compute.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_COV_MATRIX_COMPUTE_ODV ), "compute covariance matrix" )
360  (m_option_corrMatrix_compute.c_str(), po::value<bool >()->default_value(UQ_SEQUENCE_CORR_MATRIX_COMPUTE_ODV ), "compute correlation matrix" )
361  ;
362 
363  return;
364 }
365 
366 void
367 SequenceStatisticalOptions::getMyOptionValues(
368  po::options_description& optionsDesc)
369 {
370  if (m_env.allOptionsMap().count(m_option_help)) {
371  if (m_env.subDisplayFile()) {
372  *m_env.subDisplayFile() << optionsDesc
373  << std::endl;
374  }
375  }
376 
377  if (m_env.allOptionsMap().count(m_option_initialDiscardedPortions)) {
378  m_ov.m_initialDiscardedPortions.clear();
379  std::vector<double> tmpPortions(0,0.);
380  std::string inputString = m_env.allOptionsMap()[m_option_initialDiscardedPortions].as<std::string>();
381  MiscReadDoublesFromString(inputString,tmpPortions);
382  //if (m_env.subDisplayFile()) {
383  // *m_env.subDisplayFile() << "In SequenceStatisticalOptions::getMyOptionValues(): percents = ";
384  // for (unsigned int i = 0; i < tmpPortions.size(); ++i) {
385  // *m_env.subDisplayFile() << " " << tmpPortions[i];
386  // }
387  // *m_env.subDisplayFile() << std::endl;
388  //}
389 
390  if (tmpPortions.size() > 0) {
391  m_ov.m_initialDiscardedPortions.resize(tmpPortions.size(),0.);
392  for (unsigned int i = 0; i < m_ov.m_initialDiscardedPortions.size(); ++i) {
393  m_ov.m_initialDiscardedPortions[i] = tmpPortions[i];
394  }
395  }
396  }
397 
398 #ifdef QUESO_COMPUTES_EXTRA_POST_PROCESSING_STATISTICS
399  if (m_env.allOptionsMap().count(m_option_mean_monitorPeriod)) {
400  m_ov.m_meanMonitorPeriod = m_env.allOptionsMap()[m_option_mean_monitorPeriod].as<unsigned int>();
401  }
402 
403  if (m_env.allOptionsMap().count(m_option_bmm_run)) {
404  m_ov.m_bmmRun = m_env.allOptionsMap()[m_option_bmm_run].as<bool>();
405  }
406 
407  if (m_env.allOptionsMap().count(m_option_bmm_lengths)) {
408  m_ov.m_bmmLengths.clear();
409  std::vector<double> tmpLengths(0,0.);
410  std::string inputString = m_env.allOptionsMap()[m_option_bmm_lengths].as<std::string>();
411  MiscReadDoublesFromString(inputString,tmpLengths);
412  //if (m_env.subDisplayFile()) {
413  // *m_env.subDisplayFile() << "In SequenceStatisticalOptions::getMyOptionValues(): lengths for BMM = ";
414  // for (unsigned int i = 0; i < tmpLengths.size(); ++i) {
415  // *m_env.subDisplayFile() << " " << tmpLengths[i];
416  // }
417  // *m_env.subDisplayFile() << std::endl;
418  //}
419 
420  if (tmpLengths.size() > 0) {
421  m_ov.m_bmmLengths.resize(tmpLengths.size(),0);
422  for (unsigned int i = 0; i < m_ov.m_bmmLengths.size(); ++i) {
423  m_ov.m_bmmLengths[i] = (unsigned int) tmpLengths[i];
424  }
425  }
426  }
427 
428  if (m_env.allOptionsMap().count(m_option_fft_compute)) {
429  m_ov.m_fftCompute = m_env.allOptionsMap()[m_option_fft_compute].as<bool>();
430  }
431 
432  if (m_env.allOptionsMap().count(m_option_fft_paramId)) {
433  m_ov.m_fftParamId = m_env.allOptionsMap()[m_option_fft_paramId].as<unsigned int>();
434  }
435 
436  if (m_env.allOptionsMap().count(m_option_fft_size)) {
437  m_ov.m_fftSize = m_env.allOptionsMap()[m_option_fft_size].as<unsigned int>();
438  }
439 
440  if (m_env.allOptionsMap().count(m_option_fft_testInversion)) {
441  m_ov.m_fftTestInversion = m_env.allOptionsMap()[m_option_fft_testInversion].as<bool>();
442  }
443 
444  if (m_env.allOptionsMap().count(m_option_fft_write)) {
445  m_ov.m_fftWrite = m_env.allOptionsMap()[m_option_fft_write].as<bool>();
446  }
447 
448  if (m_env.allOptionsMap().count(m_option_psd_compute)) {
449  m_ov.m_psdCompute = m_env.allOptionsMap()[m_option_psd_compute].as<bool>();
450  }
451 
452  if (m_env.allOptionsMap().count(m_option_psd_numBlocks)) {
453  m_ov.m_psdNumBlocks = m_env.allOptionsMap()[m_option_psd_numBlocks].as<unsigned int>();
454  }
455 
456  if (m_env.allOptionsMap().count(m_option_psd_hopSizeRatio)) {
457  m_ov.m_psdHopSizeRatio = m_env.allOptionsMap()[m_option_psd_hopSizeRatio].as<double>();
458  }
459 
460  if (m_env.allOptionsMap().count(m_option_psd_paramId)) {
461  m_ov.m_psdParamId = m_env.allOptionsMap()[m_option_psd_paramId].as<unsigned int>();
462  }
463 
464  if (m_env.allOptionsMap().count(m_option_psd_write)) {
465  m_ov.m_psdWrite = m_env.allOptionsMap()[m_option_psd_write].as<bool>();
466  }
467 
468  if (m_env.allOptionsMap().count(m_option_psdAtZero_compute)) {
469  m_ov.m_psdAtZeroCompute = m_env.allOptionsMap()[m_option_psdAtZero_compute].as<bool>();
470  }
471 
472  if (m_env.allOptionsMap().count(m_option_psdAtZero_numBlocks)) {
473  m_ov.m_psdAtZeroNumBlocks.clear();
474  std::vector<double> tmpNumBlocks(0,0.);
475  std::string inputString = m_env.allOptionsMap()[m_option_psdAtZero_numBlocks].as<std::string>();
476  MiscReadDoublesFromString(inputString,tmpNumBlocks);
477  //if (m_env.subDisplayFile()) {
478  // *m_env.subDisplayFile() << "In SequenceStatisticalOptions::getMyOptionValues(): numBlocks for psdAtZero = ";
479  // for (unsigned int i = 0; i < tmpNumBlocks.size(); ++i) {
480  // *m_env.subDisplayFile() << " " << numBlocks[i];
481  // }
482  // *m_env.subDisplayFile() << std::endl;
483  //}
484 
485  if (tmpNumBlocks.size() > 0) {
486  m_ov.m_psdAtZeroNumBlocks.resize(tmpNumBlocks.size(),0);
487  for (unsigned int i = 0; i < m_ov.m_psdAtZeroNumBlocks.size(); ++i) {
488  m_ov.m_psdAtZeroNumBlocks[i] = (unsigned int) tmpNumBlocks[i];
489  }
490  }
491  }
492 
493  if (m_env.allOptionsMap().count(m_option_psdAtZero_hopSizeRatio)) {
494  m_ov.m_psdAtZeroHopSizeRatio = m_env.allOptionsMap()[m_option_psdAtZero_hopSizeRatio].as<double>();
495  }
496 
497  if (m_env.allOptionsMap().count(m_option_psdAtZero_display)) {
498  m_ov.m_psdAtZeroDisplay = m_env.allOptionsMap()[m_option_psdAtZero_display].as<bool>();
499  }
500 
501  if (m_env.allOptionsMap().count(m_option_psdAtZero_write)) {
502  m_ov.m_psdAtZeroWrite = m_env.allOptionsMap()[m_option_psdAtZero_write].as<bool>();
503  }
504 
505  if (m_env.allOptionsMap().count(m_option_geweke_compute)) {
506  m_ov.m_gewekeCompute = m_env.allOptionsMap()[m_option_geweke_compute].as<bool>();
507  }
508 
509  if (m_env.allOptionsMap().count(m_option_geweke_naRatio)) {
510  m_ov.m_gewekeNaRatio = m_env.allOptionsMap()[m_option_geweke_naRatio].as<double>();
511  }
512 
513  if (m_env.allOptionsMap().count(m_option_geweke_nbRatio)) {
514  m_ov.m_gewekeNbRatio = m_env.allOptionsMap()[m_option_geweke_nbRatio].as<double>();
515  }
516 
517  if (m_env.allOptionsMap().count(m_option_geweke_display)) {
518  m_ov.m_gewekeDisplay = m_env.allOptionsMap()[m_option_geweke_display].as<bool>();
519  }
520 
521  if (m_env.allOptionsMap().count(m_option_geweke_write)) {
522  m_ov.m_gewekeWrite = m_env.allOptionsMap()[m_option_geweke_write].as<bool>();
523  }
524 
525  if (m_env.allOptionsMap().count(m_option_meanStacc_compute)) {
526  m_ov.m_meanStaccCompute = m_env.allOptionsMap()[m_option_meanStacc_compute].as<bool>();
527  }
528 
529  if (m_env.allOptionsMap().count(m_option_hist_compute)) {
530  m_ov.m_histCompute = m_env.allOptionsMap()[m_option_hist_compute].as<bool>();
531  }
532 
533  if (m_env.allOptionsMap().count(m_option_hist_numInternalBins)) {
534  m_ov.m_histNumInternalBins = m_env.allOptionsMap()[m_option_hist_numInternalBins].as<unsigned int>();
535  }
536 
537  if (m_env.allOptionsMap().count(m_option_cdfStacc_compute)) {
538  m_ov.m_cdfStaccCompute = m_env.allOptionsMap()[m_option_cdfStacc_compute].as<bool>();
539  }
540 
541  if (m_env.allOptionsMap().count(m_option_cdfStacc_numEvalPositions)) {
542  m_ov.m_cdfStaccNumEvalPositions = m_env.allOptionsMap()[m_option_cdfStacc_numEvalPositions].as<unsigned int>();
543  }
544 #endif
545  if (m_env.allOptionsMap().count(m_option_autoCorr_computeViaDef)) {
546  m_ov.m_autoCorrComputeViaDef = m_env.allOptionsMap()[m_option_autoCorr_computeViaDef].as<bool>();
547  }
548 
549  if (m_env.allOptionsMap().count(m_option_autoCorr_computeViaFft)) {
550  m_ov.m_autoCorrComputeViaFft = m_env.allOptionsMap()[m_option_autoCorr_computeViaFft].as<bool>();
551  }
552 
553  if (m_env.allOptionsMap().count(m_option_autoCorr_secondLag)) {
554  m_ov.m_autoCorrSecondLag = m_env.allOptionsMap()[m_option_autoCorr_secondLag].as<unsigned int>();
555  }
556 
557  if (m_env.allOptionsMap().count(m_option_autoCorr_lagSpacing)) {
558  m_ov.m_autoCorrLagSpacing = m_env.allOptionsMap()[m_option_autoCorr_lagSpacing].as<unsigned int>();
559  }
560 
561  if (m_env.allOptionsMap().count(m_option_autoCorr_numLags)) {
562  m_ov.m_autoCorrNumLags = m_env.allOptionsMap()[m_option_autoCorr_numLags].as<unsigned int>();
563  }
564 
565  if (m_env.allOptionsMap().count(m_option_autoCorr_display)) {
566  m_ov.m_autoCorrDisplay = m_env.allOptionsMap()[m_option_autoCorr_display].as<bool>();
567  }
568 
569  if (m_env.allOptionsMap().count(m_option_autoCorr_write)) {
570  m_ov.m_autoCorrWrite = m_env.allOptionsMap()[m_option_autoCorr_write].as<bool>();
571  }
572 
573  if (m_env.allOptionsMap().count(m_option_kde_compute)) {
574  m_ov.m_kdeCompute = m_env.allOptionsMap()[m_option_kde_compute].as<bool>();
575  }
576 
577  if (m_env.allOptionsMap().count(m_option_kde_numEvalPositions)) {
578  m_ov.m_kdeNumEvalPositions = m_env.allOptionsMap()[m_option_kde_numEvalPositions].as<unsigned int>();
579  }
580 
581  if (m_env.allOptionsMap().count(m_option_covMatrix_compute)) {
582  m_ov.m_covMatrixCompute = m_env.allOptionsMap()[m_option_covMatrix_compute].as<bool>();
583  }
584 
585  if (m_env.allOptionsMap().count(m_option_corrMatrix_compute)) {
586  m_ov.m_corrMatrixCompute = m_env.allOptionsMap()[m_option_corrMatrix_compute].as<bool>();
587  }
588 
589  return;
590 }
591 
592 const std::vector<double>&
593 SequenceStatisticalOptions::initialDiscardedPortions() const
594 {
595  return m_ov.m_initialDiscardedPortions;
596 }
597 
598 #ifdef QUESO_COMPUTES_EXTRA_POST_PROCESSING_STATISTICS
599 unsigned int
600 SequenceStatisticalOptions::meanMonitorPeriod() const
601 {
602  return m_ov.m_meanMonitorPeriod;
603 }
604 
605 bool
606 SequenceStatisticalOptions::bmmRun() const
607 {
608  return m_ov.m_bmmRun;
609 }
610 
611 const std::vector<unsigned int>&
612 SequenceStatisticalOptions::bmmLengths() const
613 {
614  return m_ov.m_bmmLengths;
615 }
616 
617 bool
618 SequenceStatisticalOptions::bmmDisplay() const
619 {
620  return m_ov.m_bmmDisplay;
621 }
622 
623 bool
624 SequenceStatisticalOptions::bmmWrite() const
625 {
626  return m_ov.m_bmmWrite;
627 }
628 
629 bool
630 SequenceStatisticalOptions::fftCompute() const
631 {
632  return m_ov.m_fftCompute;
633 }
634 
635 unsigned int
636 SequenceStatisticalOptions::fftParamId() const
637 {
638  return m_ov.m_fftParamId;
639 }
640 
641 unsigned int
642 SequenceStatisticalOptions::fftSize() const
643 {
644  return m_ov.m_fftSize;
645 }
646 
647 bool
648 SequenceStatisticalOptions::fftTestInversion() const
649 {
650  return m_ov.m_fftTestInversion;
651 }
652 
653 bool
654 SequenceStatisticalOptions::fftWrite() const
655 {
656  return m_ov.m_fftWrite;
657 }
658 
659 bool
660 SequenceStatisticalOptions::psdCompute() const
661 {
662  return m_ov.m_psdCompute;
663 }
664 
665 unsigned int
666 SequenceStatisticalOptions::psdNumBlocks() const
667 {
668  return m_ov.m_psdNumBlocks;
669 }
670 
671 double
672 SequenceStatisticalOptions::psdHopSizeRatio() const
673 {
674  return m_ov.m_psdHopSizeRatio;
675 }
676 
677 unsigned int
678 SequenceStatisticalOptions::psdParamId() const
679 {
680  return m_ov.m_psdParamId;
681 }
682 
683 bool
684 SequenceStatisticalOptions::psdWrite() const
685 {
686  return m_ov.m_psdWrite;
687 }
688 
689 bool
690 SequenceStatisticalOptions::psdAtZeroCompute() const
691 {
692  return m_ov.m_psdAtZeroCompute;
693 }
694 
695 const std::vector<unsigned int>&
696 SequenceStatisticalOptions::psdAtZeroNumBlocks() const
697 {
698  return m_ov.m_psdAtZeroNumBlocks;
699 }
700 
701 double
702 SequenceStatisticalOptions::psdAtZeroHopSizeRatio() const
703 {
704  return m_ov.m_psdAtZeroHopSizeRatio;
705 }
706 
707 bool
708 SequenceStatisticalOptions::psdAtZeroDisplay() const
709 {
710  return m_ov.m_psdAtZeroDisplay;
711 }
712 
713 bool
714 SequenceStatisticalOptions::psdAtZeroWrite() const
715 {
716  return m_ov.m_psdAtZeroWrite;
717 }
718 
719 bool
720 SequenceStatisticalOptions::gewekeCompute() const
721 {
722  return m_ov.m_gewekeCompute;
723 }
724 
725 double
726 SequenceStatisticalOptions::gewekeNaRatio() const
727 {
728  return m_ov.m_gewekeNaRatio;
729 }
730 
731 double
732 SequenceStatisticalOptions::gewekeNbRatio() const
733 {
734  return m_ov.m_gewekeNbRatio;
735 }
736 
737 bool
738 SequenceStatisticalOptions::gewekeDisplay() const
739 {
740  return m_ov.m_gewekeDisplay;
741 }
742 
743 bool
744 SequenceStatisticalOptions::gewekeWrite() const
745 {
746  return m_ov.m_gewekeWrite;
747 }
748 
749 bool
750 SequenceStatisticalOptions::meanStaccCompute() const
751 {
752  return m_ov.m_meanStaccCompute;
753 }
754 
755 bool
756 SequenceStatisticalOptions::histCompute() const
757 {
758  return m_ov.m_histCompute;
759 }
760 
761 unsigned int
762 SequenceStatisticalOptions::histNumInternalBins() const
763 {
764  return m_ov.m_histNumInternalBins;
765 }
766 
767 bool
768 SequenceStatisticalOptions::cdfStaccCompute() const
769 {
770  return m_ov.m_cdfStaccCompute;
771 }
772 
773 unsigned int
774 SequenceStatisticalOptions::cdfStaccNumEvalPositions() const
775 {
776  return m_ov.m_cdfStaccNumEvalPositions;
777 }
778 #endif
779 bool
780 SequenceStatisticalOptions::autoCorrComputeViaDef() const
781 {
782  return m_ov.m_autoCorrComputeViaDef;
783 }
784 
785 bool
786 SequenceStatisticalOptions::autoCorrComputeViaFft() const
787 {
788  return m_ov.m_autoCorrComputeViaFft;
789 }
790 
791 unsigned int
792 SequenceStatisticalOptions::autoCorrSecondLag() const
793 {
794  return m_ov.m_autoCorrSecondLag;
795 }
796 
797 unsigned int
798 SequenceStatisticalOptions::autoCorrLagSpacing() const
799 {
800  return m_ov.m_autoCorrLagSpacing;
801 }
802 
803 unsigned int
804 SequenceStatisticalOptions::autoCorrNumLags() const
805 {
806  return m_ov.m_autoCorrNumLags;
807 }
808 
809 bool
810 SequenceStatisticalOptions::autoCorrDisplay() const
811 {
812  return m_ov.m_autoCorrDisplay;
813 }
814 
815 bool
816 SequenceStatisticalOptions::autoCorrWrite() const
817 {
818  return m_ov.m_autoCorrWrite;
819 }
820 
821 bool
822 SequenceStatisticalOptions::kdeCompute() const
823 {
824  return m_ov.m_kdeCompute;
825 }
826 
827 unsigned int
828 SequenceStatisticalOptions::kdeNumEvalPositions() const
829 {
830  return m_ov.m_kdeNumEvalPositions;
831 }
832 
833 bool
834 SequenceStatisticalOptions::covMatrixCompute() const
835 {
836  return m_ov.m_covMatrixCompute;
837 }
838 
839 bool
840 SequenceStatisticalOptions::corrMatrixCompute() const
841 {
842  return m_ov.m_corrMatrixCompute;
843 }
844 
845 void
846 SequenceStatisticalOptions::print(std::ostream& os) const
847 {
848  os << "\n" << m_option_initialDiscardedPortions << " = ";
849  for (unsigned int i = 0; i < m_ov.m_initialDiscardedPortions.size(); ++i) {
850  os << m_ov.m_initialDiscardedPortions[i] << " ";
851  }
852  os
853 #ifdef QUESO_COMPUTES_EXTRA_POST_PROCESSING_STATISTICS
854  << "\n" << m_option_mean_monitorPeriod << " = " << m_ov.m_meanMonitorPeriod
855  << "\n" << m_option_bmm_run << " = " << m_ov.m_bmmRun
856  << "\n" << m_option_bmm_lengths << " = ";
857  for (unsigned int i = 0; i < m_ov.m_bmmLengths.size(); ++i) {
858  os << m_ov.m_bmmLengths[i] << " ";
859  }
860  os << "\n" << m_option_fft_compute << " = " << m_ov.m_fftCompute
861  << "\n" << m_option_fft_paramId << " = " << m_ov.m_fftParamId
862  << "\n" << m_option_fft_size << " = " << m_ov.m_fftSize
863  << "\n" << m_option_fft_testInversion << " = " << m_ov.m_fftTestInversion
864  << "\n" << m_option_fft_write << " = " << m_ov.m_fftWrite
865  << "\n" << m_option_psd_compute << " = " << m_ov.m_psdCompute
866  << "\n" << m_option_psd_paramId << " = " << m_ov.m_psdParamId
867  << "\n" << m_option_psd_numBlocks << " = " << m_ov.m_psdNumBlocks
868  << "\n" << m_option_psd_hopSizeRatio << " = " << m_ov.m_psdHopSizeRatio
869  << "\n" << m_option_psd_write << " = " << m_ov.m_psdWrite
870  << "\n" << m_option_psdAtZero_compute << " = " << m_ov.m_psdAtZeroCompute
871  << "\n" << m_option_psdAtZero_numBlocks << " = ";
872  for (unsigned int i = 0; i < m_ov.m_psdAtZeroNumBlocks.size(); ++i) {
873  os << m_ov.m_psdAtZeroNumBlocks[i] << " ";
874  }
875  os << "\n" << m_option_psdAtZero_hopSizeRatio << " = " << m_ov.m_psdAtZeroHopSizeRatio
876  << "\n" << m_option_psdAtZero_display << " = " << m_ov.m_psdAtZeroDisplay
877  << "\n" << m_option_psdAtZero_write << " = " << m_ov.m_psdAtZeroWrite
878  << "\n" << m_option_geweke_compute << " = " << m_ov.m_gewekeCompute
879  << "\n" << m_option_geweke_naRatio << " = " << m_ov.m_gewekeNaRatio
880  << "\n" << m_option_geweke_nbRatio << " = " << m_ov.m_gewekeNbRatio
881  << "\n" << m_option_geweke_display << " = " << m_ov.m_gewekeDisplay
882  << "\n" << m_option_geweke_write << " = " << m_ov.m_gewekeWrite
883  << "\n" << m_option_meanStacc_compute << " = " << m_ov.m_meanStaccCompute
884  << "\n" << m_option_hist_compute << " = " << m_ov.m_histCompute
885  << "\n" << m_option_hist_numInternalBins << " = " << m_ov.m_histNumInternalBins
886  << "\n" << m_option_cdfStacc_compute << " = " << m_ov.m_cdfStaccCompute
887  << "\n" << m_option_cdfStacc_numEvalPositions << " = " << m_ov.m_cdfStaccNumEvalPositions
888 #endif
889  << "\n" << m_option_autoCorr_computeViaDef << " = " << m_ov.m_autoCorrComputeViaDef
890  << "\n" << m_option_autoCorr_computeViaFft << " = " << m_ov.m_autoCorrComputeViaFft
891  << "\n" << m_option_autoCorr_secondLag << " = " << m_ov.m_autoCorrSecondLag
892  << "\n" << m_option_autoCorr_lagSpacing << " = " << m_ov.m_autoCorrLagSpacing
893  << "\n" << m_option_autoCorr_numLags << " = " << m_ov.m_autoCorrNumLags
894  << "\n" << m_option_autoCorr_display << " = " << m_ov.m_autoCorrDisplay
895  << "\n" << m_option_autoCorr_write << " = " << m_ov.m_autoCorrWrite
896  << "\n" << m_option_kde_compute << " = " << m_ov.m_kdeCompute
897  << "\n" << m_option_kde_numEvalPositions << " = " << m_ov.m_kdeNumEvalPositions
898  << "\n" << m_option_covMatrix_compute << " = " << m_ov.m_covMatrixCompute
899  << "\n" << m_option_corrMatrix_compute << " = " << m_ov.m_corrMatrixCompute
900  << std::endl;
901 
902  return;
903 }
904 
905 std::ostream&
906 operator<<(std::ostream& os, const SequenceStatisticalOptions& obj)
907 {
908  obj.print(os);
909 
910  return os;
911 }
912 
913 } // End namespace QUESO
914 
915 #endif // ifdef QUESO_USES_SEQUENCE_STATISTICAL_OPTIONS
QUESO::BaseEnvironment::worldRank
int worldRank() const
Returns the process world rank.
Definition: Environment.C:235
QUESO::BaseEnvironment::subDisplayFile
std::ofstream * subDisplayFile() const
Access function for m_subDisplayFile (displays file on stream).
Definition: Environment.C:305
QUESO::operator<<
std::ostream & operator<<(std::ostream &os, const BaseEnvironment &obj)
QUESO::BaseVectorSequence::m_env
const BaseEnvironment & m_env
Definition: VectorSequence.h:423
QUESO::MiscReadDoublesFromString
void MiscReadDoublesFromString(const std::string &inputString, std::vector< double > &outputDoubles)
Definition: Miscellaneous.C:39
QUESO::BaseEnvironment::optionsInputFileName
std::string optionsInputFileName() const
Access to the attribute m_optionsInputFileName, which stores the name of the input file passed by the...
Definition: Environment.C:341
UQ_FATAL_TEST_MACRO
#define UQ_FATAL_TEST_MACRO(test, givenRank, where, what)
Definition: Defines.h:223
QUESO::BaseEnvironment::scanInputFileForMyOptions
void scanInputFileForMyOptions(const po::options_description &optionsDesc) const
This method scans the input file provided by the user to QUESO.
Definition: Environment.C:378
Generated on Thu Apr 23 2015 19:26:15 for queso-0.51.1 by  doxygen 1.8.5