def METHOD_NAME(self):
    no_virtual = self.system.part.select(virtual=False)
    P_kin = np.sum(
        no_virtual.mass * np.sum(no_virtual.v**2, axis=1),
        axis=0) / (3 * self.system.volume())
    np.testing.assert_allclose(
        P_kin, self.system.analysis.pressure()["kinetic"])
    analyze_pressure_tensor = np.diag(
        self.system.analysis.pressure_tensor()["kinetic"])
    expected_pressure_tensor = np.matmul(
        (no_virtual.v**2).T,
        no_virtual.mass) / (self.system.volume())
    np.testing.assert_allclose(
        expected_pressure_tensor, analyze_pressure_tensor)