Hi,
I think the changes needed to resolve this bug reside in the resize function of identity_matrices. The patch below makes the bug posted code work well with the expected results. Regression tests pass (although I did not look if there are tests for the identity_matrix in there).
Very Best
Nasos Iliopoulos
HotmailŪ is up to 70% faster. Now good news travels really fast. Try it now.