按照随机顺序对计算结果进行采样
-
目前做一个后处理算法,需要按照随机的位置对计算结果进行采样。后处理中需要用到复杂的算法,目前是基于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进行插值
-
重要的是在后处理程序中,能够使用puthon生成坐标点。这就需要 C++ 和 python 混合编程。对python不是很明白,百度了一下看起来也不想C++ FORTRAN那样简单。
既然是随机坐标,为什么不用c++生成随机坐标?
-
准备试试vtk,好像没有信息损失。。
-
@李东岳 生成的随机坐标与前一个点的值有关
-
参考这个stack overflow 问题:https://stackoverflow.com/questions/6915191/simple-ipc-between-c-and-python-cross-platform
推荐是用zeromq:http://www.zeromq.org/
这样对c++/python两边改动应该是最小的吧。
-
@程迪 厉害。确实有一些计算机方面的方法可以搞,不过用vtk解决了,多谢!