Skip to content

Commit 0eea203

Browse files
committed
copyto -> copy
1 parent c0cf0b2 commit 0eea203

File tree

2 files changed

+4
-4
lines changed

2 files changed

+4
-4
lines changed

Project.toml

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,6 @@
11
name = "SweepOperator"
22
uuid = "7522ee7d-7047-56d0-94d9-4bc626e7058d"
3-
version = "0.3.1"
3+
version = "0.3.2"
44

55
[deps]
66
LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"

src/SweepOperator.jl

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -40,15 +40,15 @@ function sweep_with_buffer!(akk::AVec{T}, A::AMat{T}, k::Integer, inv::Bool = fa
4040
p == length(akk) || throw(DimensionError("Incorrect buffer size."))
4141
@inbounds @views begin
4242
d = one(T) / A[k, k] # pivot
43-
copyto!(akk, Symmetric(A, :U)[:, k]) # akk = A[:, k]
43+
copy!(akk, Symmetric(A, :U)[:, k]) # akk = A[:, k]
4444
if A isa StridedMatrix{<:Union{LinearAlgebra.BlasFloat, LinearAlgebra.BlasComplex}}
4545
BLAS.syrk!('U', 'N', -d, akk, one(T), A) # everything not in col/row k
4646
else
4747
A .+= UpperTriangular(-d * akk * akk')
4848
end
4949
rmul!(akk, d * (-one(T)) ^ inv) # akk .* d (negated if inv=true)
50-
copyto!(A[1:k-1,k], akk[1:k-1]) # col k
51-
copyto!(A[k, k+1:end], akk[k+1:end]) # row k
50+
copy!(A[1:k-1,k], akk[1:k-1]) # col k
51+
copy!(A[k, k+1:end], akk[k+1:end]) # row k
5252
A[k, k] = -d # pivot element
5353
end
5454
return A

0 commit comments

Comments
 (0)