How can we overcome numerical optimization’s high computational complexity, while still capturing complex dynamics and providing reliable convergence for field robots? This project seeks to answer that question by co-designing new theoretically sound algorithms that are optimized to take advantage of the large-scale parallelism available on GPUs. Through support from the NSF this project seeks to go beyond developing point solutions to release a broadly applicable toolbox for the robotics and optimization communities.