Add a check for blocksquarity

Author: kshyatt

src/linalg/factorizations.jl

@@ -16,6 +16,8 @@ function MatrixAlgebraKit.zero!(A::BlockBlasMat)
 end
 
 function MatrixAlgebraKit.one!(A::BlockBlasMat)
+    mb, nb = blocksize(A)
+    mb == nb || throw(DimensionMismatch("A is not block-square. Number of row-blocks ($mb) does not match number of column-blocks ($nb)"))
     for bj in blockaxes(A, 2), bi in blockaxes(A, 1)
         a = view(A, bi, bj)
         bi == bj ? MAK.one!(a) : MAK.zero!(a)