X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?p=CreaPhase.git;a=blobdiff_plain;f=octave_packages%2Fdataframe-0.9.1%2F%40dataframe%2Ffind.m;fp=octave_packages%2Fdataframe-0.9.1%2F%40dataframe%2Ffind.m;h=46ba45e9bd0c380907562aec5463d16b536e67ce;hp=0000000000000000000000000000000000000000;hb=c880e8788dfc484bf23ce13fa2787f2c6bca4863;hpb=1705066eceaaea976f010f669ce8e972f3734b05

diff --git a/octave_packages/dataframe-0.9.1/@dataframe/find.m b/octave_packages/dataframe-0.9.1/@dataframe/find.m
new file mode 100644
index 0000000..46ba45e
--- /dev/null
+++ b/octave_packages/dataframe-0.9.1/@dataframe/find.m
@@ -0,0 +1,60 @@
+function varargout = find(df, varargin) 
+
+  %% Copyright (C) 2009-2012 Pascal Dupuis <Pascal.Dupuis@uclouvain.be>
+  %%
+  %% This file is part of Octave.
+  %%
+  %% Octave is free software; you can redistribute it and/or
+  %% modify it under the terms of the GNU General Public
+  %% License as published by the Free Software Foundation;
+  %% either version 2, or (at your option) any later version.
+  %%
+  %% Octave is distributed in the hope that it will be useful,
+  %% but WITHOUT ANY WARRANTY; without even the implied
+  %% warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+  %% PURPOSE.  See the GNU General Public License for more
+  %% details.
+  %%
+  %% You should have received a copy of the GNU General Public
+  %% License along with Octave; see the file COPYING.  If not,
+  %% write to the Free Software Foundation, 51 Franklin Street -
+  %% Fifth Floor, Boston, MA 02110-1301, USA.
+
+  %#
+  %# $Id: find.m 9585 2012-02-05 15:32:46Z cdemills $
+  %#
+
+  switch nargout
+    case {0, 1}
+      resu = []; mz = max(cellfun(@length, df._rep));
+      for indc = 1:df._cnt(2),
+	[indr, inds] = feval(@find, df._data{indc}(:, df._rep{indc}));
+	%# create a vector the same size as indr
+	dummy = indr; dummy(:) = indc;
+	resu = [resu; sub2ind([df._cnt(1:2) mz], indr, dummy, inds)];
+      endfor
+      varargout{1} = sort(resu);
+    case 2
+      nz = 0; idx_i = []; idx_j = [];
+      for indc = 1:df._cnt(2),
+	[dum1, dum2] = feval(@find, df._data{indc}(:, df._rep{indc}));
+	idx_i = [idx_i; dum1];
+	idx_j = [idx_j; nz + dum2];
+	nz = nz + df._cnt(1)*length(df._rep{indc});
+      endfor
+      varargout{1} = idx_i; varargout{2} = idx_j;
+    case 3
+      nz = 0; idx_i = []; idx_j = []; val = [];
+      for indc = 1:df._cnt(2),
+	[dum1, dum2, dum3] = feval(@find, df._data{indc}(:, df._rep{indc}));
+	idx_i = [idx_i; dum1];
+	idx_j = [idx_j; nz + dum2];
+	val = [val; dum3];
+	nz = nz + df._cnt(1)*length(df._rep{indc});
+      endfor
+      varargout{1} = idx_i; varargout{2} = idx_j; varargout{3} = val;
+    otherwise
+      print_usage('find');
+  endswitch
+
+endfunction