summaryrefslogtreecommitdiff
path: root/src/auxiliaryFunctions/find2d/dfind2da.c
blob: 1358b5d439b41d19d4b064fdcb3d3412a6d8acb7 (plain)
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
29
30
31
32
33
34
35
36
37
/*
 *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 *  Copyright (C) 2007-2008 - INRIA - Arnaud TORSET
 *
 *  This file must be used under the terms of the CeCILL.
 *  This source file is licensed as described in the file COPYING, which
 *  you should have received as part of this distribution.  The terms
 *  are also available at
 *  http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 *
 */

#include "find2d.h"
#include <malloc.h>

void dfind2da(double* x, int rows, int columns, double** out1, double** out2) {
  int i = 0, j=0;
  int indiceOut = 0;
  double* out_1 = NULL, *out_2=NULL;
  
  /* We suppose out1 and out2 are NULL on input */
  
  for (i = 0; i < rows ; ++i) {
  	for	(j = 0; j < columns ; ++j) {
	    if (x[i*columns+j] != 0) {
	    	out_1 = realloc(out_1, (uint)(indiceOut+1)*sizeof(double));
	    	out_2 = realloc(out_2, (uint)(indiceOut+1)*sizeof(double));
	    	
		out_1[indiceOut] = (double)i;
		out_2[indiceOut] = (double)j;
		indiceOut++;
	    }
	}
  }
  *out1=out_1;
  *out2=out_2;
}