34 Foam::refinementParameters::refinementParameters
40 maxGlobalCells_(readLabel(dict.
lookup(
"cellLimit"))),
41 maxLocalCells_(readLabel(dict.
lookup(
"procCellLimit"))),
42 minRefineCells_(readLabel(dict.
lookup(
"minimumRefine"))),
44 nBufferLayers_(readLabel(dict.
lookup(
"nBufferLayers"))),
45 keepPoints_(dict.
lookup(
"keepPoints")),
46 allowFreeStandingZoneFaces_
50 "allowFreeStandingZoneFaces",
58 Foam::refinementParameters::refinementParameters(
const dictionary& dict)
60 maxGlobalCells_(readLabel(dict.lookup(
"maxGlobalCells"))),
61 maxLocalCells_(readLabel(dict.lookup(
"maxLocalCells"))),
62 minRefineCells_(readLabel(dict.lookup(
"minRefinementCells"))),
63 nBufferLayers_(readLabel(dict.lookup(
"nCellsBetweenLevels"))),
64 keepPoints_(
pointField(1, dict.lookup(
"locationInMesh"))),
65 allowFreeStandingZoneFaces_
67 dict.lookupOrDefault<
Switch>
69 "allowFreeStandingZoneFaces",
73 maxLoadUnbalance_(dict.lookupOrDefault<scalar>(
"maxLoadUnbalance",0))
77 if (featAngle < 0 || featAngle > 180)
101 const point& keepPoint = keepPoints_[i];
103 label localCellI = mesh.
findCell(keepPoint);
105 label globalCellI = -1;
107 if (localCellI != -1)
109 Pout<<
"Found point " << keepPoint <<
" in cell " << localCellI
111 globalCellI = globalCells.toGlobal(localCellI);
116 if (globalCellI == -1)
120 "refinementParameters::findCells(const polyMesh&) const"
121 ) <<
"Point " << keepPoint
122 <<
" is not inside the mesh or on a face or edge." <<
nl
123 <<
"Bounding box of the mesh:" << mesh.
bounds()
127 if (globalCells.isLocal(globalCellI))
129 cellLabels[i] = localCellI;