34 lduMatrix::solver::addsymMatrixConstructorToTable<PCG>
43 const word& fieldName,
75 lduMatrix::preconditioner::getName(controlDict_) + typeName,
79 register label nCells = psi.
size();
81 scalar* __restrict__ psiPtr = psi.
begin();
84 scalar* __restrict__ pAPtr = pA.
begin();
87 scalar* __restrict__ wAPtr = wA.
begin();
89 scalar wArA = matrix_.great_;
90 scalar wArAold = wArA;
93 matrix_.Amul(wA, psi, interfaceBouCoeffs_, interfaces_, cmpt);
97 scalar* __restrict__ rAPtr = rA.
begin();
100 scalar normFactor = this->normFactor(psi, source, wA, pA);
102 if (lduMatrix::debug >= 2)
104 Info<<
" Normalisation factor = " << normFactor <<
endl;
116 lduMatrix::preconditioner::New
129 preconPtr->precondition(wA, rA, cmpt);
143 scalar
beta = wArA/wArAold;
153 matrix_.Amul(wA, pA, interfaceBouCoeffs_, interfaces_, cmpt);
164 scalar alpha = wArA/wApA;