43 scalar* __restrict__ ApsiPtr = Apsi.
begin();
46 const scalar*
const __restrict__ psiPtr = psi.
begin();
48 const scalar*
const __restrict__ diagPtr =
diag().begin();
50 const label*
const __restrict__ uPtr = lduAddr().upperAddr().begin();
51 const label*
const __restrict__ lPtr = lduAddr().lowerAddr().begin();
53 const scalar*
const __restrict__ upperPtr = upper().begin();
54 const scalar*
const __restrict__ lowerPtr = lower().begin();
66 register const label nCells =
diag().size();
73 register const label nFaces = upper().size();
77 ApsiPtr[uPtr[
face]] += lowerPtr[
face]*psiPtr[lPtr[
face]];
78 ApsiPtr[lPtr[
face]] += upperPtr[
face]*psiPtr[uPtr[
face]];
82 updateMatrixInterfaces
104 scalar* __restrict__ TpsiPtr = Tpsi.
begin();
107 const scalar*
const __restrict__ psiPtr = psi.
begin();
109 const scalar*
const __restrict__ diagPtr =
diag().begin();
111 const label*
const __restrict__ uPtr = lduAddr().upperAddr().begin();
112 const label*
const __restrict__ lPtr = lduAddr().lowerAddr().begin();
114 const scalar*
const __restrict__ lowerPtr = lower().begin();
115 const scalar*
const __restrict__ upperPtr = upper().begin();
127 register const label nCells =
diag().size();
133 register const label nFaces = upper().size();
136 TpsiPtr[uPtr[
face]] += upperPtr[
face]*psiPtr[lPtr[
face]];
137 TpsiPtr[lPtr[
face]] += lowerPtr[
face]*psiPtr[uPtr[
face]];
141 updateMatrixInterfaces
161 scalar* __restrict__ sumAPtr = sumA.
begin();
163 const scalar* __restrict__ diagPtr =
diag().begin();
165 const label* __restrict__ uPtr = lduAddr().upperAddr().begin();
166 const label* __restrict__ lPtr = lduAddr().lowerAddr().begin();
168 const scalar* __restrict__ lowerPtr = lower().begin();
169 const scalar* __restrict__ upperPtr = upper().begin();
171 register const label nCells =
diag().size();
172 register const label nFaces = upper().size();
181 sumAPtr[uPtr[
face]] += lowerPtr[
face];
182 sumAPtr[lPtr[
face]] += upperPtr[
face];
187 forAll(interfaces, patchI)
189 if (interfaces.
set(patchI))
192 const scalarField& pCoeffs = interfaceBouCoeffs[patchI];
213 scalar* __restrict__ rAPtr = rA.
begin();
215 const scalar*
const __restrict__ psiPtr = psi.
begin();
216 const scalar*
const __restrict__ diagPtr =
diag().begin();
217 const scalar*
const __restrict__ sourcePtr = source.
begin();
219 const label*
const __restrict__ uPtr = lduAddr().upperAddr().begin();
220 const label*
const __restrict__ lPtr = lduAddr().lowerAddr().begin();
222 const scalar*
const __restrict__ upperPtr = upper().begin();
223 const scalar*
const __restrict__ lowerPtr = lower().begin();
256 register const label nCells =
diag().size();
263 register const label nFaces = upper().size();
272 updateMatrixInterfaces
293 residual(trA(), psi, source, interfaceBouCoeffs, interfaces, cmpt);