The cuprate high-temperature superconductors hold great technological promise. Non-perturbative, massively parallel simulations have played an essential role in developing an understanding of these materials, from establishing the viability of the most basic cuprate model, to the inclusion of realistic effects such as lattice vibrations, and disorder. New hybrid parallel codes make it possible to calculate the phase diagram of the Hubbard model and explore the effect of higher order bond spin correlations. A fuller understanding of these realistic effects may lead to the development of better superconducting materials.