fenics 的作者之一，A. Logg，写了一篇关于存储网格的数据结构的非常好的论文。该论文是 A. Logg (2009)。计算网格的有效表示 http://arxiv.org/abs/1205.3081

事实上，它总是在存储所有拓扑信息（节点周围的节点、节点周围的面等）或必须在使用它们之前重新计算它们之间进行权衡。如果您存储它们，那么您最终可能会存储大量数据，如果您使用迭代器重新计算它们将花费一些（太多）时间。根据算法（通常对于 FEM，对于线性元素，您只需要多面体节点来组装局部刚度矩阵），您可能需要使用一些特定的遍历，然后您可以动态生成所需的内容。

无论如何，在 fenics 中提出的解决方案基于许多整数句柄数组，并且可能是最接近最有效的方法（在内存和访问方面）。

在 deal.II 中，我们基本上只使用向量。映射太慢并且数据分散在内存中，所以如果键是整数并且在给定范围内，我们通常不使用它们。例如，对于单元格之间的连通性，您可以使用数组（STL 向量）在其中存储相邻索引，以便相邻索引$4i\ldots 4i+3$对应单元格$i$（在 2d 中，因为每个单元有 4 个邻居）。这样您就可以直接进行数组查找，而不是在 STL 映射中遍历树。

您可能已经意识到这一点，但是我们已经对我们在这个库中使用的数据结构进行了大量研究，并且我们使用的方案导致缓存未命中的数量非常少。如果您想提高效率，您将花费大量时间提出我们使用的相同类型的解决方案。