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

Bugfix: sortedToc in addressing cannot be a reference. HJ

Browse source
This commit is contained in:
Hrvoje Jasak 2017-08-24 17:25:44 +01:00
parent 2f92fb6c8a
commit f646502011
1 changed files with 12 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

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

@ -1043,7 +1043,7 @@ Foam::BlockMatrixSelection<Type>::restrictMatrix() const
forAll (coarseNbrsSets, rowI) forAll (coarseNbrsSets, rowI)
{ {
const labelList& curNbrs = coarseNbrsSets[rowI].sortedToc(); const labelList curNbrs = coarseNbrsSets[rowI].sortedToc();
forAll (curNbrs, nbrI) forAll (curNbrs, nbrI)
{ {
@ -1262,6 +1262,17 @@ Foam::BlockMatrixSelection<Type>::restrictMatrix() const
{ {
coeffLabel[lowerCoarseAddr[losortCoarseAddr[indexC]]] = coeffLabel[lowerCoarseAddr[losortCoarseAddr[indexC]]] =
losortCoarseAddr[indexC]; losortCoarseAddr[indexC];
// Stupid check
if (upperCoarseAddr[losortCoarseAddr[indexC]] != ir)
{
FatalErrorIn("STUPID CHECK")
<< "FAILED: ir, upper, lower coeff: "
<< ir << " "
<< upperCoarseAddr[losortCoarseAddr[indexC]] << " "
<< lowerCoarseAddr[losortCoarseAddr[indexC]] << " "
<< losortCoarseAddr[indexC]
<< abort(FatalError);
}
} }
// Compressed row format, get indices of coeffsR in row ir // Compressed row format, get indices of coeffsR in row ir