提示:该函数属于Vector Mathematical Functions 章节
前言
Vector Mathematical Functions章节主要提供向量之间的运算函数,Intel官网说明链接如下。v?MulByConj (intel.com)
一、v?MulByConj函数
v?MulByConj函数的功能是计算一个复数向量元素与另一个复数向量元素的共轭相乘的结果,适用于MKL_Complex8和MKL_Complex16类型的向量,其中?代表了传进去的数据类型,前面的两种数据类型分别对应的字母分别为c,z。
使用语法为v?MulByConj(n, a, b, y),n为指定要计算向量的元素个数,a和b为指向输入进行计算的向量的指针,y为指向输出结果的指针。计算公式如下
注意这里是a和b向量里逐个对应位置的元素进行相乘。
较新版的MKL库还提供可以指导计算步长的函数,即v?MulByConjI(n,a,inca,b,incb,y,incy),其中inca和incy分别指代在a和y向量元素中跳跃的步长,假如a = [2+1i 4+2i 6+4i 8+2i],b = [1+2i 3+2i 5+2i 7+2i], inca和incb为2,incy为1,则计算a中[2+1i 6+4i]和b中[1+2i 5+2i]的复共轭相乘的结果,放到y中第1和第2个位置。
下面提供在C语言中应用v?MulByConj的例子。
二、使用示例
1.引入库
代码如下(示例):
#include <stdio.h>
#include "mkl.h"
2.引用v?MulByConj函数进行一个向量与另一个向量的复共轭相乘计算
代码如下(示例):
int main()
{
MKL_Complex8* a = (MKL_Complex8*)mkl_calloc(4, sizeof(MKL_Complex8), 64);
MKL_Complex8* b = (MKL_Complex8*)mkl_calloc(4, sizeof(MKL_Complex8), 64);
MKL_Complex8* y = (MKL_Complex8*)mkl_calloc(4, sizeof(MKL_Complex8), 64); // 内存分配
for (int ii = 0; ii < 4; ii++)
{
a[ii].real = ii * 2 + 2;
a[ii].imag = ii * 2 + 4;
b[ii].real = ii * 2 + 1;
b[ii].imag = ii * 2 + 3;
}
vcMulByConj(4, a, b, y);
for (int ii = 0; ii < 4; ii++)
{
printf("a[%d].real is %f, a[%d].imag is %f,", ii, a[ii].real, ii, a[ii].imag);
printf("b[%d].real is %f, b[%d].imag is %f,", ii, b[ii].real, ii, b[ii].imag);
printf("y[%d].real is %f, y[%d].imag is %f\n", ii, y[ii].real, ii, y[ii].imag); // 打印计算结果
}
vcMulByConjI(2, a, 2, b, 2, y, 1);
for (int ii = 0; ii < 4; ii += 2)
{
printf("a[%d].real is %f, a[%d].imag is %f,", ii, a[ii].real, ii, a[ii].imag);
printf("b[%d].real is %f, b[%d].imag is %f,", ii, b[ii].real, ii, b[ii].imag);
printf("y[%d].real is %f, y[%d].imag is %f\n", ii, y[ii].real, ii, y[ii].imag); // 打印计算结果
}
mkl_free(a);
mkl_free(b);
mkl_free(y);
}
运行结果如下。
总结
以上即为MKL库的v?MulByConj函数的简介和示例。