鲲鹏社区首页
中文
注册
我要评分
文档获取效率
文档正确性
内容完整性
文档易理解
在线提单
论坛求助

cvdot

复数向量点积运算。

函数接口声明如下:

float型数的操作:

vsip_cscalar_f (vsip_cvdot_f)(const vsip_cvview_f *a, const vsip_cvview_f *b);

参数

参数名

描述

取值范围

输入/输出

a

指向源复数向量的指针。

非空

输入

b

指向源复数向量的指针。

非空

输入

异常输入

空指针输入时,函数直接返回nan复数;数据长度为0时,返回0复数。

示例

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
#include <stdio.h>
#include <stdlib.h>
#include <stdint.h>
#include "kvsip.h"
#include "vsip.h"
#include "vsip_type.h"

#define BUFFER_SIZE_T 10
void CvdotExample()
{
    float src1[BUFFER_SIZE_T] = {1.64f, 1.63f, -1.09f, 0.71f, -3.20f, -0.43f, 0.41f, -4.83f, 5.36f, -4.40f};
    float src2[BUFFER_SIZE_T] = {-0.23f, 3.15f, -7.11f, 5.94f, 4.38f, 6.51f, -7.28f, -3.27f, 1.08f, 3.28f};
    int64_t stride = 1;
    int64_t cvLen = BUFFER_SIZE_T / 2;
    uint64_t offset = 0;

    vsip_cblock_f *block_a = vsip_cblockbind_f(src1, NULL, cvLen, VSIP_MEM_NONE);
    vsip_cblock_f *block_b = vsip_cblockbind_f(src2, NULL, cvLen, VSIP_MEM_NONE);

    vsip_cvview_f *a = vsip_cvbind_f(block_a, offset, stride, cvLen);
    vsip_cblockadmit_f(block_a, VSIP_TRUE);
    vsip_cvview_f *b = vsip_cvbind_f(block_b, offset, stride, cvLen);
    vsip_cblockadmit_f(block_b, VSIP_TRUE);

    vsip_cscalar_f dst = vsip_cvdot_f(a, b);

    vsip_cvalldestroy_f(a);
    vsip_cvalldestroy_f(b);

    printf("src1:\n");
    for (int32_t i = 0; i < cvLen; ++i) {
        printf("%.4f ", src1[i * 2]);
        printf("%.4f    ", src1[i * 2 + 1]);
    }
    printf("\nsrc2:\n");
    for (int32_t i = 0; i < cvLen; ++i) {
        printf("%.4f ", src2[i * 2]);
        printf("%.4f    ", src2[i * 2 + 1]);
    }
    printf("\ndst: %.4f %.4f\n", dst.r, dst.i);
}

int main(void) {
    CvdotExample();
    return 0;
}

运行结果:

1
2
3
4
5
src1:
1.6400 1.6300    -1.0900 0.7100    -3.2000 -0.4300    0.4100 -4.8300    5.3600 -4.4000
src2:
-0.2300 3.1500    -7.1100 5.9400    4.3800 6.5100    -7.2800 -3.2700    1.0800 3.2800
dst: -11.7540 17.2035