function factorize(A::SparseMatrixCSC) m, n = size(A) if m == n if istril(A) if istriu(A) return return Diagonal(A)