From d71156ce6d37958c811d079819acbb8e80187fd2 Mon Sep 17 00:00:00 2001 From: lkdvos Date: Tue, 13 Aug 2024 11:54:56 +0200 Subject: [PATCH] Remove reference to undefined tensorkit function --- src/auxiliary.jl | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/auxiliary.jl b/src/auxiliary.jl index 8f391fb..d5af541 100644 --- a/src/auxiliary.jl +++ b/src/auxiliary.jl @@ -83,7 +83,7 @@ function _polarnewton!(A::StridedMatrix; tol=10 * scalareps(A), maxiter=5) @assert m >= n A2 = copy(A) Q, R = qr!(A2) - Ri = ldiv!(UpperTriangular(R)', TensorKit._one!(similar(R))) + Ri = ldiv!(UpperTriangular(R)', TensorKit.MatrixAlgebra.one!(similar(R))) R, Ri = _avgdiff!(R, Ri) i = 1 R2 = view(A, 1:n, 1:n) @@ -94,7 +94,7 @@ function _polarnewton!(A::StridedMatrix; tol=10 * scalareps(A), maxiter=5) _polarsdd!(Ri) break end - Ri = ldiv!(lu!(R2)', TensorKit._one!(Ri)) + Ri = ldiv!(lu!(R2)', TensorKit.MatrixAlgebra.one!(Ri)) R, Ri = _avgdiff!(R, Ri) copyto!(R2, R) i += 1 @@ -163,7 +163,7 @@ function _stiefellog(Wold::StridedMatrix, Wnew::StridedMatrix; R = Q' * dW Wext = [Wold Q] F = qr!([P; R]) - U = lmul!(F.Q, TensorKit._one!(similar(P, r, r))) + U = lmul!(F.Q, TensorKit.MatrixAlgebra.one!(similar(P, r, r))) U[1:p, 1:p] .= P U[(p + 1):r, 1:p] .= R X = view(U, 1:p, (p + 1):r)