The problem is that the shapes of L and U depend on the input matrix and it need not be square. I'm very tempted to allow Triangular to be rectangular so actually the triangular matrix would be more a trapezoid. I think it could work but would require some efforts. It would also make it it possible to use the structures for blocked/paneled factorization algorithm.
2015-01-15 9:33 GMT-05:00 Tamas Papp <[email protected]>: > I am still learning about matrices in Julia so I didn't file an issue > becuase I am not sure this is a bug, but in > > julia> VERSION > v"0.3.5" > > julia> lu = lufact([1 2; 3 4]) > LU{Float64,Array{Float64,2}}(2x2 Array{Float64,2}: > 3.0 4.0 > 0.333333 0.666667,Int32[2,2],0) > > julia> typeof(lu[:U]) > Array{Float64,2} > > julia> typeof(lu[:L]) > Array{Float64,2} > > wouldn't it make sense for the :U and :L accessors to return triangular > matrices? > > Best, > > Tamas >
