We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent a20c3f4 commit 3bfe9c6Copy full SHA for 3bfe9c6
1 file changed
src/linalg/factorizations.jl
@@ -16,7 +16,6 @@ function MatrixAlgebraKit.zero!(A::BlockBlasMat)
16
end
17
18
function MatrixAlgebraKit.one!(A::BlockBlasMat)
19
- blockaxes(A, 1) == blockaxes(A, 2) || throw(DimensionMismatch("A is not block-square"))
20
for bj in blockaxes(A, 2), bi in blockaxes(A, 1)
21
a = view(A, bi, bj)
22
bi == bj ? MAK.one!(a) : MAK.zero!(a)
0 commit comments