24 R Y(dot_complex)(C *x, INT n)
29 for (k = 0, dot = K(0.0); k < n; k++)
30 dot += CONJ(x[k])*x[k];
36 R Y(dot_double)(R *x, INT n)
41 for (k = 0, dot = K(0.0); k < n; k++)
49 R Y(dot_w_complex)(C *x, R *w, INT n)
54 for (k = 0, dot = K(0.0); k < n; k++)
55 dot += w[k]*CONJ(x[k])*x[k];
61 R Y(dot_w_double)(R *x, R *w, INT n)
66 for (k = 0, dot = K(0.0); k < n; k++)
67 dot += w[k]*x[k]*x[k];
74 R Y(dot_w_w2_complex)(C *x, R *w, R *w2, INT n)
79 for (k = 0, dot = K(0.0); k < n; k++)
80 dot += w[k]*w2[k]*w2[k]*CONJ(x[k])*x[k];
86 R Y(dot_w2_complex)(C *x, R *w2, INT n)
91 for (k = 0, dot = K(0.0); k < n; k++)
92 dot+=w2[k]*w2[k]*CONJ(x[k])*x[k];