songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

BUG: GAMG: avoid inter-level interpolation. Fixes #2710

Browse Source
This commit is contained in:
mattijs 2023-02-23 13:22:19 +00:00
parent f903707e4a
commit b339be2459
1 changed files with 11 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/GAMGSolverSolve.C
View File

@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011-2017 OpenFOAM Foundation
Copyright (C) 2016-2021 OpenCFD Ltd.
Copyright (C) 2016-2021,2023 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -347,8 +347,17 @@ void Foam::GAMGSolver::Vcycle
if (interpolateCorrection_) //&& leveli < coarsestLevel - 2)
{
if (coarseCorrFields.set(leveli+1))
if
(
coarseCorrFields.set(leveli+1)
&& (
matrixLevels_[leveli].mesh().comm()
== matrixLevels_[leveli+1].mesh().comm()
)
)
{
// Normal operation : have both coarse level and fine
// level. No processor agglomeration
interpolate
(
coarseCorrFields[leveli],