We review constitutive modeling of solutions and melts of linear polymers, focusing on changes in rheological behavior in shear and extensional flow as the concentration increases from unentangled dilute, to entangled, to dense melt. The rheological changes are captured by constitutive equations, prototypes of which are the FENE-P model for unentangled solutions and the DEMG model for entangled solutions and melts. From these equations, and supporting experimental data, for dilute solutions, the extensional viscosity increases with the strain rate from the low-strain rate to the high-strain rate asymptote, but in the densely entangled state, the high-strain rate viscosity is lower than the low-shear rate value, especially when orientation-dependent friction is accounted for. In shearing flow, shear thinning increases dramatically as the entanglement density increases, which can eventually lead to a shear-banding inhomogeneity. Recent improvements in constitutive modeling are paving the way for robust and accurate numerical simulations of polymer fluid mechanics and industrial processing of polymers.