nojhan/eodev
Archived
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

correct vector initialization within Eigen implemetation; delete useless asserts

Browse source
This commit is contained in:
nojhan 2012-07-10 14:09:10 +02:00
commit b0cbdf41ba
1 changed files with 6 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

9
edo/src/edoEstimatorNormalMulti.h
View file

@ -194,16 +194,19 @@ public:
// Be sure that the symmetric matrix got the good size // Be sure that the symmetric matrix got the good size
assert(var.innerSize() == s_size); assert(var.innerSize() == s_size);
assert(var.outerSize() == s_size); assert(var.outerSize() == s_size);
assert(var.innerSize() == _varcovar.innerSize());
assert(var.outerSize() == _varcovar.outerSize());
_varcovar = var / p_size; _varcovar = var / p_size;
// assert(var.innerSize() == _varcovar.innerSize());
// assert(var.outerSize() == _varcovar.outerSize());
// _mean.resize(s_size); // FIXME: check if it is really used because of the assignation below // _mean.resize(s_size); // FIXME: check if it is really used because of the assignation below
// unit vector // unit vector
// ublas::scalar_vector< AtomType > u( p_size, 1 ); // ublas::scalar_vector< AtomType > u( p_size, 1 );
Vector u( p_size, 1); Vector u( p_size);
u = Vector::Constant(p_size, 1);
// sum over columns // sum over columns
// _mean = ublas::prod( ublas::trans( sample ), u ); // _mean = ublas::prod( ublas::trans( sample ), u );