matlab - 四面体面的大型、大部分为空的数据结构

我试图将四面体中每个面的节点 ID 与其对应的四面体 ID 联系起来。

tetras = [1 2 3 4  % Tetra 1
          5 6 7 8] % Tetra 2

对于 tetra 1，有四个面:

faces = [1 2 3; 1 2 4; 1 3 4; 2 3 4] % Notice these are sorted

然后我想将这些存储在数据结构中:

tet_for_face = cell(8,8,8) % 8 allows for the maximum node id

tet_for_face{1,2,3} = 1;
tet_for_face{1,2,4} = 1;
tet_for_face{1,3,4} = 1;
tet_for_face{2,3,4} = 1;

这意味着我可以在 O(1) 中找到任何特定人脸的 tetra ID:

tet_for_face{2,3,3}
ans = []
tet_for_face{2,3,4}
ans = 1

这种方法的问题在于它需要连续的内存。随着我的网格变大，我的内存不足:

cell(1000, 1000, 1000)
??? Error using ==> cell
Out of memory. Type HELP MEMORY for your options.

我还尝试过使用嵌套单元格:

tet = cell(num_nodes, 1);
tet2 = cellfun(@(x) cell(num_nodes, 1), tet, 'UniformOutput', 0);
tet3 = cellfun(@(x) cellfun(@(y) cell(num_nodes, 1), x, 'UniformOutput', 0), tet2, 'UniformOutput', 0);

tet3{2}{3}{4} = 1;
...

尽管这适用于小网格，并且不需要连续内存(AFAIK)，但它有一个坏习惯，即在 N=1000 时使 MATLAB 崩溃。

有任何想法吗？

最佳答案

在尝试了一些稀疏数组(只能是 1D 或 2D，而不是 3D)并且一无所获后，我决定使用 container.Map (HashMap)。

我使用了字符串键，我发现生成它们的最快方法是使用 sprintf(而不是 int2str 或 mat2str)

示例代码:

tet = containers.Map;
for tetra_id in tetras
    for face in faces
        face_key = sprintf('%d ', face);
        tet(face_key) = tetra_id;

这给了我一张像这样的 map :

tet('1 2 3') = 1

关于matlab - 四面体面的大型、大部分为空的数据结构，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/12724751/