Please don’t post code whose copyright status is not 100% certain. If you are even vaguely unsure, you are putting the community in needless danger.
As a minimal example, it may be illegal for you to have written the e-mail you just wrote. — John On Jun 30, 2014, at 7:07 AM, 'Stéphane Laurent' via julia-users <julia-users@googlegroups.com> wrote: > A couple of years ago I used this function (maybe I copied the Matlab or > Octave code, I don't remember): > > ## similar to the Matlab orth() function > function orth(A) > if (isempty (A)) > retval = [] > else > (U, S, V) = svd(A) > (rows, cols) = size(A) > tol = max (size (A)) * S[1] * eps() > r = sum (S .> tol) > if (r >0) > retval = -U[:,1:r] > else > retval = zeros(rows, 0); > end > end > return(retval) > end
