Skip to content

Commit

Permalink
Remove incorrect condition in LU test. The transpose should be tested…
Browse files Browse the repository at this point in the history
… for complex numbers, and might be skipped for real numbers; not the other way around. For clarity just test both unconditionally.
  • Loading branch information
haampie committed Jul 12, 2017
1 parent 6e87e12 commit efb8db9
Showing 1 changed file with 2 additions and 4 deletions.
6 changes: 2 additions & 4 deletions test/linalg/lu.jl
Original file line number Diff line number Diff line change
Expand Up @@ -83,10 +83,8 @@ debug && println("(Automatic) Square LU decomposition. eltya: $eltya, eltyb: $el
@test norm(a*(lua\c) - c, 1) < ε*κ*n # c is a vector
@test norm(a'*(lua'\c) - c, 1) < ε*κ*n # c is a vector
@test AbstractArray(lua) a
if eltya <: Real && eltyb <: Real
@test norm(a.'*(lua.'\b) - b,1) < ε*κ*n*2 # Two because the right hand side has two columns
@test norm(a.'*(lua.'\c) - c,1) < ε*κ*n
end
@test norm(a.'*(lua.'\b) - b,1) < ε*κ*n*2 # Two because the right hand side has two columns
@test norm(a.'*(lua.'\c) - c,1) < ε*κ*n

# Test whether Ax_ldiv_B!(y, LU, x) indeed overwrites y
resultT = typeof(oneunit(eltyb) / oneunit(eltya))
Expand Down

0 comments on commit efb8db9

Please sign in to comment.