22 R Y(dot_complex)(C *x, INT n)
27 for (k = 0, dot = K(0.0); k < n; k++)
28 dot += CONJ(x[k])*x[k];
34 R Y(dot_double)(R *x, INT n)
39 for (k = 0, dot = K(0.0); k < n; k++)
47 R Y(dot_w_complex)(C *x, R *w, INT n)
52 for (k = 0, dot = K(0.0); k < n; k++)
53 dot += w[k]*CONJ(x[k])*x[k];
59 R Y(dot_w_double)(R *x, R *w, INT n)
64 for (k = 0, dot = K(0.0); k < n; k++)
65 dot += w[k]*x[k]*x[k];
72 R Y(dot_w_w2_complex)(C *x, R *w, R *w2, INT n)
77 for (k = 0, dot = K(0.0); k < n; k++)
78 dot += w[k]*w2[k]*w2[k]*CONJ(x[k])*x[k];
84 R Y(dot_w2_complex)(C *x, R *w2, INT n)
89 for (k = 0, dot = K(0.0); k < n; k++)
90 dot+=w2[k]*w2[k]*CONJ(x[k])*x[k];