gr cpu指令加速库volk使用复共轭相乘

东方耀

gr cpu指令加速库volk使用复共轭相乘

1. #include <iostream>
   
2. #include <volk/volk.h>
   
3. #include <complex>
   
4. #include <vector>
   
5. 
   
6. using namespace std;
   
7. 
   
8. typedef std::complex<float> gr_complex;   //这个float精度与fftwf一致
   
9. 
   
10. void print_vector_complex(std::vector<gr_complex>& v){
    
11.     for(std::vector<gr_complex>::iterator it=v.begin();it!=v.end();it++){
    
12.         cout << "value of complex: real=" << (*it).real() << ",imag=" << (*it).imag() << endl;
    
13.     }
    
14. 
    
15. }
    
16. 
    
17. 
    
18. 
    
19. void main_jjj(){
    
20.     int noutput_items = 10;
    
21.     std::vector<gr_complex> out;   //这是容器对象 vector动态数组
    
22.     std::vector<gr_complex> in_1;
    
23.     std::vector<gr_complex> in_2;
    
24. 
    
25.     out.resize(noutput_items);
    
26.     in_1.resize(noutput_items);  // 这样就不能.push_back()了
    
27.     in_2.resize(noutput_items);
    
28. 
    
29. 
    
30.     // 输入的数据  赋值
    
31.     for(int i=0;i<noutput_items; i++){
    
32.         //d_hold_in.push_back(gr_complex(i, 0));
    
33.         //d_hold_in[i] = gr_complex(i,i+10);
    
34.         in_1[i] = gr_complex(i, i);
    
35.         in_2[i] = gr_complex(i+5, i);
    
36.     }
    
37. 
    
38. 
    
39.     cout << "输入数据赋初值in_1:" << endl;
    
40.     print_vector_complex(in_1);
    
41. 
    
42.     cout << "输入数据赋初值in_2:" << endl;
    
43.     print_vector_complex(in_2);
    
44. 
    
45. 
    
46.     //target_link_libraries(main -lvolk)
    
47.     //typedef std::complex<float>   lv_32fc_t;  只能是复数
    
48.     //参数类型(lv_32fc_t *, const lv_32fc_t *, const lv_32fc_t *, unsigned int)
    
49.     // 函数作用：复数数组与复数数组相乘 就是对应元素复数相乘的结果  这个可能优化了
    
50.     // volk_32fc_x2_multiply_32fc(
    
51.     //             &out[0],  //容器对象首元素的地址   只有这个变了  *out
    
52.     //             &in_1[0],  //     *in1
    
53.     //             &in_2[0],   //  没变     *in2
    
54.     //             noutput_items);
    
55. 
    
56.     //volk_32fc_x2_multiply_conjugate_32fc(out, in0, in1, noi);
    
57.     // 函数作用：in1*in2的复共轭=out
    
58.     volk_32fc_x2_multiply_conjugate_32fc(
    
59.                 &out[0],  //容器对象首元素的地址   只有这个变了  *out
    
60.                 &in_1[0],  //     *in1
    
61.                 &in_2[0],   //  没变     *in2
    
62.                 noutput_items);
    
63.     cout << "经过volk_32fc_x2_multiply_conjugate_32fc函数之后【in1*in2的复共轭=out】：" << endl;
    
64. 
    
65. 
    
66.     cout << "输出结果out【只有这个变了】:" << endl;
    
67.     print_vector_complex(out);
    
68. 
    
69. 
    
70. }
    
71. 
    
72. //gr cpu指令加速库volk使用复共轭相乘
    
73. 
    
74. int main(){
    
75.     main_jjj();    //自己改成 好理解
    
76.     return 0;
    
77. }
    
78. 
    
79. 
    

复制代码

zouqiqi

让天下人人学会人工智能！人工智能的前景一片大好！