zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Tidy-up and renaming of long variable names.

Browse Source
This commit is contained in:
graham
2008-09-15 12:16:26 +01:00
parent 735cc6bb1b
commit 3ccfd03125
1 changed files with 83 additions and 255 deletions
Download Patch File Download Diff File Expand all files Collapse all files

338
src/lagrangian/molecularDynamics/molecule/interactionLists/referredCellList/referredCellList.C
View File

@ -46,14 +46,17 @@ void Foam::referredCellList::buildReferredCellList
DynamicList<referredCell> referredInteractionList; DynamicList<referredCell> referredInteractionList;
DynamicList<label> realCellsWithinRCutMaxOfAnyReferringPatch; // realCellsWithinRCutMaxOfAnyReferringPatch
DynamicList<label> rCellsWRRP;
DynamicList<label> realFacesWithinRCutMaxOfAnyReferringPatch; // realFacesWithinRCutMaxOfAnyReferringPatch
DynamicList<label> rFacesWRRP;
DynamicList<label> realEdgesWithinRCutMaxOfAnyReferringPatch; // realEdgesWithinRCutMaxOfAnyReferringPatch
DynamicList<label> rEdgesWRRP;
DynamicList<label> realPointsWithinRCutMaxOfAnyReferringPatch;
// realPointsWithinRCutMaxOfAnyReferringPatch
DynamicList<label> rPointsWRRP;
labelListList processorPatchSegmentMapping labelListList processorPatchSegmentMapping
( (
@ -78,17 +81,6 @@ void Foam::referredCellList::buildReferredCellList
DynamicList<word> patchNames; DynamicList<word> patchNames;
// IOobject undecomposedBoundaryHeader
// (
// "undecomposedBoundary",
// mesh.time().constant(),
// polyMesh::meshSubDir,
// mesh,
// IOobject::MUST_READ,
// IOobject::NO_WRITE,
// false
// );
Time undecomposedTime Time undecomposedTime
( (
Time::controlDictName, Time::controlDictName,
@ -129,12 +121,9 @@ void Foam::referredCellList::buildReferredCellList
} }
else else
{ {
FatalErrorIn FatalErrorIn ("referredCellList.C")
( << nl << "unable to read undecomposed boundary file from "
"moleculeCloudBuildCellInteractionLists.C" << "constant/polyMesh" << nl
)
<< nl << "unable to read undecomposed boundary file from "
"constant/polyMesh" << nl
<< abort(FatalError); << abort(FatalError);
} }
@ -168,8 +157,7 @@ void Foam::referredCellList::buildReferredCellList
forAll(procPatches,pP) forAll(procPatches,pP)
{ {
const processorPolyPatch& patch = const processorPolyPatch& patch = refCast<const processorPolyPatch>
refCast<const processorPolyPatch>
( (
mesh.boundaryMesh()[procPatches[pP]] mesh.boundaryMesh()[procPatches[pP]]
); );
@ -186,12 +174,6 @@ void Foam::referredCellList::buildReferredCellList
label globalFace = abs(faceProcAdd)-1; label globalFace = abs(faceProcAdd)-1;
// Pout << "Patch index: " << pI
// << " " << patch[pI]
// << " Mesh index: " << decomposedMeshFace
// << " faceProcAdd: " << faceProcAdd
// << " globalFace:" << globalFace << endl;
label minStart = -1; label minStart = -1;
forAll(patchNames, pN) forAll(patchNames, pN)
@ -201,7 +183,6 @@ void Foam::referredCellList::buildReferredCellList
const dictionary& patchDict = const dictionary& patchDict =
patchDictionary.subDict(patchNames[pN]); patchDictionary.subDict(patchNames[pN]);
// word faceName(patchNames[pN]);
label startFace(readLabel(patchDict.lookup("startFace"))); label startFace(readLabel(patchDict.lookup("startFace")));
label nFaces(readLabel(patchDict.lookup("nFaces"))); label nFaces(readLabel(patchDict.lookup("nFaces")));
@ -213,7 +194,7 @@ void Foam::referredCellList::buildReferredCellList
if if
( (
globalFace >= startFace globalFace >= startFace
&& globalFace < startFace + nFaces/2 && globalFace < startFace + nFaces/2
) )
{ {
procPatchSegMap[pI] = pN + 1; procPatchSegMap[pI] = pN + 1;
@ -221,7 +202,7 @@ void Foam::referredCellList::buildReferredCellList
else if else if
( (
globalFace >= startFace + nFaces/2 globalFace >= startFace + nFaces/2
&& globalFace < startFace + nFaces && globalFace < startFace + nFaces
) )
{ {
procPatchSegMap[pI] = -(pN + 1); procPatchSegMap[pI] = -(pN + 1);
@ -238,8 +219,7 @@ void Foam::referredCellList::buildReferredCellList
forAll(procPatches,pP) forAll(procPatches,pP)
{ {
const processorPolyPatch& patch = const processorPolyPatch& patch = refCast<const processorPolyPatch>
refCast<const processorPolyPatch>
( (
mesh.boundaryMesh()[procPatches[pP]] mesh.boundaryMesh()[procPatches[pP]]
); );
@ -252,20 +232,13 @@ void Foam::referredCellList::buildReferredCellList
+ 2*(sizeof(label) + patch.size()*sizeof(vector)) + 2*(sizeof(label) + patch.size()*sizeof(vector))
); );
toNeighbProc toNeighbProc << patch.faceCentres() << patch.faceAreas();
<< patch.faceCentres()
<< patch.faceAreas();
// Pout << "sending face info to processor "
// << patch.neighbProcNo() << ", size = "
// << patch.size() << endl;
} }
} }
forAll(procPatches,pP) forAll(procPatches,pP)
{ {
const processorPolyPatch& patch = const processorPolyPatch& patch = refCast<const processorPolyPatch>
refCast<const processorPolyPatch>
( (
mesh.boundaryMesh()[procPatches[pP]] mesh.boundaryMesh()[procPatches[pP]]
); );
@ -285,13 +258,8 @@ void Foam::referredCellList::buildReferredCellList
patch.neighbProcNo(), patch.neighbProcNo(),
2*(sizeof(label) + patch.size()*sizeof(vector)) 2*(sizeof(label) + patch.size()*sizeof(vector))
); );
fromNeighbProc
>> neighbFaceCentres
>> neighbFaceAreas;
// Pout << "receiving face info from processor " fromNeighbProc >> neighbFaceCentres >> neighbFaceAreas;
// << patch.neighbProcNo() << ", size = "
// << patch.size() << endl;
} }
} }
@ -303,8 +271,7 @@ void Foam::referredCellList::buildReferredCellList
forAll(procPatches,pP) forAll(procPatches,pP)
{ {
const processorPolyPatch& patch = const processorPolyPatch& patch = refCast<const processorPolyPatch>
refCast<const processorPolyPatch>
( (
mesh.boundaryMesh()[procPatches[pP]] mesh.boundaryMesh()[procPatches[pP]]
); );
@ -356,15 +323,14 @@ void Foam::referredCellList::buildReferredCellList
{ {
if if
( (
sum(mag(refOff-refOff[0]))/refOff.size() > interactionLists::transTol sum(mag(refOff-refOff[0]))/refOff.size()
|| sum(mag(refTrans-refTrans[0]))/refTrans.size() > interactionLists::transTol > interactionLists::transTol
|| sum(mag(refTrans-refTrans[0]))/refTrans.size()
> interactionLists::transTol
) )
{ {
FatalErrorIn FatalErrorIn ("referredCellList.C")
( << nl << "Face pairs on patch "
"moleculeCloud::buildCellInteractionLists"
)
<< nl << "Face pairs on patch "
<< patch.name() << patch.name()
<< ", segment " << patchNames[nUP] << ", segment " << patchNames[nUP]
<< " do not give the same referring " << " do not give the same referring "
@ -398,16 +364,13 @@ void Foam::referredCellList::buildReferredCellList
if (mesh.boundaryMesh()[patchI].type() == "cyclic") if (mesh.boundaryMesh()[patchI].type() == "cyclic")
{ {
const cyclicPolyPatch& patch = const cyclicPolyPatch& patch = refCast<const cyclicPolyPatch>
refCast<const cyclicPolyPatch>
( (
mesh.boundaryMesh()[patchI] mesh.boundaryMesh()[patchI]
); );
if(patch.size()) if (patch.size())
{ {
// Pout << patch << endl;
if (iterationNo == 0) if (iterationNo == 0)
{ {
// Tests that all combinations of face pairs produce the // Tests that all combinations of face pairs produce the
@ -450,16 +413,13 @@ void Foam::referredCellList::buildReferredCellList
if if
( (
sum(mag(refOff - refOff[0]))/(patch.size()/2) sum(mag(refOff - refOff[0]))/(patch.size()/2)
> interactionLists::transTol > interactionLists::transTol
|| sum(mag(refTrans - refTrans[0]))/(patch.size()/2) || sum(mag(refTrans - refTrans[0]))/(patch.size()/2)
> interactionLists::transTol > interactionLists::transTol
) )
{ {
FatalErrorIn FatalErrorIn ("referredCellList.C")
( << nl << "Face pairs on patch "
"moleculeCloud::buildCellInteractionLists"
)
<< nl << "Face pairs on patch "
<< patch.name() << patch.name()
<< " do not give the same referring " << " do not give the same referring "
<< " transformations to within tolerance of " << " transformations to within tolerance of "
@ -518,11 +478,7 @@ void Foam::referredCellList::buildReferredCellList
{ {
const label facePoint(facePoints[fP]); const label facePoint(facePoints[fP]);
if if (findIndex(meshPointsOnThisSegment, facePoint) == -1)
(
findIndex(meshPointsOnThisSegment, facePoint)
== -1
)
{ {
meshPointsOnThisSegment.append(facePoint); meshPointsOnThisSegment.append(facePoint);
} }
@ -574,13 +530,7 @@ void Foam::referredCellList::buildReferredCellList
forAll(referredInteractionList, rIL) forAll(referredInteractionList, rIL)
{ {
if if (cellToRefer.duplicate (referredInteractionList[rIL]))
(
cellToRefer.duplicate
(
referredInteractionList[rIL]
)
)
{ {
addCellToRefer = false; addCellToRefer = false;
@ -591,14 +541,7 @@ void Foam::referredCellList::buildReferredCellList
// Check for cellToRefer being referred back // Check for cellToRefer being referred back
// ontop of a real cell // ontop of a real cell
if if (cellToRefer.duplicate (Pstream::myProcNo(), mesh.nCells()))
(
cellToRefer.duplicate
(
Pstream::myProcNo(),
mesh.nCells()
)
)
{ {
addCellToRefer = false; addCellToRefer = false;
} }
@ -611,19 +554,9 @@ void Foam::referredCellList::buildReferredCellList
// add real cells found in range of cyclic patch // add real cells found in range of cyclic patch
// to whole mesh list // to whole mesh list
if if (findIndex (rCellsWRRP, realCell) == -1)
(
findIndex
(
realCellsWithinRCutMaxOfAnyReferringPatch,
realCell
) == -1
)
{ {
realCellsWithinRCutMaxOfAnyReferringPatch.append rCellsWRRP.append(realCell);
(
realCell
);
} }
} }
} }
@ -668,10 +601,7 @@ void Foam::referredCellList::buildReferredCellList
{ {
if if
( (
cellToReRefer.duplicate cellToReRefer.duplicate(referredInteractionList[rIL])
(
referredInteractionList[rIL]
)
) )
{ {
addCellToReRefer = false; addCellToReRefer = false;
@ -685,17 +615,13 @@ void Foam::referredCellList::buildReferredCellList
if if
( (
cellToReRefer.duplicate cellToReRefer.duplicate(Pstream::myProcNo(), mesh.nCells())
(
Pstream::myProcNo(),
mesh.nCells()
)
) )
{ {
addCellToReRefer = false; addCellToReRefer = false;
} }
if(addCellToReRefer) if (addCellToReRefer)
{ {
referredInteractionList.append(cellToReRefer); referredInteractionList.append(cellToReRefer);
} }
@ -731,7 +657,7 @@ void Foam::referredCellList::buildReferredCellList
{ {
const label faceEdge(faceEdges[fE]); const label faceEdge(faceEdges[fE]);
if(findIndex(meshEdgesOnThisSegment, faceEdge) == -1) if (findIndex(meshEdgesOnThisSegment, faceEdge) == -1)
{ {
meshEdgesOnThisSegment.append(faceEdge); meshEdgesOnThisSegment.append(faceEdge);
} }
@ -743,11 +669,7 @@ void Foam::referredCellList::buildReferredCellList
{ {
const label facePoint(facePoints[fP]); const label facePoint(facePoints[fP]);
if if (findIndex(meshPointsOnThisSegment, facePoint) == -1)
(
findIndex(meshPointsOnThisSegment, facePoint)
== -1
)
{ {
meshPointsOnThisSegment.append(facePoint); meshPointsOnThisSegment.append(facePoint);
} }
@ -799,13 +721,7 @@ void Foam::referredCellList::buildReferredCellList
forAll(referredInteractionList, rIL) forAll(referredInteractionList, rIL)
{ {
if if (cellToRefer.duplicate(referredInteractionList[rIL]))
(
cellToRefer.duplicate
(
referredInteractionList[rIL]
)
)
{ {
addCellToRefer = false; addCellToRefer = false;
@ -818,11 +734,7 @@ void Foam::referredCellList::buildReferredCellList
if if
( (
cellToRefer.duplicate cellToRefer.duplicate(Pstream::myProcNo(), mesh.nCells())
(
Pstream::myProcNo(),
mesh.nCells()
)
) )
{ {
addCellToRefer = false; addCellToRefer = false;
@ -836,19 +748,9 @@ void Foam::referredCellList::buildReferredCellList
// add real cells found in range of cyclic patch // add real cells found in range of cyclic patch
// to whole mesh list // to whole mesh list
if if (findIndex (rCellsWRRP, realCell) == -1)
(
findIndex
(
realCellsWithinRCutMaxOfAnyReferringPatch,
realCell
) == -1
)
{ {
realCellsWithinRCutMaxOfAnyReferringPatch.append rCellsWRRP.append(realCell);
(
realCell
);
} }
} }
} }
@ -888,13 +790,7 @@ void Foam::referredCellList::buildReferredCellList
forAll(referredInteractionList, rIL) forAll(referredInteractionList, rIL)
{ {
if if (cellToReRefer.duplicate (referredInteractionList[rIL]))
(
cellToReRefer.duplicate
(
referredInteractionList[rIL]
)
)
{ {
addCellToReRefer = false; addCellToReRefer = false;
@ -905,14 +801,7 @@ void Foam::referredCellList::buildReferredCellList
// Check for cellToRefer being referred back // Check for cellToRefer being referred back
// ontop of a real cell // ontop of a real cell
if if (cellToReRefer.duplicate(Pstream::myProcNo(), mesh.nCells()))
(
cellToReRefer.duplicate
(
Pstream::myProcNo(),
mesh.nCells()
)
)
{ {
addCellToReRefer = false; addCellToReRefer = false;
} }
@ -932,8 +821,7 @@ void Foam::referredCellList::buildReferredCellList
forAll(procPatches,pP) forAll(procPatches,pP)
{ {
const processorPolyPatch& patch = const processorPolyPatch& patch = refCast<const processorPolyPatch>
refCast<const processorPolyPatch>
( (
mesh.boundaryMesh()[procPatches[pP]] mesh.boundaryMesh()[procPatches[pP]]
); );
@ -992,7 +880,7 @@ void Foam::referredCellList::buildReferredCellList
{ {
const label faceEdge(faceEdges[fE]); const label faceEdge(faceEdges[fE]);
if(findIndex(meshEdgesOnThisSegment, faceEdge) == -1) if (findIndex(meshEdgesOnThisSegment, faceEdge) == -1)
{ {
meshEdgesOnThisSegment.append(faceEdge); meshEdgesOnThisSegment.append(faceEdge);
} }
@ -1004,11 +892,7 @@ void Foam::referredCellList::buildReferredCellList
{ {
const label facePoint(facePoints[fP]); const label facePoint(facePoints[fP]);
if if (findIndex(meshPointsOnThisSegment, facePoint) == -1)
(
findIndex(meshPointsOnThisSegment, facePoint)
== -1
)
{ {
meshPointsOnThisSegment.append(facePoint); meshPointsOnThisSegment.append(facePoint);
} }
@ -1023,11 +907,8 @@ void Foam::referredCellList::buildReferredCellList
{ {
if (faceT == -1) if (faceT == -1)
{ {
FatalErrorIn FatalErrorIn ("referredCellList.C")
( << nl << "faceT == -1 encountered but "
"moleculeCloudBuildReferredInteractionList.H"
)
<< nl << "faceT == -1 encountered but "
<< meshFacesOnThisSegment.size() << meshFacesOnThisSegment.size()
<< " faces found on patch segment." << " faces found on patch segment."
<< abort(FatalError); << abort(FatalError);
@ -1070,20 +951,9 @@ void Foam::referredCellList::buildReferredCellList
// add real cells found in range of processor patch // add real cells found in range of processor patch
// to whole mesh list // to whole mesh list
if if (findIndex (rCellsWRRP, realCell) == -1)
(
findIndex
(
realCellsWithinRCutMaxOfAnyReferringPatch,
realCell
)
== -1
)
{ {
realCellsWithinRCutMaxOfAnyReferringPatch.append rCellsWRRP.append(realCell);
(
realCell
);
} }
} }
} }
@ -1104,7 +974,7 @@ void Foam::referredCellList::buildReferredCellList
forAll(referredCellsFoundInRange,cFIR) forAll(referredCellsFoundInRange,cFIR)
{ {
referredCell& existingRefCell referredCell& existingRefCell
= referredInteractionList[referredCellsFoundInRange[cFIR]]; = referredInteractionList[referredCellsFoundInRange[cFIR]];
referredCell cellToReRefer = existingRefCell.reRefer referredCell cellToReRefer = existingRefCell.reRefer
( (
@ -1131,23 +1001,17 @@ void Foam::referredCellList::buildReferredCellList
patch.neighbProcNo() patch.neighbProcNo()
); );
// Pout << "Processor " << Pstream::myProcNo()
// << " sending to " << patch.neighbProcNo()
// << nl << referredCellsToTransfer << endl;
toNeighbProc << referredCellsToTransfer; toNeighbProc << referredCellsToTransfer;
} }
} }
forAll(procPatches,pP) forAll(procPatches,pP)
{ {
const processorPolyPatch& patch = const processorPolyPatch& patch = refCast<const processorPolyPatch>
refCast<const processorPolyPatch>
( (
mesh.boundaryMesh()[procPatches[pP]] mesh.boundaryMesh()[procPatches[pP]]
); );
// Receive the cells from neighbour // Receive the cells from neighbour
List<referredCell> referredCellsFromNeighbour(patch.size()); List<referredCell> referredCellsFromNeighbour(patch.size());
@ -1160,10 +1024,6 @@ void Foam::referredCellList::buildReferredCellList
); );
fromNeighbProc >> referredCellsFromNeighbour; fromNeighbProc >> referredCellsFromNeighbour;
// Pout << "Processor " << Pstream::myProcNo()
// << " receiving from " << patch.neighbProcNo()
// << nl << referredCellsFromNeighbour << endl;
} }
// Check to see if they are duplicates, if not append // Check to see if they are duplicates, if not append
@ -1171,8 +1031,7 @@ void Foam::referredCellList::buildReferredCellList
forAll(referredCellsFromNeighbour,rCFN) forAll(referredCellsFromNeighbour,rCFN)
{ {
referredCell& cellToRefer referredCell& cellToRefer = referredCellsFromNeighbour[rCFN];
= referredCellsFromNeighbour[rCFN];
// Test all existing referred and real cells to check duplicates // Test all existing referred and real cells to check duplicates
// are not being made or cells aren't being referred back onto // are not being made or cells aren't being referred back onto
@ -1184,7 +1043,7 @@ void Foam::referredCellList::buildReferredCellList
forAll(referredInteractionList, rIL) forAll(referredInteractionList, rIL)
{ {
if(cellToRefer.duplicate(referredInteractionList[rIL])) if (cellToRefer.duplicate(referredInteractionList[rIL]))
{ {
addCellToRefer = false; addCellToRefer = false;
@ -1196,11 +1055,7 @@ void Foam::referredCellList::buildReferredCellList
if if
( (
cellToRefer.duplicate cellToRefer.duplicate(Pstream::myProcNo(), mesh.nCells())
(
Pstream::myProcNo(),
mesh.nCells()
)
) )
{ {
addCellToRefer = false; addCellToRefer = false;
@ -1219,16 +1074,13 @@ void Foam::referredCellList::buildReferredCellList
// record all real cells in range of any referring patch (cyclic or // record all real cells in range of any referring patch (cyclic or
// processor) on the first iteration when the real cells are evaluated. // processor) on the first iteration when the real cells are evaluated.
realCellsWithinRCutMaxOfAnyReferringPatch.shrink(); rCellsWRRP.shrink();
// construct {faces, edges, points}WithinRCutMaxOfAnyReferringPatch // construct {faces, edges, points}WithinRCutMaxOfAnyReferringPatch
forAll(realCellsWithinRCutMaxOfAnyReferringPatch, rCWR) forAll(rCellsWRRP, rCWR)
{ {
const label realCell const label realCell(rCellsWRRP[rCWR]);
(
realCellsWithinRCutMaxOfAnyReferringPatch[rCWR]
);
const labelList& rCFaces const labelList& rCFaces
( (
@ -1239,13 +1091,9 @@ void Foam::referredCellList::buildReferredCellList
{ {
const label f(rCFaces[rCF]); const label f(rCFaces[rCF]);
if if (findIndex(rFacesWRRP,f) == -1)
(
findIndex(realFacesWithinRCutMaxOfAnyReferringPatch,f)
== -1
)
{ {
realFacesWithinRCutMaxOfAnyReferringPatch.append(f); rFacesWRRP.append(f);
} }
} }
@ -1258,13 +1106,9 @@ void Foam::referredCellList::buildReferredCellList
{ {
const label e(rCEdges[rCE]); const label e(rCEdges[rCE]);
if if (findIndex(rEdgesWRRP,e) == -1)
(
findIndex(realEdgesWithinRCutMaxOfAnyReferringPatch,e)
== -1
)
{ {
realEdgesWithinRCutMaxOfAnyReferringPatch.append(e); rEdgesWRRP.append(e);
} }
} }
@ -1277,46 +1121,39 @@ void Foam::referredCellList::buildReferredCellList
{ {
const label p(rCPoints[rCP]); const label p(rCPoints[rCP]);
if if (findIndex(rPointsWRRP,p) == -1)
(
findIndex(realPointsWithinRCutMaxOfAnyReferringPatch,p)
== -1
)
{ {
realPointsWithinRCutMaxOfAnyReferringPatch.append(p); rPointsWRRP.append(p);
} }
} }
} }
realFacesWithinRCutMaxOfAnyReferringPatch.shrink(); rFacesWRRP.shrink();
realEdgesWithinRCutMaxOfAnyReferringPatch.shrink(); rEdgesWRRP.shrink();
realPointsWithinRCutMaxOfAnyReferringPatch.shrink(); rPointsWRRP.shrink();
} }
iterationNo++; iterationNo++;
cellsReferredThisIteration = cellsReferredThisIteration =
referredInteractionList.size() - refIntListStartSize; referredInteractionList.size() - refIntListStartSize;
// Pout << tab << "Cells added this iteration: "
// << cellsReferredThisIteration << endl;
reduce(cellsReferredThisIteration, sumOp<label>()); reduce(cellsReferredThisIteration, sumOp<label>());
Info << tab << "Cells added this iteration: " Info<< tab << "Cells added this iteration: "
<< cellsReferredThisIteration << endl; << cellsReferredThisIteration << endl;
} }
referredInteractionList.shrink(); referredInteractionList.shrink();
// Info << "referredInteractionList.size() = " // Info<< "referredInteractionList.size() = "
// << referredInteractionList.size() << endl; // << referredInteractionList.size() << endl;
// forAll(referredInteractionList,rIL) // forAll(referredInteractionList,rIL)
// { // {
// Info << referredInteractionList[rIL]; // Info<< referredInteractionList[rIL];
// } // }
(*this).setSize (*this).setSize
@ -1341,12 +1178,9 @@ void Foam::referredCellList::buildReferredCellList
const vectorList& refCellPoints = refCell.vertexPositions(); const vectorList& refCellPoints = refCell.vertexPositions();
forAll(realFacesWithinRCutMaxOfAnyReferringPatch, rCF) forAll(rFacesWRRP, rCF)
{ {
const label f const label f(rFacesWRRP[rCF]);
(
realFacesWithinRCutMaxOfAnyReferringPatch[rCF]
);
if (il_.testPointFaceDistance(refCellPoints,f)) if (il_.testPointFaceDistance(refCellPoints,f))
{ {
@ -1371,12 +1205,9 @@ void Foam::referredCellList::buildReferredCellList
} }
} }
forAll(realPointsWithinRCutMaxOfAnyReferringPatch, rCP) forAll(rPointsWRRP, rCP)
{ {
const label p const label p(rPointsWRRP[rCP]);
(
realPointsWithinRCutMaxOfAnyReferringPatch[rCP]
);
if (il_.testPointFaceDistance(p,refCell)) if (il_.testPointFaceDistance(p,refCell))
{ {
@ -1397,12 +1228,9 @@ void Foam::referredCellList::buildReferredCellList
const edgeList& refCellEdges = refCell.edges(); const edgeList& refCellEdges = refCell.edges();
forAll(realEdgesWithinRCutMaxOfAnyReferringPatch, rCE) forAll(rEdgesWRRP, rCE)
{ {
const label edgeIIndex const label edgeIIndex(rEdgesWRRP[rCE]);
(
realEdgesWithinRCutMaxOfAnyReferringPatch[rCE]
);
const edge& eI(mesh.edges()[edgeIIndex]); const edge& eI(mesh.edges()[edgeIIndex]);
@ -1457,7 +1285,7 @@ void Foam::referredCellList::buildReferredCellList
// { // {
// const label cellI(ptICells[pIC]); // const label cellI(ptICells[pIC]);
// //
// if(findIndex(realCellsFoundInRange, cellI) == -1) // if (findIndex(realCellsFoundInRange, cellI) == -1)
// { // {
// realCellsFoundInRange.append(cellI); // realCellsFoundInRange.append(cellI);
// } // }