def test_galaxy0030(): ds = data_dir_load(g30) yield check_color_conservation(ds) yield assert_equal, str(ds), "galaxy0030" for test in big_patch_amr(g30, _fields): test_galaxy0030.__name__ = test.description yield test
def test_galaxy0030(): ds = data_dir_load(g30) yield check_color_conservation(ds) assert_equal(str(ds), "galaxy0030") for test in big_patch_amr(ds, _fields): test_galaxy0030.__name__ = test.description yield test assert_equal(ds.particle_type_counts, {'io': 1124453})
def test_d9p(): ds = data_dir_load(d9p) yield assert_equal, str(ds), "10MpcBox_HartGal_csf_a0.500.d" for test in big_patch_amr(d9p, _fields): test_d9p.__name__ = test.description yield test dso = [None, ("sphere", ("max", (0.1, 'unitary')))] for field in _fields: for axis in [0, 1, 2]: for dobj_name in dso: for weight_field in [None, "density"]: yield PixelizedProjectionValuesTest( d9p, axis, field, weight_field, dobj_name) ad = ds.all_data() # 'Ana' variable values output from the ART Fortran 'ANA' analysis code AnaNStars = 6255 yield assert_equal, ad[('stars','particle_type')].size, AnaNStars yield assert_equal, ad[('specie4', 'particle_type')].size, AnaNStars AnaNDM = 2833405 yield assert_equal, ad[('darkmatter','particle_type')].size, AnaNDM yield assert_equal, ad[('specie0', 'particle_type')].size + \ ad[('specie1', 'particle_type')].size + \ ad[('specie2', 'particle_type')].size + \ ad[('specie3', 'particle_type')].size, AnaNDM AnaBoxSize = yt.units.yt_array.YTQuantity(7.1442196564,'Mpc') AnaVolume = yt.units.yt_array.YTQuantity(364.640074656,'Mpc**3') Volume = 1 for i in ds.domain_width.in_units('Mpc'): yield assert_almost_equal, i, AnaBoxSize Volume *= i yield assert_almost_equal, Volume, AnaVolume AnaNCells = 4087490 yield assert_equal, len(ad[('index','cell_volume')]), AnaNCells AnaTotDMMass = yt.units.yt_array.YTQuantity(1.01191786811e+14,'Msun') yield assert_almost_equal, ad[('darkmatter','particle_mass')].sum()\ .in_units('Msun'), AnaTotDMMass AnaTotStarMass = yt.units.yt_array.YTQuantity(1776251.,'Msun') yield assert_almost_equal, ad[('stars','particle_mass')].sum()\ .in_units('Msun'), AnaTotStarMass AnaTotStarMassInitial = yt.units.yt_array.YTQuantity(2422854.,'Msun') yield assert_almost_equal, ad[('stars','particle_mass_initial')].sum()\ .in_units('Msun'), AnaTotStarMass AnaTotGasMass = yt.units.yt_array.YTQuantity(1.781994e+13,'Msun') yield assert_almost_equal, ad[('gas','cell_mass')].sum()\ .in_units('Msun'), AnaTotGasMass AnaTotTemp = yt.units.yt_array.YTQuantity(1.5019e11, 'K') #just leaves yield assert_equal, ad[('gas','temperature')].sum(), AnaTotTemp
def test_d9p(): ds = data_dir_load(d9p) yield assert_equal, str(ds), "10MpcBox_HartGal_csf_a0.500.d" for test in big_patch_amr(d9p, _fields): test_d9p.__name__ = test.description yield test dso = [None, ("sphere", ("max", (0.1, 'unitary')))] for field in _fields: for axis in [0, 1, 2]: for dobj_name in dso: for weight_field in [None, "density"]: yield PixelizedProjectionValuesTest( d9p, axis, field, weight_field, dobj_name) ad = ds.all_data() # 'Ana' variable values output from the ART Fortran 'ANA' analysis code AnaNStars = 6255 yield assert_equal, ad[('stars', 'particle_type')].size, AnaNStars yield assert_equal, ad[('specie4', 'particle_type')].size, AnaNStars AnaNDM = 2833405 yield assert_equal, ad[('darkmatter', 'particle_type')].size, AnaNDM yield assert_equal, ad[('specie0', 'particle_type')].size + \ ad[('specie1', 'particle_type')].size + \ ad[('specie2', 'particle_type')].size + \ ad[('specie3', 'particle_type')].size, AnaNDM AnaBoxSize = yt.units.yt_array.YTQuantity(7.1442196564, 'Mpc') AnaVolume = yt.units.yt_array.YTQuantity(364.640074656, 'Mpc**3') Volume = 1 for i in ds.domain_width.in_units('Mpc'): yield assert_almost_equal, i, AnaBoxSize Volume *= i yield assert_almost_equal, Volume, AnaVolume AnaNCells = 4087490 yield assert_equal, len(ad[('index', 'cell_volume')]), AnaNCells AnaTotDMMass = yt.units.yt_array.YTQuantity(1.01191786811e+14, 'Msun') yield assert_almost_equal, ad[('darkmatter','particle_mass')].sum()\ .in_units('Msun'), AnaTotDMMass AnaTotStarMass = yt.units.yt_array.YTQuantity(1776251., 'Msun') yield assert_almost_equal, ad[('stars','particle_mass')].sum()\ .in_units('Msun'), AnaTotStarMass AnaTotStarMassInitial = yt.units.yt_array.YTQuantity(2422854., 'Msun') yield assert_almost_equal, ad[('stars','particle_mass_initial')].sum()\ .in_units('Msun'), AnaTotStarMass AnaTotGasMass = yt.units.yt_array.YTQuantity(1.781994e+13, 'Msun') yield assert_almost_equal, ad[('gas','cell_mass')].sum()\ .in_units('Msun'), AnaTotGasMass AnaTotTemp = yt.units.yt_array.YTQuantity(1.5019e11, 'K') #just leaves yield assert_equal, ad[('gas', 'temperature')].sum(), AnaTotTemp
def test_d9p(): ds = data_dir_load(d9p) yield assert_equal, str(ds), "10MpcBox_HartGal_csf_a0.500.d" for test in big_patch_amr(d9p, _fields): test_d9p.__name__ = test.description yield test dso = [None, ("sphere", ("max", (0.1, 'unitary')))] for field in _fields: for axis in [0, 1, 2]: for dobj_name in dso: for weight_field in [None, "density"]: yield PixelizedProjectionValuesTest( d9p, axis, field, weight_field, dobj_name)