Dakota  Version 6.4
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Pages
Public Member Functions | Protected Member Functions | Static Protected Member Functions | Private Attributes | Static Private Attributes | List of all members
NonDLHSSampling Class Reference

Performs LHS and Monte Carlo sampling for uncertainty quantification. More...

Inheritance diagram for NonDLHSSampling:
NonDSampling NonD Analyzer Iterator

Public Member Functions

 NonDLHSSampling (ProblemDescDB &problem_db, Model &model)
 standard constructor More...
 
 NonDLHSSampling (Model &model, unsigned short sample_type, int samples, int seed, const String &rng, bool vary_pattern=true, short sampling_vars_mode=ACTIVE)
 alternate constructor for sample generation and evaluation "on the fly" More...
 
 NonDLHSSampling (unsigned short sample_type, int samples, int seed, const String &rng, const RealVector &lower_bnds, const RealVector &upper_bnds)
 alternate constructor for uniform sample generation "on the fly" More...
 
 NonDLHSSampling (unsigned short sample_type, int samples, int seed, const String &rng, const RealVector &means, const RealVector &std_devs, const RealVector &lower_bnds, const RealVector &upper_bnds, RealSymMatrix &correl)
 alternate constructor for sample generation of correlated normals "on the fly" More...
 
 ~NonDLHSSampling ()
 destructor
 
- Public Member Functions inherited from NonDSampling
 NonDSampling (Model &model, const RealMatrix &sample_matrix)
 alternate constructor for evaluating and computing statistics for the provided set of samples More...
 
void compute_statistics (const RealMatrix &vars_samples, const IntResponseMap &resp_samples)
 For the input sample set, computes mean, standard deviation, and probability/reliability/response levels (aleatory uncertainties) or intervals (epsitemic or mixed uncertainties)
 
void compute_intervals (RealRealPairArray &extreme_fns)
 called by compute_statistics() to calculate min/max intervals using allResponses
 
void compute_intervals (const IntResponseMap &samples)
 called by compute_statistics() to calculate extremeValues from samples
 
void compute_intervals (RealRealPairArray &extreme_fns, const IntResponseMap &samples)
 called by compute_statistics() to calculate min/max intervals using samples
 
void compute_moments (const IntResponseMap &samples)
 called by compute_statistics() to calculate sample moments and confidence intervals
 
void compute_level_mappings (const IntResponseMap &samples)
 called by compute_statistics() to calculate CDF/CCDF mappings of z to p/beta and of p/beta to z as well as PDFs More...
 
void print_statistics (std::ostream &s) const
 prints the statistics computed in compute_statistics()
 
void print_intervals (std::ostream &s) const
 prints the intervals computed in compute_intervals() with default qoi_type and moment_labels
 
void print_intervals (std::ostream &s, String qoi_type, const StringArray &interval_labels) const
 prints the intervals computed in compute_intervals()
 
void print_moments (std::ostream &s) const
 prints the moments computed in compute_moments() with default qoi_type and moment_labels
 
void print_moments (std::ostream &s, String qoi_type, const StringArray &moment_labels) const
 prints the moments computed in compute_moments()
 
void update_final_statistics ()
 update finalStatistics from minValues/maxValues, momentStats, and computedProbLevels/computedRelLevels/computedRespLevels
 
void compute_moments (const RealMatrix &samples)
 calculates sample moments for an array of observations for a set of QoI
 
void transform_samples (bool x_to_u=true)
 transform allSamples imported by alternate constructor. This is needed since random variable distribution parameters are not updated until run time and an imported sample_matrix is typically in x-space.
 
void transform_samples (RealMatrix &sample_matrix, bool x_to_u, int num_samples=0)
 transform the specified samples matrix from x to u or u to x
 
- Public Member Functions inherited from NonD
void initialize_random_variables (short u_space_type)
 initialize natafTransform based on distribution data from iteratedModel More...
 
void initialize_random_variables (const Pecos::ProbabilityTransformation &transform, bool deep_copy=false)
 alternate form: initialize natafTransform based on incoming data More...
 
void initialize_random_variable_transformation ()
 instantiate natafTransform
 
void initialize_random_variable_types ()
 initializes ranVarTypesX within natafTransform (u-space not needed) More...
 
void initialize_random_variable_types (short u_space_type)
 initializes ranVarTypesX and ranVarTypesU within natafTransform More...
 
void initialize_random_variable_parameters ()
 initializes ranVarMeansX, ranVarStdDevsX, ranVarLowerBndsX, ranVarUpperBndsX, and ranVarAddtlParamsX within natafTransform More...
 
void initialize_random_variable_correlations ()
 propagate iteratedModel correlations to natafTransform
 
void verify_correlation_support (short u_space_type)
 verify that correlation warping is supported by Nataf for given variable types
 
void transform_correlations ()
 perform correlation warping for variable types supported by Nataf
 
void requested_levels (const RealVectorArray &req_resp_levels, const RealVectorArray &req_prob_levels, const RealVectorArray &req_rel_levels, const RealVectorArray &req_gen_rel_levels, short resp_lev_tgt, short resp_lev_tgt_reduce, bool cdf_flag, bool pdf_output)
 set requestedRespLevels, requestedProbLevels, requestedRelLevels, requestedGenRelLevels, respLevelTarget, cdfFlag, and pdfOutput (used in combination with alternate ctors)
 
void distribution_parameter_derivatives (bool dist_param_derivs)
 set distParamDerivs
 
void print_level_mappings (std::ostream &s) const
 prints the z/p/beta/beta* mappings reflected in {requested,computed}{Resp,Prob,Rel,GenRel}Levels for default qoi_type and qoi_labels
 
void print_level_mappings (std::ostream &s, String qoi_type, const StringArray &qoi_labels) const
 prints the z/p/beta/beta* mappings reflected in {requested,computed}{Resp,Prob,Rel,GenRel}Levels More...
 
bool resize ()
 reinitializes iterator based on new variable size
 
bool pdf_output () const
 get pdfOutput
 
void pdf_output (bool output)
 set pdfOutput
 
Pecos::ProbabilityTransformation & variable_transformation ()
 return natafTransform
 
- Public Member Functions inherited from Analyzer
const VariablesArray & all_variables ()
 return the complete set of evaluated variables
 
const RealMatrix & all_samples ()
 return the complete set of evaluated samples
 
const IntResponseMap & all_responses () const
 return the complete set of computed responses
 
- Public Member Functions inherited from Iterator
 Iterator ()
 default constructor More...
 
 Iterator (ProblemDescDB &problem_db)
 standard envelope constructor, which constructs its own model(s) More...
 
 Iterator (ProblemDescDB &problem_db, Model &model)
 alternate envelope constructor which uses the ProblemDescDB but accepts a model from a higher level (meta-iterator) context, instead of constructing its own More...
 
 Iterator (const String &method_string, Model &model)
 alternate envelope constructor for instantiations by name without the ProblemDescDB More...
 
 Iterator (const Iterator &iterator)
 copy constructor More...
 
virtual ~Iterator ()
 destructor More...
 
Iterator operator= (const Iterator &iterator)
 assignment operator More...
 
virtual void derived_free_communicators (ParLevLIter pl_iter)
 derived class contributions to freeing the communicators associated with this Iterator instance
 
virtual void reset ()
 restore initial state for repeated sub-iterator executions
 
virtual void initialize_iterator (int job_index)
 used by IteratorScheduler to set the starting data for a run
 
virtual void pack_parameters_buffer (MPIPackBuffer &send_buffer, int job_index)
 used by IteratorScheduler to pack starting data for an iterator run
 
virtual void unpack_parameters_buffer (MPIUnpackBuffer &recv_buffer)
 used by IteratorScheduler to unpack starting data for an iterator run
 
virtual void unpack_parameters_initialize (MPIUnpackBuffer &recv_buffer)
 used by IteratorScheduler to unpack starting data and initialize an iterator run
 
virtual void pack_results_buffer (MPIPackBuffer &send_buffer, int job_index)
 used by IteratorScheduler to pack results data from an iterator run
 
virtual void unpack_results_buffer (MPIUnpackBuffer &recv_buffer, int job_index)
 used by IteratorScheduler to unpack results data from an iterator run
 
virtual void update_local_results (int job_index)
 used by IteratorScheduler to update local results arrays
 
virtual bool accepts_multiple_points () const
 indicates if this iterator accepts multiple initial points. Default return is false. Override to return true if appropriate.
 
virtual void initial_points (const VariablesArray &pts)
 sets the multiple initial points for this iterator. This should only be used if accepts_multiple_points() returns true.
 
virtual void initialize_graphics (int iterator_server_id=1)
 initialize the 2D graphics window and the tabular graphics data More...
 
virtual unsigned short uses_method () const
 return name of any enabling iterator used by this iterator More...
 
virtual void method_recourse ()
 perform a method switch, if possible, due to a detected conflict
 
virtual IntIntPair estimate_partition_bounds ()
 estimate the minimum and maximum partition sizes that can be utilized by this Iterator
 
void init_communicators (ParLevLIter pl_iter)
 initialize the communicators associated with this Iterator instance
 
void set_communicators (ParLevLIter pl_iter)
 set the communicators associated with this Iterator instance
 
void free_communicators (ParLevLIter pl_iter)
 free the communicators associated with this Iterator instance
 
void resize_communicators (ParLevLIter pl_iter, bool reinit_comms)
 Resize the communicators. This is called from the letter's resize()
 
void parallel_configuration_iterator (ParConfigLIter pc_iter)
 set methodPCIter
 
ParConfigLIter parallel_configuration_iterator () const
 return methodPCIter
 
void run (ParLevLIter pl_iter)
 invoke set_communicators(pl_iter) prior to run()
 
void run ()
 orchestrate initialize/pre/core/post/finalize phases More...
 
void assign_rep (Iterator *iterator_rep, bool ref_count_incr=true)
 replaces existing letter with a new one More...
 
void iterated_model (const Model &model)
 set the iteratedModel (iterators and meta-iterators using a single model instance)
 
Modeliterated_model ()
 return the iteratedModel (iterators & meta-iterators using a single model instance)
 
ProblemDescDBproblem_description_db () const
 return the problem description database (probDescDB)
 
ParallelLibraryparallel_library () const
 return the parallel library (parallelLib)
 
void method_name (unsigned short m_name)
 set the method name to an enumeration value
 
unsigned short method_name () const
 return the method name via its native enumeration value
 
void method_string (const String &m_str)
 set the method name by string
 
String method_string () const
 return the method name by string
 
String method_enum_to_string (unsigned short method_name) const
 convert a method name enumeration value to a string
 
unsigned short method_string_to_enum (const String &method_name) const
 convert a method name string to an enumeration value
 
String submethod_enum_to_string (unsigned short submethod_name) const
 convert a method name enumeration value to a string
 
const String & method_id () const
 return the method identifier (methodId)
 
int maximum_evaluation_concurrency () const
 return the maximum evaluation concurrency supported by the iterator
 
void maximum_evaluation_concurrency (int max_conc)
 set the maximum evaluation concurrency supported by the iterator
 
void convergence_tolerance (Real conv_tol)
 set the method convergence tolerance (convergenceTol)
 
Real convergence_tolerance () const
 return the method convergence tolerance (convergenceTol)
 
void output_level (short out_lev)
 set the method output level (outputLevel)
 
short output_level () const
 return the method output level (outputLevel)
 
void summary_output (bool summary_output_flag)
 Set summary output control; true enables evaluation/results summary.
 
size_t num_final_solutions () const
 return the number of solutions to retain in best variables/response arrays
 
void num_final_solutions (size_t num_final)
 set the number of solutions to retain in best variables/response arrays
 
void active_set (const ActiveSet &set)
 set the default active set vector (for use with iterators that employ evaluate_parameter_sets())
 
const ActiveSetactive_set () const
 return the default active set vector (used by iterators that employ evaluate_parameter_sets())
 
void sub_iterator_flag (bool si_flag)
 set subIteratorFlag (and update summaryOutputFlag if needed)
 
void active_variable_mappings (const SizetArray &c_index1, const SizetArray &di_index1, const SizetArray &ds_index1, const SizetArray &dr_index1, const ShortArray &c_target2, const ShortArray &di_target2, const ShortArray &ds_target2, const ShortArray &dr_target2)
 set primaryA{CV,DIV,DRV}MapIndices, secondaryA{CV,DIV,DRV}MapTargets
 
bool is_null () const
 function to check iteratorRep (does this envelope contain a letter?)
 
Iteratoriterator_rep () const
 returns iteratorRep for access to derived class member functions that are not mapped to the top Iterator level
 
virtual void eval_tag_prefix (const String &eval_id_str)
 set the hierarchical eval ID tag prefix More...
 

Protected Member Functions

void sampling_increment ()
 increment to next in sequence of refinement samples
 
void pre_run ()
 generate LHS samples in non-VBD cases
 
void core_run ()
 perform the evaluate parameter sets portion of run More...
 
void post_run (std::ostream &s)
 generate statistics for LHS runs in non-VBD cases
 
void post_input ()
 read tabular data for post-run mode
 
void compute_pca (std::ostream &s)
 compute a principal components analysis on the sample set
 
void print_results (std::ostream &s)
 print the final statistics
 
void d_optimal_parameter_set (int previous_samples, int new_samples, RealMatrix &full_samples)
 generate a d-optimal parameter set, leaving the first previous_samples columns intact and adding new_samples new columns following them More...
 
void initial_increm_lhs_set (int new_samples, RealMatrix &full_samples, IntMatrix &full_ranks)
 Populate the first new_samples columns of allSamples with an LHS design and update the stored ranks.
 
void increm_lhs_parameter_set (int previous_samples, int new_samples, RealMatrix &full_samples, IntMatrix &all_ranks)
 generate a new batch that is Latin w.r.t. the previous samples
 
void store_ranks (const RealMatrix &sample_values, IntMatrix &sample_ranks)
 store the ranks of the last generated sample for continuous (based on sampleRanks) and calculate/store discrete ranks
 
void store_ranks (IntMatrix &full_ranks)
 store the combined ranks from sampleRanks to leading submatrix local cached ranks matrix
 
void combine_discrete_ranks (const RealMatrix &initial_values, const RealMatrix &increm_values)
 merge the discrete ranks into a submatrix of sampleRanks
 
- Protected Member Functions inherited from NonDSampling
 NonDSampling (ProblemDescDB &problem_db, Model &model)
 constructor More...
 
 NonDSampling (unsigned short method_name, Model &model, unsigned short sample_type, int samples, int seed, const String &rng, bool vary_pattern, short sampling_vars_mode)
 alternate constructor for sample generation and evaluation "on the fly" More...
 
 NonDSampling (unsigned short sample_type, int samples, int seed, const String &rng, const RealVector &lower_bnds, const RealVector &upper_bnds)
 alternate constructor for sample generation "on the fly" More...
 
 NonDSampling (unsigned short sample_type, int samples, int seed, const String &rng, const RealVector &means, const RealVector &std_devs, const RealVector &lower_bnds, const RealVector &upper_bnds, RealSymMatrix &correl)
 alternate constructor for sample generation of correlated normals "on the fly" More...
 
 ~NonDSampling ()
 destructor
 
void core_run ()
 
int num_samples () const
 
void sampling_reset (int min_samples, bool all_data_flag, bool stats_flag)
 resets number of samples and sampling flags More...
 
void sampling_reference (int samples_ref)
 set reference number of samples, which is a lower bound during reset
 
unsigned short sampling_scheme () const
 return sampleType
 
void vary_pattern (bool pattern_flag)
 set varyPattern
 
void get_parameter_sets (Model &model)
 Uses lhsDriver to generate a set of samples from the distributions/bounds defined in the incoming model. More...
 
void get_parameter_sets (const RealVector &lower_bnds, const RealVector &upper_bnds)
 Uses lhsDriver to generate a set of uniform samples over lower_bnds/upper_bnds. More...
 
void get_parameter_sets (const RealVector &means, const RealVector &std_devs, const RealVector &lower_bnds, const RealVector &upper_bnds, RealSymMatrix &correl)
 Uses lhsDriver to generate a set of normal samples. More...
 
void update_model_from_sample (Model &model, const Real *sample_vars)
 Override default update of continuous vars only.
 
void sample_to_variables (const Real *sample_vars, Variables &vars)
 override default mapping of continuous variables only
 
void variables_to_sample (const Variables &vars, Real *sample_vars)
 
void initialize_lhs (bool write_message)
 increments numLHSRuns, sets random seed, and initializes lhsDriver
 
void view_design_counts (const Model &model, size_t &num_cdv, size_t &num_ddiv, size_t &num_ddsv, size_t &num_ddrv) const
 compute sampled subsets (all, active, uncertain) within all variables (acv/adiv/adrv) from samplingVarsMode and model More...
 
void view_aleatory_uncertain_counts (const Model &model, size_t &num_cauv, size_t &num_dauiv, size_t &num_dausv, size_t &num_daurv) const
 compute sampled subsets (all, active, uncertain) within all variables (acv/adiv/adrv) from samplingVarsMode and model More...
 
void view_epistemic_uncertain_counts (const Model &model, size_t &num_ceuv, size_t &num_deuiv, size_t &num_deusv, size_t &num_deurv) const
 compute sampled subsets (all, active, uncertain) within all variables (acv/adiv/adrv) from samplingVarsMode and model More...
 
void view_uncertain_counts (const Model &model, size_t &num_cuv, size_t &num_duiv, size_t &num_dusv, size_t &num_durv) const
 compute sampled subsets (all, active, uncertain) within all variables (acv/adiv/adrv) from samplingVarsMode and model More...
 
void view_state_counts (const Model &model, size_t &num_csv, size_t &num_dsiv, size_t &num_dssv, size_t &num_dsrv) const
 compute sampled subsets (all, active, uncertain) within all variables (acv/adiv/adrv) from samplingVarsMode and model
 
void mode_counts (const Model &model, size_t &cv_start, size_t &num_cv, size_t &div_start, size_t &num_div, size_t &dsv_start, size_t &num_dsv, size_t &drv_start, size_t &num_drv) const
 compute sampled subsets (all, active, uncertain) within all variables (acv/adiv/adrv) from samplingVarsMode and model More...
 
void get_lhs_samples (const Model &model, int num_samples, RealMatrix &design_matrix)
 Uses lhsDriver to generate a set of samples from the distributions/bounds defined in the incoming model and populates the specified design matrix. More...
 
- Protected Member Functions inherited from NonD
 NonD (ProblemDescDB &problem_db, Model &model)
 constructor
 
 NonD (unsigned short method_name, Model &model)
 alternate constructor for sample generation and evaluation "on the fly"
 
 NonD (unsigned short method_name, const RealVector &lower_bnds, const RealVector &upper_bnds)
 alternate constructor for sample generation "on the fly"
 
 ~NonD ()
 destructor
 
void derived_set_communicators (ParLevLIter pl_iter)
 derived class contributions to setting the communicators associated with this Iterator instance
 
void initialize_run ()
 utility function to perform common operations prior to pre_run(); typically memory initialization; setting of instance pointers More...
 
void finalize_run ()
 utility function to perform common operations following post_run(); deallocation and resetting of instance pointers More...
 
const Responseresponse_results () const
 return the final statistics from the nondeterministic iteration
 
void response_results_active_set (const ActiveSet &set)
 set the active set within finalStatistics
 
virtual void initialize_response_covariance ()
 initializes respCovariance
 
virtual void initialize_final_statistics ()
 initializes finalStatistics for storing NonD final results More...
 
void size ()
 Size local variables.
 
int generate_system_seed ()
 create a system-generated unique seed (when a seed is unspecified)
 
void initialize_final_statistics_gradients ()
 initializes finalStatistics::functionGradients
 
void update_aleatory_final_statistics ()
 update finalStatistics::functionValues from momentStats and computed{Prob,Rel,GenRel,Resp}Levels
 
void update_system_final_statistics ()
 update system metrics from component metrics within finalStatistics
 
void update_system_final_statistics_gradients ()
 update finalStatistics::functionGradients
 
void initialize_level_mappings ()
 size computed{Resp,Prob,Rel,GenRel}Levels
 
void compute_densities (const RealRealPairArray &min_max_fns, bool prob_refinement=false, bool all_levels_computed=false)
 compute the PDF bins from the CDF/CCDF values and store in computedPDF{Abscissas,Ordinates} More...
 
void print_densities (std::ostream &s) const
 output the PDFs reflected in computedPDF{Abscissas,Ordinates} using default qoi_type and pdf_labels
 
void print_densities (std::ostream &s, String qoi_type, const StringArray &pdf_labels) const
 output the PDFs reflected in computedPDF{Abscissas,Ordinates}
 
void print_system_mappings (std::ostream &s) const
 print system series/parallel mappings for response levels
 
void transform_model (Model &x_model, Model &u_model, bool truncated_bounds=false, Real bound=10.)
 recast x_model from x-space to u-space to create u_model
 
void construct_lhs (Iterator &u_space_sampler, Model &u_model, unsigned short sample_type, int num_samples, int seed, const String &rng, bool vary_pattern, short sampling_vars_mode=ACTIVE)
 assign a NonDLHSSampling instance within u_space_sampler
 
void archive_allocate_mappings ()
 allocate results array storage for distribution mappings
 
void archive_from_resp (size_t fn_index)
 archive the mappings from specified response levels for specified fn
 
void archive_to_resp (size_t fn_index)
 archive the mappings to computed response levels for specified fn
 
void archive_allocate_pdf ()
 allocate results array storage for pdf histograms
 
void archive_pdf (size_t fn_index)
 archive a single pdf histogram for specified function
 
- Protected Member Functions inherited from Analyzer
 Analyzer ()
 default constructor
 
 Analyzer (ProblemDescDB &problem_db, Model &model)
 standard constructor
 
 Analyzer (unsigned short method_name, Model &model)
 alternate constructor for instantiations "on the fly" with a Model
 
 Analyzer (unsigned short method_name)
 alternate constructor for instantiations "on the fly" without a Model
 
 ~Analyzer ()
 destructor
 
virtual void update_model_from_variables (Model &model, const Variables &vars)
 update model's current variables with data from vars
 
void update_from_model (const Model &model)
 set inherited data attributes based on extractions from incoming model
 
void pre_output ()
 
const Modelalgorithm_space_model () const
 
const Variablesvariables_results () const
 return a single final iterator solution (variables)
 
const VariablesArray & variables_array_results ()
 return multiple final iterator solutions (variables). This should only be used if returns_multiple_points() returns true.
 
const ResponseArray & response_array_results ()
 return multiple final iterator solutions (response). This should only be used if returns_multiple_points() returns true.
 
bool compact_mode () const
 returns Analyzer::compactMode
 
bool returns_multiple_points () const
 indicates if this iterator returns multiple final points. Default return is false. Override to return true if appropriate.
 
void evaluate_parameter_sets (Model &model, bool log_resp_flag, bool log_best_flag)
 perform function evaluations to map parameter sets (allVariables) into response sets (allResponses) More...
 
void variance_based_decomp (int ncont, int ndiscint, int ndiscreal, int num_samples)
 
void read_variables_responses (int num_evals, size_t num_vars)
 convenience function for reading variables/responses (used in derived classes post_input) More...
 
void print_sobol_indices (std::ostream &s) const
 Printing of VBD results. More...
 
void samples_to_variables_array (const RealMatrix &sample_matrix, VariablesArray &vars_array)
 convert samples array to variables array; e.g., allSamples to allVariables
 
void variables_array_to_samples (const VariablesArray &vars_array, RealMatrix &sample_matrix)
 convert variables array to samples array; e.g., allVariables to allSamples
 
- Protected Member Functions inherited from Iterator
 Iterator (BaseConstructor, ProblemDescDB &problem_db)
 constructor initializes the base class part of letter classes (BaseConstructor overloading avoids infinite recursion in the derived class constructors - Coplien, p. 139) More...
 
 Iterator (NoDBBaseConstructor, unsigned short method_name, Model &model)
 alternate constructor for base iterator classes constructed on the fly More...
 
 Iterator (NoDBBaseConstructor, unsigned short method_name)
 alternate constructor for base iterator classes constructed on the fly More...
 
virtual void derived_init_communicators (ParLevLIter pl_iter)
 derived class contributions to initializing the communicators associated with this Iterator instance
 
virtual const VariablesArray & initial_points () const
 gets the multiple initial points for this iterator. This will only be meaningful after a call to initial_points mutator.
 
StrStrSizet run_identifier () const
 get the unique run identifier based on method name, id, and number of executions
 

Static Protected Member Functions

static bool rank_sort (const int &x, const int &y)
 sort algorithm to compute ranks for rank correlations
 

Private Attributes

size_t numResponseFunctions
 number of response functions; used to distinguish NonD from opt/NLS usage
 
IntVector refineSamples
 list of refinement sample batch sizes
 
bool dOptimal
 whether to generate d-optimal point sets
 
Real oversampleRatio
 oversampling ratio for d-optimal candidate set generation (default 10.0)
 
bool varBasedDecompFlag
 flags computation of variance-based decomposition indices
 
bool pcaFlag
 flag to specify the calculation of principal components
 
Real percentVarianceExplained
 Threshold to keep number of principal components that explain this much variance.
 

Static Private Attributes

static RealArray rawData
 static data used by static rank_sort() fn
 

Additional Inherited Members

- Static Public Member Functions inherited from NonDSampling
static void print_moments (std::ostream &s, const RealMatrix &moment_stats, const RealMatrix moment_cis, String qoi_type, const StringArray &moment_labels, bool print_cis)
 core print moments that can be called without object
 
static void compute_moments (const RealMatrix &samples, RealMatrix &moment_stats)
 core compute moments that can be called without object
 
static int compute_wilks_sample_size (unsigned short order, Real alpha, Real beta, bool twosided=false)
 calculates the number of samples using the Wilks formula Static for now so I can test without instantiating a NonDSampling object - RWH
 
- Protected Attributes inherited from NonDSampling
const int seedSpec
 the user seed specification (default is 0)
 
int randomSeed
 the current seed
 
const int samplesSpec
 initial specification of number of samples
 
int samplesRef
 reference number of samples updated for refinement
 
int numSamples
 the current number of samples to evaluate
 
String rngName
 name of the random number generator
 
unsigned short sampleType
 
 the sample type: default, random, lhs,

< incremental random, or incremental lhs

 
bool wilksFlag
 
int samplesIncrement
 flags use of Wilks formula to calculate num samples More...
 
Pecos::LHSDriver lhsDriver
 the C++ wrapper for the F90 LHS library
 
bool statsFlag
 flags computation/output of statistics
 
bool allDataFlag
 
 flags update of allResponses

< (allVariables or allSamples already defined)

 
short samplingVarsMode
 the sampling mode: ALEATORY_UNCERTAIN{,_UNIFORM}, EPISTEMIC_UNCERTAIN{,_UNIFORM}, UNCERTAIN{,_UNIFORM}, ACTIVE{,_UNIFORM}, or ALL{,_UNIFORM}. This is a secondary control on top of the variables view that allows sampling over subsets of variables that may differ from the view.
 
short sampleRanksMode
 mode for input/output of LHS sample ranks: IGNORE_RANKS, GET_RANKS, SET_RANKS, or SET_GET_RANKS
 
bool varyPattern
 flag for generating a sequence of seed values within multiple get_parameter_sets() calls so that these executions (e.g., for SBO/SBNLS) are not repeated, but are still repeatable
 
RealMatrix sampleRanks
 data structure to hold the sample ranks
 
SensAnalysisGlobal nonDSampCorr
 initialize statistical post processing
 
bool backfillFlag
 flags whether to use backfill to enforce uniqueness of discrete LHS samples
 
RealRealPairArray extremeValues
 Minimum and maximum values of response functions for epistemic calculations (calculated in compute_intervals()),.
 
- Static Protected Attributes inherited from NonD
static NonDnondInstance
 pointer to the active object instance used within static evaluator functions in order to avoid the need for static data
 

Detailed Description

Performs LHS and Monte Carlo sampling for uncertainty quantification.

The Latin Hypercube Sampling (LHS) package from Sandia Albuquerque's Risk and Reliability organization provides comprehensive capabilities for Monte Carlo and Latin Hypercube sampling within a broad array of user-specified probabilistic parameter distributions. It enforces user-specified rank correlations through use of a mixing routine. The NonDLHSSampling class provides a C++ wrapper for the LHS library and is used for performing forward propagations of parameter uncertainties into response statistics.

Batch generation options, including D-Optimal and incremental LHS are provided.

The incremental LHS sampling capability allows one to supplement an initial sample of size n to size 2n while maintaining the correct stratification of the 2n samples and also maintaining the specified correlation structure. The incremental version of LHS will return a sample of size n, which when combined with the original sample of size n, allows one to double the size of the sample.

Constructor & Destructor Documentation

NonDLHSSampling ( ProblemDescDB problem_db,
Model model 
)

standard constructor

This constructor is called for a standard letter-envelope iterator instantiation. In this case, set_db_list_nodes has been called and probDescDB can be queried for settings from the method specification.

References NonDLHSSampling::dOptimal, Model::num_primary_fns(), NonDLHSSampling::numResponseFunctions, Iterator::outputLevel, Model::primary_fn_type(), NonDLHSSampling::refineSamples, and NonDSampling::sampleType.

NonDLHSSampling ( Model model,
unsigned short  sample_type,
int  samples,
int  seed,
const String &  rng,
bool  vary_pattern = true,
short  sampling_vars_mode = ACTIVE 
)

alternate constructor for sample generation and evaluation "on the fly"

This alternate constructor is used for generation and evaluation of Model-based sample sets. A set_db_list_nodes has not been performed so required data must be passed through the constructor. It's purpose is to avoid the need for a separate LHS specification within methods that use LHS sampling.

NonDLHSSampling ( unsigned short  sample_type,
int  samples,
int  seed,
const String &  rng,
const RealVector &  lower_bnds,
const RealVector &  upper_bnds 
)

alternate constructor for uniform sample generation "on the fly"

This alternate constructor is used by ConcurrentStrategy for generation of uniform, uncorrelated sample sets. It is not a letter-envelope instantiation and a set_db_list_nodes has not been performed. It is called with all needed data passed through the constructor and is designed to allow more flexibility in variables set definition (i.e., relax connection to a variables specification and allow sampling over parameter sets such as multiobjective weights). In this case, a Model is not used and the object must only be used for sample generation (no evaluation).

References NonDSampling::get_parameter_sets().

NonDLHSSampling ( unsigned short  sample_type,
int  samples,
int  seed,
const String &  rng,
const RealVector &  means,
const RealVector &  std_devs,
const RealVector &  lower_bnds,
const RealVector &  upper_bnds,
RealSymMatrix &  correl 
)

alternate constructor for sample generation of correlated normals "on the fly"

This alternate constructor is used to generate correlated normal sample sets. It is not a letter-envelope instantiation and a set_db_list_nodes has not been performed. It is called with all needed data passed through the constructor. In this case, a Model is not used and the object must only be used for sample generation (no evaluation).

References NonDSampling::get_parameter_sets().

Member Function Documentation

void core_run ( )
protectedvirtual

perform the evaluate parameter sets portion of run

Loop over the set of samples and compute responses. Compute statistics on the set of responses if statsFlag is set.

Reimplemented from Iterator.

References NonDSampling::allDataFlag, Analyzer::evaluate_parameter_sets(), Iterator::iteratedModel, Analyzer::numContinuousVars, Analyzer::numDiscreteIntVars, Analyzer::numDiscreteRealVars, NonDLHSSampling::numResponseFunctions, NonDSampling::numSamples, NonDSampling::statsFlag, NonDLHSSampling::varBasedDecompFlag, and Analyzer::variance_based_decomp().

void d_optimal_parameter_set ( int  previous_samples,
int  new_samples,
RealMatrix &  full_samples 
)
protected

generate a d-optimal parameter set, leaving the first previous_samples columns intact and adding new_samples new columns following them

For now, when this function is called, numSamples is the number of new samples to generate.

References Model::aleatory_distribution_parameters(), NonDSampling::get_lhs_samples(), NonD::initialize_random_variables(), Iterator::iteratedModel, NonDSampling::mode_counts(), NonD::natafTransform, NonDLHSSampling::oversampleRatio, and NonDSampling::transform_samples().

Referenced by NonDLHSSampling::pre_run().


The documentation for this class was generated from the following files: