Example #1
0
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)
    yield assert_equal, str(ds), "galaxy0030"
    for test in big_patch_amr(g30, _fields):
        test_galaxy0030.__name__ = test.description
        yield test
Example #3
0
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
Example #6
0
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)