My project uses QR factorization on large matrices to quantify expansion and contraction rates of phase space volumes of dynamical systems.  I am hoping the speedup on for loops and matrix factorization will greatly speed up runs.