Toward Efficient Simulations of Highly Magnetized Plasma