exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Fix handling for solo rows

Browse source
This commit is contained in:
Hrvoje Jasak 2017-08-01 23:15:57 +01:00
parent cc34210c87
commit 0326f7259f
1 changed files with 11 additions and 13 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

24
src/foam/matrices/blockLduMatrix/BlockAMG/BlockMatrixCoarsening/BlockMatrixSelection/BlockMatrixSelection.C
View file

@ -355,19 +355,17 @@ void Foam::BlockMatrixSelection<Type>::calcCoarsening()
equationWeight.set(i, tRow[i + 1] - tRow[i]);
}
// // HJ, REMOVE THIS: it is not allowed to have empty prolongation rows
// // because it destroys the diagonal dominance on the coarse level
// // HJ, 29/Jul/2017
// // Mark highly diagonally dominant rows as fine. This also removes solo
// // equations from coarsening
// for (label i = 0; i < nRows; i++)
// {
// // Label rows without strong connections as FINE
// if (strongRow[i + 1] == strongRow[i])
// {
// rowLabel_[i] = FINE;
// }
// }
// Mark disconnected rows as fine. This also removes solo
// equations from coarsening
// HJ, 29/Jul/2017
for (label i = 0; i < nRows; i++)
{
// Label rows without strong connections as FINE
if (strongRow[i + 1] == strongRow[i])
{
rowLabel_[i] = FINE;
}
}
// Start counting coarse equations
nCoarseEqns_ = 0;