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 subsequent sections. The review concludes with some recent developments and future directions.