38 label n = tmpMatrix.
n();
41 for (
register label i=0; i<n; i++)
44 scalar largestCoeff =
mag(tmpMatrix[iMax][i]);
47 for (
register label j=i+1; j<n; j++)
49 if (
mag(tmpMatrix[j][i]) > largestCoeff)
52 largestCoeff =
mag(tmpMatrix[iMax][i]);
60 for (
register label
k=i;
k<n;
k++)
62 Swap(tmpMatrix[i][
k], tmpMatrix[iMax][k]);
64 Swap(sourceSol[i], sourceSol[iMax]);
68 if (
mag(tmpMatrix[i][i]) < 1
e-20)
70 FatalErrorIn(
"solve(scalarSquareMatrix&, Field<Type>& sourceSol)")
76 for (
register label j=i+1; j<n; j++)
78 sourceSol[j] -= sourceSol[i]*(tmpMatrix[j][i]/tmpMatrix[i][i]);
80 for (
register label
k=n-1;
k>=i;
k--)
83 tmpMatrix[i][
k]*tmpMatrix[j][i]/tmpMatrix[i][i];
89 for (
register label j=n-1; j>=0; j--)
93 for (
register label
k=j+1;
k<n;
k++)
95 ntempvec += tmpMatrix[j][
k]*sourceSol[
k];
98 sourceSol[j] = (sourceSol[j] - ntempvec)/tmpMatrix[j][j];
113 solve(tmpMatrix, psi);
125 label n = luMatrix.
n();
129 for (
register label i=0; i<n; i++)
131 label ip = pivotIndices[i];
132 Type
sum = sourceSol[ip];
133 sourceSol[ip] = sourceSol[i];
134 const scalar* __restrict__ luMatrixi = luMatrix[i];
138 for (label j=ii-1; j<i; j++)
140 sum -= luMatrixi[j]*sourceSol[j];
151 for (
register label i=n-1; i>=0; i--)
153 Type
sum = sourceSol[i];
154 const scalar* __restrict__ luMatrixi = luMatrix[i];
156 for (
register label j=i+1; j<n; j++)
158 sum -= luMatrixi[j]*sourceSol[j];
161 sourceSol[i] = sum/luMatrixi[i];