diff options
Diffstat (limited to 'sourcecodes/bnt-master/KPMtools/ind2subv.m')
| -rw-r--r-- | sourcecodes/bnt-master/KPMtools/ind2subv.m | 22 |
1 files changed, 22 insertions, 0 deletions
diff --git a/sourcecodes/bnt-master/KPMtools/ind2subv.m b/sourcecodes/bnt-master/KPMtools/ind2subv.m new file mode 100644 index 00000000..6d9c4308 --- /dev/null +++ b/sourcecodes/bnt-master/KPMtools/ind2subv.m @@ -0,0 +1,22 @@ +function sub = ind2subv(siz,index) +%IND2SUBV Subscript vector from linear index. +% IND2SUBV(SIZ,IND) returns a vector of the equivalent subscript values +% corresponding to a single index into an array of size SIZ. +% If IND is a vector, then the result is a matrix, with subscript vectors +% as rows. + +%sub = ind2subvTest(siz, index); +n = length(siz); +if n==0, sub = index; return; end % added by KPM 17 Nov 07 +cum_size = cumprod(siz(:)'); +prev_cum_size = [1 cum_size(1:end-1)]; +index = index(:) - 1; +sub = rem(repmat(index,1,n),repmat(cum_size,length(index),1)); +sub = floor(sub ./ repmat(prev_cum_size,length(index),1))+1; + +% slow way +%for dim = n:-1:1 +% sub(:,dim) = floor(index/cum_size(dim))+1; +% index = rem(index,cum_size(dim)); +%end + |