@pytest.fixture
def mock_galaxy(disc_particles_all, halo_particles):
    'Mock galaxy.'
    (mass_s, pos_s, vel_s, mass_g, pos_g, vel_g) = disc_particles_all
    (mass_dm, pos_dm, vel_dm) = halo_particles(N_part=100, seed=42)
    g = core.Galaxy(m_s=(mass_s * u.M_sun), x_s=(pos_s[:, 0] * u.kpc), y_s=(pos_s[:, 1] * u.kpc), z_s=(pos_s[:, 2] * u.kpc), vx_s=(vel_s[:, 0] * (u.km / u.s)), vy_s=(vel_s[:, 1] * (u.km / u.s)), vz_s=(vel_s[:, 2] * (u.km / u.s)), m_dm=(mass_dm * u.M_sun), x_dm=(pos_dm[:, 0] * u.kpc), y_dm=(pos_dm[:, 1] * u.kpc), z_dm=(pos_dm[:, 2] * u.kpc), vx_dm=(vel_dm[:, 0] * (u.km / u.s)), vy_dm=(vel_dm[:, 1] * (u.km / u.s)), vz_dm=(vel_dm[:, 2] * (u.km / u.s)), m_g=(mass_g * u.M_sun), x_g=(pos_g[:, 0] * u.kpc), y_g=(pos_g[:, 1] * u.kpc), z_g=(pos_g[:, 2] * u.kpc), vx_g=(vel_g[:, 0] * (u.km / u.s)), vy_g=(vel_g[:, 1] * (u.km / u.s)), vz_g=(vel_g[:, 2] * (u.km / u.s)))
    return g
Mock galaxy.
@pytest.fixture
def mock_galaxy(disc_particles_all, halo_particles):
    
    (mass_s, pos_s, vel_s, mass_g, pos_g, vel_g) = disc_particles_all
    (mass_dm, pos_dm, vel_dm) = halo_particles(N_part=100, seed=42)
    g = core.Galaxy(m_s=(mass_s * u.M_sun), x_s=(pos_s[:, 0] * u.kpc), y_s=(pos_s[:, 1] * u.kpc), z_s=(pos_s[:, 2] * u.kpc), vx_s=(vel_s[:, 0] * (u.km / u.s)), vy_s=(vel_s[:, 1] * (u.km / u.s)), vz_s=(vel_s[:, 2] * (u.km / u.s)), m_dm=(mass_dm * u.M_sun), x_dm=(pos_dm[:, 0] * u.kpc), y_dm=(pos_dm[:, 1] * u.kpc), z_dm=(pos_dm[:, 2] * u.kpc), vx_dm=(vel_dm[:, 0] * (u.km / u.s)), vy_dm=(vel_dm[:, 1] * (u.km / u.s)), vz_dm=(vel_dm[:, 2] * (u.km / u.s)), m_g=(mass_g * u.M_sun), x_g=(pos_g[:, 0] * u.kpc), y_g=(pos_g[:, 1] * u.kpc), z_g=(pos_g[:, 2] * u.kpc), vx_g=(vel_g[:, 0] * (u.km / u.s)), vy_g=(vel_g[:, 1] * (u.km / u.s)), vz_g=(vel_g[:, 2] * (u.km / u.s)))
    return g