QDP++
qdp_globalfuncs.h
Go to the documentation of this file.
1// -*- C++ -*-
2
6
7#ifndef QDP_GLOBALFUNCS_H
8#define QDP_GLOBALFUNCS_H
9
10namespace QDP
11{
12
18 //-----------------------------------------------
19 // Global sums
21
25 template<class T, class C>
26 inline typename UnaryReturn<C, FnSum>::Type_t
27 sum(const QDPType<T,C>& s1)
28 {
29 return sum(PETE_identity(s1));
30 }
31
32
34
38 template<class T, class C>
39 inline typename UnaryReturn<C, FnSum>::Type_t
40 sum(const QDPType<T,C>& s1, const Subset& s)
41 {
42 return sum(PETE_identity(s1),s);
43 }
44
45
47
53 template<class T, class C>
54 inline typename UnaryReturn<C, FnNorm2>::Type_t
55 norm2(const QDPType<T,C>& s1)
56 {
57 return sum(localNorm2(s1));
58 }
59
60 template<class T, class C>
61 inline typename UnaryReturn<C, FnNorm2>::Type_t
62 norm2(const QDPExpr<T,C>& s1)
63 {
64 return sum(localNorm2(s1));
65 }
66
67
69
75 template<class T, class C>
76 inline typename UnaryReturn<C, FnNorm2>::Type_t
77 norm2(const QDPType<T,C>& s1, const Subset& s)
78 {
79 return sum(localNorm2(s1),s);
80 }
81
82 template<class T, class C>
83 inline typename UnaryReturn<C, FnNorm2>::Type_t
84 norm2(const QDPExpr<T,C>& s1, const Subset& s)
85 {
86 return sum(localNorm2(s1),s);
87 }
88
89
91
96 template<class T1, class C1, class T2, class C2>
97 inline typename BinaryReturn<C1, C2, FnInnerProduct>::Type_t
98 innerProduct(const QDPType<T1,C1>& s1, const QDPType<T2,C2>& s2)
99 {
100 return sum(localInnerProduct(s1,s2));
101 }
102
103 template<class T1, class C1, class T2, class C2>
104 inline typename BinaryReturn<C1, C2, FnInnerProduct>::Type_t
105 innerProduct(const QDPType<T1,C1>& s1, const QDPExpr<T2,C2>& s2)
106 {
107 return sum(localInnerProduct(s1,s2));
108 }
109
110 template<class T1, class C1, class T2, class C2>
111 inline typename BinaryReturn<C1, C2, FnInnerProduct>::Type_t
112 innerProduct(const QDPExpr<T1,C1>& s1, const QDPType<T2,C2>& s2)
113 {
114 return sum(localInnerProduct(s1,s2));
115 }
116
117 template<class T1, class C1, class T2, class C2>
118 inline typename BinaryReturn<C1, C2, FnInnerProduct>::Type_t
119 innerProduct(const QDPExpr<T1,C1>& s1, const QDPExpr<T2,C2>& s2)
120 {
121 return sum(localInnerProduct(s1,s2));
122 }
123
124
126
131 template<class T1, class C1, class T2, class C2>
132 inline typename BinaryReturn<C1, C2, FnInnerProduct>::Type_t
133 innerProduct(const QDPType<T1,C1>& s1, const QDPType<T2,C2>& s2,
134 const Subset& s)
135 {
136 return sum(localInnerProduct(s1,s2),s);
137 }
138
139 template<class T1, class C1, class T2, class C2>
140 inline typename BinaryReturn<C1, C2, FnInnerProduct>::Type_t
141 innerProduct(const QDPType<T1,C1>& s1, const QDPExpr<T2,C2>& s2,
142 const Subset& s)
143 {
144 return sum(localInnerProduct(s1,s2),s);
145 }
146
147 template<class T1, class C1, class T2, class C2>
148 inline typename BinaryReturn<C1, C2, FnInnerProduct>::Type_t
149 innerProduct(const QDPExpr<T1,C1>& s1, const QDPType<T2,C2>& s2,
150 const Subset& s)
151 {
152 return sum(localInnerProduct(s1,s2),s);
153 }
154
155 template<class T1, class C1, class T2, class C2>
156 inline typename BinaryReturn<C1, C2, FnInnerProduct>::Type_t
157 innerProduct(const QDPExpr<T1,C1>& s1, const QDPExpr<T2,C2>& s2,
158 const Subset& s)
159 {
160 return sum(localInnerProduct(s1,s2),s);
161 }
162
163
165
170 template<class T1, class C1, class T2, class C2>
171 inline typename BinaryReturn<C1, C2, FnInnerProductReal>::Type_t
172 innerProductReal(const QDPType<T1,C1>& s1, const QDPType<T2,C2>& s2)
173 {
174 return sum(localInnerProductReal(s1,s2));
175 }
176
177 template<class T1, class C1, class T2, class C2>
178 inline typename BinaryReturn<C1, C2, FnInnerProductReal>::Type_t
179 innerProductReal(const QDPType<T1,C1>& s1, const QDPExpr<T2,C2>& s2)
180 {
181 return sum(localInnerProductReal(s1,s2));
182 }
183
184 template<class T1, class C1, class T2, class C2>
185 inline typename BinaryReturn<C1, C2, FnInnerProductReal>::Type_t
186 innerProductReal(const QDPExpr<T1,C1>& s1, const QDPType<T2,C2>& s2)
187 {
188 return sum(localInnerProductReal(s1,s2));
189 }
190
191 template<class T1, class C1, class T2, class C2>
192 inline typename BinaryReturn<C1, C2, FnInnerProductReal>::Type_t
193 innerProductReal(const QDPExpr<T1,C1>& s1, const QDPExpr<T2,C2>& s2)
194 {
195 return sum(localInnerProductReal(s1,s2));
196 }
197
198
200
205 template<class T1, class C1, class T2, class C2>
206 inline typename BinaryReturn<C1, C2, FnInnerProductReal>::Type_t
207 innerProductReal(const QDPType<T1,C1>& s1, const QDPType<T2,C2>& s2,
208 const Subset& s)
209 {
210 return sum(localInnerProductReal(s1,s2),s);
211 }
212
213 template<class T1, class C1, class T2, class C2>
214 inline typename BinaryReturn<C1, C2, FnInnerProductReal>::Type_t
215 innerProductReal(const QDPType<T1,C1>& s1, const QDPExpr<T2,C2>& s2,
216 const Subset& s)
217 {
218 return sum(localInnerProductReal(s1,s2),s);
219 }
220
221 template<class T1, class C1, class T2, class C2>
222 inline typename BinaryReturn<C1, C2, FnInnerProductReal>::Type_t
223 innerProductReal(const QDPExpr<T1,C1>& s1, const QDPType<T2,C2>& s2,
224 const Subset& s)
225 {
226 return sum(localInnerProductReal(s1,s2),s);
227 }
228
229 template<class T1, class C1, class T2, class C2>
230 inline typename BinaryReturn<C1, C2, FnInnerProductReal>::Type_t
231 innerProductReal(const QDPExpr<T1,C1>& s1, const QDPExpr<T2,C2>& s2,
232 const Subset& s)
233 {
234 return sum(localInnerProductReal(s1,s2),s);
235 }
236
237
238 //-----------------------------------------------------------------------------
239 // Multiple global sums
241
249 template<class T, class C>
250 inline typename UnaryReturn<C, FnSumMulti>::Type_t
251 sumMulti(const QDPType<T,C>& s1, const Set& ss)
252 {
253 return sumMulti(PETE_identity(s1), ss);
254 }
255
256
257 //-----------------------------------------------
258 // Global max and min
260
263 template<class T, class C>
264 inline typename UnaryReturn<C, FnGlobalMax>::Type_t
265 globalMax(const QDPType<T,C>& s1)
266 {
267 return globalMax(PETE_identity(s1));
268 }
269
270
272
275 template<class T, class C>
276 inline typename UnaryReturn<C, FnGlobalMin>::Type_t
277 globalMin(const QDPType<T,C>& s1)
278 {
279 return globalMin(PETE_identity(s1));
280 }
281
282
283 //-----------------------------------------------------------------------------
284 // These functions always return bool
286
289 template<class T, class C>
290 inline bool
291 isnan(const QDPExpr<T,C>& s1)
292 {
293 return isnan(C(s1));
294 }
295
297
300 template<class T, class C>
301 inline bool
302 isinf(const QDPExpr<T,C>& s1)
303 {
304 return isinf(C(s1));
305 }
306
308
311 template<class T, class C>
312 inline bool
313 isnormal(const QDPExpr<T,C>& s1)
314 {
315 return isnormal(C(s1));
316 }
317
319
322 template<class T, class C>
323 inline bool
324 isfinite(const QDPExpr<T,C>& s1)
325 {
326 return isfinite(C(s1));
327 }
328
329 // end of group3
331
332
338
339 //-----------------------------------------------
340 // Spin projection
342
343 template<class T, class C>
344 inline typename UnaryReturn<C, FnSpinProject>::Type_t
345 spinProject(const QDPType<T,C>& s1, int mu, int isign)
346 {
347 typedef typename UnaryReturn<C, FnSpinProject>::Type_t Ret_t;
348 Ret_t d;
349
350 switch (isign)
351 {
352 case -1:
353 switch (mu)
354 {
355 case 0:
356 d = spinProjectDir0Minus(s1);
357 break;
358 case 1:
359 d = spinProjectDir1Minus(s1);
360 break;
361 case 2:
362 d = spinProjectDir2Minus(s1);
363 break;
364 case 3:
365 d = spinProjectDir3Minus(s1);
366 break;
367 default:
368 std::cerr << "Spin project: illegal direction\n";
369 exit(1);
370 }
371 break;
372
373 case +1:
374 switch (mu)
375 {
376 case 0:
377 d = spinProjectDir0Plus(s1);
378 break;
379 case 1:
380 d = spinProjectDir1Plus(s1);
381 break;
382 case 2:
383 d = spinProjectDir2Plus(s1);
384 break;
385 case 3:
386 d = spinProjectDir3Plus(s1);
387 break;
388 default:
389 std::cerr << "Spin project: illegal direction\n";
390 exit(1);
391 }
392 break;
393
394 default:
395 std::cerr << "Spin project: isign must be pos or neg.\n";
396 exit(1);
397 }
398
399 return d;
400 }
401
402
403 //-----------------------------------------------
404 // Spin reconstruction
406
407 template<class T, class C>
408 inline typename UnaryReturn<C, FnSpinReconstruct>::Type_t
409 spinReconstruct(const QDPType<T,C>& s1, int mu, int isign)
410 {
411 // typedef typename UnaryReturn<C, FnSpinReconstruct>::Type_t Ret_t;
412 // Ret_t d;
413
414 switch (isign)
415 {
416 case -1:
417 switch (mu)
418 {
419 case 0:
420 return spinReconstructDir0Minus(s1);
421 case 1:
422 return spinReconstructDir1Minus(s1);
423 case 2:
424 return spinReconstructDir2Minus(s1);
425 case 3:
426 return spinReconstructDir3Minus(s1);
427 default:
428 std::cerr << "Spin reconstruct: illegal direction\n";
429 exit(1);
430 }
431 break;
432
433 case +1:
434 switch (mu)
435 {
436 case 0:
437 return spinReconstructDir0Plus(s1);
438 case 1:
439 return spinReconstructDir1Plus(s1);
440 case 2:
441 return spinReconstructDir2Plus(s1);
442 case 3:
443 return spinReconstructDir3Plus(s1);
444 default:
445 std::cerr << "Spin reconstruct: illegal direction\n";
446 exit(1);
447 }
448 break;
449
450 default:
451 std::cerr << "Spin reconstruct: isign must be pos or neg.\n";
452 exit(1);
453 }
454
455 // return d;
456 }
457
458 // end of group5
460
461} // namespace QDP
462
463#endif
QDP::QDPExpr
Expression class for QDP.
Definition qdp_qdpexpr.h:16
QDP::QDPType
QDPType - major type class/container for all QDP objects.
Definition qdp_qdptype.h:29
QDP::Set
Set - collection of subsets controlling which sites are involved in an operation.
Definition qdp_subset.h:96
QDP::Subset
Subsets - controls how lattices are looped.
Definition qdp_subset.h:39
QDP::innerProductReal
BinaryReturn< C1, C2, FnInnerProductReal >::Type_t innerProductReal(const QDPType< T1, C1 > &s1, const QDPType< T2, C2 > &s2)
OScalar = innerProductReal(adj(source1)*source2).
Definition qdp_globalfuncs.h:172
QDP::globalMin
UnaryReturn< C, FnGlobalMin >::Type_t globalMin(const QDPType< T, C > &s1)
OScalar = globalMin(source).
Definition qdp_globalfuncs.h:277
QDP::norm2
UnaryReturn< C, FnNorm2 >::Type_t norm2(const QDPType< T, C > &s1)
OScalar = norm2(trace(adj(source)*source)).
Definition qdp_globalfuncs.h:55
QDP::isnan
bool isnan(const QDPExpr< T, C > &s1)
bool = isnan(source)
Definition qdp_globalfuncs.h:291
QDP::isfinite
bool isfinite(const QDPExpr< T, C > &s1)
bool = isfinite(source)
Definition qdp_globalfuncs.h:324
QDP::innerProduct
BinaryReturn< C1, C2, FnInnerProduct >::Type_t innerProduct(const QDPType< T1, C1 > &s1, const QDPType< T2, C2 > &s2)
OScalar = innerProduct(adj(source1)*source2).
Definition qdp_globalfuncs.h:98
QDP::isnormal
bool isnormal(const QDPExpr< T, C > &s1)
bool = isnormal(source)
Definition qdp_globalfuncs.h:313
QDP::isinf
bool isinf(const QDPExpr< T, C > &s1)
bool = isinf(source)
Definition qdp_globalfuncs.h:302
QDP::globalMax
UnaryReturn< C, FnGlobalMax >::Type_t globalMax(const QDPType< T, C > &s1)
OScalar = globalMax(source).
Definition qdp_globalfuncs.h:265
QDP::sumMulti
UnaryReturn< C, FnSumMulti >::Type_t sumMulti(const QDPType< T, C > &s1, const Set &ss)
dest = sumMulti(source1,Set)
Definition qdp_globalfuncs.h:251
QDP::sum
UnaryReturn< C, FnSum >::Type_t sum(const QDPType< T, C > &s1)
OScalar = sum(source).
Definition qdp_globalfuncs.h:27
QDP::spinProject
UnaryReturn< C, FnSpinProject >::Type_t spinProject(const QDPType< T, C > &s1, int mu, int isign)
dest = spinProject(source1)
Definition qdp_globalfuncs.h:345
QDP::spinReconstruct
UnaryReturn< C, FnSpinReconstruct >::Type_t spinReconstruct(const QDPType< T, C > &s1, int mu, int isign)
dest = spinReconstruct(source1)
Definition qdp_globalfuncs.h:409
QDP::C
OLattice< PScalar< PColorMatrix< RComplexFloat, 3 > > > C
Definition qdp_scalarsite_generic_linalg.h:12
QDP
Yet another random number generator.
Definition qdp_bgq_threadbind.cc:16
QDP::spinProjectDir2Minus
MakeReturn< UnaryNode< FnSpinProjectDir2Minus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinProjectDir2Minus >::Type_t >::Expression_t spinProjectDir2Minus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5128
QDP::spinProjectDir1Minus
MakeReturn< UnaryNode< FnSpinProjectDir1Minus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinProjectDir1Minus >::Type_t >::Expression_t spinProjectDir1Minus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5115
QDP::spinProjectDir0Plus
MakeReturn< UnaryNode< FnSpinProjectDir0Plus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinProjectDir0Plus >::Type_t >::Expression_t spinProjectDir0Plus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5050
QDP::spinReconstructDir2Minus
MakeReturn< UnaryNode< FnSpinReconstructDir2Minus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinReconstructDir2Minus >::Type_t >::Expression_t spinReconstructDir2Minus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5232
QDP::localInnerProductReal
MakeReturn< BinaryNode< FnLocalInnerProductReal, typenameCreateLeaf< QDPType< T1, C1 > >::Leaf_t, typenameCreateLeaf< QDPExpr< T2, C2 > >::Leaf_t >, typenameBinaryReturn< C1, C2, FnLocalInnerProductReal >::Type_t >::Expression_t localInnerProductReal(const QDPType< T1, C1 > &l, const QDPExpr< T2, C2 > &r)
Definition qdp.h:2428
QDP::PETE_identity
MakeReturn< UnaryNode< OpIdentity, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, OpIdentity >::Type_t >::Expression_t PETE_identity(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5531
QDP::spinReconstructDir2Plus
MakeReturn< UnaryNode< FnSpinReconstructDir2Plus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinReconstructDir2Plus >::Type_t >::Expression_t spinReconstructDir2Plus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5180
QDP::spinProjectDir3Minus
MakeReturn< UnaryNode< FnSpinProjectDir3Minus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinProjectDir3Minus >::Type_t >::Expression_t spinProjectDir3Minus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5141
QDP::spinReconstructDir0Plus
MakeReturn< UnaryNode< FnSpinReconstructDir0Plus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinReconstructDir0Plus >::Type_t >::Expression_t spinReconstructDir0Plus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5154
QDP::spinReconstructDir1Plus
MakeReturn< UnaryNode< FnSpinReconstructDir1Plus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinReconstructDir1Plus >::Type_t >::Expression_t spinReconstructDir1Plus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5167
QDP::spinReconstructDir1Minus
MakeReturn< UnaryNode< FnSpinReconstructDir1Minus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinReconstructDir1Minus >::Type_t >::Expression_t spinReconstructDir1Minus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5219
QDP::spinProjectDir2Plus
MakeReturn< UnaryNode< FnSpinProjectDir2Plus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinProjectDir2Plus >::Type_t >::Expression_t spinProjectDir2Plus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5076
QDP::localInnerProduct
MakeReturn< BinaryNode< FnLocalInnerProduct, typenameCreateLeaf< QDPType< T1, C1 > >::Leaf_t, typenameCreateLeaf< QDPExpr< T2, C2 > >::Leaf_t >, typenameBinaryReturn< C1, C2, FnLocalInnerProduct >::Type_t >::Expression_t localInnerProduct(const QDPType< T1, C1 > &l, const QDPExpr< T2, C2 > &r)
Definition qdp.h:2412
QDP::spinProjectDir3Plus
MakeReturn< UnaryNode< FnSpinProjectDir3Plus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinProjectDir3Plus >::Type_t >::Expression_t spinProjectDir3Plus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5089
QDP::spinReconstructDir3Minus
MakeReturn< UnaryNode< FnSpinReconstructDir3Minus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinReconstructDir3Minus >::Type_t >::Expression_t spinReconstructDir3Minus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5245
QDP::spinReconstructDir0Minus
MakeReturn< UnaryNode< FnSpinReconstructDir0Minus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinReconstructDir0Minus >::Type_t >::Expression_t spinReconstructDir0Minus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5206
QDP::spinProjectDir0Minus
MakeReturn< UnaryNode< FnSpinProjectDir0Minus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinProjectDir0Minus >::Type_t >::Expression_t spinProjectDir0Minus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5102
QDP::localNorm2
MakeReturn< UnaryNode< FnLocalNorm2, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnLocalNorm2 >::Type_t >::Expression_t localNorm2(const QDPExpr< T1, C1 > &l)
Definition qdp.h:4998
QDP::spinReconstructDir3Plus
MakeReturn< UnaryNode< FnSpinReconstructDir3Plus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinReconstructDir3Plus >::Type_t >::Expression_t spinReconstructDir3Plus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5193
QDP::spinProjectDir1Plus
MakeReturn< UnaryNode< FnSpinProjectDir1Plus, typenameCreateLeaf< QDPExpr< T1, C1 > >::Leaf_t >, typenameUnaryReturn< C1, FnSpinProjectDir1Plus >::Type_t >::Expression_t spinProjectDir1Plus(const QDPExpr< T1, C1 > &l)
Definition qdp.h:5063
BinaryReturn::Type_t
Promote< T1, T2 >::Type_t Type_t
Definition TypeComputations.h:376
Generated by doxygen 1.16.1