Many physical problems are governed by constrained hyperbolic systems including Maxwell's equations of electromagnetism and Einstein's equations of gravitation. Numerically solving these systems accurately requires the development of special methods to guarantee constraint satisfaction. This talk will discuss various ways to augment standard and least-squares finite element methods to aid in this goal.