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



  • 目前做一个后处理算法,需要按照随机的位置对计算结果进行采样。后处理中需要用到复杂的算法,目前是基于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功能,问题就简单了。

    各位有何想法?多谢!



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


  • CORE 网格教授 OpenFOAM教授 管理员

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

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



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



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


  • OpenFOAM讲师

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



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


 

Forest
Mountains