Fixed loop ranges and array indexing for backward substitution routine.
Work on backward/forward substitution -- still seems incorrect.

Check on Finite Cholesky routine -- do not shift on last iteration?
Check on range of rows which need to use finite Cholesky piece (bw or bw-1?)