QDP++
qdp_scalarsite_pabasm.h
Go to the documentation of this file.
1// -*- C++ -*-
8
9#ifndef QDP_SCALARSITE_GENERIC_H
10#define QDP_SCALARSITE_GENERIC_H
11
12
13#include "qdp_diagnostics.h"
14
15QDPXX_MESSAGE("Using PABASM Scalarsite")
16
17namespace QDP {
18
25
26// Use this def just to safe some typing later on in the file
27#define RComplexFloat RComplex<float>
28
29
30#include "scalarsite_generic/generic_mult_nn.h"
31#include "scalarsite_generic/generic_mult_na.h"
32#include "scalarsite_generic/generic_mult_an.h"
33#include "scalarsite_generic/generic_mult_aa.h"
34#include "scalarsite_generic/generic_mat_vec.h"
35#include "scalarsite_generic/generic_adj_mat_vec.h"
36#include "scalarsite_generic/generic_addvec.h"
37
38
39// #define QDP_SCALARSITE_DEBUG
40
41// Optimized version of
42// PColorMatrix<RComplexFloat,3> <- PColorMatrix<RComplexFloat,3> * PColorMatrix<RComplexFloat,3>
43template<>
44inline BinaryReturn<PMatrix<RComplexFloat,3,PColorMatrix>,
45 PMatrix<RComplexFloat,3,PColorMatrix>, OpMultiply>::Type_t
46operator*(const PMatrix<RComplexFloat,3,PColorMatrix>& l,
47 const PMatrix<RComplexFloat,3,PColorMatrix>& r)
48{
49 BinaryReturn<PMatrix<RComplexFloat,3,PColorMatrix>,
50 PMatrix<RComplexFloat,3,PColorMatrix>, OpMultiply>::Type_t d;
51
52#if defined(QDP_SCALARSITE_DEBUG)
53 cout << "M*M" << endl;
54#endif
55
56 _inline_generic_mult_su3_nn(l,r,d);
57
58 return d;
59}
60
61
62// Optimized version of
63// PScalar<PColorMatrix<RComplexFloat,3>> <- PScalar<PColorMatrix<RComplexFloat,3>> *
64// PScalar<PColorMatrix<RComplexFloat,3>>
65template<>
66inline BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
67 PScalar<PColorMatrix<RComplexFloat,3> >, OpMultiply>::Type_t
68operator*(const PScalar<PColorMatrix<RComplexFloat,3> >& l,
69 const PScalar<PColorMatrix<RComplexFloat,3> >& r)
70{
71 BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
72 PScalar<PColorMatrix<RComplexFloat,3> >, OpMultiply>::Type_t d;
73
74#if defined(QDP_SCALARSITE_DEBUG)
75 cout << "PSc<M>*PSc<M>" << endl;
76#endif
77
78 _inline_generic_mult_su3_nn(l.elem(),r.elem(),d.elem());
79
80 return d;
81}
82
83
84// Optimized version of
85// PColorMatrix<RComplexFloat,3> <- adj(PColorMatrix<RComplexFloat,3>) * PColorMatrix<RComplexFloat,3>
86template<>
87inline BinaryReturn<PMatrix<RComplexFloat,3,PColorMatrix>,
88 PMatrix<RComplexFloat,3,PColorMatrix>, OpAdjMultiply>::Type_t
89adjMultiply(const PMatrix<RComplexFloat,3,PColorMatrix>& l,
90 const PMatrix<RComplexFloat,3,PColorMatrix>& r)
91{
92 BinaryReturn<PMatrix<RComplexFloat,3,PColorMatrix>,
93 PMatrix<RComplexFloat,3,PColorMatrix>, OpAdjMultiply>::Type_t d;
94
95#if defined(QDP_SCALARSITE_DEBUG)
96 cout << "adj(M)*M" << endl;
97#endif
98
99 _inline_generic_mult_su3_an(l,r,d);
100
101 return d;
102}
103
104
105// Optimized version of
106// PScalar<PColorMatrix<RComplexFloat,3>> <- adj(PScalar<PColorMatrix<RComplexFloat,3>>) * PScalar<PColorMatrix<RComplexFloat,3>>
107template<>
108inline BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
109 PScalar<PColorMatrix<RComplexFloat,3> >, OpAdjMultiply>::Type_t
110adjMultiply(const PScalar<PColorMatrix<RComplexFloat,3> >& l,
111 const PScalar<PColorMatrix<RComplexFloat,3> >& r)
112{
113 BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
114 PScalar<PColorMatrix<RComplexFloat,3> >, OpAdjMultiply>::Type_t d;
115
116#if defined(QDP_SCALARSITE_DEBUG)
117 cout << "adj(PSc<M>)*PSc<M>" << endl;
118#endif
119
120 _inline_generic_mult_su3_an(l.elem(),r.elem(),d.elem());
121
122 return d;
123}
124
125
126// Optimized version of
127// PColorMatrix<RComplexFloat,3> <- PColorMatrix<RComplexFloat,3> * adj(PColorMatrix<RComplexFloat,3>)
128template<>
129inline BinaryReturn<PMatrix<RComplexFloat,3,PColorMatrix>,
130 PMatrix<RComplexFloat,3,PColorMatrix>, OpMultiplyAdj>::Type_t
131multiplyAdj(const PMatrix<RComplexFloat,3,PColorMatrix>& l,
132 const PMatrix<RComplexFloat,3,PColorMatrix>& r)
133{
134 BinaryReturn<PMatrix<RComplexFloat,3,PColorMatrix>,
135 PMatrix<RComplexFloat,3,PColorMatrix>, OpMultiplyAdj>::Type_t d;
136
137#if defined(QDP_SCALARSITE_DEBUG)
138 cout << "M*adj(M)" << endl;
139#endif
140
141 _inline_generic_mult_su3_na(l,r,d);
142
143 return d;
144}
145
146
147// Optimized version of
148// PScalar<PColorMatrix<RComplexFloat,3>> <- PScalar<PColorMatrix<RComplexFloat,3>> * adj(PScalar<PColorMatrix<RComplexFloat,3>>)
149template<>
150inline BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
151 PScalar<PColorMatrix<RComplexFloat,3> >, OpMultiplyAdj>::Type_t
152multiplyAdj(const PScalar<PColorMatrix<RComplexFloat,3> >& l,
153 const PScalar<PColorMatrix<RComplexFloat,3> >& r)
154{
155 BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
156 PScalar<PColorMatrix<RComplexFloat,3> >, OpMultiplyAdj>::Type_t d;
157
158#if defined(QDP_SCALARSITE_DEBUG)
159 cout << "PSc<M>*adj(PSc<M>)" << endl;
160#endif
161
162 _inline_generic_mult_su3_na(l.elem(),r.elem(),d.elem());
163
164 return d;
165}
166
167
168#if 0
169// Do not use this routine - for some unknown reason on P4 it is slower than
170// the default code!!
171
172// Optimized version of
173// PColorMatrix<RComplexFloat,3> <- adj(PColorMatrix<RComplexFloat,3>) * adj(PColorMatrix<RComplexFloat,3>)
174template<>
175inline BinaryReturn<PMatrix<RComplexFloat,3,PColorMatrix>,
176 PMatrix<RComplexFloat,3,PColorMatrix>, OpAdjMultiplyAdj>::Type_t
177adjMultiplyAdj(const PMatrix<RComplexFloat,3,PColorMatrix>& l,
178 const PMatrix<RComplexFloat,3,PColorMatrix>& r)
179{
180 BinaryReturn<PMatrix<RComplexFloat,3,PColorMatrix>,
181 PMatrix<RComplexFloat,3,PColorMatrix>, OpAdjMultiplyAdj>::Type_t d;
182
183#if defined(QDP_SCALARSITE_DEBUG)
184 cout << "adj(PSc<M>)*adj(PSc<M>)" << endl;
185#endif
186
187 _inline_generic_mult_su3_aa(l,r,d);
188
189 return d;
190}
191#endif
192
193
194// Optimized version of
195// PColorVector<RComplexFloat,3> <- PColorMatrix<RComplexFloat,3> * PColorVector<RComplexFloat,3>
196template<>
197inline BinaryReturn<PMatrix<RComplexFloat,3,PColorMatrix>,
198 PVector<RComplexFloat,3,PColorVector>, OpMultiply>::Type_t
199operator*(const PMatrix<RComplexFloat,3,PColorMatrix>& l,
200 const PVector<RComplexFloat,3,PColorVector>& r)
201{
202 BinaryReturn<PMatrix<RComplexFloat,3,PColorMatrix>,
203 PVector<RComplexFloat,3,PColorVector>, OpMultiply>::Type_t d;
204
205#if defined(QDP_SCALARSITE_DEBUG)
206 cout << "M*V" << endl;
207#endif
208
209 _inline_generic_mult_su3_mat_vec(l,r,d);
210
211 return d;
212}
213
214
215// Optimized version of
216// PScalar<PColorVector<RComplexFloat,3>> <- PScalar<PColorMatrix<RComplexFloat,3>> * PScalar<PColorVector<RComplexFloat,3>>
217template<>
218inline BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
219 PScalar<PColorVector<RComplexFloat,3> >, OpMultiply>::Type_t
220operator*(const PScalar<PColorMatrix<RComplexFloat,3> >& l,
221 const PScalar<PColorVector<RComplexFloat,3> >& r)
222{
223 BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
224 PScalar<PColorVector<RComplexFloat,3> >, OpMultiply>::Type_t d;
225
226#if defined(QDP_SCALARSITE_DEBUG)
227 cout << "PSc<M>*PSc<V>" << endl;
228#endif
229
230 _inline_generic_mult_su3_mat_vec(l.elem(),r.elem(),d.elem());
231
232 return d;
233}
234
235
236// Optimized version of
237// PColorVector<RComplexFloat,3> <- adj(PColorMatrix<RComplexFloat,3>) * PColorVector<RComplexFloat,3>
238template<>
239inline BinaryReturn<PMatrix<RComplexFloat,3,PColorMatrix>,
240 PVector<RComplexFloat,3,PColorVector>, OpAdjMultiply>::Type_t
241adjMultiply(const PMatrix<RComplexFloat,3,PColorMatrix>& l,
242 const PVector<RComplexFloat,3,PColorVector>& r)
243{
244 BinaryReturn<PMatrix<RComplexFloat,3,PColorMatrix>,
245 PVector<RComplexFloat,3,PColorVector>, OpAdjMultiply>::Type_t d;
246
247#if defined(QDP_SCALARSITE_DEBUG)
248 cout << "adj(M)*V" << endl;
249#endif
250
251 _inline_generic_mult_adj_su3_mat_vec(l,r,d);
252
253 return d;
254}
255
256
257// Optimized version of StaggeredFermion <- ColorMatrix*StaggeredFermion
258// PSpinVector<PColorVector<RComplexFloat,3>,1> <- PScalar<PColorMatrix<RComplexFloat,3>> * PSpinVector<PColorVector<RComplexFloat,3>,1>
259template<>
260inline BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
261 PVector<PColorVector<RComplexFloat,3>,1,PSpinVector>, OpMultiply>::Type_t
262operator*(const PScalar<PColorMatrix<RComplexFloat,3> >& l,
263 const PVector<PColorVector<RComplexFloat,3>,1,PSpinVector>& r)
264{
265 BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
266 PVector<PColorVector<RComplexFloat,3>,1,PSpinVector>, OpMultiply>::Type_t d;
267
268#if defined(QDP_SCALARSITE_DEBUG)
269 cout << "PSc<M>*S" << endl;
270#endif
271
272 _inline_generic_mult_su3_mat_vec(l.elem(),r.elem(0),d.elem(0));
273
274 return d;
275}
276
277
278// Optimized version of
279// PScalar<PColorVector<RComplexFloat,3>> <- adj(PScalar<PColorMatrix<RComplexFloat,3>>) * PScalar<PColorVector<RComplexFloat,3>>
280template<>
281inline BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
282 PScalar<PColorVector<RComplexFloat,3> >, OpAdjMultiply>::Type_t
283adjMultiply(const PScalar<PColorMatrix<RComplexFloat,3> >& l,
284 const PScalar<PColorVector<RComplexFloat,3> >& r)
285{
286 BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
287 PScalar<PColorVector<RComplexFloat,3> >, OpAdjMultiply>::Type_t d;
288
289#if defined(QDP_SCALARSITE_DEBUG)
290 cout << "adj(PSc<M>)*PSc<V>" << endl;
291#endif
292
293 _inline_generic_mult_adj_su3_mat_vec(l.elem(),r.elem(),d.elem());
294
295 return d;
296}
297
298
299// Optimized version of StaggeredFermion <- adj(ColorMatrix)*StaggeredFermion
300// PSpinVector<PColorVector<RComplexFloat,3>,1> <- adj(PScalar<PColorMatrix<RComplexFloat,3>>) * PSpinVector<PColorVector<RComplexFloat,3>,1>
301template<>
302inline BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
303 PVector<PColorVector<RComplexFloat,3>,1,PSpinVector>, OpAdjMultiply>::Type_t
304adjMultiply(const PScalar<PColorMatrix<RComplexFloat,3> >& l,
305 const PVector<PColorVector<RComplexFloat,3>,1,PSpinVector>& r)
306{
307 BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
308 PVector<PColorVector<RComplexFloat,3>,1,PSpinVector>, OpAdjMultiply>::Type_t d;
309
310#if defined(QDP_SCALARSITE_DEBUG)
311 cout << "adj(PSc<M>)*S" << endl;
312#endif
313
314 _inline_generic_mult_adj_su3_mat_vec(l.elem(),r.elem(0),d.elem(0));
315
316 return d;
317}
318
319
320// Optimized version of HalfFermion <- ColorMatrix*HalfFermion
321// PSpinVector<PColorVector<RComplexFloat,3>,2> <- PScalar<PColorMatrix<RComplexFloat,3>> *
322// PSpinVector<ColorVector<RComplexFloat,3>,2>
323template<>
324inline BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
325 PVector<PColorVector<RComplexFloat,3>,2,PSpinVector>, OpMultiply>::Type_t
326operator*(const PScalar<PColorMatrix<RComplexFloat,3> >& l,
327 const PVector<PColorVector<RComplexFloat,3>,2,PSpinVector>& r)
328{
329 BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
330 PVector<PColorVector<RComplexFloat,3>,2,PSpinVector>, OpMultiply>::Type_t d;
331
332#if defined(QDP_SCALARSITE_DEBUG)
333 cout << "PSc<M>*H" << endl;
334#endif
335
336 _inline_generic_mult_su3_mat_vec(l.elem(),r.elem(0),d.elem(0));
337 _inline_generic_mult_su3_mat_vec(l.elem(),r.elem(1),d.elem(1));
338
339 return d;
340}
341
342
343// Optimized version of HalfFermion <- ColorMatrix*HalfFermion
344// PSpinVector<PColorVector<RComplexFloat,3>,2> <- adj(PScalar<PColorMatrix<RComplexFloat,3>>) *
345// PSpinVector<ColorVector<RComplexFloat,3>,2>
346template<>
347inline BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
348 PVector<PColorVector<RComplexFloat,3>,2,PSpinVector>, OpAdjMultiply>::Type_t
349adjMultiply(const PScalar<PColorMatrix<RComplexFloat,3> >& l,
350 const PVector<PColorVector<RComplexFloat,3>,2,PSpinVector>& r)
351{
352 BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
353 PVector<PColorVector<RComplexFloat,3>,2,PSpinVector>, OpAdjMultiply>::Type_t d;
354
355#if defined(QDP_SCALARSITE_DEBUG)
356 cout << "adj(PSc<M>)*H" << endl;
357#endif
358
359 _inline_generic_mult_adj_su3_mat_vec(l.elem(),r.elem(0),d.elem(0));
360 _inline_generic_mult_adj_su3_mat_vec(l.elem(),r.elem(1),d.elem(1));
361
362 return d;
363}
364
365
366// Optimized version of
367// PColorVector<RComplexFloat,3> <- PColorVector<RComplexFloat,3> + PColorVector<RComplexFloat,3>
368template<>
369inline BinaryReturn<PVector<RComplexFloat,3,PColorVector>,
370 PVector<RComplexFloat,3,PColorVector>, OpAdd>::Type_t
371operator+(const PVector<RComplexFloat,3,PColorVector>& l,
372 const PVector<RComplexFloat,3,PColorVector>& r)
373{
374 BinaryReturn<PVector<RComplexFloat,3,PColorVector>,
375 PVector<RComplexFloat,3,PColorVector>, OpAdd>::Type_t d;
376
377#if defined(QDP_SCALARSITE_DEBUG)
378 cout << "V+V" << endl;
379#endif
380
381 _inline_generic_add_su3_vector(l,r,d);
382
383 return d;
384}
385
386
387// Optimized version of DiracFermion <- ColorMatrix*DiracFermion
388// PSpinVector<PColorVector<RComplexFloat,3>,4> <- PScalar<PColorMatrix<RComplexFloat,3>>
389// * PSpinVector<PColorVector<RComplexFloat,3>,4>
390template<>
391inline BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
392 PVector<PColorVector<RComplexFloat,3>,4,PSpinVector>, OpMultiply>::Type_t
393operator*(const PScalar<PColorMatrix<RComplexFloat,3> >& l,
394 const PVector<PColorVector<RComplexFloat,3>,4,PSpinVector>& r)
395{
396 BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
397 PVector<PColorVector<RComplexFloat,3>,4,PSpinVector>, OpMultiply>::Type_t d;
398
399#if defined(QDP_SCALARSITE_DEBUG)
400 cout << "PSc<M>*D" << endl;
401#endif
402
403 _inline_generic_mult_su3_mat_vec(l.elem(),r.elem(0),d.elem(0));
404 _inline_generic_mult_su3_mat_vec(l.elem(),r.elem(1),d.elem(1));
405 _inline_generic_mult_su3_mat_vec(l.elem(),r.elem(2),d.elem(2));
406 _inline_generic_mult_su3_mat_vec(l.elem(),r.elem(3),d.elem(3));
407
408 return d;
409}
410
411
412// Optimized version of DiracFermion <- adj(ColorMatrix)*DiracFermion
413// PSpinVector<PColorVector<RComplexFloat,3>,4> <- adj(PScalar<PColorMatrix<RComplexFloat,3>>)
414// * PSpinVector<PColorVector<RComplexFloat,3>,4>
415template<>
416inline BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
417 PVector<PColorVector<RComplexFloat,3>,4,PSpinVector>, OpAdjMultiply>::Type_t
418adjMultiply(const PScalar<PColorMatrix<RComplexFloat,3> >& l,
419 const PVector<PColorVector<RComplexFloat,3>,4,PSpinVector>& r)
420{
421 BinaryReturn<PScalar<PColorMatrix<RComplexFloat,3> >,
422 PVector<PColorVector<RComplexFloat,3>,4,PSpinVector>, OpAdjMultiply>::Type_t d;
423
424#if defined(QDP_SCALARSITE_DEBUG)
425 cout << "adj(PSc<M>)*D" << endl;
426#endif
427
428 _inline_generic_mult_adj_su3_mat_vec(l.elem(),r.elem(0),d.elem(0));
429 _inline_generic_mult_adj_su3_mat_vec(l.elem(),r.elem(1),d.elem(1));
430 _inline_generic_mult_adj_su3_mat_vec(l.elem(),r.elem(2),d.elem(2));
431 _inline_generic_mult_adj_su3_mat_vec(l.elem(),r.elem(3),d.elem(3));
432
433 return d;
434}
435
436
437
438// Optimized version of
439// PScalar<PColorVector<RComplexFloat,3>> <- PScalar<PColorVector<RComplexFloat,3>>
440// + PScalar<PColorVector<RComplexFloat,3>>
441template<>
442inline BinaryReturn<PScalar<PColorVector<RComplexFloat,3> >,
443 PScalar<PColorVector<RComplexFloat,3> >, OpAdd>::Type_t
444operator+(const PScalar<PColorVector<RComplexFloat,3> >& l,
445 const PScalar<PColorVector<RComplexFloat,3> >& r)
446{
447 BinaryReturn<PScalar<PColorVector<RComplexFloat,3> >,
448 PScalar<PColorVector<RComplexFloat,3> >, OpAdd>::Type_t d;
449
450#if defined(QDP_SCALARSITE_DEBUG)
451 cout << "PSc<V>+PSc<V>" << endl;
452#endif
453
454 _inline_generic_add_su3_vector(l.elem(),r.elem(),d.elem());
455
456 return d;
457}
458
459
460#if 1
461// Specialization to optimize the case
462// LatticeHalfFermion = LatticeColorMatrix * LatticeHalfFermion
463// NOTE: let this be a subroutine to save space
464template<>
465inline
466void evaluate(OLattice<PSpinVector<PColorVector<RComplexFloat, 3>, 2> >& d,
467 const OpAssign& op,
468 const QDPExpr<BinaryNode<OpMultiply,
469 Reference<QDPType<PScalar<PColorMatrix<RComplexFloat, 3> >,
470 OLattice<PScalar<PColorMatrix<RComplexFloat, 3> > > > >,
471 Reference<QDPType<PSpinVector<PColorVector<RComplexFloat, 3>, 2>,
472 OLattice<PSpinVector<PColorVector<RComplexFloat, 3>, 2> > > > >,
473 OLattice<PSpinVector<PColorVector<RComplexFloat, 3>, 2> > >& rhs,
474 const Subset& s)
475{
476#if defined(QDP_SCALARSITE_DEBUG)
477 cout << "specialized QDP_H_M_times_H" << endl;
478#endif
479
480 typedef OLattice<PScalar<PColorMatrix<RComplexFloat, 3> > > C;
481 typedef OLattice<PSpinVector<PColorVector<RComplexFloat, 3>, 2> > H;
482
483 const C& l = static_cast<const C&>(rhs.expression().left());
484 const H& r = static_cast<const H&>(rhs.expression().right());
485
486 const int *tab = s.siteTable().slice();
487
488 for(int j=0; j <= s.numSiteTable(); ++j)
489 {
490 int i = tab[j];
491
492 _inline_generic_mult_su3_mat_vec(l.elem(i).elem(),
493 r.elem(i).elem(0),
494 d.elem(i).elem(0));
495 _inline_generic_mult_su3_mat_vec(l.elem(i).elem(),
496 r.elem(i).elem(1),
497 d.elem(i).elem(1));
498 }
499}
500
501#endif
502 // end of group optimizations
504
505#if defined(QDP_SCALARSITE_DEBUG)
506#undef QDP_SCALARSITE_DEBUG
507#endif
508
509} // namespace QDP
510
511#endif
operator*
MakeReturn< BinaryNode< OpMultiply, typenameCreateLeaf< QDPType< T1, C1 > >::Leaf_t, typenameCreateLeaf< QDPExpr< T2, C2 > >::Leaf_t >, typenameBinaryReturn< C1, C2, OpMultiply >::Type_t >::Expression_t operator*(const QDPType< T1, C1 > &l, const QDPExpr< T2, C2 > &r)
Definition QDPOperators.h:2587
operator+
MakeReturn< BinaryNode< OpAdd, typenameCreateLeaf< QDPType< T1, C1 > >::Leaf_t, typenameCreateLeaf< QDPExpr< T2, C2 > >::Leaf_t >, typenameBinaryReturn< C1, C2, OpAdd >::Type_t >::Expression_t operator+(const QDPType< T1, C1 > &l, const QDPExpr< T2, C2 > &r)
Definition QDPOperators.h:2555
QDP::OLattice
Outer grid Lattice type.
Definition qdp_outer.h:264
QDP::OLattice::elem
T & elem(int i)
Definition qdp_outer.h:400
QDP::PColorMatrix
Primitive color Matrix class.
Definition qdp_primcolormat.h:26
QDP::PColorVector
Primitive color Vector class.
Definition qdp_primcolorvec.h:24
QDP::PMatrix
Primitive Matrix class.
Definition qdp_primmatrix.h:31
QDP::PScalar
Primitive Scalar.
Definition qdp_primscalar.h:25
QDP::PSpinVector
Primitive spin Vector class.
Definition qdp_primspinvec.h:33
QDP::PVector
Primitive Vector class.
Definition qdp_primvector.h:31
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::Subset
Subsets - controls how lattices are looped.
Definition qdp_subset.h:39
QDP::Subset::siteTable
const multi1d< int > & siteTable() const
Definition qdp_subset.h:83
QDP::Subset::numSiteTable
int numSiteTable() const
Definition qdp_subset.h:84
QDP::multi1d::slice
const T * slice() const
Return ref to a column slice.
Definition qdp_multi.h:225
generic_addvec.h
generic_adj_mat_vec.h
generic_mat_vec.h
generic_mult_aa.h
generic_mult_an.h
generic_mult_na.h
generic_mult_nn.h
QDP::multiplyAdj
BinaryReturn< IScalar< T1 >, IScalar< T2 >, OpMultiplyAdj >::Type_t multiplyAdj(const IScalar< T1 > &l, const IScalar< T2 > &r)
Definition qdp_inner.h:1189
QDP::adjMultiply
BinaryReturn< IScalar< T1 >, IScalar< T2 >, OpAdjMultiply >::Type_t adjMultiply(const IScalar< T1 > &l, const IScalar< T2 > &r)
Definition qdp_inner.h:1180
QDP::H
OLattice< PSpinVector< PColorVector< RComplexFloat, 3 >, 2 > > H
Definition qdp_scalarsite_generic_linalg.h:13
QDP::evaluate
void evaluate(OLattice< DCol > &d, const OpAssign &op, const QDPExpr< BinaryNode< OpMultiply, Reference< QDPType< DCol, OLattice< DCol > > >, Reference< QDPType< DCol, OLattice< DCol > > > >, OLattice< DCol > > &rhs, const Subset &s)
Definition qdp_scalarsite_linalg_double.cc:78
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_diagnostics.h
QDPXX_MESSAGE
#define QDPXX_MESSAGE(s)
Definition qdp_diagnostics.h:20
_inline_generic_mult_adj_su3_mat_vec
#define _inline_generic_mult_adj_su3_mat_vec(aa, bb, cc)
_inline_generic_mult_su3_an
#define _inline_generic_mult_su3_an(aa, bb, cc)
_inline_generic_mult_su3_aa
#define _inline_generic_mult_su3_aa(aa, bb, cc)
_inline_generic_mult_su3_na
#define _inline_generic_mult_su3_na(aa, bb, cc)
_inline_generic_add_su3_vector
#define _inline_generic_add_su3_vector(aa, bb, cc)
_inline_generic_mult_su3_nn
#define _inline_generic_mult_su3_nn(aa, bb, cc)
_inline_generic_mult_su3_mat_vec
#define _inline_generic_mult_su3_mat_vec(aa, bb, cc)
Generated by doxygen 1.16.1