applications/solvers/solidMechanics/elasticSolidFoam/aitkenRelaxation.H

   1 // aitken acceleration
   2
   3 aitkenDelta.storePrevIter();
   4
   5 // update delta
   6 aitkenDelta = (U - U.prevIter()) / aitkenInitialRes;
   7
   8 // update relaxation factor
   9 if(iCorr == 0)
  10 {
  11     aitkenTheta = 0.1;
  12 }
  13 else
  14 {
  15     vectorField b = aitkenDelta.internalField()
  16         - aitkenDelta.prevIter().internalField();
  17
  18     scalar sumMagB = gSum(magSqr(b));
  19     if(sumMagB < SMALL)
  20     {
  21         // Warning << "Aitken under-relaxation: denominator less then SMALL"
  22         // << endl;
  23         sumMagB += SMALL;
  24     }
  25
  26     aitkenTheta = -aitkenTheta*
  27         gSum(aitkenDelta.prevIter().internalField() & b)/
  28         sumMagB;
  29 }
  30
  31 // correction to the latest U
  32 U += aitkenTheta*aitkenDelta*aitkenInitialRes;