问Matlab中顶点面隶属度矩阵的矢量化

EN

我正在Matlab上做一些网格的实验。以下是我所拥有的：

• 包含所有顶点的3X#矩阵
• 3X # faces 矩阵，包含我的网格中的所有面。

以下是我正在努力高效地创造的东西：

• 一个#Verts X#面对二进制矩阵Mem。一个条目Mem(i，j) = 1，如果顶点i是三角形j的顶点之一。

下面是我以一种非常低效的方式(我对Matlab矢量化不太了解，所以请原谅我的无知)：

numVerts = size(verts, 2);
numFaces = size(faces, 2);
mem = sparse(numVerts, numFaces);

for i = 1:size(verts, 2)
    mem(i,:) = any(faces == i);
end

下面是一些能够测试这段代码的数值数据：

verts = [0 0 1 1; 0 1 0 1; 0 0 0 0];
faces = [1 3; 2 2; 3 4];

矩阵的每一列都是顶点的坐标。例如，这里的4个顶点位于位置(0,0,0)，(0,1,0)，(1,0,0)和(1,1,0)。

面矩阵的每一列都是属于一个面的3个顶点的指数。例如，上面定义的网格的面是两个三角形，它们被定义为包含顶点(1,2,3)和(3,2,4)。

这将创建一个简单的网格，包含4个顶点和2个三角形。本例中的mem矩阵是：

mem = [1 0; 1 1; 1 1; 0 1]

例如，在这里，mem( 2 ,1) =1，这意味着面1包含顶点数2，任何关于我如何更有效地完成它的指针？