#@ s*: Label=FastTest
#@ *: DakotaConfig=HAVE_NPSOL
environment,
	method_pointer = 'SBLO'

method,
	id_method = 'SBLO'
	surrogate_based_local
	model_pointer = 'OPTIM_M'
	  approx_method_pointer = 'OPTIM'
	  max_iterations = 50
	  soft_convergence_limit = 2
	  trust_region
	    initial_size = 1.0
	    contraction_factor = 0.5
	    expansion_factor   = 1.50
	  output verbose

###########################
# begin opt specification #
###########################
method,
	id_method = 'OPTIM'
	npsol_sqp
#	dot_sqp
	  convergence_tolerance = 1.e-6
	  output verbose

model,
	id_model = 'OPTIM_M'
	variables_pointer  = 'OPTIM_V'
	responses_pointer  = 'OPTIM_R'
	surrogate hierarchical
	  ordered_model_fidelities = 'LO_FI' 'HI_FI'
	  correction additive first_order	#s0
#	  correction additive second_order	#s1

variables,
	id_variables = 'OPTIM_V'
	continuous_design = 3
	  initial_point   300.     20.      300.
	  lower_bounds    200.     10.      100.
	  upper_bounds    400.     30.      500.
	  descriptors     'b'      'd'      'h'

responses,
# NOTE: This specifies the TOTAL RESPONSE for the optimization,
#       which is a combination of nested & interface responses.
	id_responses = 'OPTIM_R'
	objective_functions = 1
	nonlinear_inequality_constraints = 1
	nonlinear_inequality_upper_bounds = 4000.
	analytic_gradients
	no_hessians				#s0
#	quasi_hessians sr1			#s1

responses,
# NOTE: This specifies the TOTAL RESPONSE for the optimization,
#       which is a combination of nested & interface responses.
	id_responses = 'OPTIM_R_MIXED'
	objective_functions = 1
	nonlinear_inequality_constraints = 1
	nonlinear_inequality_upper_bounds = 4000.
	mixed_gradients
	  id_analytic_gradients = 2
	  id_numerical_gradients = 1
	  method_source dakota
	  interval_type central
	  fd_gradient_step_size = 1.e-4
	no_hessians				#s0
#	quasi_hessians sr1			#s1

##########################
# begin LF specification #
##########################
model,
	id_model = 'LO_FI'
	nested
	  variables_pointer  = 'OPTIM_V'
	  sub_method_pointer = 'UQ_MV'
	  optional_interface_pointer  = 'OPTIONAL_I'
	  optional_interface_responses_pointer = 'OPTIONAL_I_R'
	  responses_pointer  = 'OPTIM_R_MIXED'
	  primary_variable_mapping   = 'B'    'D'    'H'
	  secondary_variable_mapping = 'mean' 'mean' 'mean'
# maximize beta
# s.t.     Cost <= 4000.
	  primary_response_mapping = 0. 0. -1.

method,
	id_method = 'UQ_MV'
	model_pointer = 'UQ_M_LO_FI'
	local_reliability
	  response_levels = 0.0
	  compute reliabilities
	  cumulative distribution

##########################
# begin HF specification #
##########################
model,
	id_model = 'HI_FI'
	nested
	  variables_pointer  = 'OPTIM_V'
	  sub_method_pointer = 'UQ_UNCERTAIN_EXP'
	  optional_interface_pointer  = 'OPTIONAL_I'
	  optional_interface_responses_pointer = 'OPTIONAL_I_R'
	  responses_pointer  = 'OPTIM_R'
	  primary_variable_mapping   = 'B'    'D'    'H'
	  secondary_variable_mapping = 'mean' 'mean' 'mean'
# maximize beta
# s.t.     Cost <= 4000.
	  primary_response_mapping = 0. 0. -1.

method,
	id_method = 'UQ_UNCERTAIN_EXP'
	model_pointer = 'UQ_M_HI_FI'
	stoch_collocation askey
	  sparse_grid_level = 2 non_nested
	  response_levels = 0.0
	  compute reliabilities
	  cumulative distribution
	  output silent

##########################################
# begin optional interface specification #
##########################################
interface,
	id_interface = 'OPTIONAL_I'
	direct
	  analysis_driver = 'steel_column_cost'

responses,
# NOTE: This is the response set from the optional interface.
	id_responses = 'OPTIONAL_I_R'
	objective_functions = 0
	nonlinear_inequality_constraints = 1
	analytic_gradients
	no_hessians

#################################
# begin UQ model specifications #
#################################
# NOTE: must have different models to enable different variables views
model,
	id_model = 'UQ_M_LO_FI'
	single
	  variables_pointer = 'UQ_V'
	  interface_pointer = 'UQ_I'
	  responses_pointer = 'UQ_R_LO_FI'

model,
	id_model = 'UQ_M_HI_FI'
	single
	  variables_pointer = 'UQ_V'
	  interface_pointer = 'UQ_I'
	  responses_pointer = 'UQ_R_HI_FI'

variables,
	id_variables = 'UQ_V'
	normal_uncertain = 2
	  means           =  30. 500000.
	  std_deviations  =  10.  50000.
	  descriptors     = 'F0'  'P1'
	lognormal_uncertain = 4
	  means          = 300.  20. 300. 400.
	  std_deviations =   3.   2.   5.  35.
	  descriptors    =  'B'  'D'  'H' 'Fs'
	gumbel_uncertain = 2
	  alphas         = 1.4250554e-5 1.4250554e-5
	  betas          = 559496.31    559496.31
	  descriptors    =  'P2'    'P3'
	weibull_uncertain = 1
	  alphas          =     5.7974
	  betas           =	22679.4777
	  descriptors     =    'E'

interface,
	id_interface = 'UQ_I'
	direct
	  analysis_driver = 'steel_column_perf'
#	  deactivate evaluation_cache restart_file

responses,
	id_responses = 'UQ_R_LO_FI'
	response_functions = 1
	analytic_gradients
	no_hessians

responses,
	id_responses = 'UQ_R_HI_FI'
	response_functions = 1
	analytic_gradients
	no_hessians