I've submitted some changes on the matrix class (and dependencies) today, so it now uses jagged arrays, as you suggested. Indeed, I observed nearly 50% perf improvement on the included small Solve perf analysis project. You may want to have a look at the changes first before updating your local SVN...