def test_models_all_xspec():
"""Check all models are returned"""
import sherpa.astro.xspec
s = Session()
s._add_model_types(sherpa.models.basic)
s._add_model_types(sherpa.astro.xspec)
# Actually, it's not guaranteed that these models are disjoint
# but they are at the moment.
#
mall = s.list_models('all')
m1d = s.list_models('1d')
m2d = s.list_models('2d')
mxs = s.list_models('xspec')
nall = len(mall)
assert nall == len(m1d) + len(m2d) + len(mxs)
mall = set(mall)
assert nall == len(mall)
mcomb = set(m1d).union(set(m2d)).union(set(mxs))
assert nall == len(mcomb)
def test_models_all():
"""We can get a list of all models"""
s = Session()
s._add_model_types(sherpa.models.basic)
models = s.list_models()
assert type(models) == list
assert len(models) > 1
assert all([type(m) == str for m in models])
def test_models_filtered(mtype):
"""We can get a list of a subset of models
Check the fix for issue #749
"""
import sherpa.astro.xspec
s = Session()
s._add_model_types(sherpa.models.basic)
s._add_model_types(sherpa.astro.xspec)
models = s.list_models(mtype)
assert type(models) == list
assert len(models) > 1
assert all([type(m) == str for m in models])
def test_models_models():
"""There are no models available by default"""
s = Session()
assert s.list_models() == []
def test_paramprompt():
"""Does paramprompt work?"""
s = Session()
assert s.list_model_ids() == []
assert s.list_model_components() == []
# Add in some models
s._add_model_types(sherpa.models.basic)
models = s.list_models()
assert models != []
assert 'const1d' in models
s.create_model_component('const1d', 'm1')
assert s.list_model_ids() == []
expected = ['m1']
assert s.list_model_components() == expected
# Now there are multiple components do not rely on any ordering
# provided by list_model_components()
#
s.paramprompt(True)
# paramprompt doesn't affect create_model_component (as shown
# by the fact that we don't have to mock stdin)
#
s.create_model_component('const1d', 'm2')
assert s.list_model_ids() == []
expected = set(expected)
expected.add('m2')
assert set(s.list_model_components()) == expected
# it does affect set_model
#
# pytest errors out if you try to read from stdin in a test, so
# try to work around this here using
# https://stackoverflow.com/questions/13238566/python-equivalent-of-input-using-sys-stdin
# An alternative would be just to check that the error is
# raised, and the text, but then we don't get to test the
# behaviour of the paramprompt code.
#
with patch("sys.stdin", StringIO(u"2.1")):
s.set_model('const1d.mx')
assert s.list_model_ids() == [1]
expected.add('mx')
assert set(s.list_model_components()) == expected
mx = s.get_model_component('mx')
assert mx.c0.val == pytest.approx(2.1)
with patch("sys.stdin", StringIO(u"2.1e-3 , 2.0e-3, 1.2e-2")):
s.set_model('x', 'const1d.my')
assert set(s.list_model_ids()) == set([1, 'x'])
expected.add('my')
assert set(s.list_model_components()) == expected
my = s.get_model_component('my')
assert my.c0.val == pytest.approx(2.1e-3)
assert my.c0.min == pytest.approx(2.0e-3)
assert my.c0.max == pytest.approx(1.2e-2)
with patch("sys.stdin", StringIO(u"2.1e-3 ,, 1.2e-2")):
s.set_model(2, 'const1d.mz')
assert set(s.list_model_ids()) == set([1, 2, 'x'])
expected.add('mz')
assert set(s.list_model_components()) == expected
mz = s.get_model_component('mz')
assert mz.c0.val == pytest.approx(2.1e-3)
assert mz.c0.min == pytest.approx(- parameter.hugeval)
assert mz.c0.max == pytest.approx(1.2e-2)
with patch("sys.stdin", StringIO(u"-12.1,-12.1")):
s.set_model('const1d.f1')
assert set(s.list_model_ids()) == set([1, 2, 'x'])
expected.add('f1')
assert set(s.list_model_components()) == expected
f1 = s.get_model_component('f1')
assert f1.c0.val == pytest.approx(-12.1)
assert f1.c0.min == pytest.approx(-12.1)
assert f1.c0.max == pytest.approx(parameter.hugeval)
with patch("sys.stdin", StringIO(u" ,-12.1")):
s.set_model('const1d.f2')
# stop checking list_model_ids
expected.add('f2')
assert set(s.list_model_components()) == expected
f2 = s.get_model_component('f2')
assert f2.c0.val == pytest.approx(1.0)
assert f2.c0.min == pytest.approx(-12.1)
assert f2.c0.max == pytest.approx(parameter.hugeval)
with patch("sys.stdin", StringIO(u" ,")):
s.set_model('const1d.f3')
f3 = s.get_model_component('f3')
assert f3.c0.val == pytest.approx(1.0)
assert f3.c0.min == pytest.approx(- parameter.hugeval)
assert f3.c0.max == pytest.approx(parameter.hugeval)
with patch("sys.stdin", StringIO(u" ,, ")):
s.set_model('const1d.f4')
f4 = s.get_model_component('f4')
assert f4.c0.val == pytest.approx(1.0)
assert f4.c0.min == pytest.approx(- parameter.hugeval)
assert f4.c0.max == pytest.approx(parameter.hugeval)
def test_default_models():
"""There are no models available by default"""
s = Session()
assert s.list_models() == []
def test_paramprompt():
"""Does paramprompt work?"""
s = Session()
assert s.list_model_ids() == []
assert s.list_model_components() == []
# Add in some models
import sherpa.models.basic
s._add_model_types(sherpa.models.basic)
assert s.list_models() != []
s.create_model_component('const1d', 'm1')
assert s.list_model_ids() == []
expected = ['m1']
assert s.list_model_components() == expected
# Now there are multiple components do not rely on any ordering
# provided by list_model_components()
#
s.paramprompt(True)
# paramprompt doesn't affect create_model_component (as shown
# by the fact that we don't have to mock stdin)
#
s.create_model_component('const1d', 'm2')
assert s.list_model_ids() == []
expected = set(expected)
expected.add('m2')
assert set(s.list_model_components()) == expected
# it does affect set_model
#
# pytest errors out if you try to read from stdin in a test, so
# try to work around this here using
# https://stackoverflow.com/questions/13238566/python-equivalent-of-input-using-sys-stdin
# An alternative would be just to check that the error is
# raised, and the text, but then we don't get to test the
# behaviour of the paramprompt code.
#
with patch("sys.stdin", StringIO(u"2.1")):
s.set_model('const1d.mx')
assert s.list_model_ids() == [1]
expected.add('mx')
assert set(s.list_model_components()) == expected
mx = s.get_model_component('mx')
assert mx.c0.val == pytest.approx(2.1)
with patch("sys.stdin", StringIO(u"2.1e-3 , 2.0e-3, 1.2e-2")):
s.set_model('x', 'const1d.my')
assert set(s.list_model_ids()) == set([1, 'x'])
expected.add('my')
assert set(s.list_model_components()) == expected
my = s.get_model_component('my')
assert my.c0.val == pytest.approx(2.1e-3)
assert my.c0.min == pytest.approx(2.0e-3)
assert my.c0.max == pytest.approx(1.2e-2)
with patch("sys.stdin", StringIO(u"2.1e-3 ,, 1.2e-2")):
s.set_model(2, 'const1d.mz')
assert set(s.list_model_ids()) == set([1, 2, 'x'])
expected.add('mz')
assert set(s.list_model_components()) == expected
mz = s.get_model_component('mz')
assert mz.c0.val == pytest.approx(2.1e-3)
assert mz.c0.min == pytest.approx(- parameter.hugeval)
assert mz.c0.max == pytest.approx(1.2e-2)
with patch("sys.stdin", StringIO(u"-12.1,-12.1")):
s.set_model('const1d.f1')
assert set(s.list_model_ids()) == set([1, 2, 'x'])
expected.add('f1')
assert set(s.list_model_components()) == expected
f1 = s.get_model_component('f1')
assert f1.c0.val == pytest.approx(-12.1)
assert f1.c0.min == pytest.approx(-12.1)
assert f1.c0.max == pytest.approx(parameter.hugeval)
with patch("sys.stdin", StringIO(u" ,-12.1")):
s.set_model('const1d.f2')
# stop checking list_model_ids
expected.add('f2')
assert set(s.list_model_components()) == expected
f2 = s.get_model_component('f2')
assert f2.c0.val == pytest.approx(1.0)
assert f2.c0.min == pytest.approx(-12.1)
assert f2.c0.max == pytest.approx(parameter.hugeval)
with patch("sys.stdin", StringIO(u" ,")):
s.set_model('const1d.f3')
f3 = s.get_model_component('f3')
assert f3.c0.val == pytest.approx(1.0)
assert f3.c0.min == pytest.approx(- parameter.hugeval)
assert f3.c0.max == pytest.approx(parameter.hugeval)
with patch("sys.stdin", StringIO(u" ,, ")):
s.set_model('const1d.f4')
f4 = s.get_model_component('f4')
assert f4.c0.val == pytest.approx(1.0)
assert f4.c0.min == pytest.approx(- parameter.hugeval)
assert f4.c0.max == pytest.approx(parameter.hugeval)
