research/OpenFOAM-2.4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

polyMesh::findCell: Add support for octree search and cell tet-decomposition

Browse source 
Make this new method the default so that probes find the locations more reliably
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1544
This commit is contained in:
Henry 2015-02-24 16:52:51 +00:00
parent ab7c7a7d1c
commit 3321553039
2 changed files with 17 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

13
src/OpenFOAM/meshes/polyMesh/polyMesh.C
View file

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -1456,6 +1456,17 @@ bool Foam::polyMesh::pointInCell
return true;
}
break;
case CELL_TETS:
{
label tetFacei;
label tetPti;
findTetFacePt(cellI, p, tetFacei, tetPti);
return tetFacei != -1;
}
break;
}
return false;
}

9
src/OpenFOAM/meshes/polyMesh/polyMesh.H
View file

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -100,7 +100,8 @@ public:
{
FACEPLANES, // cell bound by planes of faces
FACECENTRETETS, // tet decomposition using facectr and cellctr
FACEDIAGTETS // tet decomposition using face diagonal and cellctr
FACEDIAGTETS, // tet decomposition using face diagonal and cellctr
CELL_TETS // Cell decomposed into tets
};
@ -693,14 +694,14 @@ public:
(
const point&,
label cellI,
const cellRepresentation = FACEDIAGTETS
const cellRepresentation = CELL_TETS
) const;
//- Find cell enclosing this location (-1 if not in mesh)
label findCell
(
const point&,
const cellRepresentation = FACEDIAGTETS
const cellRepresentation = CELL_TETS
) const;
};