1#@ *: DakotaConfig=HAVE_NOWPAC
2
3environment,
4	  method_pointer = 'OPTIM'
5
6###########################
7# begin opt specification #
8###########################
9method,
10	id_method = 'OPTIM'
11	model_pointer = 'OPTIM_M'
12	snowpac
13	  seed = 25041981
14#	  max_iterations = 500
15#	  convergence_tolerance = 1e-4
16	  max_function_evaluations = 250
17	  trust_region
18	    initial_size = 0.10
19	    minimum_size = 1.0e-6
20	    contract_threshold = 0.25
21	    expand_threshold   = 0.75
22	    contraction_factor = 0.50
23	    expansion_factor   = 1.50
24#	  output debug
25
26model,
27	id_model = 'OPTIM_M'
28	nested
29	  variables_pointer  = 'OPTIM_V'
30	  sub_method_pointer = 'UQ'
31	  responses_pointer  = 'OPTIM_R'
32	  primary_response_mapping   = 0. 1.e+6 0. 1.e+6 0. 1.e+6
33	  secondary_response_mapping = 1. 0. 0. 0. 0. 0.
34	  			       0. 0. 1. 0. 0. 0.
35                                       0. 0. 0. 0. 1. 0.
36
37variables,
38	id_variables = 'OPTIM_V'
39	continuous_design = 3
40	  initial_point    3*.5
41	  lower_bounds     3*-1.
42	  upper_bounds     3* 1.
43
44responses,
45	id_responses = 'OPTIM_R'
46	objective_functions = 1
47	nonlinear_inequality_constraints = 3
48	nonlinear_inequality_upper_bounds = .02 .2 .02
49	no_gradients
50	no_hessians
51
52
53##########################
54# begin UQ specification #
55##########################
56method,
57	id_method = 'UQ'
58	model_pointer = 'HIERARCH'
59        multilevel_sampling
60#	  sample_type lhs
61	  pilot_samples = 20 seed = 1237
62	  max_iterations = 2
63#	  convergence_tolerance =    .1
64	  convergence_tolerance =   .01
65#	  convergence_tolerance =  .001
66#	  convergence_tolerance = .0001
67	output silent
68	final_moments central
69
70model,
71	id_model = 'HIERARCH'
72	surrogate hierarchical
73	  ordered_model_fidelities = 'SIM1'
74
75model,
76	id_model = 'SIM1'
77	variables_pointer  = 'UQ_V'
78	interface_pointer  = 'UQ_I'
79	responses_pointer  = 'UQ_R'
80	simulation
81	  solution_level_control = 'mesh_size'
82# dense matrices --> direct solvers --> cubic growth in expense
83	  solution_level_cost = 1. 8. 64. 512. 4096.
84
85variables,
86	id_variables = 'UQ_V'
87	continuous_design = 3	 # number of cosine terms in soln expansion
88	uniform_uncertain = 6	 # number of cosine terms in soln expansion
89	  lower_bounds      =  6*-1.
90	  upper_bounds      =  6* 1.
91	discrete_state_set
92	  integer = 1 # even poly degree for odd mesh pts including boundary
93      	  set_values = 4 8 16 32 64
94      	  descriptors = 'mesh_size'
95    	real = 4
96      	  elements_per_variable = 2 2 1 1
97      	  set_values = 0.1 1 # field mean
98                       0.5 4 # field std dev
99                       1     # kernel order
100                       0.2   # kernel length
101      	  descriptors = 'field_mean' 'field_std_dev' 'kernel_order' 'kernel_length'
102      	  initial_state = 1 4 1 0.2
103	string = 2  # form of kernel (refer to Yang,Lei,Baker,Lin)
104      	  elements_per_variable = 2 2
105	  initial_state = 'cosine' 'off'
106      	  set_values = 'cosine' 'exponential'
107                       'off' 'on'
108      	  descriptors = 'kernel_type' 'positivity'
109
110interface,
111	id_interface = 'UQ_I'
112	direct
113	  analysis_driver = 'steady_state_diffusion_1d'
114	  deactivate evaluation_cache restart_file
115
116responses,
117	id_responses = 'UQ_R'
118	response_functions = 3	 # number of soln pts excl. +/-5% (symmetric)
119	no_gradients
120	no_hessians
121