template <class T>

boost::numeric::ublas::matrix<T> SubstractMeans( boost::numeric::ublas::matrix<T> &m, boost::numeric::ublas::vector<T> &vMeans)

{

ublas::matrix<T> mout(m.size1(), m.size2());

for(size_t i=0;i<m.size1();i++)

for(size_t j=0;j<m.size2();j++)

{

mout(i,j)=m(i, j)-vMeans(i);

}

//cout<<mout<<endl;

return mout;

}

thanks,

Arman.