Commit 0fcade3f authored by Peter Gottschling's avatar Peter Gottschling
Browse files

Used accumulated slot size.

parent d341b941
......@@ -72,7 +72,14 @@ namespace AMDiS {
matrix.change_dim(block_starts[ns], block_starts[ns]);
set_to_zero(matrix);
DOFMatrix::inserter_type ins(matrix);
int nnz= 0;
for (int rb= 0; rb < ns; ++rb)
for (int cb= 0; cb < ns; ++cb)
if (A[rb][cb])
nnz+= A[rb][cb]->getBaseMatrix().nnz();
DOFMatrix::inserter_type ins(matrix, int(1.2 * nnz / matrix.dim1()));
for (int rb= 0; rb < ns; ++rb)
for (int cb= 0; cb < ns; ++cb)
if (A[rb][cb])
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment