Time domain localized-density-matrix method