#@ s1: TimeoutAbsolute=3600
#@ s1: TimeoutDelay=1200
#@ s2: TimeoutDelay=180

environment,
	output_precision = 16

method,
	model_pointer = 'HIERARCH'
        multifidelity_polynomial_chaos
	  allocation_control greedy
	  p_refinement dimension_adaptive generalized		#s0
	    sparse_grid_level_sequence = 0 unrestricted		#s0
	    convergence_tolerance 1.e-8				#s0
#	  p_refinement uniform	       	 			#s1,#s2
#	    expansion_order_sequence = 0			#s1
#	    collocation_ratio = .9  seed = 160415		#s1
#	    orthogonal_matching_pursuit				#s1
#	    convergence_tolerance 1.e-2				#s1
#	    sparse_grid_level_sequence = 0 unrestricted		#s2
#	    convergence_tolerance 1.e-3				#s2
#	    max_refinement_iterations 100
#	output quiet						#s1,#s2
	output silent						#s0
#	variance_based_decomp

model,
	id_model = 'HIERARCH'
	surrogate hierarchical
	  ordered_model_fidelities = 'SIM1'
	  correction additive zeroth_order

model,
	id_model = 'SIM1'
	simulation
	  solution_level_control = 'mesh_size'
	  solution_level_cost = 1. 8. 64. 512. 4096.

variables,
	uniform_uncertain = 9	 # number of cosine terms in soln expansion
	  lower_bounds      =  9*-1.
	  upper_bounds      =  9* 1.
	discrete_state_set
	  integer = 1 # even poly degree for odd mesh pts including boundary
      	  set_values = 4 8 16 32 64
      	  descriptors = 'mesh_size'
    	real = 4
      	  elements_per_variable = 2 2 1 1
      	  set_values = 0.1 1 # field mean
                       0.5 4 # field std dev
                       1     # kernel order
                       0.2   # kernel length
      	  descriptors = 'field_mean' 'field_std_dev' 'kernel_order' 'kernel_length'
      	  initial_state = 1 4 1 0.2
	string = 2  # form of kernel (refer to Yang,Lei,Baker,Lin)
      	  elements_per_variable = 2 2
	  initial_state = 'cosine' 'off'
      	  set_values = 'cosine' 'exponential'
                       'off' 'on'
      	  descriptors = 'kernel_type' 'positivity'

interface,
	direct
	  analysis_driver = 'steady_state_diffusion_1d'
	  deactivate restart_file

responses,
	response_functions = 3	 # number of soln pts excl. +/-5% (symmetric)
	no_gradients
	no_hessians