1#@ *: Label=Experimental
2#@ *: DakotaConfig=HAVE_QUESO
3#@ s*: TimeoutAbsolute=3600
4#@ s1: ReqFiles=dakota_bayes_expdesign.dat
5#@ s2: ReqFiles=dakota_bayes_expdesign.dat,dakota_bayes_expdesign.cand.dat
6#@ s2: UserMan=bayes_experimental_design
7
8# This example tests the Bayesian experimental_design algorithm. The
9# high-fidelity model is a steady-state heat equation example, which the
10# low-fidelity model attempts to approximate with a quadratic function,
11# Ax^2+Bx+C. The coefficients A, B, and C are the calibration parameters,
12# while x is the configuration variable. See Section 4.1 of the JCP paper
13# "An information theoretic approach to use high-fidelity codes to
14# calibration low-fidelity codes" by Lewis, et al, published in 2016.
15#
16# s3 and s4: Second-order diffs due to changes in Boost versions are
17# expected. These changes occur in boost/random/normal_distribution.hpp.
18# Baseline uses Boost 1.53.0. Platforms using 1.58.0 are known to diff.
19# Moments of posterior parameter and response distributions in s3 and s4
20# should differ from those in s0; baselines for these tests can be
21# updated as the variate generator is udpated.
22
23environment
24  tabular_data
25
26method,
27#          #output debug						#s1
28       bayes_calibration queso
29          dram
30	  seed = 34785
31          chain_samples = 500
32          experimental_design
33          initial_samples = 2
34	  num_candidates = 2				#s0,#s1,#s2,#s3,#s4
35#	  num_candidates = 5				#s5
36# 	  batch_size = 2				#s5
37#	  import_candidate_points_file = 'dakota_bayes_expdesign.cand.dat' #s2
38#		freeform						#s2
39#	  #ksg2						#s2
40	  max_hifi_evaluations = 1			#s0,#s1,#s2,#s3,#s4
41#	  max_hifi_evaluations = 3			#s5
42          model_pointer = 'HIERARCH'
43
44model,
45	id_model = 'HIERARCH'
46	variables_pointer = 'ALL_VARS'
47	surrogate hierarchical
48	  ordered_model_fidelities = 'LF' 'HF'
49
50model,
51        id_model = 'LF'
52        single
53          interface_pointer = 'lofi_IF'
54        variables_pointer = 'ALL_VARS'
55        responses_pointer = 'lofi_resp'
56
57model,
58        id_model = 'HF'
59        single
60          interface_pointer = 'hifi_IF'
61        variables_pointer = 'CONFIG_VARS'
62        responses_pointer = 'hifi_resp'
63
64variables,
65	id_variables = 'ALL_VARS'
66        continuous_state = 1
67	  upper_bounds = 70
68          lower_bounds = 10
69        uniform_uncertain = 3
70	  upper_bounds   0.06  	0   	260
71	  lower_bounds   0	-8	0
72
73variables,
74	id_variables = 'CONFIG_VARS'
75# MSE: for hierarchical model self-consistency, would prefer to leave state
76#      as inactive and have zero active vars
77	active state
78	continuous_state = 1
79          upper_bounds = 70
80          lower_bounds = 10
81
82interface,
83        id_interface = 'hifi_IF'
84        analysis_drivers = 'expdesign_high'
85          fork
86        asynchronous evaluation_concurrency=2              #s0,#s3,#s4,#s5
87
88interface,
89        id_interface = 'lofi_IF'
90        analysis_drivers = 'expdesign_low'
91          fork
92        asynchronous evaluation_concurrency=4              #s0,#s3,#s4,#s5
93
94responses,
95 	id_responses = 'lofi_resp'
96        calibration_terms = 1
97	no_gradients
98	no_hessians
99#	simulation_variance = 130			#s3
100
101responses,
102 	id_responses = 'hifi_resp'
103 	calibration_terms = 1
104#	calibration_data_file = 'dakota_bayes_expdesign.dat'	#s1,#s2
105#	  freeform						#s1,#s2
106#	  num_config_variables = 1				#s1,#s2
107#	  num_experiments = 2				        #s2
108#	  num_experiments = 1					#s1
109#	  experiment_variance_type = 'none'			#s1,#s2
110	no_gradients
111	no_hessians
112#	simulation_variance = 0.5			#s4
113
114
115