def test_material_mu_components2():
mu = material_mu('TiO2', 10000, density=4.23)
assert_allclose(mu, 290.7, rtol=0.001)
mu = material_mu('TiO2', 10000)
assert_allclose(mu, 290.7, rtol=0.001)
mu = material_mu('TiO2', 10000, density=4.5)
assert_allclose(mu, 309.26, rtol=0.001)
comps = material_mu_components('TiO2', 10000, density=4.23)
known_comps = {
'mass': 79.88,
'density': 4.23,
'elements': ['Ti', 'O'],
'Ti': (1, 47.88, 110.676),
'O': (2.0, 15.9994, 5.953)
}
assert 'Ti' in comps['elements']
assert 'O' in comps['elements']
for attr in ('mass', 'density'):
assert_allclose(comps[attr], known_comps[attr], rtol=0.01)
for attr in ('Ti', 'O'):
assert_allclose(comps[attr][0], known_comps[attr][0], rtol=0.01)
assert_allclose(comps[attr][1], known_comps[attr][1], rtol=0.01)
assert_allclose(comps[attr][2], known_comps[attr][2], rtol=0.01)
with pytest.raises(Warning):
c = material_mu_components('TiO2', 10000)
def test_material_mu_components1():
mu = material_mu('quartz', 10000)
assert_allclose(mu, 50.368, rtol=0.001)
comps = material_mu_components('quartz', 10000)
known_comps = {'mass': 60.08, 'density': 2.65, 'elements': ['Si', 'O'],
'Si': (1, 28.1, 33.879), 'O': (2.0, 16.0, 5.953)}
assert 'Si'in comps['elements']
assert 'O'in comps['elements']
for attr in ('mass', 'density'):
assert_allclose(comps[attr], known_comps[attr], rtol=0.01)
for attr in ('Si', 'O'):
assert_allclose(comps[attr][0], known_comps[attr][0], rtol=0.01)
assert_allclose(comps[attr][1], known_comps[attr][1], rtol=0.01)
assert_allclose(comps[attr][2], known_comps[attr][2], rtol=0.01)