We propose a general and highly efficient method for solving and estimating general equilibrium heterogeneous-agent models with aggregate shocks in discrete time. Our approach relies on the rapid computation and composition of sequence-space Jacobians--the derivatives of perfect-foresight...