Guard against non-progressing blocking size estimation.
New functions for QL factored matrices.
Fix wrong flag bit in WY-update. Added some comments.
New tests.