applications/solvers/solidMechanics/elasticAcpSolidFoam/aitkenRelaxation.H

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