def METHOD_NAME(self):
    # Particles of type 0
    I0 = self.i_tensor(self.system.part.select(
        lambda p: (not p.virtual) and p.type == 0).id)
    np.testing.assert_allclose(
        I0, self.system.analysis.moment_of_inertia_matrix(p_type=0), atol=1E-9)
    # Particles of type 1
    I1 = self.i_tensor(self.system.part.select(type=1).id)
    np.testing.assert_allclose(
        I1, self.system.analysis.moment_of_inertia_matrix(p_type=1), atol=1E-9)