Merge pull request #37 from YingboMa/cy/fix

No `ivdep` when it's not ivdep

Author: YingboMa

src/lu.jl

@@ -93,7 +93,7 @@ end
 function apply_permutation_threaded!(P, A)
     batchsize = cld(2000, length(P))
     @batch minbatch=batchsize for j in axes(A, 2)
-        @inbounds @simd ivdep for i in axes(P, 1)
+        @inbounds for i in axes(P, 1)
             i′ = P[i]
             tmp = A[i, j]
             A[i, j] = A[i′, j]

test/runtests.jl

@@ -11,7 +11,7 @@ const mylu = RecursiveFactorization.lu
 
 function testlu(A, MF, BF)
     @test MF.info == BF.info
-    @test norm(MF.L*MF.U - A[MF.p, :], Inf) < length(A)*sqrt(eps(real(one(float(first(A))))))/16
+    @test norm(MF.L*MF.U - A[MF.p, :], Inf) < 200sqrt(eps(real(one(float(first(A))))))
     nothing
 end
 testlu(A::Adjoint, MF::Adjoint, BF) = testlu(parent(A), parent(MF), BF)