按照随机顺序对计算结果进行采样

Wayne

目前做一个后处理算法，需要按照随机的位置对计算结果进行采样。后处理中需要用到复杂的算法，目前是基于python的类库实现。
最简单的想法是在python中调用probe，但是这样每采样一次就需要重新从硬盘读取整个网格的结果，太慢。
另一个想法是写一个of的后处理程序，首先读入整个网格的结果进内存，然后调用python代码，有个问题是此时python怎么访问读进来的volScalarField？即：

int main(int argc, char *argv[])
{
    #include "setRootCase.H"

    #include "createTime.H"
    #include "createMesh.H"

    Info<< "Reading field T\n" << endl;
    volScalarField T
    (
        IOobject
        (
            "T",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh
    );

    调用一段python，里面有一个函数的功能是：对于python随机生成的坐标(x,y,z)，能够返回T[mesh.findCell(point(x,y,z))]
}

如果能在python中实现读入mesh和field，以及mesh.findCell功能，问题就简单了。

各位有何想法？多谢！

Wayne

发现paraview导出csv然后处理比较简单，虽然这样没有网格信息了，只能基于网格点上的值对任意点的T进行插值

李东岳

重要的是在后处理程序中，能够使用puthon生成坐标点。这就需要 C++ 和 python 混合编程。对python不是很明白，百度了一下看起来也不想C++ FORTRAN那样简单。

既然是随机坐标，为什么不用c++生成随机坐标？

Wayne

准备试试vtk，好像没有信息损失。。

Wayne

@李东岳生成的随机坐标与前一个点的值有关

程迪

参考这个stack overflow 问题：https://stackoverflow.com/questions/6915191/simple-ipc-between-c-and-python-cross-platform
推荐是用zeromq:http://www.zeromq.org/
这样对c++/python两边改动应该是最小的吧。

Wayne

@程迪厉害。确实有一些计算机方面的方法可以搞，不过用vtk解决了，多谢！

CFD中文网

按照随机顺序对计算结果进行采样