1#@ *: DakotaConfig=HAVE_C3
2#@ s*: TimeoutAbsolute=3600
3#@ s*: TimeoutDelay=1200
4
5#environment,
6#	graphics
7
8method,
9	model_pointer = 'HIERARCH'
10	multilevel_function_train
11	  pilot_samples = 100 60 50 40 30			#s0,#s1
12#	  collocation_ratio = 1.75				#s2,#s3,#s4
13	  allocation_control rank_sampling			#s0,#s2,#s3,#s4
14#	  allocation_control estimator_variance			#s1
15#	  start_rank_sequence = 3 2 2 1 1			#s1,#s3
16	  start_rank_sequence = 2 2 2 2 2			#s0,#s2,#s4
17	  adapt_rank   kick_rank  = 1  max_rank  = 5		#s0,#s2,#s4
18#	  adapt_order  kick_order = 1  max_order = 4		#s3,#s4
19#
20	  start_order_sequence = 2 2 2 2 2
21#	  convergence_tolerance = 1e-4				#s1
22          solver_tolerance = 1e-10  rounding_tolerance = 1e-8
23	  seed = 1237						#s0,#s1,#s2
24#	  seed_sequence = 12345 23456 34567 45678 56789   	#s3,#s4
25	  output quiet
26
27model,
28	id_model = 'HIERARCH'
29	surrogate hierarchical
30	  ordered_model_fidelities = 'SIM1'
31	  correction additive zeroth_order
32
33model,
34	id_model = 'SIM1'
35	simulation
36	  solution_level_control = 'mesh_size'
37# dense matrices --> direct solvers --> cubic growth in expense
38	  solution_level_cost = 1. 8. 64. 512. 4096.
39
40variables,
41	uniform_uncertain = 9	 # number of cosine terms in soln expansion
42	  lower_bounds      =  9*-1.
43	  upper_bounds      =  9* 1.
44	discrete_state_set
45	  integer = 1 # even poly degree for odd mesh pts including boundary
46      	  set_values = 4 8 16 32 64
47      	  descriptors = 'mesh_size'
48    	real = 4
49      	  elements_per_variable = 2 2 1 1
50      	  set_values = 0.1 1 # field mean
51                       0.5 4 # field std dev
52                       1     # kernel order
53                       0.2   # kernel length
54      	  descriptors = 'field_mean' 'field_std_dev' 'kernel_order' 'kernel_length'
55      	  initial_state = 1 4 1 0.2
56	string = 2  # form of kernel (refer to Yang,Lei,Baker,Lin)
57      	  elements_per_variable = 2 2
58	  initial_state = 'cosine' 'off'
59      	  set_values = 'cosine' 'exponential'
60                       'off' 'on'
61      	  descriptors = 'kernel_type' 'positivity'
62
63interface,
64	direct
65	  analysis_driver = 'steady_state_diffusion_1d'
66
67responses,
68	response_functions = 3	 # number of soln pts excl. +/-5% (symmetric)
69	no_gradients
70	no_hessians
71