modified_integralImage2D2D.c

/********************************
Author: Sravanthi Kota Venkata
********************************/

#include <stdio.h>
#include <stdlib.h>
#include "disparity.h"
//#include "m5op.h"
void integralImage2D2D(F2D* SAD, F2D* integralImg)
{
    int nr, nc, i, j;

    nr = SAD->height;
    nc = SAD->width;
    int togo,size;
    for(i=0; i<nc; i++)
        subsref(integralImg,0,i) = subsref(SAD,0,i);

      /* original code
      for(i=1; i<nr; i++)
            for(j=0; j<nc; j++)
            {
                subsref(integralImg,i,j) = subsref(integralImg, (i-1), j) + subsref(SAD,i,j);
            }
       */


      int N=nr*nc;
      //int N=(nr*nc)-nc;
      for(i=nr;i<N; i+=8)
      {
        togo = N - i;
        size = (togo > 8) ? 8 : togo;

        for(j=0; j<size; j++)
        {
            asubsref(integralImg,i+j) = asubsref(integralImg, i-nc +j) + asubsref(SAD,i+j);
        }

      }



     /*original code
      for(i=0; i<nr; i++)
        for(j=1; j<nc; j++)
            subsref(integralImg,i,j) = subsref(integralImg, i, (j-1)) + subsref(integralImg,i,j);
       */

      for(i=0; i<N; i+=8)
            {

              togo = N - i;
              size = (togo > 8) ? 8 : togo;
              for(j=0; j<size; j++)
              {
                    if (((i+j) % nc) != 0)
                       {
                          asubsref(integralImg,i+j) = asubsref(integralImg, i+j-1) + asubsref(integralImg,i+j);

                       }
              }

            }

    return;

}