Subversion Repositories gelsvn

Rev

Rev 12 | Rev 431 | Go to most recent revision | Show entire file | Regard whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 12 Rev 45
Line 10... Line 10...
10 
{
 10 
{
11 
  template<class VT, class MT>
 11 
		template<class VT, class MT>
12 
  VT covariance(const vector<VT>& vec, MT& C_out)
 12 
		VT covariance(const vector<VT>& vec, MT& C_out)
13 
  {
 13 
		{
14 
    VT m = mean(vec);
 14 
				VT m = mean(vec);
15 
    int n = vec.size();
 15 
				unsigned int n = vec.size();
16 
 
 16 
 
17 
    MT C(0);
 17 
				MT C(0);
18 
    for(int i=0;i<n;++i)
 18 
				for(unsigned int i=0;i<n;++i)
19 
      {
 19 
				{
20 
	MT B;
 20 
						MT B;
21 
	VT v = vec[i]-m;
 21 
						VT v = vec[i]-m;
22 
	outer_product(v,v,B);
 22 
						outer_product(v,v,B);
23 
 
 23