STL like template based coding with the MMX/SSE extension

#include "stllcv/iplvhistogram.hxx"
#include "stllcv/iplvbasicimage.hxx"
#include "highgui.h"

#define PEXORGIMG unsigned char

#if !DebugFileRead
#define Filename argv[1]
#define Filename2 argv[2] 
#define OutImageName argv[3] 
#endif

float iplvMutualInformation(
   CiplvHistogram<1> &hist1
   ,CiplvHistogram<1> &hist2
   ,CiplvHistogram<2> &hist12)
{
	float sum=0;
	float N=hist1.size(0);
	vigra::MultiIterator< 1, float> i1=hist1.traverser_begin() ,iend1=hist1.traverser_end();
		//vigra::MultiIterator< 2, float> i;
        //i3.template dim<1>()++;  // increment outer dimension
        //i3.template dim<0>()++;  // increment inner dimension
	vigra::MultiIterator<2, float> k2 = hist12.traverser_begin(),kend2 =hist12.traverser_end();
    for(; k2 != kend2; ++k2, ++i1)
    {
		vigra::MultiIterator< 1, float> j1=hist2.traverser_begin() ,jend1=hist2.traverser_end();
        vigra::MultiIterator<2, float>::next_type k1 = k2.begin(), kend1 = k2.end();
        for(; k1 != kend1; ++k1, ++j1)
        {
			if(k1[0]>0 && i1[0]>0 && j1[0]>0){
				sum += k1[0]/N * log( k1[0]*N/(i1[0])*(j1[0]) ); // do something with the current element
			}
        }
    }
	return sum;
}


int main( int argc, char** argv )
{


	
	CiplvBasicImage<PEXORGIMG> image1("lena.jpg");
	CiplvBasicImage<PEXORGIMG> image2("lena2.jpg");


	int hist_size = 256; 
	float range_0[]={0,255};
	float* ranges[]={ range_0 };

	CiplvHistogram<1> hist1(&hist_size,ranges);
	CiplvHistogram<1> hist2(&hist_size,ranges);

	cvCalcHist( &image1.pIplImage, hist1.pCvHistogram, 0, 0 );
	cvCalcHist( &image2.pIplImage, hist2.pCvHistogram, 0, 0 );



	float range_1[]={0,255};
	float* ranges2[]={ range_0 ,range_1 };
	int hist_sizes[] = {256,256};
	CiplvHistogram<2> hist12(hist_sizes,ranges2);
	float mutInf= iplvMutualInformation(hist1,hist2,hist12);

	return 0;
}