Merge pull request #5 from andreasnoack/master

Update syntax and test setup

Author: simonbyrne

.travis.yml

@@ -4,9 +4,9 @@ os:
   - linux
   - osx
 julia:
-  - 0.4
   - 0.5
-#  - nightly
+  - 0.6
+  - nightly
 notifications:
   email: false
 # uncomment the following lines to override the default test script

appveyor.yml

@@ -2,6 +2,8 @@ environment:
   matrix:
   - JULIA_URL: "https://julialang-s3.julialang.org/bin/winnt/x86/0.5/julia-0.5-latest-win32.exe"
   - JULIA_URL: "https://julialang-s3.julialang.org/bin/winnt/x64/0.5/julia-0.5-latest-win64.exe"
+  - JULIA_URL: "https://julialang-s3.julialang.org/bin/winnt/x86/0.6/julia-0.6-latest-win32.exe"
+  - JULIA_URL: "https://julialang-s3.julialang.org/bin/winnt/x64/0.6/julia-0.6-latest-win64.exe"
   - JULIA_URL: "https://julialangnightlies-s3.julialang.org/bin/winnt/x86/julia-latest-win32.exe"
   - JULIA_URL: "https://julialangnightlies-s3.julialang.org/bin/winnt/x64/julia-latest-win64.exe"
 

src/GenericSVD.jl

@@ -100,7 +100,7 @@ function svd!{T<:Real}(B::Bidiagonal{T}, U=nothing, Vt=nothing, ɛ=eps(T))
     n = size(B, 1)
     n₂ = n
 
-    maxB = max(maxabs(B.dv),maxabs(B.ev))
+    maxB = max(maximum(abs, B.dv), maximum(abs, B.ev))
 
     if istriu(B)
         while true

src/bidiagonalize.jl

@@ -46,8 +46,8 @@ end
 function bidiagonalize_tall!{T}(A::Matrix{T})
     m,n = size(A)
     R = real(T)
-    B = Bidiagonal{R}(Array(R,n),Array(R,n-1),true)
-    bidiagonalize_tall!(A,B)
+    B = Bidiagonal{R}(Vector{R}(n), Vector{R}(n - 1), true)
+    bidiagonalize_tall!(A, B)
 end
 
 function Base.full{T}(P::PackedUVt{T};thin=true)

test/bigfloat.jl

@@ -6,13 +6,13 @@ srand(1)
 n,m = 100,20
 
 X = randn(n,m)
-bX = big(X)
+bX = big.(X)
 bS = svdfact(bX)
 @test isapprox(full(bS), bX, rtol=1e3*eps(BigFloat))
 @test isapprox(svdvals(bX), svdvals(X), rtol=1e3*eps())
 @test bX == X # check we didn't modify the input
 
-bY = big(randn(n))
+bY = big.(randn(n))
 @test isapprox(qrfact(bX,Val{false}) \ bY, bS \ bY, rtol=1e3*eps(BigFloat))
 @test bX == X # check we didn't modify the input
 
@@ -23,32 +23,32 @@ bSt = svdfact(bXt)
 @test bXt == X' # check we didn't modify the input
 
 X = Float64[1 2 0; 0 1 2; 0 0 0]
-bX = big(X)
+bX = big.(X)
 bS = svdfact(bX)
 @test isapprox(full(bS), bX, rtol=1e3*eps(BigFloat))
 @test isapprox(svdvals(bX), svdvals(X), rtol=1e3*eps())
 @test bX == X # check we didn't modify the input
 
 X = Float64[0 2 0; 0 1 2; 0 0 1]
-bX = big(X)
+bX = big.(X)
 bS = svdfact(bX)
 @test isapprox(full(bS), bX, rtol=1e3*eps(BigFloat))
 @test isapprox(svdvals(bX), svdvals(X), rtol=1e3*eps())
 @test bX == X # check we didn't modify the input
 
 
-bD = big(randn(m))
+bD = big.(randn(m))
 bX = diagm(bD)
 bS = svdfact(bX)
 
 @test isapprox(full(bS), bX, rtol=1e3*eps(BigFloat))
-@test bS.S == sort(abs(bD),rev=true)
+@test bS.S == sort(abs.(bD),rev=true)
 
 
 
 
 X = randn(n,m)+im*randn(n,m)
-bX = big(X)
+bX = big.(X)
 bS = svdfact(bX)
 @test isapprox(full(bS), bX, rtol=1e3*eps(BigFloat))
 @test isapprox(svdvals(bX), svdvals(X), rtol=1e3*eps())
@@ -62,14 +62,14 @@ bSt = svdfact(bXt)
 
 
 X = Complex128[1 2 0; 0 1 2; 0 0 0]
-bX = big(X)
+bX = big.(X)
 bS = svdfact(bX)
 @test isapprox(full(bS), bX, rtol=1e3*eps(BigFloat))
 @test isapprox(svdvals(bX), svdvals(X), rtol=1e3*eps())
 @test bX == X # check we didn't modify the input
 
 X = Complex128[0 2 0; 0 1 2; 0 0 1]
-bX = big(X)
+bX = big.(X)
 bS = svdfact(bX)
 @test isapprox(full(bS), bX, rtol=1e3*eps(BigFloat))
 @test isapprox(svdvals(bX), svdvals(X), rtol=1e3*eps())

test/misc.jl

@@ -9,7 +9,7 @@ c = randn()
 U = [a c; 0 b]
 x,y = GenericSVD.svdvals2x2(a,b,c)
 
-@test sort(sqrt(eigvals(U'*U))) ≈ [x,y]
+@test sort(sqrt.(eigvals(U'*U))) ≈ [x,y]
 @test sort(svdvals(U)) ≈ [x,y]
 
 U = eye(3)