Subversion Repositories gelsvn

Rev

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

Rev 2 Rev 12
Line 12... Line 12...
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 
    int n = vec.size();
16 
 
 16 
 
17 
		MT C;
 17 
    MT C(0);
18 
		for(int i=0;i<n;++i)
 18 
    for(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);