By the way. The same thing is true for R in the QR. 2015-01-15 9:39 GMT-05:00 Andreas Noack <[email protected]>:
> 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 >> > >
