def test_array_inputs(self):
mm = MetaModel()
mm.add_input('x', np.zeros((2,2)))
mm.add_output('y1', 0.)
mm.add_output('y2', 0.)
mm.default_surrogate = FloatKrigingSurrogate()
prob = Problem()
prob.model.add_subsystem('mm', mm)
prob.setup(check=False)
mm.metadata['train:x'] = [
[[1.0, 1.0], [1.0, 1.0]],
[[2.0, 1.0], [1.0, 1.0]],
[[1.0, 2.0], [1.0, 1.0]],
[[1.0, 1.0], [2.0, 1.0]],
[[1.0, 1.0], [1.0, 2.0]]
]
mm.metadata['train:y1'] = [3.0, 2.0, 1.0, 6.0, -2.0]
mm.metadata['train:y2'] = [1.0, 4.0, 7.0, -3.0, 3.0]
prob['mm.x'] = [[1.0, 2.0], [1.0, 1.0]]
prob.run_model()
assert_rel_error(self, prob['mm.y1'], 1.0, .00001)
assert_rel_error(self, prob['mm.y2'], 7.0, .00001)
def test_unequal_training_outputs(self):
mm = MetaModel()
mm.add_input('x', 0.)
mm.add_input('y', 0.)
mm.add_output('f', 0.)
mm.default_surrogate = FloatKrigingSurrogate()
prob = Problem()
prob.model.add_subsystem('mm', mm)
prob.setup(check=False)
mm.metadata['train:x'] = [1.0, 1.0, 1.0, 1.0]
mm.metadata['train:y'] = [1.0, 2.0, 3.0, 4.0]
mm.metadata['train:f'] = [1.0, 1.0]
prob['mm.x'] = 1.0
prob['mm.y'] = 1.0
with self.assertRaises(RuntimeError) as cm:
prob.run_model()
expected = ("MetaModel: Each variable must have the same number"
" of training points. Expected 4 but found"
" 2 points for 'f'.")
self.assertEqual(str(cm.exception), expected)
def test_derivatives(self):
mm = MetaModel()
mm.add_input('x', 0.)
mm.add_output('f', 0.)
mm.default_surrogate = FloatKrigingSurrogate()
prob = Problem()
prob.model.add_subsystem('p',
IndepVarComp('x', 0.),
promotes_outputs=['x'])
prob.model.add_subsystem('mm', mm, promotes_inputs=['x'])
prob.setup()
mm.metadata['train:x'] = [0., .25, .5, .75, 1.]
mm.metadata['train:f'] = [1., .75, .5, .25, 0.]
prob['x'] = 0.125
prob.run_model()
data = prob.check_partials()
Jf = data['mm'][('f', 'x')]['J_fwd']
Jr = data['mm'][('f', 'x')]['J_rev']
assert_rel_error(self, Jf[0][0], -1., 1.e-3)
assert_rel_error(self, Jr[0][0], -1., 1.e-3)
# TODO: complex step not currently supported in check_partial_derivs
# data = prob.check_partials(global_options={'method': 'cs'})
abs_errors = data['mm'][('f', 'x')]['abs error']
self.assertTrue(len(abs_errors) > 0)
for match in abs_errors:
abs_error = float(match)
self.assertTrue(abs_error < 1.e-6)
def test_metamodel_feature_vector(self):
# Like simple sine example, but with input of length n instead of scalar
# The expected behavior is that the output is also of length n, with
# each one being an independent prediction.
# Its as if you stamped out n copies of metamodel, ran n scalars
# through its input, then muxed all those outputs into one contiguous
# array but you skip all the n-copies thing and do it all as an array
size = 3
# create a vectorized MetaModel for sine
trig = MetaModel(vectorize=size,
default_surrogate=FloatKrigingSurrogate())
trig.add_input('x', np.zeros(size))
trig.add_output('y', np.zeros(size))
# add it to a Problem
prob = Problem()
prob.model.add_subsystem('trig', trig)
prob.setup(check=False)
# provide training data
trig.metadata['train:x'] = np.linspace(0, 10, 20)
trig.metadata['train:y'] = .5 * np.sin(trig.metadata['train:x'])
# train the surrogate and check predicted value
prob['trig.x'] = np.array([2.1, 3.2, 4.3])
prob.run_model()
assert_rel_error(self, prob['trig.y'],
np.array(.5 * np.sin(prob['trig.x'])), 1e-4)
def test_metamodel_feature_vector2d(self):
# similar to previous example, but processes 3 inputs/outputs at a time
size = 3
# create a vectorized MetaModel for sine and cosine
trig = MetaModel(vectorize=size,
default_surrogate=FloatKrigingSurrogate())
trig.add_input('x', np.zeros(size))
trig.add_output('y', np.zeros((size, 2)))
# add it to a Problem
prob = Problem()
prob.model.add_subsystem('trig', trig)
prob.setup(check=False)
# provide training data
trig.metadata['train:x'] = np.linspace(0, 10, 20)
trig.metadata['train:y'] = np.column_stack(
(.5 * np.sin(trig.metadata['train:x']),
.5 * np.cos(trig.metadata['train:x'])))
# train the surrogate and check predicted value
prob['trig.x'] = np.array([2.1, 3.2, 4.3])
prob.run_model()
assert_rel_error(
self, prob['trig.y'],
np.column_stack(
(.5 * np.sin(prob['trig.x']), .5 * np.cos(prob['trig.x']))),
1e-4)
def test_metamodel_feature2d(self):
# similar to previous example, but output is 2d
# create a MetaModel that predicts sine and cosine as an array
trig = MetaModel(default_surrogate=FloatKrigingSurrogate())
trig.add_input('x', 0)
trig.add_output('y', np.zeros(2))
# add it to a Problem
prob = Problem()
prob.model.add_subsystem('trig', trig)
prob.setup(check=False)
# provide training data
trig.metadata['train:x'] = np.linspace(0, 10, 20)
trig.metadata['train:y'] = np.column_stack(
(.5 * np.sin(trig.metadata['train:x']),
.5 * np.cos(trig.metadata['train:x'])))
# train the surrogate and check predicted value
prob['trig.x'] = 2.1
prob.run_model()
assert_rel_error(
self, prob['trig.y'],
np.append(.5 * np.sin(prob['trig.x']),
.5 * np.cos(prob['trig.x'])), 1e-4)
def test_sin_metamodel(self):
# create a MetaModel for sine and add it to a Problem
sin_mm = MetaModel()
sin_mm.add_input('x', 0.)
sin_mm.add_output('f_x', 0.)
prob = Problem()
prob.model.add_subsystem('sin_mm', sin_mm)
# check that missing surrogate is detected in check_config
testlogger = TestLogger()
prob.setup(logger=testlogger)
# Conclude setup but don't run model.
prob.final_setup()
msg = ("No default surrogate model is defined and the "
"following outputs do not have a surrogate model:\n"
"['f_x']\n"
"Either specify a default_surrogate, or specify a "
"surrogate model for all outputs.")
self.assertEqual(len(testlogger.get('error')), 1)
self.assertTrue(msg in testlogger.get('error')[0])
# check that output with no specified surrogate gets the default
sin_mm.default_surrogate = FloatKrigingSurrogate()
prob.setup(check=False)
surrogate = sin_mm._metadata('f_x').get('surrogate')
self.assertTrue(isinstance(surrogate, FloatKrigingSurrogate),
'sin_mm.f_x should get the default surrogate')
# check error message when no training data is provided
with self.assertRaises(RuntimeError) as cm:
prob.run_model()
msg = ("MetaModel: The following training data sets must be "
"provided as metadata for sin_mm: ['train:x', 'train:f_x']")
self.assertEqual(str(cm.exception), msg)
# train the surrogate and check predicted value
sin_mm.metadata['train:x'] = np.linspace(0, 10, 20)
sin_mm.metadata['train:f_x'] = .5 * np.sin(sin_mm.metadata['train:x'])
prob['sin_mm.x'] = 2.1
prob.run_model()
assert_rel_error(self, prob['sin_mm.f_x'],
.5 * np.sin(prob['sin_mm.x']), 1e-4)
def test_array_outputs(self):
meta = MetaModel()
meta.add_param('x', np.zeros((2, 2)))
meta.add_output('y', np.zeros(2,))
meta.default_surrogate = FloatKrigingSurrogate()
prob = Problem(Group())
prob.root.add('meta', meta)
prob.setup(check=False)
prob['meta.train:x'] = [
[[1.0, 1.0], [1.0, 1.0]],
[[2.0, 1.0], [1.0, 1.0]],
[[1.0, 2.0], [1.0, 1.0]],
[[1.0, 1.0], [2.0, 1.0]],
[[1.0, 1.0], [1.0, 2.0]]
]
prob['meta.train:y'] = [[3.0, 1.0],
[2.0, 4.0],
[1.0, 7.0],
[6.0, -3.0],
[-2.0, 3.0]]
prob['meta.x'] = [[1.0, 2.0], [1.0, 1.0]]
prob.run()
assert_rel_error(self, prob['meta.y'], np.array([1.0, 7.0]), .00001)
def test_derivatives(self):
meta = MetaModel()
meta.add_param('x', 0.)
meta.add_output('f', 0.)
meta.default_surrogate = FloatKrigingSurrogate()
prob = Problem(Group())
prob.root.add('meta', meta, promotes=['x'])
prob.root.add('p', IndepVarComp('x', 0.), promotes=['x'])
prob.setup(check=False)
prob['meta.train:x'] = [0., .25, .5, .75, 1.]
prob['meta.train:f'] = [1., .75, .5, .25, 0.]
prob['x'] = 0.125
prob.run()
Jf = prob.calc_gradient(['x'], ['meta.f'], mode='fwd')
Jr = prob.calc_gradient(['x'], ['meta.f'], mode='rev')
assert_rel_error(self, Jf[0][0], -1., 1.e-5)
assert_rel_error(self, Jr[0][0], -1., 1.e-5)
stream = cStringIO()
prob.check_partial_derivatives(out_stream=stream,
global_options={'check_type': 'cs'})
abs_errors = findall('Absolute Error \(.+\) : (.+)', stream.getvalue())
self.assertTrue(len(abs_errors) > 0)
for match in abs_errors:
abs_error = float(match)
self.assertTrue(abs_error < 1.e-6)
def test_unequal_training_outputs(self):
meta = MetaModel()
meta.add_param('x', 0.)
meta.add_param('y', 0.)
meta.add_output('f', 0.)
meta.default_surrogate = FloatKrigingSurrogate()
prob = Problem(Group())
prob.root.add('meta', meta)
prob.setup(check=False)
prob['meta.train:x'] = [1.0, 1.0, 1.0, 1.0]
prob['meta.train:y'] = [1.0, 2.0, 3.0, 4.0]
prob['meta.train:f'] = [1.0, 1.0]
prob['meta.x'] = 1.0
prob['meta.y'] = 1.0
with self.assertRaises(RuntimeError) as cm:
prob.run()
expected = "MetaModel: Each variable must have the same number" \
" of training points. Expected 4 but found" \
" 2 points for 'f'."
self.assertEqual(str(cm.exception), expected)
def test_vector_inputs(self):
meta = MetaModel()
meta.add_param('x', np.zeros(4))
meta.add_output('y1', 0.)
meta.add_output('y2', 0.)
meta.default_surrogate = FloatKrigingSurrogate()
prob = Problem(Group())
prob.root.add('meta', meta)
prob.setup(check=False)
prob['meta.train:x'] = [
[1.0, 1.0, 1.0, 1.0],
[2.0, 1.0, 1.0, 1.0],
[1.0, 2.0, 1.0, 1.0],
[1.0, 1.0, 2.0, 1.0],
[1.0, 1.0, 1.0, 2.0]
]
prob['meta.train:y1'] = [3.0, 2.0, 1.0, 6.0, -2.0]
prob['meta.train:y2'] = [1.0, 4.0, 7.0, -3.0, 3.0]
prob['meta.x'] = [1.0, 2.0, 1.0, 1.0]
prob.run()
assert_rel_error(self, prob['meta.y1'], 1.0, .00001)
assert_rel_error(self, prob['meta.y2'], 7.0, .00001)
def __init__(self):
super(TrigMM, self).__init__()
# Create meta_model for f_x as the response
sin_mm = self.add("sin_mm", MetaModel())
sin_mm.add_param('x', val=0.)
sin_mm.add_output('f_x:float',
val=0.,
surrogate=FloatKrigingSurrogate())
sin_mm.add_output('f_x:norm_dist',
val=(0., 0.),
surrogate=KrigingSurrogate())
def test_derivatives(self):
meta = MetaModel()
meta.add_param('x', 0.)
meta.add_output('f', 0.)
meta.default_surrogate = FloatKrigingSurrogate()
prob = Problem(Group())
prob.root.add('meta', meta, promotes=['x'])
prob.root.add('p', IndepVarComp('x', 0.), promotes=['x'])
prob.setup(check=False)
prob['meta.train:x'] = [0., .25, .5, .75, 1.]
prob['meta.train:f'] = [1., .75, .5, .25, 0.]
prob['x'] = 0.125
prob.run()
Jf = prob.calc_gradient(['x'], ['meta.f'], mode='fwd')
Jr = prob.calc_gradient(['x'], ['meta.f'], mode='rev')
assert_rel_error(self, Jf[0][0], -1., 1.e-5)
assert_rel_error(self, Jr[0][0], -1., 1.e-5)
stream = cStringIO()
prob.check_partial_derivatives(out_stream=stream, global_options={'check_type': 'cs'})
abs_errors = findall('Absolute Error \(.+\) : (.+)', stream.getvalue())
self.assertTrue(len(abs_errors) > 0)
for match in abs_errors:
abs_error = float(match)
self.assertTrue(abs_error < 1.e-6)
def test_unequal_training_outputs(self):
meta = MetaModel()
meta.add_param('x', 0.)
meta.add_param('y', 0.)
meta.add_output('f', 0.)
meta.default_surrogate = FloatKrigingSurrogate()
prob = Problem(Group())
prob.root.add('meta', meta)
prob.setup(check=False)
prob['meta.train:x'] = [1.0, 1.0, 1.0, 1.0]
prob['meta.train:y'] = [1.0, 2.0, 3.0, 4.0]
prob['meta.train:f'] = [1.0, 1.0]
prob['meta.x'] = 1.0
prob['meta.y'] = 1.0
with self.assertRaises(RuntimeError) as cm:
prob.run()
expected = "MetaModel: Each variable must have the same number" \
" of training points. Expected 4 but found" \
" 2 points for 'f'."
self.assertEqual(str(cm.exception), expected)
def test_array_outputs(self):
meta = MetaModel()
meta.add_param('x', np.zeros((2, 2)))
meta.add_output('y', np.zeros(2,))
meta.default_surrogate = FloatKrigingSurrogate()
prob = Problem(Group())
prob.root.add('meta', meta)
prob.setup(check=False)
prob['meta.train:x'] = [
[[1.0, 1.0], [1.0, 1.0]],
[[2.0, 1.0], [1.0, 1.0]],
[[1.0, 2.0], [1.0, 1.0]],
[[1.0, 1.0], [2.0, 1.0]],
[[1.0, 1.0], [1.0, 2.0]]
]
prob['meta.train:y'] = [[3.0, 1.0],
[2.0, 4.0],
[1.0, 7.0],
[6.0, -3.0],
[-2.0, 3.0]]
prob['meta.x'] = [[1.0, 2.0], [1.0, 1.0]]
prob.run()
assert_rel_error(self, prob['meta.y'], np.array([1.0, 7.0]), .00001)
def test_vector_inputs(self):
meta = MetaModel()
meta.add_param('x', np.zeros(4))
meta.add_output('y1', 0.)
meta.add_output('y2', 0.)
meta.default_surrogate = FloatKrigingSurrogate()
prob = Problem(Group())
prob.root.add('meta', meta)
prob.setup(check=False)
prob['meta.train:x'] = [
[1.0, 1.0, 1.0, 1.0],
[2.0, 1.0, 1.0, 1.0],
[1.0, 2.0, 1.0, 1.0],
[1.0, 1.0, 2.0, 1.0],
[1.0, 1.0, 1.0, 2.0]
]
prob['meta.train:y1'] = [3.0, 2.0, 1.0, 6.0, -2.0]
prob['meta.train:y2'] = [1.0, 4.0, 7.0, -3.0, 3.0]
prob['meta.x'] = [1.0, 2.0, 1.0, 1.0]
prob.run()
assert_rel_error(self, prob['meta.y1'], 1.0, .00001)
assert_rel_error(self, prob['meta.y2'], 7.0, .00001)
def test_sin_metamodel_rmse(self):
# create MetaModel with Kriging, using the rmse option
sin_mm = MetaModel()
sin_mm.add_input('x', 0.)
sin_mm.add_output('f_x', 0.)
sin_mm.default_surrogate = KrigingSurrogate(eval_rmse=True)
# add it to a Problem
prob = Problem()
prob.model.add_subsystem('sin_mm', sin_mm)
prob.setup(check=False)
# train the surrogate and check predicted value
sin_mm.metadata['train:x'] = np.linspace(0,10,20)
sin_mm.metadata['train:f_x'] = np.sin(sin_mm.metadata['train:x'])
prob['sin_mm.x'] = 2.1
prob.run_model()
assert_rel_error(self, prob['sin_mm.f_x'], np.sin(2.1), 1e-4) # mean
self.assertTrue(self, sin_mm._metadata('f_x')['rmse'] < 1e-5) # std deviation
def test_warm_start(self):
# create metamodel with warm_restart = True
meta = MetaModel()
meta.add_param('x1', 0.)
meta.add_param('x2', 0.)
meta.add_output('y1', 0.)
meta.add_output('y2', 0.)
meta.default_surrogate = ResponseSurface()
meta.warm_restart = True
# add to problem
prob = Problem(Group())
prob.root.add('meta', meta)
prob.setup(check=False)
# provide initial training data
prob['meta.train:x1'] = [1.0, 3.0]
prob['meta.train:x2'] = [1.0, 4.0]
prob['meta.train:y1'] = [3.0, 1.0]
prob['meta.train:y2'] = [1.0, 7.0]
# run against a data point and check result
prob['meta.x1'] = 2.0
prob['meta.x2'] = 3.0
prob.run()
assert_rel_error(self, prob['meta.y1'], 1.9085, .001)
assert_rel_error(self, prob['meta.y2'], 3.9203, .001)
# Add 3rd training point, moves the estimate for that point
# back to where it should be.
prob['meta.train:x1'] = [2.0]
prob['meta.train:x2'] = [3.0]
prob['meta.train:y1'] = [2.0]
prob['meta.train:y2'] = [4.0]
meta.train = True # currently need to tell meta to re-train
prob.run()
assert_rel_error(self, prob['meta.y1'], 2.0, .00001)
assert_rel_error(self, prob['meta.y2'], 4.0, .00001)
def test_sin_metamodel_preset_data(self):
# preset training data
x = np.linspace(0, 10, 200)
f_x = .5 * np.sin(x)
# create a MetaModel for Sin and add it to a Problem
sin_mm = MetaModel()
sin_mm.add_input('x', 0., training_data=np.linspace(0, 10, 200))
sin_mm.add_output('f_x', 0., training_data=f_x)
prob = Problem()
prob.model.add_subsystem('sin_mm', sin_mm)
# check that missing surrogate is detected in check_setup
testlogger = TestLogger()
prob.setup(logger=testlogger)
# Conclude setup but don't run model.
prob.final_setup()
msg = ("No default surrogate model is defined and the "
"following outputs do not have a surrogate model:\n"
"['f_x']\n"
"Either specify a default_surrogate, or specify a "
"surrogate model for all outputs.")
self.assertEqual(len(testlogger.get('error')), 1)
self.assertTrue(msg in testlogger.get('error')[0])
# check that output with no specified surrogate gets the default
sin_mm.default_surrogate = FloatKrigingSurrogate()
prob.setup(check=False)
surrogate = sin_mm._metadata('f_x').get('surrogate')
self.assertTrue(isinstance(surrogate, FloatKrigingSurrogate),
'sin_mm.f_x should get the default surrogate')
prob['sin_mm.x'] = 2.22
prob.run_model()
assert_rel_error(self, prob['sin_mm.f_x'],
.5 * np.sin(prob['sin_mm.x']), 1e-4)
def test_warm_start(self):
# create metamodel with warm_restart = True
meta = MetaModel()
meta.add_param('x1', 0.)
meta.add_param('x2', 0.)
meta.add_output('y1', 0.)
meta.add_output('y2', 0.)
meta.default_surrogate = ResponseSurface()
meta.warm_restart = True
# add to problem
prob = Problem(Group())
prob.root.add('meta', meta)
prob.setup(check=False)
# provide initial training data
prob['meta.train:x1'] = [1.0, 3.0]
prob['meta.train:x2'] = [1.0, 4.0]
prob['meta.train:y1'] = [3.0, 1.0]
prob['meta.train:y2'] = [1.0, 7.0]
# run against a data point and check result
prob['meta.x1'] = 2.0
prob['meta.x2'] = 3.0
prob.run()
assert_rel_error(self, prob['meta.y1'], 1.9085, .001)
assert_rel_error(self, prob['meta.y2'], 3.9203, .001)
# Add 3rd training point, moves the estimate for that point
# back to where it should be.
prob['meta.train:x1'] = [2.0]
prob['meta.train:x2'] = [3.0]
prob['meta.train:y1'] = [2.0]
prob['meta.train:y2'] = [4.0]
meta.train = True # currently need to tell meta to re-train
prob.run()
assert_rel_error(self, prob['meta.y1'], 2.0, .00001)
assert_rel_error(self, prob['meta.y2'], 4.0, .00001)
def __init__(self):
super(CompressorMap, self).__init__()
compmap = self.add('compmap', MetaModel())
compmap.add_param('Nc', val=1.0)
compmap.add_param('Rline', val=2.0)
compmap.add_param('alpha', val=0.0)
compmap.add_output(
'PR',
val=1.0,
surrogate=NearestNeighbor(interpolant_type='linear'))
compmap.add_output(
'eff',
val=1.0,
surrogate=NearestNeighbor(interpolant_type='linear'))
compmap.add_output(
'Wc',
val=1.0,
surrogate=NearestNeighbor(interpolant_type='linear'))
def test_metamodel_vector_errors(self):
# invalid values for vectorize argument. Bad.
for bad_value in [True, -1, 0, 1, 1.5]:
with self.assertRaises(RuntimeError) as cm:
MetaModel(vectorize=True)
self.assertEqual(str(cm.exception),
"Metamodel: The value of the 'vectorize' "
"argument must be an integer greater than "
"one, found '%s'." % str(bad_value))
# first dimension of all inputs/outputs must be 3
mm = MetaModel(vectorize=3)
with self.assertRaises(RuntimeError) as cm:
mm.add_input('x', np.zeros(2))
self.assertEqual(str(cm.exception),
"Metamodel: First dimension of input 'x' must be 3")
with self.assertRaises(RuntimeError) as cm:
mm.add_output('y', np.zeros(4))
self.assertEqual(str(cm.exception),
"Metamodel: First dimension of output 'y' must be 3")
def test_sin_metamodel(self):
# create a MetaModel for Sin and add it to a Problem
sin_mm = MetaModel()
sin_mm.add_param('x', 0.)
sin_mm.add_output('f_x', 0.)
prob = Problem(Group())
prob.root.add('sin_mm', sin_mm)
# check that missing surrogate is detected in check_setup
stream = cStringIO()
prob.setup(out_stream=stream)
msg = ("No default surrogate model is defined and the "
"following outputs do not have a surrogate model:\n"
"['f_x']\n"
"Either specify a default_surrogate, or specify a "
"surrogate model for all outputs.")
self.assertTrue(msg in stream.getvalue())
# check that output with no specified surrogate gets the default
sin_mm.default_surrogate = FloatKrigingSurrogate()
prob.setup(check=False)
surrogate = prob.root.unknowns.metadata('sin_mm.f_x').get('surrogate')
self.assertTrue(isinstance(surrogate, FloatKrigingSurrogate),
'sin_mm.f_x should get the default surrogate')
# train the surrogate and check predicted value
prob['sin_mm.train:x'] = np.linspace(0, 10, 200)
prob['sin_mm.train:f_x'] = .5 * np.sin(prob['sin_mm.train:x'])
prob['sin_mm.x'] = 2.22
prob.run()
self.assertAlmostEqual(prob['sin_mm.f_x'],
.5 * np.sin(prob['sin_mm.x']),
places=5)
def test_sin_metamodel_obj_return(self):
# create a MetaModel for Sin and add it to a Problem
sin_mm = MetaModel()
sin_mm.add_param('x', 0.)
sin_mm.add_output('f_x', (0., 0.))
prob = Problem(Group())
prob.root.add('sin_mm', sin_mm)
# check that missing surrogate is detected in check_setup
stream = cStringIO()
prob.setup(out_stream=stream)
msg = ("No default surrogate model is defined and the "
"following outputs do not have a surrogate model:\n"
"['f_x']\n"
"Either specify a default_surrogate, or specify a "
"surrogate model for all outputs.")
self.assertTrue(msg in stream.getvalue())
# check that output with no specified surrogate gets the default
sin_mm.default_surrogate = KrigingSurrogate()
prob.setup(check=False)
surrogate = prob.root.unknowns.metadata('sin_mm.f_x').get('surrogate')
self.assertTrue(isinstance(surrogate, KrigingSurrogate),
'sin_mm.f_x should get the default surrogate')
# train the surrogate and check predicted value
prob['sin_mm.train:x'] = np.linspace(0, 10, 20)
prob['sin_mm.train:f_x'] = np.sin(prob['sin_mm.train:x'])
prob['sin_mm.x'] = 2.1
prob.run()
assert_rel_error(self, prob['sin_mm.f_x'][0], np.sin(2.1),
1e-4) # mean
self.assertTrue(self, prob['sin_mm.f_x'][1] < 1e-5) #std deviation
def test_sin_metamodel(self):
# create a MetaModel for Sin and add it to a Problem
sin_mm = MetaModel()
sin_mm.add_param('x', 0.)
sin_mm.add_output('f_x', 0.)
prob = Problem(Group())
prob.root.add('sin_mm', sin_mm)
# check that missing surrogate is detected in check_setup
stream = cStringIO()
prob.setup(out_stream=stream)
msg = ("No default surrogate model is defined and the "
"following outputs do not have a surrogate model:\n"
"['f_x']\n"
"Either specify a default_surrogate, or specify a "
"surrogate model for all outputs.")
self.assertTrue(msg in stream.getvalue())
# check that output with no specified surrogate gets the default
sin_mm.default_surrogate = FloatKrigingSurrogate()
prob.setup(check=False)
surrogate = prob.root.unknowns.metadata('sin_mm.f_x').get('surrogate')
self.assertTrue(isinstance(surrogate, FloatKrigingSurrogate),
'sin_mm.f_x should get the default surrogate')
# train the surrogate and check predicted value
prob['sin_mm.train:x'] = np.linspace(0,10,200)
prob['sin_mm.train:f_x'] = .5*np.sin(prob['sin_mm.train:x'])
prob['sin_mm.x'] = 2.22
prob.run()
self.assertAlmostEqual(prob['sin_mm.f_x'],
.5*np.sin(prob['sin_mm.x']),
places=5)
def test_sin_metamodel_obj_return(self):
# create a MetaModel for Sin and add it to a Problem
sin_mm = MetaModel()
sin_mm.add_param('x', 0.)
sin_mm.add_output('f_x', (0.,0.))
prob = Problem(Group())
prob.root.add('sin_mm', sin_mm)
# check that missing surrogate is detected in check_setup
stream = cStringIO()
prob.setup(out_stream=stream)
msg = ("No default surrogate model is defined and the "
"following outputs do not have a surrogate model:\n"
"['f_x']\n"
"Either specify a default_surrogate, or specify a "
"surrogate model for all outputs.")
self.assertTrue(msg in stream.getvalue())
# check that output with no specified surrogate gets the default
sin_mm.default_surrogate = KrigingSurrogate()
prob.setup(check=False)
surrogate = prob.root.unknowns.metadata('sin_mm.f_x').get('surrogate')
self.assertTrue(isinstance(surrogate, KrigingSurrogate),
'sin_mm.f_x should get the default surrogate')
# train the surrogate and check predicted value
prob['sin_mm.train:x'] = np.linspace(0,10,20)
prob['sin_mm.train:f_x'] = np.sin(prob['sin_mm.train:x'])
prob['sin_mm.x'] = 2.1
prob.run()
assert_rel_error(self, prob['sin_mm.f_x'][0], 0.86323233, 1e-4) # mean
self.assertTrue(self, prob['sin_mm.f_x'][1] < 1e-5) #std deviation
def test_metamodel_feature(self):
# create a MetaModel, specifying surrogates for the outputs
trig = MetaModel()
trig.add_input('x', 0.)
trig.add_output('sin_x', 0., surrogate=FloatKrigingSurrogate())
trig.add_output('cos_x', 0.)
trig.default_surrogate = FloatKrigingSurrogate()
# provide training data
trig.metadata['train:x'] = np.linspace(0, 10, 20)
trig.metadata['train:sin_x'] = .5 * np.sin(trig.metadata['train:x'])
trig.metadata['train:cos_x'] = .5 * np.cos(trig.metadata['train:x'])
# add it to a Problem, run and check the predicted values
prob = Problem()
prob.model.add_subsystem('trig', trig)
prob.setup(check=False)
prob['trig.x'] = 2.1
prob.run_model()
assert_rel_error(self, prob['trig.sin_x'], .5 * np.sin(prob['trig.x']),
1e-4)
assert_rel_error(self, prob['trig.cos_x'], .5 * np.cos(prob['trig.x']),
1e-4)
def test_basics(self):
# create a metamodel component
mm = MetaModel()
mm.add_param('x1', 0.)
mm.add_param('x2', 0.)
mm.add_output('y1', 0.)
mm.add_output('y2', 0., surrogate=FloatKrigingSurrogate())
mm.default_surrogate = ResponseSurface()
# add metamodel to a problem
prob = Problem(root=Group())
prob.root.add('mm', mm)
prob.setup(check=False)
# check that surrogates were properly assigned
surrogate = prob.root.unknowns.metadata('mm.y1').get('surrogate')
self.assertTrue(isinstance(surrogate, ResponseSurface))
surrogate = prob.root.unknowns.metadata('mm.y2').get('surrogate')
self.assertTrue(isinstance(surrogate, FloatKrigingSurrogate))
# populate training data
prob['mm.train:x1'] = [1.0, 2.0, 3.0]
prob['mm.train:x2'] = [1.0, 3.0, 4.0]
prob['mm.train:y1'] = [3.0, 2.0, 1.0]
prob['mm.train:y2'] = [1.0, 4.0, 7.0]
# run problem for provided data point and check prediction
prob['mm.x1'] = 2.0
prob['mm.x2'] = 3.0
self.assertTrue(mm.train) # training will occur before 1st run
prob.run()
assert_rel_error(self, prob['mm.y1'], 2.0, .00001)
assert_rel_error(self, prob['mm.y2'], 4.0, .00001)
# run problem for interpolated data point and check prediction
prob['mm.x1'] = 2.5
prob['mm.x2'] = 3.5
self.assertFalse(mm.train) # training will not occur before 2nd run
prob.run()
assert_rel_error(self, prob['mm.y1'], 1.5934, .001)
# change default surrogate, re-setup and check that metamodel re-trains
mm.default_surrogate = FloatKrigingSurrogate()
prob.setup(check=False)
surrogate = prob.root.unknowns.metadata('mm.y1').get('surrogate')
self.assertTrue(isinstance(surrogate, FloatKrigingSurrogate))
self.assertTrue(mm.train) # training will occur after re-setup
mm.warm_restart = True # use existing training data
prob['mm.x1'] = 2.5
prob['mm.x2'] = 3.5
prob.run()
assert_rel_error(self, prob['mm.y1'], 2., 1e-2)
def test_basics(self):
# create a metamodel component
mm = MetaModel()
mm.add_param('x1', 0.)
mm.add_param('x2', 0.)
mm.add_output('y1', 0.)
mm.add_output('y2', 0., surrogate=FloatKrigingSurrogate())
mm.default_surrogate = ResponseSurface()
# add metamodel to a problem
prob = Problem(root=Group())
prob.root.add('mm', mm)
prob.setup(check=False)
# check that surrogates were properly assigned
surrogate = prob.root.unknowns.metadata('mm.y1').get('surrogate')
self.assertTrue(isinstance(surrogate, ResponseSurface))
surrogate = prob.root.unknowns.metadata('mm.y2').get('surrogate')
self.assertTrue(isinstance(surrogate, FloatKrigingSurrogate))
# populate training data
prob['mm.train:x1'] = [1.0, 2.0, 3.0]
prob['mm.train:x2'] = [1.0, 3.0, 4.0]
prob['mm.train:y1'] = [3.0, 2.0, 1.0]
prob['mm.train:y2'] = [1.0, 4.0, 7.0]
# run problem for provided data point and check prediction
prob['mm.x1'] = 2.0
prob['mm.x2'] = 3.0
self.assertTrue(mm.train) # training will occur before 1st run
prob.run()
assert_rel_error(self, prob['mm.y1'], 2.0, .00001)
assert_rel_error(self, prob['mm.y2'], 4.0, .00001)
# run problem for interpolated data point and check prediction
prob['mm.x1'] = 2.5
prob['mm.x2'] = 3.5
self.assertFalse(mm.train) # training will not occur before 2nd run
prob.run()
assert_rel_error(self, prob['mm.y1'], 1.5934, .001)
# change default surrogate, re-setup and check that metamodel re-trains
mm.default_surrogate = FloatKrigingSurrogate()
prob.setup(check=False)
surrogate = prob.root.unknowns.metadata('mm.y1').get('surrogate')
self.assertTrue(isinstance(surrogate, FloatKrigingSurrogate))
self.assertTrue(mm.train) # training will occur after re-setup
mm.warm_restart = True # use existing training data
prob['mm.x1'] = 2.5
prob['mm.x2'] = 3.5
prob.run()
assert_rel_error(self, prob['mm.y1'], 2., 1e-2)
