Update src/linalg/factorizations.jl

Co-authored-by: Lukas Devos <ldevos98@gmail.com>

Author: kshyatt

src/linalg/factorizations.jl

@@ -13,7 +13,7 @@ function MatrixAlgebraKit.one!(A::BlockBlasMat)
     n_blocks = blocksize(A)[1]
     # awful workaround to BlockArrays indexing interface
     for bi in 1:n_blocks
-        A[Block(bi), Block(bi)] .= diagm(fill(_one, blocksizes(A)[bi, bi][1]))
+        MAK.one!(A[Block(bi), Block(bi)])
     end
     return A
 end