From 1594a1c8d2c4887d0a6cfa445aed8d0cfef93e0a Mon Sep 17 00:00:00 2001 From: Rafael Fourquet Date: Mon, 25 Dec 2017 16:21:04 +0100 Subject: [PATCH] LinearAlgebra: delete _rand_pm1! in favor of rand! --- stdlib/SparseArrays/src/linalg.jl | 10 ++-------- 1 file changed, 2 insertions(+), 8 deletions(-) diff --git a/stdlib/SparseArrays/src/linalg.jl b/stdlib/SparseArrays/src/linalg.jl index 6e9a751dd4306..48f7e00ae54fb 100644 --- a/stdlib/SparseArrays/src/linalg.jl +++ b/stdlib/SparseArrays/src/linalg.jl @@ -735,12 +735,6 @@ function opnormestinv(A::SparseMatrixCSC{T}, t::Integer = min(2,maximum(size(A)) S = zeros(T <: Real ? Int : Ti, n, t) - function _rand_pm1!(v) - for i in eachindex(v) - v[i] = rand()<0.5 ? 1 : -1 - end - end - function _any_abs_eq(v,n::Int) for vv in v if abs(vv)==n @@ -755,7 +749,7 @@ function opnormestinv(A::SparseMatrixCSC{T}, t::Integer = min(2,maximum(size(A)) X[1:n,1] .= 1 for j = 2:t while true - _rand_pm1!(view(X,1:n,j)) + rand!(view(X,1:n,j), (-1, 1)) yaux = X[1:n,j]' * X[1:n,1:j-1] if !_any_abs_eq(yaux,n) break @@ -816,7 +810,7 @@ function opnormestinv(A::SparseMatrixCSC{T}, t::Integer = min(2,maximum(size(A)) end end if repeated - _rand_pm1!(view(S,1:n,j)) + rand!(view(S,1:n,j), (-1, 1)) else break end