背景
我正在处理大的几何形状,这些几何形状将被网格划分或分解成更小的部分然后重新分组。因此,例如,如果一个块被分成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
编辑 - gpmetisfile是我用来将几何体分解成碎片的方法。它具有与第n个元素对应的第n行的格式,该第n个元素只具有单个值,即它所属的补丁。 ElementsList是几何中元素的列表。