def METHOD_NAME(self):
    for es in self.es_list:
        dos = Dos(es_obj=es, n_bins=100)
        if es.grand_dos_matrix is not None:
            self.assertIsInstance(es.grand_dos_matrix, np.ndarray)
            _, _, _, _, n_orbitals = np.shape(dos.es_obj.grand_dos_matrix)
            total_rdos = np.zeros_like(dos.t_dos)
            for i in range(n_orbitals):
                r_dos = dos.get_orbital_resolved_dos([i])
                total_rdos += r_dos
            self.assertTrue(np.allclose(dos.t_dos, total_rdos))