Skip to content
Snippets Groups Projects
Commit 2c75608a authored by Rebecca Neumann's avatar Rebecca Neumann
Browse files

use correct initialization for BlockSizes > 1

[[Imported from SVN: r1429]]
parent 9f38781f
No related branches found
No related tags found
No related merge requests found
......@@ -638,7 +638,6 @@ namespace Dune
for(RowIterator row = fine.begin(); row != endRow; ++row)
if(aggregates[row.index()] != AggregatesMap<V>::ISOLATED) {
assert(aggregates[row.index()]!=AggregatesMap<V>::UNAGGREGATED);
//typedef typename RowIterator::Iterator ColIterator;
typedef typename M::ConstColIterator ColIterator;
ColIterator endCol = row->end();
......@@ -650,8 +649,8 @@ namespace Dune
}
// get the right diagonal matrix values on copy lines from owner processes
std::vector<typename M::block_type> rowsize(coarse.N(),0);
typedef typename M::block_type BlockType;
std::vector<BlockType> rowsize(coarse.N(),BlockType(0));
for (RowIterator row = coarse.begin(); row != coarse.end(); ++row)
rowsize[row.index()]=coarse[row.index()][row.index()];
pinfo.copyOwnerToAll(rowsize,rowsize);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment