Error in JD found in prototype: on step 3 t*u is not orthogonal, and it causes residual to blow up. 5x5 matrix has shown itself to be the ideal testing size.