Density-matrix based numerical methods for discovering order and correlations in interacting systems