MKL 实现复数向量点乘 - 两种方法详解

MKL 提供了两种方法实现复数向量点乘：

1. 使用 MKL 的 BLAS 函数库中的 cblas_zdotu_sub 函数实现复数向量点乘。

该函数的原型为：

void cblas_zdotu_sub(const MKL_INT N, const void *X, const MKL_INT incX, const void *Y, const MKL_INT incY, void *dotc);

其中，N 为向量长度，X 和 Y 分别为长度为 N 的复数向量，incX 和 incY 为 X 和 Y 的步长，dotc 为结果。

示例代码：

#include <mkl.h>

void complex_dot_product()
{
    // 定义两个长度为 3 的复数向量
    double complex x[3] = {1.0 + 2.0*I, 2.0 + 3.0*I, 3.0 + 4.0*I};
    double complex y[3] = {2.0 + 1.0*I, 3.0 + 2.0*I, 4.0 + 3.0*I};

    // 计算复数向量点乘
    double complex dotc;
    cblas_zdotu_sub(3, x, 1, y, 1, &dotc);

    // 输出结果
    printf('Complex dot product: %g + %gi
', creal(dotc), cimag(dotc));
}

2. 使用 MKL 的 Vector Mathematics 函数库中的 vdDotc 函数实现复数向量点乘。

该函数的原型为：

void vdDotc(const MKL_INT N, const MKL_Complex16 *x, const MKL_INT incx, const MKL_Complex16 *y, const MKL_INT incy, MKL_Complex16 *result);

其中，N 为向量长度，x 和 y 分别为长度为 N 的复数向量，incx 和 incy 为 x 和 y 的步长，result 为结果。

示例代码：

#include <mkl.h>

void complex_dot_product()
{
    // 定义两个长度为 3 的复数向量
    double complex x[3] = {1.0 + 2.0*I, 2.0 + 3.0*I, 3.0 + 4.0*I};
    double complex y[3] = {2.0 + 1.0*I, 3.0 + 2.0*I, 4.0 + 3.0*I};

    // 计算复数向量点乘
    double complex dotc;
    vdDotc(3, x, 1, y, 1, &dotc);

    // 输出结果
    printf('Complex dot product: %g + %gi
', creal(dotc), cimag(dotc));
}