octave/Mat2Hrows.m


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28

function [Hrows Hcols] = Mat2Hrows(HRA, RAmode); 

H = full(HRA); 
[Nr Nc] = size(H);

if nargin==1, RAmode = 1; end; 

if RAmode==1
   H= H(:,1:Nc-Nr); 
   [Nr Nc] = size(H); 
end 

Max_colwt = max(sum(H));
Max_rowwt = max(sum(H')); 
Hcols = zeros(Nc, Max_colwt); 
Hrows = zeros(Nr, Max_rowwt); 

for i = 1:Nr
    nz = find(H(i,:)); 
    Hrows(i,1:length(nz)) = nz; 
end

H = H'; 
for i = 1:Nc
    nz = find(H(i,:)); 
    Hcols(i,1:length(nz)) = nz; 
end