Hi Jed, I think there are two ways to improve the performance of the code. (1) Change G-S to Jacobi. (2) Continue chaining Jacobi to a block version. Fande,