NFFT  3.3.0
vector1.c
1 /*
2  * Copyright (c) 2002, 2015 Jens Keiner, Stefan Kunis, Daniel Potts
3  *
4  * This program is free software; you can redistribute it and/or modify it under
5  * the terms of the GNU General Public License as published by the Free Software
6  * Foundation; either version 2 of the License, or (at your option) any later
7  * version.
8  *
9  * This program is distributed in the hope that it will be useful, but WITHOUT
10  * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
11  * FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
12  * details.
13  *
14  * You should have received a copy of the GNU General Public License along with
15  * this program; if not, write to the Free Software Foundation, Inc., 51
16  * Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
17  */
18 
19 /* $Id: util.c 3483 2010-04-23 19:02:34Z keiner $ */
20 
21 #include "infft.h"
22 
24 R Y(dot_complex)(C *x, INT n)
25 {
26  INT k;
27  R dot;
28 
29  for (k = 0, dot = K(0.0); k < n; k++)
30  dot += CONJ(x[k])*x[k];
31 
32  return dot;
33 }
34 
36 R Y(dot_double)(R *x, INT n)
37 {
38  INT k;
39  R dot;
40 
41  for (k = 0, dot = K(0.0); k < n; k++)
42  dot += x[k]*x[k];
43 
44  return dot;
45 }
46 
47 
49 R Y(dot_w_complex)(C *x, R *w, INT n)
50 {
51  INT k;
52  R dot;
53 
54  for (k = 0, dot = K(0.0); k < n; k++)
55  dot += w[k]*CONJ(x[k])*x[k];
56 
57  return dot;
58 }
59 
61 R Y(dot_w_double)(R *x, R *w, INT n)
62 {
63  INT k;
64  R dot;
65 
66  for (k = 0, dot = K(0.0); k < n; k++)
67  dot += w[k]*x[k]*x[k];
68 
69  return dot;
70 }
71 
72 
74 R Y(dot_w_w2_complex)(C *x, R *w, R *w2, INT n)
75 {
76  INT k;
77  R dot;
78 
79  for (k = 0, dot = K(0.0); k < n; k++)
80  dot += w[k]*w2[k]*w2[k]*CONJ(x[k])*x[k];
81 
82  return dot;
83 }
84 
86 R Y(dot_w2_complex)(C *x, R *w2, INT n)
87 {
88  INT k;
89  R dot;
90 
91  for (k = 0, dot = K(0.0); k < n; k++)
92  dot+=w2[k]*w2[k]*CONJ(x[k])*x[k];
93 
94  return dot;
95 }