1#@ *: DakotaConfig=HAVE_QUESO
2#@ s*: TimeoutDelay=360
3
4method,
5	model_pointer = 'HIERARCH'
6	bayes_calibration queso
7	  chain_samples = 40000 seed = 348
8	  delayed_rejection
9	  emulator
10	    pce expansion_order = 2				#s0
11#	    mf_pce   expansion_order_sequence = 3 2 2 1 1	#s1
12#	    mf_sc  sparse_grid_level_sequence = 3 2 2 1 1	#s2
13#	    ml_pce   expansion_order_sequence = 2		#s3
14#	      pilot_samples = 100 80 60 40 20			#s3
15#	      convergence_tolerance = .005			#s3
16#	      max_iterations = 5				#s3
17	      collocation_ratio = 0.9				#s0,#s1,#s3
18#	    mf_pce sparse_grid_level_sequence = 0     		#s4,#s5
19#	    mf_sc  sparse_grid_level_sequence = 0     		#s6,#s7,#s8,#s9
20#	      hierarchical		      			#s8,#s9
21#	      allocation_control greedy				#s4,#s6,#s8
22#	      p_refinement dimension_adaptive generalized	#s4,#s5,#s6,#s7,#s8,#s9
23#	      convergence_tolerance 1.e-4			#s4,#s5,#s6,#s7,#s8,#s9
24	  proposal_covariance
25	    derivatives
26	      update_period = 40
27#	  logit_transform         # default is off
28	  # TODO: debug negative information gain on tests s0--s3
29##	  posterior_stats kl_divergence
30
31model,
32	id_model = 'HIERARCH'
33	surrogate hierarchical
34	  ordered_model_fidelities = 'SIM1'
35	  correction additive zeroth_order
36
37model,
38	id_model = 'SIM1'
39	simulation
40	  solution_level_control = 'mesh_size'
41# dense matrices --> direct solvers --> cubic growth in expense
42	  solution_level_cost = 1. 8. 64. 512. 4096.
43
44variables,
45	uniform_uncertain = 9	 # number of cosine terms in soln expansion
46	  lower_bounds      =  9*-1.
47	  upper_bounds      =  9* 1.
48	discrete_state_set
49	  integer = 1 # even poly degree for odd mesh pts including boundary
50      	  set_values = 4 8 16 32 64
51      	  descriptors = 'mesh_size'
52    	real = 4
53      	  elements_per_variable = 2 2 1 1
54      	  set_values = 0.1 1 # field mean
55                       0.5 4 # field std dev
56                       1     # kernel order
57                       0.2   # kernel length
58      	  descriptors = 'field_mean' 'field_std_dev' 'kernel_order' 'kernel_length'
59      	  initial_state = 1 4 1 0.2
60	string = 2  # form of kernel (refer to Yang,Lei,Baker,Lin)
61      	  elements_per_variable = 2 2
62	  initial_state = 'cosine' 'off'
63      	  set_values = 'cosine' 'exponential'
64                       'off' 'on'
65      	  descriptors = 'kernel_type' 'positivity'
66
67interface,
68	direct
69	  analysis_driver = 'steady_state_diffusion_1d'
70
71responses,
72	calibration_terms = 3	 # number of soln pts excl. +/-5% (symmetric)
73        calibration_data_file = 'dakota_diffusion_queso_mlpce.withsigma.dat'
74          freeform num_experiments = 1 variance_type = 'scalar'
75 	no_gradients
76	no_hessians
77