Go to the documentation of this file.
29 double c = 2.0 / (len-1.0);
31 x86_reg i = -n2*
sizeof(int32_t);
34 "movsd %4, %%xmm7 \n\t"
37 "movlhps %%xmm7, %%xmm7 \n\t"
38 "subpd %%xmm5, %%xmm7 \n\t"
39 "addsd %%xmm6, %%xmm7 \n\t"
42 #define WELCH(MOVPD, offset)\
44 "movapd %%xmm7, %%xmm1 \n\t"\
45 "mulpd %%xmm1, %%xmm1 \n\t"\
46 "movapd %%xmm6, %%xmm0 \n\t"\
47 "subpd %%xmm1, %%xmm0 \n\t"\
48 "pshufd $0x4e, %%xmm0, %%xmm1 \n\t"\
49 "cvtpi2pd (%3,%0), %%xmm2 \n\t"\
50 "cvtpi2pd "#offset"*4(%3,%1), %%xmm3 \n\t"\
51 "mulpd %%xmm0, %%xmm2 \n\t"\
52 "mulpd %%xmm1, %%xmm3 \n\t"\
53 "movapd %%xmm2, (%2,%0,2) \n\t"\
54 MOVPD" %%xmm3, "#offset"*8(%2,%1,2) \n\t"\
55 "subpd %%xmm5, %%xmm7 \n\t"\
66 :
"r"(w_data+n2),
"r"(data+n2),
"m"(c),
"r"(len)
68 "%xmm5",
"%xmm6",
"%xmm7")
81 for(j=0; j<lag; j+=2){
82 x86_reg i = -len*
sizeof(double);
89 "movapd (%2,%0), %%xmm3 \n\t"
90 "movupd -8(%3,%0), %%xmm4 \n\t"
91 "movapd (%3,%0), %%xmm5 \n\t"
92 "mulpd %%xmm3, %%xmm4 \n\t"
93 "mulpd %%xmm3, %%xmm5 \n\t"
94 "mulpd -16(%3,%0), %%xmm3 \n\t"
95 "addpd %%xmm4, %%xmm1 \n\t"
96 "addpd %%xmm5, %%xmm0 \n\t"
97 "addpd %%xmm3, %%xmm2 \n\t"
100 "movhlps %%xmm0, %%xmm3 \n\t"
101 "movhlps %%xmm1, %%xmm4 \n\t"
102 "movhlps %%xmm2, %%xmm5 \n\t"
103 "addsd %%xmm3, %%xmm0 \n\t"
104 "addsd %%xmm4, %%xmm1 \n\t"
105 "addsd %%xmm5, %%xmm2 \n\t"
106 "movsd %%xmm0, (%1) \n\t"
107 "movsd %%xmm1, 8(%1) \n\t"
108 "movsd %%xmm2, 16(%1) \n\t"
110 :
"r"(autoc+j),
"r"(data+
len),
"r"(data+len-j)
118 "movapd (%3,%0), %%xmm3 \n\t"
119 "movupd -8(%4,%0), %%xmm4 \n\t"
120 "mulpd %%xmm3, %%xmm4 \n\t"
121 "mulpd (%4,%0), %%xmm3 \n\t"
122 "addpd %%xmm4, %%xmm1 \n\t"
123 "addpd %%xmm3, %%xmm0 \n\t"
126 "movhlps %%xmm0, %%xmm3 \n\t"
127 "movhlps %%xmm1, %%xmm4 \n\t"
128 "addsd %%xmm3, %%xmm0 \n\t"
129 "addsd %%xmm4, %%xmm1 \n\t"
130 "movsd %%xmm0, %1 \n\t"
131 "movsd %%xmm1, %2 \n\t"
132 :
"+&r"(i),
"=m"(autoc[j]),
"=m"(autoc[j+1])
133 :
"r"(data+len),
"r"(data+len-j)