#@ s*: Label=FastTest
#@ *: DakotaConfig=HAVE_NPSOL
# DAKOTA Input File: dakota_pcbdo_short_column.in
# Optimization under uncertainty using polynomial chaos methods within a
# nested OUU environment.  The test problem is the "short_column" problem from
# Kuschel and Rackwitz, 1997.

environment,
	  method_pointer = 'OPTIM'

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

model,
	id_model = 'OPTIM_M'
	nested
	  variables_pointer  = 'OPTIM_V'
	  sub_method_pointer = 'UQ'
	  responses_pointer  = 'OPTIM_R'
# use projection of analytic PCE moments: constrain beta
	  primary_response_mapping   = 1. 0. 0. 0. 0.
	  secondary_response_mapping = 0. 0. 0. 0. 1.

variables,
	id_variables = 'OPTIM_V'
	continuous_design = 2
	  initial_point    10.    15.
	  lower_bounds      5.    15.
	  upper_bounds     15.    25.
	  descriptors      'b'    'h'

responses,
# minimize b*h
# s.t.     beta >= 2.5
# 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_lower_bounds = 2.5
	nonlinear_inequality_upper_bounds = 1.e+50
	analytic_gradients			#s0,#s1,#s4,#s5
#	numerical_gradients			#s2,#s3,#s6,#s7
#	  method_source dakota			#s2,#s3,#s6,#s7
#	  interval_type central			#s2,#s3,#s6,#s7
#	  fd_gradient_step_size = 1.e-6		#s2,#s3,#s6,#s7
	no_hessians

##########################
# begin UQ specification #
##########################
method,
	id_method = 'UQ'
	model_pointer = 'UQ_M'
	polynomial_chaos askey
 	  expansion_order = 2			#s0,#s2
# 	  expansion_order = 3			#s1,#s3
	  collocation_ratio = 2.		#s0,#s1,#s2,#s3
	  seed = 1234567 fixed_seed 		#s0,#s1,#s2,#s3
#	  quadrature_order = 3
#	  quadrature_order = 4
#	  sparse_grid_level = 2 non_nested	#s4,#s6
#	  sparse_grid_level = 3 non_nested	#s5,#s7
	  num_response_levels = 0 1
	  response_levels = 0.0
	  compute reliabilities
	  cumulative distribution

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

variables,
	id_variables = 'UQ_V'
#	active all				#s1,#s3,#s5,#s7
	continuous_design = 2
	normal_uncertain = 2
	  means              =  500.0   2000.0
	  std_deviations     =  100.0    400.0
	  descriptors        =   'P'      'M'
	lognormal_uncertain = 1
       	  means             =  5.0
	  std_deviations    =  0.5
	  descriptors       =  'Y'
	uncertain_correlation_matrix =  1   0.5 0
					0.5 1   0
					0   0   1

interface,
	id_interface = 'UQ_I'
	direct
	  analysis_driver = 'short_column'

responses,
	id_responses = 'UQ_R'
	response_functions = 2
	analytic_gradients			#s0,#s4
#	no_gradients				#s1,#s2,#s3,#s5,#s6,#s7
	no_hessians