MKL库---cblas_?copy

提示：该函数属于BLAS Level 1 Routines and Functions章节

文章目录

• 前言
• 一、cblas_?copy函数
• 二、使用示例
  • 1.引入库
  • 2.使用cblas_?copy函数进行向量复制操作
• 总结

前言

BLAS Level 1 Routines and Functions章节主要提供向量与向量之间的操作函数，该函数在Intel官网说明链接如下。

cblas_?copy (intel.com)

一、cblas_?copy函数

cblas_?copy函数的功能是将一个向量的元素复制给另一个向量，适用于float，double，MKL_Complex8，MKL_Complex16类型的向量，其中？代表了传进去的数据类型，前面的四种数据类型分别对应的字母分别为s，d，c，z。

使用语法为cblas_?copy(n, a, inca, y, incy)，n为指定要计算向量的元素个数，a为指向输入进行计算的向量的指针，y为指向输出结果的指针。inca和incy分别指代在a和y向量元素中跳跃的步长，假如a = [2 4 6 8]，inca为2，incy为1，则将a中[2 6]复制到y中第1和第2个位置。

下面提供在C语言中应用cblas_?copy的例子。

二、使用示例

1.引入库

代码如下（示例）：

#include <stdio.h>
#include "mkl.h"

2.使用cblas_?copy函数进行向量复制操作

代码如下（示例）：

int main()
{
    MKL_Complex8* a = (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;
    }
    cblas_ccopy(4, a, 1, y, 1);
    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("y[%d].real is %f, y[%d].imag is %f\n", ii, y[ii].real, ii, y[ii].imag); // 打印复制结果
    }
    
    mkl_free(a);
    mkl_free(y);
}

代码运行结果如下

总结

以上即为MKL库的cblas_?copy函数的简介和示例。