Neel order in square and triangular lattice Heisenberg models

We show that the density matrix renormalization group can be used to study magnetic ordering in two-dimensional spin models. Local quantities should be extrapolated with the truncation error, not with its square root. We introduce sequences of clusters, using cylindrical boundary conditions with pinning fields, which provide for rapidly converging finite-size scaling. We determine the magnetization for both the square and triangular Heisenberg lattices with errors comparable to the best alternative approaches.