背景:
我正在与大的几何图形,将网格或分解成更小的部分,然后重新组合。因此,举例来说,如果一个块被分成16块,它可能会被重新组织成4个“补丁”,每个补丁有4个块或元素。在这个过程之后,我仍然需要跟踪元素。我指定patchIndex(我有一个补丁数量的列表)作为“键”,它将返回补丁中的所有元素及其1)局部索引(补丁中元素的索引)和2)全局索引(整个几何体中元素的索引)。在
问题:
我究竟如何才能将这些信息放入HDF5文件中?在
我的代码:
下面是我如何设置字典,如果这有助于了解:def readAscii(ElementsList,gpmetisfile):
f = open(gpmetisfile, 'r')
indexer={}
i=0
for line in gpmetisfile:
patchIndex = eval(line)
if patchIndex in indexer:
localIndex=indexer[patchIndex]
else:
indexer[patchIndex]=0
test = ElementsList[i].setLocalIndex(patchIndex,localIndex)
if test:
indexer[patchIndex] +=1
编辑-gpmetis文件是我用来将几何体分解成碎片的。它的格式是第n行,对应于第n个元素,该元素只有一个值,即它所属的面片。ElementsList是几何体中元素的列表。在