research/OpenFOAM-5.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

scalarMatrices: update [i][j] -> (i, j)

Browse source
This commit is contained in:
Henry Weller 2016-07-18 13:12:26 +01:00
parent 3420fee771
commit af6b67c719
1 changed files with 3 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
src/OpenFOAM/matrices/scalarMatrices/scalarMatricesTemplates.C
View file

@ -42,7 +42,7 @@ void Foam::solve
for (label i=0; i<m; i++) for (label i=0; i<m; i++)
{ {
label iMax = i; label iMax = i;
scalar largestCoeff = mag(tmpMatrix[iMax][i]); scalar largestCoeff = mag(tmpMatrix(iMax, i));
// Swap elements around to find a good pivot // Swap elements around to find a good pivot
for (label j=i+1; j<m; j++) for (label j=i+1; j<m; j++)
@ -50,7 +50,7 @@ void Foam::solve
if (mag(tmpMatrix(j, i)) > largestCoeff) if (mag(tmpMatrix(j, i)) > largestCoeff)
{ {
iMax = j; iMax = j;
largestCoeff = mag(tmpMatrix[iMax][i]); largestCoeff = mag(tmpMatrix(iMax, i));
} }
} }
@ -58,7 +58,7 @@ void Foam::solve
{ {
for (label k=i; k<m; k++) for (label k=i; k<m; k++)
{ {
Swap(tmpMatrix(i, k), tmpMatrix[iMax][k]); Swap(tmpMatrix(i, k), tmpMatrix(iMax, k));
} }
Swap(sourceSol[i], sourceSol[iMax]); Swap(sourceSol[i], sourceSol[iMax]);
} }