def _get_covmat(self):
"returns covariance matrix of vector"
if self.cluster == None or self.cluster.covmat == None:
# raise ObsVectorError, "No covariance matrix - no cluster"
return None
if self.index == None:
# raise ObsVectorError, "No row index in covariance matrix"
return None
if len(self.cluster) == 1: # just one vector in cluster
if self.cluster.covmat.dim == 3:
return self.cluster.covmat
cm = CovMat(dim=self._dim, band=self._dim - 1)
cm.var = self.var
cov = self.cov
if self._dim == 2:
cm.set_cov(0, 1, cov)
else:
cm.set_cov(0, 1, cov[0])
cm.set_cov(0, 2, cov[1])
cm.set_cov(1, 2, cov[2])
return cm
