输入如下:点数组 (x,y,z,f)。都是 np.arrays。这些 (x,y,z) 只是单元格的网格,f 是每个网格点的能量值。我正在尝试解决以下任务:构造渗透此单元格的等值面(包含 (x,y,z) 点的表面,其 f 值与某些容差值相同),即您可以绘制路径(在这个表面中)它是连续的并且接触单元的相对(不相交)平面(或相对的顶点,边缘)。无限互连是指我有周期性的边界条件(但那是另一回事)。我想到了以下伪代码:
input (x, y, z, f)
take minimum f value
determine iso-surfaces at this f value
check connectivity of iso-surfaces
are they interconnected inside the cell&and they intersect opossite vertexes ( or edges, plane surfaces) of the cell?
yes
good work
break
no
repeat for bigger f value
那么,你能推荐一些工具来解决它吗?