wolffd@0: function cind = som_ind2cod(msize,ind)
wolffd@0: 
wolffd@0: %SOM_IND2COD SOM_PAK style linear indeces from Matlab linear index.
wolffd@0: %
wolffd@0: % Cind = som_ind2cod(msize,inds)
wolffd@0: %
wolffd@0: %  cind = som_ind2cod([10 15],44);
wolffd@0: %  cind = som_ind2cod(sMap,44);
wolffd@0: %  cind = som_ind2cod(sMap.msize,44);
wolffd@0: %  Cind = som_ind2cod([10 15],[44 13 91]');
wolffd@0: %
wolffd@0: %  Input and output arguments: 
wolffd@0: %   msize  (struct) map or topology struct
wolffd@0: %          (vector) size 1 x m, specifies the map grid size
wolffd@0: %   ind    (vector) size n x 1, linear indeces of n map units
wolffd@0: % 
wolffd@0: %   cind   (matrix) size n x 1, SOM_PAK style linear indeces
wolffd@0: %                   (row first, then column)
wolffd@0: %
wolffd@0: % See also SOM_COD2IND.
wolffd@0: 
wolffd@0: % Contributed to SOM Toolbox vs2, January 14th, 2002 by Juha Vesanto
wolffd@0: % http://www.cis.hut.fi/projects/somtoolbox/
wolffd@0: 
wolffd@0: % Version 2.0beta juuso 140102
wolffd@0: 
wolffd@0: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
wolffd@0: 
wolffd@0: if isstruct(msize), 
wolffd@0:   if strcmp(msize.type,'som_map'), msize = msize.topol.msize; 
wolffd@0:   elseif strcmp(msize.type,'som_topol'), msize = msize.msize;
wolffd@0:   else error('Invalid first argument.'); end
wolffd@0: end
wolffd@0: 
wolffd@0: if nargin<2, ind = 1:prod(msize); end
wolffd@0: 
wolffd@0: Co = som_unit_coords(msize,'rect','sheet');
wolffd@0: 
wolffd@0: switch size(Co,2),
wolffd@0: case 1, I2C = [1:prod(msize)]; 
wolffd@0: case 2, I2C = 1 + Co(:,1) + Co(:,2)*msize(2); 
wolffd@0: case 3, I2C = 1 + Co(:,1) + Co(:,2)*msize(2) + Co(:,3)*msize(1)*msize(2); % ?????
wolffd@0: end
wolffd@0: 
wolffd@0: cind = I2C(ind); 
wolffd@0: return; 
wolffd@0: 
wolffd@0: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%