Moved asym and apd in test/linalg/svd.jl.

JuliaLang · Mar 4, 2016 · 81ded47 · 81ded47
1 parent ef3264a
commit 81ded47
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/test/linalg/svd.jl b/test/linalg/svd.jl
@@ -21,16 +21,18 @@ a2img  = randn(n,n)/2
 for eltya in (Float32, Float64, Complex64, Complex128, Int)
     aa = eltya == Int ? rand(1:7, n, n) : convert(Matrix{eltya}, eltya <: Complex ? complex(areal, aimg) : areal)
     aa2 = eltya == Int ? rand(1:7, n, n) : convert(Matrix{eltya}, eltya <: Complex ? complex(a2real, a2img) : a2real)
+    asym = aa'+aa                 # symmetric indefinite
+    apd  = aa'*aa                 # symmetric positive-definite
     for atype in ("Array", "SubArray")
         if atype == "Array"
             a = aa
             a2 = aa2
         else
             a = sub(aa, 1:n, 1:n)
             a2 = sub(aa2, 1:n, 1:n)
+            asym = sub(asym, 1:n, 1:n)
+            apd = sub(apd, 1:n, 1:n)
         end
-        asym = a'+a                  # symmetric indefinite
-        apd  = a'*a                 # symmetric positive-definite
         ε = εa = eps(abs(float(one(eltya))))
 
     debug && println("\ntype of a: ", eltya, "\n")