#@ s*: Label=FastTest
#@ *: DakotaConfig=HAVE_NPSOL
# DAKOTA INPUT FILE: dakota_sbouu4_tbch.in
#
# OUU with surrogates: surrogate containing nested containing surrogate
#

environment,
	  method_pointer = 'OPTIM'

###########################
# begin opt specification #
###########################
method,
	id_method = 'OPTIM'
	model_pointer = 'OPTIM_M'
        npsol_sqp
	  convergence_tolerance = 1.e-8

model,
	id_model = 'OPTIM_M'
	surrogate global
	  variables_pointer = 'OPTIM_V'
	  responses_pointer = 'OPTIM_R'
	  polynomial quadratic
	  dace_method_pointer = 'OPTIM_DACE'

variables,
	id_variables = 'OPTIM_V'
	continuous_design = 2
	  initial_point    1.8    1.0
	  upper_bounds     2.164  4.0
	  lower_bounds     1.5    0.0
	  descriptors      'd1'   'd2'

responses,
# minimize f + p_fail_r1 + p_fail_r3
# s.t.     g_1/2/3 <= 0
#          mean_r2 + 3*sigma_r2 <= 1.6e+05
# 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 = 4
	nonlinear_inequality_upper_bounds = 0.0 0.0 0.0 1.6e+5
	numerical_gradients
	  method_source dakota
	  interval_type central
	  fd_gradient_step_size = 1.e-5
	no_hessians

############################
# begin DACE specification #
############################
method,
	id_method = 'OPTIM_DACE'
	model_pointer = 'OPTIM_DACE_M'
	sampling
	  sample_type lhs seed = 12347 
#	dace lhs seed = 5

model,
	id_model = 'OPTIM_DACE_M'
	nested
	  variables_pointer  = 'OPTIM_V'
	  sub_method_pointer = 'UQ'
	  optional_interface_pointer  = 'OPTIONAL_I'
	  optional_interface_responses_pointer = 'OPTIONAL_I_R'
	  responses_pointer  = 'OPTIM_R_NO_GRAD'
	  primary_response_mapping   = 0. 0. 1. 0. 0. 0. 0. 0. 1.
	  secondary_response_mapping = 0. 0. 0. 1. 3. 0. 0. 0. 0.

responses,
	id_responses = 'OPTIM_R_NO_GRAD'
	objective_functions = 1
	nonlinear_inequality_constraints = 4
	no_gradients
	no_hessians

##########################################
# begin optional interface specification #
##########################################
interface,
	id_interface = 'OPTIONAL_I'
	system
	  analysis_driver = 'cyl_head'

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

##########################
# begin UQ specification #
##########################
method,
	id_method = 'UQ'
	model_pointer = 'UQ_M'
        sampling
	  samples = 5000 sample_type lhs
	  seed = 12347 
#	  fixed_seed					#s1
	  response_levels = 4.2e+11 1.6e+05 7.5e+05
	  complementary distribution
	  output quiet

model,
	id_model = 'UQ_M'
	surrogate global
	  variables_pointer = 'UQ_V'
	  responses_pointer = 'UQ_R'
	  dace_method_pointer = 'UQ_DACE'
	  polynomial quadratic

variables,
	id_variables = 'UQ_V'
	continuous_design = 2
	normal_uncertain = 2
	  means             =  248.89, 593.33
	  std_deviations    =   12.4,   29.7
# 3 sigma bounds for DDACE:
#	  lower_bounds =  211.7,  504.2
#	  upper_bounds =  286.1,  682.4
	  descriptors  =  'nuv1'  'nuv2'
	uniform_uncertain = 2
	  lower_bounds =  199.3,  474.63
	  upper_bounds =  298.5,  712.
	  descriptors  =  'uuv1'  'uuv2'
	weibull_uncertain = 2
	  alphas       =   12.,    30.
	  betas        =  250.,   590.
#	  lower_bounds =  225.,   540.
#	  upper_bounds =  275.,   640.
	  descriptors  =  'wuv1'  'wuv2'

responses,
	id_responses = 'UQ_R'
	response_functions = 3
	no_gradients
	no_hessians

############################
# begin DACE specification #
############################
method,
	id_method = 'UQ_DACE'
	model_pointer = 'UQ_DACE_M'
        sampling
	  seed = 12347 
	  sample_type lhs
#	dace lhs seed = 5
#	  samples = 50

model,
	id_model = 'UQ_DACE_M'
	single
	  variables_pointer = 'UQ_DACE_V'
	  interface_pointer = 'UQ_ACTUAL'
	  responses_pointer = 'UQ_R'

variables,
	id_variables = 'UQ_DACE_V'
	active all
	continuous_design = 2
	normal_uncertain = 2
	  means             =  248.89, 593.33
	  std_deviations    =   12.4,   29.7
# 3 sigma bounds for DDACE:
#	  lower_bounds =  211.7,  504.2
#	  upper_bounds =  286.1,  682.4
	  descriptors  =  'nuv1'  'nuv2'
	uniform_uncertain = 2
	  lower_bounds =  199.3,  474.63
	  upper_bounds =  298.5,  712.
	  descriptors  =  'uuv1'  'uuv2'
	weibull_uncertain = 2
	  alphas       =   12.,    30.
	  betas        =  250.,   590.
#	  lower_bounds =  225.,   540.
#	  upper_bounds =  275.,   640.
	  descriptors  =  'wuv1'  'wuv2'

interface,
	id_interface = 'UQ_ACTUAL'
	system asynch evaluation_concurrency = 5
	  analysis_driver = 'text_book_ouu'