# Tests for D-optimal sampling
# Currently only support aleatory variables
# With Leja option, only support continuous aleatory variables
# test 0: continuous aleatory
# test 1: discrete aleatory
# test 2: continuous aleatory with Leja option
environment,				
	  tabular_data		

method, 				 
  sampling
    sample_type random
    samples = 100                     #s0,#s1
#   samples = 20                      #s2
      d_optimal
# Default candidate designs = 100
#       candidate_designs 50          #s1
#       leja_oversample_ratio = 10.0  #s2
    seed 8272

## ERROR: Want samples greater, but get this error on s1:
#  terminate called after throwing an instance of 'std::runtime_error'
#    what():  enrich_samples: The basis matrix has rank 21 less than num_total_samples 30. Try increasing the degree of the basis.


variables,								 
# (default is active uncertain)
# active design                                   
# active state                                    

## DESIGN
#	continuous_design = 1                           
#          upper_bounds = 20.			                
#          lower_bounds = 10.			                
## ERROR:
#   Distribution required, not found for DiscDesRange2   
#   Error: code 1 returned from lhs_prep in LHSDriver.
#	discrete_design_range = 1         
#          upper_bounds =  7			 
#          lower_bounds = -3			 
## ERROR (no discrete histogram):
##   Error: RandomVariable type 0 not available.
#	discrete_design_set     		 
#	  integer = 1				
#            elements = 3 4 5	
#	  string = 1				
#            elements = 'dog' 'fleas' 'has' 'my'
#	  real = 1				
#            elements = 30. 40. 50.       

## ALEATORY UNCERTAIN
  normal_uncertain = 1                            #s0,#s2
    means = 69.0       	                          #s0,#s2
    std_deviations = 1.0  	                      #s0,#s2
	lognormal_uncertain = 1		                      #s0,#s2
    means             =  1.	                      #s0,#s2
    std_deviations    =  0.5                      #s0,#s2
	uniform_uncertain = 1                           #s0,#s2     
	  lower_bounds      = 200.	                    #s0,#s2
    upper_bounds = 300.		                        #s0,#s2
 	loguniform_uncertain = 1	                      #s0,#s2
    lower_bounds      =  0.1	                    #s0,#s2
    upper_bounds      =  10.	                    #s0,#s2
 	triangular_uncertain = 1	                      #s0,#s2
    modes             =  0.  	                    #s0,#s2
    lower_bounds      = -2.	                      #s0,#s2
	  upper_bounds      =  2.  	                    #s0,#s2
	exponential_uncertain = 1	                      #s0,#s2
	  betas             =  2.	                      #s0,#s2
	beta_uncertain = 1		                          #s0,#s2
	  alphas            =  1.5  	                  #s0,#s2
	  betas             =  2. 	                    #s0,#s2
	  lower_bounds      = -2.	                      #s0,#s2
	  upper_bounds      =  2.	                      #s0,#s2
	gamma_uncertain = 1                             #s0,#s2
	  alphas            =  2.5  	                  #s0,#s2
	  betas             =  2.  	                    #s0,#s2
	gumbel_uncertain = 1		                        #s0,#s2
	  alphas            = 1.5 	                    #s0,#s2
	  betas             = 2.5	                      #s0,#s2
	frechet_uncertain = 1		                        #s0,#s2
	  alphas            = 2.5 	                    #s0,#s2
	  betas             = 3.5 	                    #s0,#s2
	weibull_uncertain = 1		                        #s0,#s2
	  alphas            =   12.                     #s0,#s2
	  betas             =  250.	                    #s0,#s2
	histogram_bin_uncertain = 1	                    #s0,#s2
	  abscissas   =  5  8 10 	                      #s0,#s2
	  counts      = 17 21  0	                      #s0,#s2
# poisson_uncertain = 1		                        #s1
#   lambdas           =  0.2   	                  #s1
# binomial_uncertain = 1		                      #s1
#   prob_per_trial    =  0.3  	                  #s1
#   num_trials        =  10  	                    #s1
# negative_binomial_uncertain = 1	                #s1
#   prob_per_trial    =  0.3  	                  #s1
#   num_trials        =  10    	                  #s1
# geometric_uncertain = 1		                      #s1
#   prob_per_trial    = 0.1                       #s1
#	hypergeometric_uncertain = 1	                  #s1
#         total_population    = 50 	              #s1
#         selected_population = 20	              #s1
#         num_drawn           =  5 	              #s1
#	histogram_point_uncertain                       #s1
#   integer = 1                                   #s1
#	    abscissas   = 3 4 5                         #s1
#	    counts      = 1 2 1                         #s1
#   string = 1                                    #s1
#	    abscissas   = 's3' 's4' 's5'                #s1
#	    counts      = 1 2 1                         #s1
#   real = 1                                      #s1
#	    abscissas   = 3.2 4.3 5.4                   #s1
#	    counts      = 1 2 1                         #s1

## EPISTEMIC UNCERTAIN
#	continuous_interval_uncertain = 1               
#	  interval_probabilities = .4 .6	              
#	  lower_bounds = 2. 4.	      		              
#	  upper_bounds = 5. 8.			                    
## ERROR (no discrete histogram)
#	discrete_interval_uncertain = 1		
#	  interval_probabilities = .4 .6	
#	  lower_bounds = 2 4	      		
#	  upper_bounds = 5 8			
#	discrete_uncertain_set		
#	  integer = 1			
#	    set_probabilities = .4 .3 .2 .1	
#	    set_values        =  2  4  5  8	
## TODO: string
#	  real = 1		       	  
#	    set_probabilities =  .1  .2  .3  .4	
#	    set_values        = 3.1 5.2 6.3 7.4	

## STATE
#	continuous_state = 1                            
#          upper_bounds = 20.			                
#          lower_bounds = 10.			                
## ERROR:
#   Distribution required, not found for DiscDesRange2   
#   Error: code 1 returned from lhs_prep in LHSDriver.
#	discrete_state_range = 1		
#          upper_bounds =  7	
#          lower_bounds = -3	
## ERROR (no discrete histogram):
##   Error: RandomVariable type 0 not available.
#	discrete_state_set		
#	  integer = 1			
#            set_values = 3 4 5		
## TODO: string
#	  real = 1		
#            set_values = 30. 40. 50.		

interface,								 
	direct					 
	  analysis_driver = 'text_book'
	  deactivate evaluation_cache restart_file

responses,                                      			 
	response_functions = 1					 
	no_gradients
	no_hessians