按列主顺序重新排序 3D 矢量三元组的速度很慢

Question

我有很多 (x1,y1,z1),(x2,y2,z2),(x3,y3,z3) 单精度向量三元组，我想对它们重新排序，所以 (x1,y1,z1),(x2,y2,z2),(x3,y3,z3) 变成 (x1,x2,x3,0,y1,y2,y3,0,z1,z2,z3,0)

目标是为基于 SSE 的计算准备数据集。我有以下代码来执行此操作：

for (int i=0;i<count;i++)
{
    Vect3F p0 = get_first_point(i);
    Vect3F p1 = get_second_point(i);
    Vect3F p2 = get_third_point(i);
    int idx = i*3;
    scratch[idx] = Vec4F(p0.x, p1.x, p2.x, 0); // These 3 rows are the slowest
    scratch[idx+1] = Vec4F(p0.y, p1.y, p2.y, 0);
    scratch[idx+2] = Vec4F(p0.z, p1.z, p2.z, 0);
}

循环的最后 3 行非常慢，它们占用了整个算法 90% 的时间！

正常吗？我可以让这种洗牌更快吗？ （scratch是一个静态变量，并且是16对齐的。该函数被频繁调用，所以我认为scratch的块不应该从缓存中消失。）

Answer 1

首先，您不应该创建 3 个临时矢量对象。
而不是：

tri = triangles[i];
Vect3F p0 = points[indices[tri]];
Vect3F p1 = points[indices[tri+1]];
Vect3F p2 = points[indices[tri+2]];

您应该只使用 memcpy() 复制数据；创建一个适用于整个集合的循环并复制原始数据。这是我能想到的最快的方法。

使用 3 个变量会运行很多构造函数，速度非常慢。出于同样的原因，第二种方法（来自评论）也好不了多少。