1#@ *: DakotaConfig=HAVE_NPSOL
2# DAKOTA INPUT FILE: dakota_trsbouu4_tbch.in
3#
4# OUU with surrogates: surrogate containing nested containing surrogate
5#
6
7environment,
8	method_pointer = 'SBLO'
9
10method,
11	id_method = 'SBLO'
12	surrogate_based_local
13	model_pointer = 'OPTIM_M'
14	  approx_method_pointer = 'OPTIM'
15	  max_iterations = 50
16	  soft_convergence_limit = 2
17	  trust_region
18	    initial_size = 0.4
19	    contraction_factor = 0.5
20	    expansion_factor   = 1.50
21
22###########################
23# begin opt specification #
24###########################
25method,
26	id_method = 'OPTIM'
27        npsol_sqp
28	  convergence_tolerance = 1.e-8
29
30model,
31	id_model = 'OPTIM_M'
32	surrogate global
33	  variables_pointer = 'OPTIM_V'
34	  responses_pointer = 'OPTIM_R'
35	  polynomial quadratic
36	  dace_method_pointer = 'OPTIM_DACE'
37
38variables,
39	id_variables = 'OPTIM_V'
40	continuous_design = 2
41	  initial_point    1.8    1.0
42	  upper_bounds     2.164  4.0
43	  lower_bounds     1.5    0.0
44	  descriptors      'd1'   'd2'
45
46responses,
47# minimize f + p_fail_r1 + p_fail_r2
48# s.t.     g_1/2/3 <= 0
49#          p_fail_r3 <= 0.1
50# NOTE: This specifies the TOTAL RESPONSE for the optimization,
51#       which is a combination of nested & interface responses.
52	id_responses = 'OPTIM_R'
53	objective_functions = 1
54	nonlinear_inequality_constraints = 4
55	nonlinear_inequality_upper_bounds = 0.0 0.0 0.0 1.6e+5
56	analytic_gradients
57#	numerical_gradients
58#	  method_source dakota
59#	  interval_type central
60#	  fd_gradient_step_size = 1.e-5
61	no_hessians
62
63############################
64# begin DACE specification #
65############################
66method,
67	id_method = 'OPTIM_DACE'
68	model_pointer = 'OPTIM_DACE_M'
69	sampling
70	  sample_type lhs seed = 12347
71#	  fixed_seed					#s1
72#	dace lhs seed = 5
73#	  samples = 50 symbols = 50
74
75model,
76	id_model = 'OPTIM_DACE_M'
77	nested
78	  variables_pointer  = 'OPTIM_V'
79	  sub_method_pointer = 'UQ'
80	  optional_interface_pointer  = 'OPTIONAL_I'
81	  optional_interface_responses_pointer = 'OPTIONAL_I_R'
82	  responses_pointer  = 'OPTIM_R_NO_GRAD'
83	  primary_response_mapping   = 0. 0. 1. 0. 0. 0. 0. 0. 1.
84	  secondary_response_mapping = 0. 0. 0. 1. 3. 0. 0. 0. 0.
85
86responses,
87	id_responses = 'OPTIM_R_NO_GRAD'
88	objective_functions = 1
89	nonlinear_inequality_constraints = 4
90	no_gradients
91	no_hessians
92
93##########################################
94# begin optional interface specification #
95##########################################
96interface,
97	id_interface = 'OPTIONAL_I'
98	direct
99	  analysis_driver = 'cyl_head'
100
101responses,
102# NOTE: This is the response set from the optional interface.
103	id_responses = 'OPTIONAL_I_R'
104	objective_functions = 1
105	nonlinear_inequality_constraints = 3
106	no_gradients
107	no_hessians
108
109##########################
110# begin UQ specification #
111##########################
112method,
113	id_method = 'UQ'
114	model_pointer = 'UQ_M'
115        sampling
116	  output quiet
117	  samples = 5000 sample_type lhs
118	  seed = 12347
119#	  fixed_seed					#s1
120	  response_levels = 4.2e+11 1.6e+05 7.5e+05
121	  complementary distribution
122
123model,
124	id_model = 'UQ_M'
125	surrogate global
126	  variables_pointer = 'UQ_V'
127	  responses_pointer = 'UQ_R'
128	  dace_method_pointer = 'UQ_DACE'
129	  polynomial quadratic
130
131variables,
132	id_variables = 'UQ_V'
133	continuous_design = 2
134	normal_uncertain = 2
135	  means             =  248.89, 593.33
136	  std_deviations    =   12.4,   29.7
137# 3 sigma bounds for DDACE:
138#	  lower_bounds =  211.7,  504.2
139#	  upper_bounds =  286.1,  682.4
140	uniform_uncertain = 2
141	  lower_bounds =  199.3,  474.63
142	  upper_bounds =  298.5,  712.
143	weibull_uncertain = 2
144	  alphas       =   12.,    30.
145	  betas        =  250.,   590.
146#	  lower_bounds =  225.,   540.
147#	  upper_bounds =  275.,   640.
148
149responses,
150	id_responses = 'UQ_R'
151	response_functions = 3
152	no_gradients
153	no_hessians
154
155############################
156# begin DACE specification #
157############################
158method,
159	id_method = 'UQ_DACE'
160	model_pointer = 'UQ_DACE_M'
161        sampling
162	  sample_type lhs #samples = 50
163	  seed = 12347
164#	  fixed_seed					#s1
165
166model,
167	id_model = 'UQ_DACE_M'
168	single
169	  variables_pointer = 'UQ_DACE_V'
170	  interface_pointer = 'UQ_ACTUAL'
171	  responses_pointer = 'UQ_R'
172
173variables,
174	id_variables = 'UQ_DACE_V'
175	active all
176	continuous_design = 2
177	normal_uncertain = 2
178	  means             =  248.89, 593.33
179	  std_deviations    =   12.4,   29.7
180# 3 sigma bounds for DDACE:
181#	  lower_bounds =  211.7,  504.2
182#	  upper_bounds =  286.1,  682.4
183	uniform_uncertain = 2
184	  lower_bounds =  199.3,  474.63
185	  upper_bounds =  298.5,  712.
186	weibull_uncertain = 2
187	  alphas       =   12.,    30.
188	  betas        =  250.,   590.
189#	  lower_bounds =  225.,   540.
190#	  upper_bounds =  275.,   640.
191
192interface,
193	id_interface = 'UQ_ACTUAL'
194	direct
195	  analysis_driver=        'text_book_ouu'
196	  deactivate evaluation_cache restart_file
197