wolffd@0: 
wolffd@0: #include "mex.h"
wolffd@0: #include <stdio.h>
wolffd@0: 
wolffd@0: #define MAX(x,y)  ((x)>(y) ? (x) : (y))
wolffd@0: #define MIN(x,y)  ((x)<(y) ? (x) : (y))
wolffd@0: 
wolffd@0: void mexFunction(
wolffd@0:         int nlhs,       mxArray *plhs[],
wolffd@0:         int nrhs, const mxArray *prhs[]
wolffd@0:         )
wolffd@0: {
wolffd@0:   int j,k,m,n,nzmax,*irs,*jcs, *irs2, *jcs2;
wolffd@0:   double *overlap, *overlap2, tmp, areaA, areaB;
wolffd@0:   double *leftA, *rightA, *topA, *bottomA;
wolffd@0:   double *leftB, *rightB, *topB, *bottomB;
wolffd@0:   double *verbose;
wolffd@0: 
wolffd@0:   m = MAX(mxGetM(prhs[0]), mxGetN(prhs[0]));
wolffd@0:   n = MAX(mxGetM(prhs[4]), mxGetN(prhs[4]));
wolffd@0:   /* printf("A=%d, B=%d\n", m, n); */
wolffd@0: 
wolffd@0:   leftA = mxGetPr(prhs[0]);
wolffd@0:   rightA = mxGetPr(prhs[1]);
wolffd@0:   topA = mxGetPr(prhs[2]);
wolffd@0:   bottomA = mxGetPr(prhs[3]);
wolffd@0: 
wolffd@0:   leftB = mxGetPr(prhs[4]);
wolffd@0:   rightB = mxGetPr(prhs[5]);
wolffd@0:   topB = mxGetPr(prhs[6]);
wolffd@0:   bottomB = mxGetPr(prhs[7]);
wolffd@0: 
wolffd@0:   verbose = mxGetPr(prhs[8]);
wolffd@0: 
wolffd@0:   plhs[0] = mxCreateDoubleMatrix(m,n, mxREAL);
wolffd@0:   overlap  = mxGetPr(plhs[0]);
wolffd@0: 
wolffd@0:   plhs[1] = mxCreateDoubleMatrix(m,n, mxREAL);
wolffd@0:   overlap2  = mxGetPr(plhs[1]);
wolffd@0: 
wolffd@0:   k = 0; 
wolffd@0:   for (j = 0; (j < n); j++) {
wolffd@0:     int i;
wolffd@0:     for (i = 0; (i < m); i++) {
wolffd@0:       tmp = (MAX(0, MIN(rightA[i], rightB[j]) - MAX(leftA[i], leftB[j]) )) * 
wolffd@0: 	(MAX(0, MIN(topA[i], topB[j]) - MAX(bottomA[i], bottomB[j]) ));
wolffd@0: 
wolffd@0:       if (tmp > 0) {
wolffd@0: 	overlap[k] = tmp;
wolffd@0: 	
wolffd@0: 	areaA = (rightA[i]-leftA[i])*(topA[i]-bottomA[i]);
wolffd@0: 	areaB = (rightB[j]-leftB[j])*(topB[j]-bottomB[j]);
wolffd@0: 	overlap2[k] = tmp/MIN(areaA, areaB);
wolffd@0: 	
wolffd@0: 	if (*verbose) {
wolffd@0: 	  printf("j=%d,i=%d,overlap=%5.3f, norm=%5.3f\n", j,i, overlap[k], overlap2[k]);
wolffd@0: 	}
wolffd@0:       }
wolffd@0: 
wolffd@0:       k++;
wolffd@0:     }
wolffd@0:   }
wolffd@0: }
wolffd@0: 
wolffd@0: 
wolffd@0: 
wolffd@0: 
wolffd@0: 
wolffd@0: 
wolffd@0: 
wolffd@0: