@马乔 您好，这么讲的话，也许是可能是截断误差影响。然后我之所以用findCell，是因为我用到一个delta function，做一个点的物理量分布到周围网格范围内流体网格的distribution的操作，和点周围流体网格物理量插值到点的interpolation操作。就是grid to point and point to grid。需要找到点所在cell及其周围一圈两倍网格间距范围内的cell，不知道想实现这种插值分布函数，还有没有其他更简单的方法吗？请问您有什么看法吗？谢谢！

@李东岳 在 OpenFOAM编程findCell的诡异问题，对同一坐标寻找cell结果不一致 中说：

另外你没用过primitiveMesh的findCell()试试么，这种

point a(1,2,3);
mesh.findCell(a);

我几年前记得在某个SCI里面看到过说meshSearch好像效率低，不过具体细节我忘了

老师，mesh search和mesh的findCell我都是用了的，一样的结果，都是循环外单独找那两个点是对的，循环内就不对，是重复的CellID

@李东岳 在 OpenFOAM编程findCell的诡异问题，对同一坐标寻找cell结果不一致 中说：

代码上看起来是没问题。

label cellId = searchEngine_.findCell(tmpOffSetPoint);

这句话之前赋值为0然后再赋值呢？

老师，这样做也是一样的结果，还是重复的CellID

我把代码和测试算例精简后打包发在了网盘，如果有兴趣的老师和朋友可以帮忙看一下，谢谢！blockMesh后运行
链接：https://pan.baidu.com/s/1HEy6w8MgbeXu9EsMcayung
提取码：1111
--来自百度网盘超级会员V4的分享

各位老师同学好，最近写的一个代码的一个小部分用到了findcell去找点所在的cell，但是发现了一个奇怪的问题。对于两个点，我单独findcell出来结果是对的，但是在循环里面算出来这两个点后再findcell，结果不对。首先，代码我已经把无用的都删了后如下：

#include "fvCFD.H"
#include "meshSearch.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
{
    argList::addBoolOption
    (
        "noOverwrite",
        "increment time to avoid overwriting initial fields"
    );
    #include "setRootCase.H"
    #include "createTime.H"

    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

    #include "createMesh.H"
IOdictionary blockMeshDict
(
    IOobject
    (
        "blockMeshDict",
        runTime.system(),
        mesh,
        IOobject::MUST_READ,
        IOobject::NO_WRITE
    )
);  
IOdictionary NetProperties
(
    IOobject
    (
        "NetProperties",
        runTime.constant(),
        mesh,
        IOobject::MUST_READ,
        IOobject::NO_WRITE
    )
);
const dictionary& NetDict(NetProperties.subDict("deltaFunction"));
const label  nx(blockMeshDict.lookup<scalar>("nx"));
const label  ny(blockMeshDict.lookup<scalar>("ny"));
const label  nz(blockMeshDict.lookup<scalar>("nz"));
const scalar lx(blockMeshDict.lookup<scalar>("lx"));
const scalar ly(blockMeshDict.lookup<scalar>("ly"));
const scalar lz(blockMeshDict.lookup<scalar>("lz"));
const scalar dx = lx/nx;
const scalar dy = ly/ny;
const scalar dz = lz/nz;
const List<label> offSet_
(
    NetDict.lookup("a")
);

meshSearch searchEngine_(mesh);

vector tmpKnotPoint(1.5,0.72,1.359);
vector point1(1.47727,0.694583,1.30181);
vector point2(1.5,0.694583,1.30181);
Pout << "point1 (1.47727 0.694583 1.30181) should be inside cell " << searchEngine_.findCell(point1) << endl;
Pout << "point2(1.5 0.694583 1.30181) should be inside cell " << searchEngine_.findCell(point2) << endl;

vector tmpOffSetPoint;
label num = 0;
forAll(offSet_, zi)
{
    tmpOffSetPoint.z() = tmpKnotPoint.z() + offSet_[zi]*dz;
    forAll(offSet_, yi)
    {
        tmpOffSetPoint.y() = tmpKnotPoint.y() + offSet_[yi]*dy; 
        forAll(offSet_, xi)
        {
            tmpOffSetPoint.x() = tmpKnotPoint.x() + offSet_[xi]*dx;
            label cellId = searchEngine_.findCell(tmpOffSetPoint);
            Pout << num << "  " << cellId << "  " << tmpOffSetPoint  << endl;
        }                    
    }
}

    return 0;
}

输出如下：

point1 (1.47727 0.694583 1.30181) should be inside cell 423586
point2(1.5 0.694583 1.30181) should be inside cell 423587
0  423587  (1.47727 0.694583 1.30181)
0  423587  (1.5 0.694583 1.30181)
0  423785  (1.47727 0.72 1.30181)
0  423785  (1.5 0.72 1.30181)
0  442595  (1.47727 0.694583 1.359)
0  442595  (1.5 0.694583 1.359)
0  442793  (1.47727 0.72 1.359)
0  442793  (1.5 0.72 1.359)

循环里输出的Cell 的Id对于这两个点(1.47727 0.694583 1.30181)和(1.5 0.694583 1.30181)都是423587，但是上面的输出直接输出这两个坐标findcell的结果，那就是一个423586，另一个是423587.

难道是findcell用的方法不对，还是什么其他原因呢？请各位同学老师帮下忙。

@anubis 那就用位压缩，只需要建立一个m大小的一维数组。每个位置只存储一个该网格对应的粒子有无信息。比如一共十个粒子，网格23对应的状态就是0000100011，转换成十进制为35，只需要把35存进数组位置，因为一个整数35就代表了1，2，6三个位置处状态为1也就是这三个编号粒子存在。下次需要该网格23的粒子编号信息时候，取出35，再写个小程序转换进制判断一下，时间复杂度很低的，一般就是一点简单的位运算，就能知道换成二进制后第1，2，6个位置处是1.

@joann 我是这样的 版本不同可能情况也有所区别

我用tec读openfoam数据二进制格式无法读取，你要输出十进制。而且就算十进制，一些复杂网格的openfoam数据我用tecplot读进来就闪退，所以还是用paraview了

@桎梏 我都是在求解器main函数写点东西输出文件，写在哪都行吧，如果跟时间无关，就卸载pimple的while循环外面就行吧

@桎梏 直接openfoam代码里mesh.v()乘以U[cellI]就行吧 然后输出到一个文件里