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

ENH: add findCell for PDRblockMesh

Browse Source
This commit is contained in:
Mark Olesen
2019-03-14 16:24:20 +01:00
parent 2336b46cd6
commit 659be87944
6 changed files with 97 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

30
applications/test/IjkField/Test-IjkField.C
View File

@ -49,6 +49,36 @@ Ostream& print(const IjkField<T>& fld)
int main(int argc, char *argv[])
{
// Basic addressing checks
#if 0
{
ijkAddressing addr1(3, 4, 5);
Info<< "addressing: " << addr1.sizes() << nl;
Info<< "index of (2,2,2) " << addr1.index(2,2,2) << nl;
for (const label idx : labelRange(0, addr1.size()))
{
Info<< "index of " << idx << " => " << addr1.index(idx) << nl;
}
for (label k=0; k < addr1.sizes().z(); ++k)
{
for (label j=0; j < addr1.sizes().y(); ++j)
{
for (label i=0; i < addr1.sizes().x(); ++i)
{
labelVector ijk(i,j,k);
Info<< "index of " << addr1.index(ijk)
<< " <= " << ijk << nl;
}
}
}
}
#endif
// Create with inconsistent sizes
IjkField<label> field0({1, 2, 3}, identity(10));

22
applications/test/PDRblockMesh/Test-PDRblockMesh.C
View File

@ -157,6 +157,28 @@ int main(int argc, char *argv[])
<< " length = " << mesh.width(mid) << nl;
}
// Test findCell
{
Info<< nl << "findCell:" << nl;
for
(
const point& pt
: {
mesh.bounds().centre(),
mesh.bounds().min() - 0.1 * mesh.bounds().span(),
mesh.bounds().max() + 0.1 * mesh.bounds().span()
}
)
{
labelVector ijk = mesh.findCell(pt);
Info<< " " << pt << " = " << ijk;
if (cmptMin(ijk) < 0) Info<< " [not found]";
Info<< nl;
}
}
Info<< nl;
// Fatal with FULLDEBUG