STL like template based coding with the MMX/SSE extension

/************************************************************************/
/*                                                                      */
/*               Copyright 2004-2005 by Hirotaka Niitsuma               */
/*                                                                      */
/*                                                                      */
/*    This file is part of the STLLCV computer vision library.          */
/*    ( Version 0.1, Sep 28 2005 )                                      */
/*    You may use, modify, and distribute this software according       */
/*    to the terms stated in the LICENSE file included in               */
/*    the STLLCV distribution.                                          */
/*                                                                      */
/*    The STLLCV Website is                                             */
/*        http://www2s.biglobe.ne.jp/~niitsuma/STLLCV/                  */
/*    Please direct questions, bug reports, and contributions to        */
/*        niitsuma@mub.biglobe.ne.jp�@�@�@                              */
/*                                                                      */
/*  THIS SOFTWARE IS PROVIDED AS IS AND WITHOUT ANY EXPRESS OR          */
/*  IMPLIED WARRANTIES, INCLUDING, WITHOUT LIMITATION, THE IMPLIED      */
/*  WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE. */
/*                                                                      */
/************************************************************************/

#include "stllcv/ublascvmatrix.hxx"

using namespace stllcv; 

int main(int argc, char * argv[])
{
	//init as identity_matrix
	CWrapCvMat<boost::numeric::ublas::matrix<float> > mat((boost::numeric::ublas::matrix<float>)boost::numeric::ublas::identity_matrix<float>(4));
	std::cout << "mat= " << mat  << std::endl;
	
	

	CWrapCvMat<boost::numeric::ublas::matrix<float> > matI((boost::numeric::ublas::matrix<float>)boost::numeric::ublas::identity_matrix<float>(4)); 

	//multiple as uBLAS
	mat *= 2 ; 
	std::cout << "mat= " << mat  << std::endl;

	cvmInvert(&(CvMat)mat, &(CvMat)matI );
	std::cout << "matI= " << matI  << std::endl;


	cvInvert(mat.pCvMat, matI.pCvMat );
	std::cout << "matI= " << matI  << std::endl;




	return 0;
}