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

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 = 'DACE'

variables,
	id_variables = 'OPTIM_V'
	continuous_design = 2
	  initial_point    2.5    2.5
	  upper_bounds    10.0   10.0
	  lower_bounds     1.0    1.0
	  descriptors 'w' 't'

responses,
# minimize mean Weight
# s.t.     mean_S + 3 sigma_S <= 0
#          mean_D + 3 sigma_D <= 0
#
# 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 = 2
	analytic_gradients
	no_hessians

############################
# begin DACE specification #
############################
method,
	id_method = 'DACE'
	model_pointer = 'DACE_M'
	sampling
	  sample_type lhs #samples = 9
	  seed = 12347 
#	dace lhs seed = 5
#	  samples = 50 symbols = 50

model,
	id_model = 'DACE_M'
	nested
	  variables_pointer  = 'OPTIM_V'
	  sub_method_pointer = 'UQ'
	  responses_pointer  = 'OPTIM_R_NO_GRAD'
	  primary_response_mapping   = 1. 0. 0. 0. 0. 0.
	  secondary_response_mapping = 0. 0. 1. 3. 0. 0.
                                       0. 0. 0. 0. 1. 3.

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

##########################
# begin UQ specification #
##########################
method,
	id_method = 'UQ'
	model_pointer = 'UQ_M'
        sampling
	  sample_type lhs samples = 50
	  seed = 12347 
#	  fixed_seed					#s1
	  output quiet

model,
	id_model = 'UQ_M'
	single
	  variables_pointer = 'UQ_V'
	  interface_pointer = 'UQ_I'
	  responses_pointer = 'UQ_R'

variables,
	id_variables = 'UQ_V'
	continuous_design = 2
	normal_uncertain = 4
	  means             =  40000. 29.E+6 500. 1000.
	  std_deviations    =  2000. 1.45E+6 100. 100.
	  descriptors       =  'R' 'E' 'X' 'Y'

interface,
	id_interface = 'UQ_I'
	direct
	  analysis_driver = 'cantilever'
	  deactivate evaluation_cache restart_file

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