Steboss

cpdef rectangle(np.ndarray[double, ndim=3, mode="c"] boxes,
               np.ndarray[double, ndim=2, mode="c"] points):
  #as input we have 2 arrays, one with the boxes  cartesian coords
  #the other with points cartesian coords (2D)
  #mode C ensure you have a one D array
  #boxes and points generated in python
  #now we need C to compute the number of points within the rectangle
  #transform into memory view the arrays
  B, M = boxes.shape[0], boxes.shape[1]*2  #B number of rectangles, M rect coords
  N, y = points.shape[0], points.shape[1]

Steboss

all:
	cython rect_in.pyx
	gcc -g -O3 -fpic -c rect_in.c `python-config --cflags`  -std=c99
	gcc rect_in.o -o rect_in.so -shared `python-config --ldflags` -std=c99

Steboss

double kendalT(double *x, int len_x, double * y, int len_y)
{
    //the two array MUST be have the same length
    int n2=0, n1=0, k, j;
    int is =0;
    double aa, a2, a1;

    for(j=1;j<len_x;j++){
      for(k=(j+1);k<=len_x;k++){
        a1=x[j]-x[k];

Steboss

Steboss

Steboss

Steboss

Steboss

Steboss

Steboss