The density matrix renormalization group theory is reviewed as a numerical
variational method. The variational state, expressed as a product of local
tensors, is improved through locally tuning each tensor. The first section
is a tutorial with simplified discussions. Details are discussed in the sub
sequent sections. The review concludes with some recent developments and fu
ture directions.