/** Check the filling of EigenMatrix per rows of blocks. Return true if the test succeeds.*/
bool checkEigenMatrixBlockRowFilling()
{
EigenBlockSparseMatrix mb;
FullMatrix ma;
unsigned br=3, bc=3;
ma.resize(br*BROWS,bc*BCOLS);
mb.resizeBlocks(br,bc);
for( unsigned i=0; i<br; i++ )
{
mb.beginBlockRow(i);
if( i%2==0 ) // leave some rows empty
{
for( int j=bc-1; j>=0; j--) // set the blocs in reverse order, for fun.
{
// create a block and give it some value
BlockMN b;
for( unsigned k=0; k<BROWS && k<BCOLS; k++ ){
b[k][k] = i+j;
ma.set(i*BROWS+k, j*BCOLS+k, i+j);
}
// insert the block in the matrix
mb.createBlock(j,b);
}
}
mb.endBlockRow();
}
mb.compress();
// serr()<<"MatrixTest<Real,RN,CN>::checkEigenMatrixBlockRowFilling, ma = " << ma << endl;
// serr()<<"MatrixTest<Real,RN,CN>::checkEigenMatrixBlockRowFilling, mb = " << mb << endl;
return Sofa_test<_Real>::matricesAreEqual(ma,mb);
}
