applications/solvers/solidMechanics/elasticOrthoAcpSolidFoam/setReference.H

   1 {
   2     // Set reference
   3     for (label regI = 0; regI < mesh.regions().nRegions(); regI++)
   4     {
   5         if (needRef[regI])
   6         {
   7             UEqn.source()[refCells[regI]] +=
   8                 UEqn.diag()[refCells[regI]]*vector::zero;
   9
  10             UEqn.diag()[refCells[regI]] +=
  11                 UEqn.diag()[refCells[regI]];
  12         }
  13     }
  14
  15     // Set value for dead cells (small particles or regions)
  16 //     label size = mesh.nCells() - label(sum(gamma.internalField()));
  17
  18     if (deadCells.size())
  19     {
  20         labelList curRefCell(deadCells.size(), -1);
  21         label index = 0;
  22         forAll(deadCells, cellI)
  23         {
  24             curRefCell[index++] = deadCells[cellI];
  25
  26 //             if (gamma.internalField()[cellI] < 1 - SMALL)
  27 //             {
  28 //                 curRefCell[index++] = cellI;
  29 //             }
  30         }
  31
  32         vectorField curRefValue(curRefCell.size(), vector::zero);
  33
  34         UEqn.setValues(curRefCell, curRefValue);
  35     }
  36 }