C# 计算两两坐标之间的距离(SIMD加速)

问题描述:对于Point[] points坐标数组(Point是坐标点,包含X与Y),请计算两两坐标点之间的距离。

思路:既然是算距离,没有别的要求,只需要最简单的曼哈顿距离就行。

测试数据:

cs 复制代码
Random random = new Random();
int n = 3000;
Point[] points = new Point[n];
for (int i = 0; i < n; i++)
{
    points[i] = new Point(random.Next(10), random.Next(10));
}

方法一:采用普通遍历法

cs 复制代码
Stopwatch sw = Stopwatch.StartNew();
int m = (n * n - n) / 2;
List<int> list2 = new List<int>(m);
for (int i = 0; i < n; i++)
{
    Point p= points[i];
    for (int j = i + 1; j < n; j++)
    {
        Point p2 = points[j];
        int dif = Math.Abs(p.X - p2.X) + Math.Abs(p.Y - p2.Y);
        list2.Add(dif);
    }
}
sw.Stop();
Console.WriteLine($"普通耗时:{sw.ElapsedMilliseconds}ms");

方法二:采用SIMD加速

cs 复制代码
 public static int[] CalculateDistance(Point[] points)
        {
            int n = points.Length;
            int m = (n * n - n) >> 1;
            int[] ints = new int[m];
            int k = 0;
            if (n <= 4)// 如果数据量小于等于4采用指针处理
            {
                fixed (Point* ptr = points)
                {
                    for (int i = 0; i < n; i++)
                    {
                        Point* p1 = ptr + i;
                        for (int j = i + 1; j < n; j++)
                        {
                            Point* p2 = ptr + j;
                            int dif = Math.Abs(p1->X - p2->X) + Math.Abs(p1->Y - p2->Y);
                            ints[k++] = dif;
                        }
                    }
                    return ints;
                }
            }
            int vsize = Vector128<int>.Count;
            int t = 0;
            int[] xs = new int[n];
            int[] ys = new int[n];
            fixed (Point* ptr = points)
            {
                fixed (int* xptr = xs, yptr = ys, rptr = ints)
                {
                    // 分离x与y
                    for (int i = 0; i < n; i++)
                    {
                        Point* p = ptr + i;
                        *(xptr + i) = p->X;
                        *(yptr + i) += p->Y;
                    }
                    // 开始向量计算
                    for (int i = 0; i < n; i += vsize)
                    {
                        Vector128<int> vx1 = *(Vector128<int>*)(xptr + i);
                        Vector128<int> vy1 = *(Vector128<int>*)(yptr + i);
                        for (int j = i + 1; j < n; j++)
                        {
                            Vector128<int> vx2 = *(Vector128<int>*)(xptr + j);
                            Vector128<int> vy2 = *(Vector128<int>*)(yptr + j);
                            // 向量运算
                            var vDifference = Sse2.Add(Vector128.Abs(Sse2.Subtract(vx1, vx2)), Vector128.Abs(Sse2.Subtract(vy1, vy2)));
                            int* maskPtr = (int*)&vDifference;
                            int h = n - j;
                            for (int a = 0; a < vsize; a++)
                            {
                                if (a >= h) break;
                                *(rptr + t++) = *(maskPtr + a);
                            }
                        }
                    }
                }
            }
            return ints;
        }

注:根据CPU选择合适的Vector,比如Vector256、Vector512等等。本文只用Vector128.

测试结果(Release环境):

|----------|------|------|------|-------|-------|
| 数据量 | 3000 | 5000 | 7000 | 10000 | 30000 |
| 普通耗时(ms) | 43 | 127 | 211 | 538 | 3872 |
| 向量耗时(ms) | 19 | 56 | 68 | 190 | 1092 |

结论:明显向量运算性能优于普通计算。

相关推荐
Yhame.15 分钟前
【使用层次序列构建二叉树(数据结构C)】
c语言·开发语言·数据结构
言之。21 分钟前
【Go语言】RPC 使用指南(初学者版)
开发语言·rpc·golang
投笔丶从戎1 小时前
Kotlin Multiplatform--01:项目结构基础
android·开发语言·kotlin
杜小暑2 小时前
动态内存管理
c语言·开发语言·动态内存管理
想不明白的过度思考者2 小时前
Java从入门到“放弃”(精通)之旅——JavaSE终篇(异常)
java·开发语言
我真的不会C2 小时前
QT窗口相关控件及其属性
开发语言·qt
CodeCraft Studio2 小时前
Excel处理控件Aspose.Cells教程:使用 Python 在 Excel 中进行数据验
开发语言·python·excel
火柴盒zhang2 小时前
websheet之 编辑器
开发语言·前端·javascript·编辑器·spreadsheet·websheet
景天科技苑2 小时前
【Rust】Rust中的枚举与模式匹配,原理解析与应用实战
开发语言·后端·rust·match·enum·枚举与模式匹配·rust枚举与模式匹配
阿让啊2 小时前
C语言中操作字节的某一位
c语言·开发语言·数据结构·单片机·算法