zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: [fa|fv]Matrix setValuesFromList - corrected loops. Fixes #2642

Browse Source
This commit is contained in:
Andrew Heather
2022-11-25 09:23:50 +00:00
parent adf95d483c
commit 360af221fe
2 changed files with 14 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
src/finiteArea/faMatrices/faMatrix/faMatrix.C
View File

@ -335,9 +335,15 @@ void Foam::faMatrix<Type>::setValuesFromList
psi[facei] = value; psi[facei] = value;
source_[facei] = value*Diag[facei]; source_[facei] = value*Diag[facei];
}
if (symmetric() || asymmetric()) if (symmetric() || asymmetric())
{
forAll(faceLabels, i)
{ {
const label facei = faceLabels[i];
const Type& value = values[i];
for (const label edgei : edges[facei]) for (const label edgei : edges[facei])
{ {
if (mesh.isInternalEdge(edgei)) if (mesh.isInternalEdge(edgei))

8
src/finiteVolume/fvMatrices/fvMatrix/fvMatrix.C
View File

@ -249,9 +249,15 @@ void Foam::fvMatrix<Type>::setValuesFromList
psi[celli] = value; psi[celli] = value;
source_[celli] = value*Diag[celli]; source_[celli] = value*Diag[celli];
}
if (symmetric() || asymmetric()) if (symmetric() || asymmetric())
{
forAll(cellLabels, i)
{ {
const label celli = cellLabels[i];
const Type& value = values[i];
for (const label facei : cells[celli]) for (const label facei : cells[celli])
{ {
if (mesh.isInternalFace(facei)) if (mesh.isInternalFace(facei))