如果
S=稀疏（I，J，V）；
那么这个函数是稀疏的“逆”，从某种意义上说
V=聚集（S，I，J）；
对于小矩阵，等效（且更有效）的版本是
S(sub2ind(size(S),I,J));
此版本分解为 prod(size) 大于 2^31 左右的矩阵。 对于方阵，这会将大小限制为大约 30,000 x 30,000。 对于有限元等用途，30,000 x 30,000 的稀疏矩阵实际上很小。 除了在网格细分例程中，大多数情况下不需要收集操作，它需要能够在稀疏矩阵中查找顶点对以恢复边标识符（以便可以细分这些边。）
在这种情况下，对于大网格，我发现不会溢出 31 位整数的唯一方法是使用这种收集操作的实现。
性能相当不错。 因为从稀疏矩阵中读取数据，所以每次访问预计需要 log(n) 时间 (n=size(S,1))。 对于 m 次访问，运行时间因此应为 mlogn。 相关的屏幕截图
