81 sumOff[l[face]] +=
mag(Upper[face]);
93 <<
"lduMatrix::operator=(const lduMatrix&) : "
94 <<
"attempted assignment to self"
153 upper() += A.
upper();
166 upper() += A.
upper();
167 lower() += A.
lower();
173 lower() += A.
upper();
174 upper() += A.
upper();
178 lower() += A.
lower();
179 upper() += A.
lower();
185 lower() += A.
lower();
186 upper() += A.
upper();
205 FatalErrorIn(
"lduMatrix::operator+=(const lduMatrix& A)")
206 <<
"Unknown matrix type combination"
221 upper() -= A.
upper();
234 upper() -= A.
upper();
235 lower() -= A.
lower();
241 lower() -= A.
upper();
242 upper() -= A.
upper();
246 lower() -= A.
lower();
247 upper() -= A.
lower();
253 lower() -= A.
lower();
254 upper() -= A.
upper();
260 upper() = -A.
upper();
265 lower() = -A.
lower();
273 FatalErrorIn(
"lduMatrix::operator-=(const lduMatrix& A)")
274 <<
"Unknown matrix type combination"
339 if (lowerPtr_ || upperPtr_)
343 scalar* __restrict__ H1Ptr = H1_.
begin();
345 const label* __restrict__ uPtr = lduAddr().upperAddr().begin();
346 const label* __restrict__ lPtr = lduAddr().lowerAddr().begin();
348 const scalar* __restrict__ lowerPtr = lower().begin();
349 const scalar* __restrict__ upperPtr = upper().begin();
351 register const label nFaces = upper().size();
355 H1Ptr[uPtr[
face]] -= lowerPtr[
face];
356 H1Ptr[lPtr[
face]] -= upperPtr[
face];