Skip to content

Commit

Permalink
Remove reference to undefined tensorkit function
Browse files Browse the repository at this point in the history
  • Loading branch information
lkdvos committed Aug 13, 2024
1 parent 7391fe0 commit d71156c
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions src/auxiliary.jl
Original file line number Diff line number Diff line change
Expand Up @@ -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)
Expand All @@ -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
Expand Down Expand Up @@ -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)
Expand Down

0 comments on commit d71156c

Please sign in to comment.