1# Tests for D-optimal sampling
2# Currently only support aleatory variables
3# With Leja option, only support continuous aleatory variables
4# test 0: continuous aleatory
5# test 1: discrete aleatory
6# test 2: continuous aleatory with Leja option
7environment,
8	  tabular_data
9
10method,
11  sampling
12    sample_type random
13    samples = 100                     #s0,#s1
14#   samples = 20                      #s2
15      d_optimal
16# Default candidate designs = 100
17#       candidate_designs 50          #s1
18#       leja_oversample_ratio = 10.0  #s2
19    seed 8272
20
21## ERROR: Want samples greater, but get this error on s1:
22#  terminate called after throwing an instance of 'std::runtime_error'
23#    what():  enrich_samples: The basis matrix has rank 21 less than num_total_samples 30. Try increasing the degree of the basis.
24
25
26variables,
27# (default is active uncertain)
28# active design
29# active state
30
31## DESIGN
32#	continuous_design = 1
33#          upper_bounds = 20.
34#          lower_bounds = 10.
35## ERROR:
36#   Distribution required, not found for DiscDesRange2
37#   Error: code 1 returned from lhs_prep in LHSDriver.
38#	discrete_design_range = 1
39#          upper_bounds =  7
40#          lower_bounds = -3
41## ERROR (no discrete histogram):
42##   Error: RandomVariable type 0 not available.
43#	discrete_design_set
44#	  integer = 1
45#            elements = 3 4 5
46#	  string = 1
47#            elements = 'dog' 'fleas' 'has' 'my'
48#	  real = 1
49#            elements = 30. 40. 50.
50
51## ALEATORY UNCERTAIN
52  normal_uncertain = 1                            #s0,#s2
53    means = 69.0       	                          #s0,#s2
54    std_deviations = 1.0  	                      #s0,#s2
55	lognormal_uncertain = 1		                      #s0,#s2
56    means             =  1.	                      #s0,#s2
57    std_deviations    =  0.5                      #s0,#s2
58	uniform_uncertain = 1                           #s0,#s2
59	  lower_bounds      = 200.	                    #s0,#s2
60    upper_bounds = 300.		                        #s0,#s2
61 	loguniform_uncertain = 1	                      #s0,#s2
62    lower_bounds      =  0.1	                    #s0,#s2
63    upper_bounds      =  10.	                    #s0,#s2
64 	triangular_uncertain = 1	                      #s0,#s2
65    modes             =  0.  	                    #s0,#s2
66    lower_bounds      = -2.	                      #s0,#s2
67	  upper_bounds      =  2.  	                    #s0,#s2
68	exponential_uncertain = 1	                      #s0,#s2
69	  betas             =  2.	                      #s0,#s2
70	beta_uncertain = 1		                          #s0,#s2
71	  alphas            =  1.5  	                  #s0,#s2
72	  betas             =  2. 	                    #s0,#s2
73	  lower_bounds      = -2.	                      #s0,#s2
74	  upper_bounds      =  2.	                      #s0,#s2
75	gamma_uncertain = 1                             #s0,#s2
76	  alphas            =  2.5  	                  #s0,#s2
77	  betas             =  2.  	                    #s0,#s2
78	gumbel_uncertain = 1		                        #s0,#s2
79	  alphas            = 1.5 	                    #s0,#s2
80	  betas             = 2.5	                      #s0,#s2
81	frechet_uncertain = 1		                        #s0,#s2
82	  alphas            = 2.5 	                    #s0,#s2
83	  betas             = 3.5 	                    #s0,#s2
84	weibull_uncertain = 1		                        #s0,#s2
85	  alphas            =   12.                     #s0,#s2
86	  betas             =  250.	                    #s0,#s2
87	histogram_bin_uncertain = 1	                    #s0,#s2
88	  abscissas   =  5  8 10 	                      #s0,#s2
89	  counts      = 17 21  0	                      #s0,#s2
90# poisson_uncertain = 1		                        #s1
91#   lambdas           =  0.2   	                  #s1
92# binomial_uncertain = 1		                      #s1
93#   prob_per_trial    =  0.3  	                  #s1
94#   num_trials        =  10  	                    #s1
95# negative_binomial_uncertain = 1	                #s1
96#   prob_per_trial    =  0.3  	                  #s1
97#   num_trials        =  10    	                  #s1
98# geometric_uncertain = 1		                      #s1
99#   prob_per_trial    = 0.1                       #s1
100#	hypergeometric_uncertain = 1	                  #s1
101#         total_population    = 50 	              #s1
102#         selected_population = 20	              #s1
103#         num_drawn           =  5 	              #s1
104#	histogram_point_uncertain                       #s1
105#   integer = 1                                   #s1
106#	    abscissas   = 3 4 5                         #s1
107#	    counts      = 1 2 1                         #s1
108#   string = 1                                    #s1
109#	    abscissas   = 's3' 's4' 's5'                #s1
110#	    counts      = 1 2 1                         #s1
111#   real = 1                                      #s1
112#	    abscissas   = 3.2 4.3 5.4                   #s1
113#	    counts      = 1 2 1                         #s1
114
115## EPISTEMIC UNCERTAIN
116#	continuous_interval_uncertain = 1
117#	  interval_probabilities = .4 .6
118#	  lower_bounds = 2. 4.
119#	  upper_bounds = 5. 8.
120## ERROR (no discrete histogram)
121#	discrete_interval_uncertain = 1
122#	  interval_probabilities = .4 .6
123#	  lower_bounds = 2 4
124#	  upper_bounds = 5 8
125#	discrete_uncertain_set
126#	  integer = 1
127#	    set_probabilities = .4 .3 .2 .1
128#	    set_values        =  2  4  5  8
129## TODO: string
130#	  real = 1
131#	    set_probabilities =  .1  .2  .3  .4
132#	    set_values        = 3.1 5.2 6.3 7.4
133
134## STATE
135#	continuous_state = 1
136#          upper_bounds = 20.
137#          lower_bounds = 10.
138## ERROR:
139#   Distribution required, not found for DiscDesRange2
140#   Error: code 1 returned from lhs_prep in LHSDriver.
141#	discrete_state_range = 1
142#          upper_bounds =  7
143#          lower_bounds = -3
144## ERROR (no discrete histogram):
145##   Error: RandomVariable type 0 not available.
146#	discrete_state_set
147#	  integer = 1
148#            set_values = 3 4 5
149## TODO: string
150#	  real = 1
151#            set_values = 30. 40. 50.
152
153interface,
154	direct
155	  analysis_driver = 'text_book'
156	  deactivate evaluation_cache restart_file
157
158responses,
159	response_functions = 1
160	no_gradients
161	no_hessians
162