On Tuesday, April 1, 2014 10:44:03 PM UTC-5, Gustavo Lacerda wrote:
>
> yes, I think that would be useful.
>
Do you know what Matlab's orth function returns for a rank deficient
matrix? Suppose you ask for an orthogonal basis of the column space of a
matrix like
julia> hcat(fill(1.,(20,)), eye(4)[iceil([1:20] ./ 5),:])
20x5 Array{Float64,2}:
1.0 1.0 0.0 0.0 0.0
1.0 1.0 0.0 0.0 0.0
1.0 1.0 0.0 0.0 0.0
1.0 1.0 0.0 0.0 0.0
1.0 1.0 0.0 0.0 0.0
1.0 0.0 1.0 0.0 0.0
1.0 0.0 1.0 0.0 0.0
1.0 0.0 1.0 0.0 0.0
1.0 0.0 1.0 0.0 0.0
1.0 0.0 1.0 0.0 0.0
1.0 0.0 0.0 1.0 0.0
1.0 0.0 0.0 1.0 0.0
1.0 0.0 0.0 1.0 0.0
1.0 0.0 0.0 1.0 0.0
1.0 0.0 0.0 1.0 0.0
1.0 0.0 0.0 0.0 1.0
1.0 0.0 0.0 0.0 1.0
1.0 0.0 0.0 0.0 1.0
1.0 0.0 0.0 0.0 1.0
1.0 0.0 0.0 0.0 1.0
Do you get 4 or 5 columns?
The simple calculation gives you 5 columns but there should only be 4
columns.
