QDP++
generic_fused_spin_recon.h
Go to the documentation of this file.
1#ifndef GENERIC_FUSED_SPIN_RECON_H
2#define GENERIC_FUSED_SPIN_RECON_H
3
4namespace QDP {
5
6// Convenience Types
7typedef PColorVector<RComplex<REAL>, 3> ColVec;
8typedef PSpinVector< ColVec, 4> Spin4;
9typedef PSpinVector< ColVec, 2> Spin2;
10typedef PColorMatrix<RComplex<REAL>, 3> ColMat;
11
12
13// Tell the compiler what the result of the multiply
14// half spinor with gauge field operations are
15// OScalar
16template<>
23
24// Tell the compiler what the result of opreations are
25// For OLattice
26template<>
33
37struct FnSReconDir0MinusProd
38{
39 PETE_EMPTY_CONSTRUCTORS(FnSReconDir0MinusProd)
40
41 template<class T1, class T2, int N>
42 inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N>, FnSReconDir0MinusProd>::Type_t
43 operator()(const PScalar<T1>& a, const PSpinVector<T2,N>& b) const
44 {
45 return sreconDir0Minus(a, b);
46 }
47};
48
49
50// Proclaim that the return type of the fused op is the same as the return type
51// of the SpinReconstruction on a normal half vector
52template<class T1, class T2, int N>
57
58// For Generic Subtypes of OLattice<> and OScalar<>
59// Do the generic thing
60template<class T1, class T2, int N>
61inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N> , FnSReconDir0MinusProd>::Type_t
69
70// This rewrites spinReconstructDir0Minus( Matrix * Vec )
71// as a fused reconstruct: FnSReconDir0MinusProd< Matrix, Vec >
72template<class T1, class T2, template<class> class C, int N >
73inline typename MakeReturn<BinaryNode< FnSReconDir0MinusProd,
74 typename CreateLeaf< QDPType< PScalar<T1> , C<PScalar<T1> > > >::Leaf_t,
75 typename CreateLeaf< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >::Leaf_t>,
76 typename BinaryReturn< C<PScalar<T1> >, C<PSpinVector<T2,N> >, FnSReconDir0MinusProd>::Type_t >::Expression_t
102
106struct FnSReconDir0PlusProd
107{
108 PETE_EMPTY_CONSTRUCTORS(FnSReconDir0PlusProd)
109
110 template<class T1, class T2, int N>
111 inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N>, FnSReconDir0PlusProd>::Type_t
112 operator()(const PScalar<T1>& a, const PSpinVector<T2,N>& b) const
113 {
114 return sreconDir0Plus(a, b);
115 }
116};
117
118
119// Proclaim that the return type of the fused op is the same as the return type
120// of the SpinReconstruction on a normal half vector
121template<class T1, class T2, int N>
126
127// For Generic Subtypes of OLattice<> and OScalar<>
128// Do the generic thing
129template<class T1, class T2, int N>
130inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N> , FnSReconDir0PlusProd>::Type_t
138
139// This rewrites spinReconstructDir0Plus( Matrix * Vec )
140// as a fused reconstruct: FnSReconDir0PlusProd< Matrix, Vec >
141template<class T1, class T2, template<class> class C, int N >
142inline typename MakeReturn<BinaryNode< FnSReconDir0PlusProd,
143 typename CreateLeaf< QDPType< PScalar<T1> , C<PScalar<T1> > > >::Leaf_t,
144 typename CreateLeaf< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >::Leaf_t>,
145 typename BinaryReturn< C<PScalar<T1> >, C<PSpinVector<T2,N> >, FnSReconDir0PlusProd>::Type_t >::Expression_t
146spinReconstructDir0Plus(const QDPExpr<
147 BinaryNode< OpMultiply,
148 Reference< QDPType< PScalar<T1> , C<PScalar<T1> > > >,
149 Reference< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >
150 >,
151 typename BinaryReturn< QDPType< PScalar<T1> , C<PScalar<T1> > >,
152 QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > >,
153 OpMultiply >::Type_t >& l)
154{
155
156
157 typedef BinaryNode< FnSReconDir0PlusProd,
158 typename CreateLeaf< QDPType< PScalar<T1> , C<PScalar<T1> > > >::Leaf_t,
159 typename CreateLeaf< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >::Leaf_t> Tree_t;
160
161
162 return MakeReturn<Tree_t,typename BinaryReturn< C<PScalar<T1> >, C<PSpinVector<T2,N> >, FnSReconDir0PlusProd>::Type_t>::make(
163 Tree_t(
164 CreateLeaf<QDPType<PScalar<T1> ,C<PScalar<T1> > > >::make(l.expression().left()),
165 CreateLeaf<QDPType<PSpinVector<T2,N>,C<PSpinVector<T2,N> > > >::make(l.expression().right())
166 )
167 );
168
169
170}
171
175struct FnSReconDir1MinusProd
176{
177 PETE_EMPTY_CONSTRUCTORS(FnSReconDir1MinusProd)
178
179 template<class T1, class T2, int N>
180 inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N>, FnSReconDir1MinusProd>::Type_t
181 operator()(const PScalar<T1>& a, const PSpinVector<T2,N>& b) const
182 {
183 return sreconDir1Minus(a, b);
184 }
185};
186
187
188// Proclaim that the return type of the fused op is the same as the return type
189// of the SpinReconstruction on a normal half vector
190template<class T1, class T2, int N>
195
196// For Generic Subtypes of OLattice<> and OScalar<>
197// Do the generic thing
198template<class T1, class T2, int N>
199inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N> , FnSReconDir1MinusProd>::Type_t
207
208// This rewrites spinReconstructDir1Minus( Matrix * Vec )
209// as a fused reconstruct: FnSReconDir1MinusProd< Matrix, Vec >
210template<class T1, class T2, template<class> class C, int N >
211inline typename MakeReturn<BinaryNode< FnSReconDir1MinusProd,
212 typename CreateLeaf< QDPType< PScalar<T1> , C<PScalar<T1> > > >::Leaf_t,
213 typename CreateLeaf< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >::Leaf_t>,
214 typename BinaryReturn< C<PScalar<T1> >, C<PSpinVector<T2,N> >, FnSReconDir1MinusProd>::Type_t >::Expression_t
215spinReconstructDir1Minus(const QDPExpr<
216 BinaryNode< OpMultiply,
217 Reference< QDPType< PScalar<T1> , C<PScalar<T1> > > >,
218 Reference< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >
219 >,
220 typename BinaryReturn< QDPType< PScalar<T1> , C<PScalar<T1> > >,
221 QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > >,
222 OpMultiply >::Type_t >& l)
223{
224
225
226 typedef BinaryNode< FnSReconDir1MinusProd,
227 typename CreateLeaf< QDPType< PScalar<T1> , C<PScalar<T1> > > >::Leaf_t,
228 typename CreateLeaf< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >::Leaf_t> Tree_t;
229
230
231 return MakeReturn<Tree_t,typename BinaryReturn< C<PScalar<T1> >, C<PSpinVector<T2,N> >, FnSReconDir1MinusProd>::Type_t>::make(
232 Tree_t(
233 CreateLeaf<QDPType<PScalar<T1> ,C<PScalar<T1> > > >::make(l.expression().left()),
234 CreateLeaf<QDPType<PSpinVector<T2,N>,C<PSpinVector<T2,N> > > >::make(l.expression().right())
235 )
236 );
237
238
239}
240
244struct FnSReconDir1PlusProd
245{
246 PETE_EMPTY_CONSTRUCTORS(FnSReconDir1PlusProd)
247
248 template<class T1, class T2, int N>
249 inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N>, FnSReconDir1PlusProd>::Type_t
250 operator()(const PScalar<T1>& a, const PSpinVector<T2,N>& b) const
251 {
252 return sreconDir1Plus(a, b);
253 }
254};
255
256
257// Proclaim that the return type of the fused op is the same as the return type
258// of the SpinReconstruction on a normal half vector
259template<class T1, class T2, int N>
264
265// For Generic Subtypes of OLattice<> and OScalar<>
266// Do the generic thing
267template<class T1, class T2, int N>
268inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N> , FnSReconDir1PlusProd>::Type_t
276
277// This rewrites spinReconstructDir1Plus( Matrix * Vec )
278// as a fused reconstruct: FnSReconDir1PlusProd< Matrix, Vec >
279template<class T1, class T2, template<class> class C, int N >
280inline typename MakeReturn<BinaryNode< FnSReconDir1PlusProd,
281 typename CreateLeaf< QDPType< PScalar<T1> , C<PScalar<T1> > > >::Leaf_t,
282 typename CreateLeaf< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >::Leaf_t>,
283 typename BinaryReturn< C<PScalar<T1> >, C<PSpinVector<T2,N> >, FnSReconDir1PlusProd>::Type_t >::Expression_t
284spinReconstructDir1Plus(const QDPExpr<
285 BinaryNode< OpMultiply,
286 Reference< QDPType< PScalar<T1> , C<PScalar<T1> > > >,
287 Reference< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >
288 >,
289 typename BinaryReturn< QDPType< PScalar<T1> , C<PScalar<T1> > >,
290 QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > >,
291 OpMultiply >::Type_t >& l)
292{
293
294
295 typedef BinaryNode< FnSReconDir1PlusProd,
296 typename CreateLeaf< QDPType< PScalar<T1> , C<PScalar<T1> > > >::Leaf_t,
297 typename CreateLeaf< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >::Leaf_t> Tree_t;
298
299
300 return MakeReturn<Tree_t,typename BinaryReturn< C<PScalar<T1> >, C<PSpinVector<T2,N> >, FnSReconDir1PlusProd>::Type_t>::make(
301 Tree_t(
302 CreateLeaf<QDPType<PScalar<T1> ,C<PScalar<T1> > > >::make(l.expression().left()),
303 CreateLeaf<QDPType<PSpinVector<T2,N>,C<PSpinVector<T2,N> > > >::make(l.expression().right())
304 )
305 );
306
307
308}
309
310
314struct FnSReconDir2MinusProd
315{
316 PETE_EMPTY_CONSTRUCTORS(FnSReconDir2MinusProd)
317
318 template<class T1, class T2, int N>
319 inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N>, FnSReconDir2MinusProd>::Type_t
320 operator()(const PScalar<T1>& a, const PSpinVector<T2,N>& b) const
321 {
322 return sreconDir2Minus(a, b);
323 }
324};
325
326
327// Proclaim that the return type of the fused op is the same as the return type
328// of the SpinReconstruction on a normal half vector
329template<class T1, class T2, int N>
334
335// For Generic Subtypes of OLattice<> and OScalar<>
336// Do the generic thing
337template<class T1, class T2, int N>
338inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N> , FnSReconDir2MinusProd>::Type_t
346
347// This rewrites spinReconstructDir2Minus( Matrix * Vec )
348// as a fused reconstruct: FnSReconDir2MinusProd< Matrix, Vec >
349template<class T1, class T2, template<class> class C, int N >
350inline typename MakeReturn<BinaryNode< FnSReconDir2MinusProd,
351 typename CreateLeaf< QDPType< PScalar<T1> , C<PScalar<T1> > > >::Leaf_t,
352 typename CreateLeaf< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >::Leaf_t>,
353 typename BinaryReturn< C<PScalar<T1> >, C<PSpinVector<T2,N> >, FnSReconDir2MinusProd>::Type_t >::Expression_t
354spinReconstructDir2Minus(const QDPExpr<
355 BinaryNode< OpMultiply,
356 Reference< QDPType< PScalar<T1> , C<PScalar<T1> > > >,
357 Reference< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >
358 >,
359 typename BinaryReturn< QDPType< PScalar<T1> , C<PScalar<T1> > >,
360 QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > >,
361 OpMultiply >::Type_t >& l)
362{
363
364
365 typedef BinaryNode< FnSReconDir2MinusProd,
366 typename CreateLeaf< QDPType< PScalar<T1> , C<PScalar<T1> > > >::Leaf_t,
367 typename CreateLeaf< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >::Leaf_t> Tree_t;
368
369
370 return MakeReturn<Tree_t,typename BinaryReturn< C<PScalar<T1> >, C<PSpinVector<T2,N> >, FnSReconDir2MinusProd>::Type_t>::make(
371 Tree_t(
372 CreateLeaf<QDPType<PScalar<T1> ,C<PScalar<T1> > > >::make(l.expression().left()),
373 CreateLeaf<QDPType<PSpinVector<T2,N>,C<PSpinVector<T2,N> > > >::make(l.expression().right())
374 )
375 );
376
377
378}
379
383struct FnSReconDir2PlusProd
384{
385 PETE_EMPTY_CONSTRUCTORS(FnSReconDir2PlusProd)
386
387 template<class T1, class T2, int N>
388 inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N>, FnSReconDir2PlusProd>::Type_t
389 operator()(const PScalar<T1>& a, const PSpinVector<T2,N>& b) const
390 {
391 return sreconDir2Plus(a, b);
392 }
393};
394
395
396// Proclaim that the return type of the fused op is the same as the return type
397// of the SpinReconstruction on a normal half vector
398template<class T1, class T2, int N>
403
404// For Generic Subtypes of OLattice<> and OScalar<>
405// Do the generic thing
406template<class T1, class T2, int N>
407inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N> , FnSReconDir2PlusProd>::Type_t
415
416// This rewrites spinReconstructDir2Plus( Matrix * Vec )
417// as a fused reconstruct: FnSReconDir2PlusProd< Matrix, Vec >
418template<class T1, class T2, template<class> class C, int N >
419inline typename MakeReturn<BinaryNode< FnSReconDir2PlusProd,
420 typename CreateLeaf< QDPType< PScalar<T1> , C<PScalar<T1> > > >::Leaf_t,
421 typename CreateLeaf< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >::Leaf_t>,
422 typename BinaryReturn< C<PScalar<T1> >, C<PSpinVector<T2,N> >, FnSReconDir2PlusProd>::Type_t >::Expression_t
423spinReconstructDir2Plus(const QDPExpr<
424 BinaryNode< OpMultiply,
425 Reference< QDPType< PScalar<T1> , C<PScalar<T1> > > >,
426 Reference< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >
427 >,
428 typename BinaryReturn< QDPType< PScalar<T1> , C<PScalar<T1> > >,
429 QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > >,
430 OpMultiply >::Type_t >& l)
431{
432
433
434 typedef BinaryNode< FnSReconDir2PlusProd,
435 typename CreateLeaf< QDPType< PScalar<T1> , C<PScalar<T1> > > >::Leaf_t,
436 typename CreateLeaf< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >::Leaf_t> Tree_t;
437
438
439 return MakeReturn<Tree_t,typename BinaryReturn< C<PScalar<T1> >, C<PSpinVector<T2,N> >, FnSReconDir2PlusProd>::Type_t>::make(
440 Tree_t(
441 CreateLeaf<QDPType<PScalar<T1> ,C<PScalar<T1> > > >::make(l.expression().left()),
442 CreateLeaf<QDPType<PSpinVector<T2,N>,C<PSpinVector<T2,N> > > >::make(l.expression().right())
443 )
444 );
445
446
447}
448
449
453struct FnSReconDir3MinusProd
454{
455 PETE_EMPTY_CONSTRUCTORS(FnSReconDir3MinusProd)
456
457 template<class T1, class T2, int N>
458 inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N>, FnSReconDir3MinusProd>::Type_t
459 operator()(const PScalar<T1>& a, const PSpinVector<T2,N>& b) const
460 {
461 return sreconDir3Minus(a, b);
462 }
463};
464
465
466// Proclaim that the return type of the fused op is the same as the return type
467// of the SpinReconstruction on a normal half vector
468template<class T1, class T2, int N>
473
474// For Generic Subtypes of OLattice<> and OScalar<>
475// Do the generic thing
476template<class T1, class T2, int N>
477inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N> , FnSReconDir3MinusProd>::Type_t
485
486// This rewrites spinReconstructDir3Minus( Matrix * Vec )
487// as a fused reconstruct: FnSReconDir3MinusProd< Matrix, Vec >
488template<class T1, class T2, template<class> class C, int N >
489inline typename MakeReturn<BinaryNode< FnSReconDir3MinusProd,
490 typename CreateLeaf< QDPType< PScalar<T1> , C<PScalar<T1> > > >::Leaf_t,
491 typename CreateLeaf< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >::Leaf_t>,
492 typename BinaryReturn< C<PScalar<T1> >, C<PSpinVector<T2,N> >, FnSReconDir3MinusProd>::Type_t >::Expression_t
517
521struct FnSReconDir3PlusProd
522{
523 PETE_EMPTY_CONSTRUCTORS(FnSReconDir3PlusProd)
524
525 template<class T1, class T2, int N>
526 inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N>, FnSReconDir3PlusProd>::Type_t
527 operator()(const PScalar<T1>& a, const PSpinVector<T2,N>& b) const
528 {
529 return sreconDir3Plus(a, b);
530 }
531};
532
533
534// Proclaim that the return type of the fused op is the same as the return type
535// of the SpinReconstruction on a normal half vector
536template<class T1, class T2, int N>
541
542// For Generic Subtypes of OLattice<> and OScalar<>
543// Do the generic thing
544template<class T1, class T2, int N>
545inline typename BinaryReturn<PScalar<T1>, PSpinVector<T2,N> , FnSReconDir3PlusProd>::Type_t
553
554// This rewrites spinReconstructDir3Plus( Matrix * Vec )
555// as a fused reconstruct: FnSReconDir3PlusProd< Matrix, Vec >
556template<class T1, class T2, template<class> class C, int N >
557inline typename MakeReturn<BinaryNode< FnSReconDir3PlusProd,
558 typename CreateLeaf< QDPType< PScalar<T1> , C<PScalar<T1> > > >::Leaf_t,
559 typename CreateLeaf< QDPType< PSpinVector<T2,N>, C<PSpinVector<T2,N> > > >::Leaf_t>,
560 typename BinaryReturn< C<PScalar<T1> >, C<PSpinVector<T2,N> >, FnSReconDir3PlusProd>::Type_t >::Expression_t
585
586
587
591
592// The actual fused op - Replace this with specialist code
593template<>
594inline BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir0MinusProd>::Type_t
595sreconDir0Minus(const PScalar<ColMat>& a, const PSpinVector<ColVec,2>& b)
596{
597 PSpinVector<ColVec,2> d;
598 BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir0MinusProd>::Type_t ret;
599
600 _inline_mult_su3_mat_vec(a.elem(),
601 b.elem(0),
602 d.elem(0));
603
604 _inline_mult_su3_mat_vec(a.elem(),
605 b.elem(1),
606 d.elem(1));
607
608 inlineSpinReconDir0Minus(&(d.elem(0).elem(0).real()),
609 &(ret.elem(0).elem(0).real()),
610 1);
611
612 return ret;
613}
614
615
616// The actual fused op - Replace this with specialist code
617template<>
618inline BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir0PlusProd>::Type_t
619sreconDir0Plus(const PScalar<ColMat>& a, const PSpinVector<ColVec,2>& b)
620{
621 PSpinVector<ColVec,2> d;
622 BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir0MinusProd>::Type_t ret;
623
624 _inline_mult_su3_mat_vec(a.elem(),
625 b.elem(0),
626 d.elem(0));
627
628 _inline_mult_su3_mat_vec(a.elem(),
629 b.elem(1),
630 d.elem(1));
631
632 inlineSpinReconDir0Plus(&(d.elem(0).elem(0).real()),
633 &(ret.elem(0).elem(0).real()),
634 1);
635
636
637 return ret;
638}
639
640// The actual fused op - Replace this with specialist code
641template<>
642inline BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir1MinusProd>::Type_t
643sreconDir1Minus(const PScalar<ColMat>& a, const PSpinVector<ColVec,2>& b)
644{
645 PSpinVector<ColVec,2> d;
646 BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir0MinusProd>::Type_t ret;
647
648 _inline_mult_su3_mat_vec(a.elem(),
649 b.elem(0),
650 d.elem(0));
651
652 _inline_mult_su3_mat_vec(a.elem(),
653 b.elem(1),
654 d.elem(1));
655
656 inlineSpinReconDir1Minus(&(d.elem(0).elem(0).real()),
657 &(ret.elem(0).elem(0).real()),
658 1);
659
660 return ret;
661}
662
663// The actual fused op - Replace this with specialist code
664template<>
665inline BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir1PlusProd>::Type_t
666sreconDir1Plus(const PScalar<ColMat>& a, const PSpinVector<ColVec,2>& b)
667{
668 PSpinVector<ColVec,2> d;
669 BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir0MinusProd>::Type_t ret;
670
671 _inline_mult_su3_mat_vec(a.elem(),
672 b.elem(0),
673 d.elem(0));
674
675 _inline_mult_su3_mat_vec(a.elem(),
676 b.elem(1),
677 d.elem(1));
678
679 inlineSpinReconDir1Plus(&(d.elem(0).elem(0).real()),
680 &(ret.elem(0).elem(0).real()),
681 1);
682
683 return ret;
684}
685
686// The actual fused op - Replace this with specialist code
687template<>
688inline BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir2MinusProd>::Type_t
689sreconDir2Minus(const PScalar<ColMat>& a, const PSpinVector<ColVec,2>& b)
690{
691 PSpinVector<ColVec,2> d;
692 BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir0MinusProd>::Type_t ret;
693
694 _inline_mult_su3_mat_vec(a.elem(),
695 b.elem(0),
696 d.elem(0));
697
698 _inline_mult_su3_mat_vec(a.elem(),
699 b.elem(1),
700 d.elem(1));
701
702 inlineSpinReconDir2Minus(&(d.elem(0).elem(0).real()),
703 &(ret.elem(0).elem(0).real()),
704 1);
705
706 return ret;
707}
708
709// The actual fused op - Replace this with specialist code
710template<>
711inline BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir2PlusProd>::Type_t
712sreconDir2Plus(const PScalar<ColMat>& a, const PSpinVector<ColVec,2>& b)
713{
714 PSpinVector<ColVec,2> d;
715 BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir0MinusProd>::Type_t ret;
716
717 _inline_mult_su3_mat_vec(a.elem(),
718 b.elem(0),
719 d.elem(0));
720
721 _inline_mult_su3_mat_vec(a.elem(),
722 b.elem(1),
723 d.elem(1));
724
725 inlineSpinReconDir2Plus(&(d.elem(0).elem(0).real()),
726 &(ret.elem(0).elem(0).real()),
727 1);
728
729 return ret;
730
731}
732
733// The actual fused op - Replace this with specialist code
734template<>
735inline BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir3MinusProd>::Type_t
736sreconDir3Minus(const PScalar<ColMat>& a, const PSpinVector<ColVec,2>& b)
737{
738 PSpinVector<ColVec,2> d;
739 BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir0MinusProd>::Type_t ret;
740
741 _inline_mult_su3_mat_vec(a.elem(),
742 b.elem(0),
743 d.elem(0));
744
745 _inline_mult_su3_mat_vec(a.elem(),
746 b.elem(1),
747 d.elem(1));
748
749 inlineSpinReconDir3Minus(&(d.elem(0).elem(0).real()),
750 &(ret.elem(0).elem(0).real()),
751 1);
752
753 return ret;
754}
755
756
757// The actual fused op - Replace this with specialist code
758template<>
759inline BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir3PlusProd>::Type_t
760sreconDir3Plus(const PScalar<ColMat>& a, const PSpinVector<ColVec,2>& b)
761{
762 PSpinVector<ColVec,2> d;
763 BinaryReturn< PScalar<ColMat>, PSpinVector<ColVec,2>, FnSReconDir0PlusProd>::Type_t ret;
764
765 _inline_mult_su3_mat_vec(a.elem(),
766 b.elem(0),
767 d.elem(0));
768
769 _inline_mult_su3_mat_vec(a.elem(),
770 b.elem(1),
771 d.elem(1));
772
773 inlineSpinReconDir3Plus(&(d.elem(0).elem(0).real()),
774 &(ret.elem(0).elem(0).real()),
775 1);
776
777 return ret;
778}
779
780
781} // namespace QDP
782
783#endif
PETE_EMPTY_CONSTRUCTORS
#define PETE_EMPTY_CONSTRUCTORS(CLASS)
Definition PETE.h:58
QDP::OLattice
Outer grid Lattice type.
Definition qdp_outer.h:264
QDP::OScalar
Outer grid Scalar class *‍/.
Definition qdp_outer.h:37
QDP::PColorMatrix
Primitive color Matrix class.
Definition qdp_primcolormat.h:26
QDP::PColorVector
Primitive color Vector class.
Definition qdp_primcolorvec.h:24
QDP::PScalar
Primitive Scalar.
Definition qdp_primscalar.h:25
QDP::PSpinVector
Primitive spin Vector class.
Definition qdp_primspinvec.h:33
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
_inline_mult_su3_mat_vec
#define _inline_mult_su3_mat_vec(aa, bb, cc)
Definition generic_mv_switchbox.h:7
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::sreconDir2Plus
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir2PlusProd >::Type_t sreconDir2Plus(const PScalar< T1 > &a, const PSpinVector< T2, N > &b)
Definition generic_fused_spin_recon.h:408
QDP::inlineSpinReconDir0Minus
void inlineSpinReconDir0Minus(const REAL *src, REAL *dst, unsigned int n_vec)
Spin recon (1/2)(1-\gamma_0).
Definition generic_spin_recon_inlines.h:88
QDP::ColVec
PColorVector< RComplex< REAL >, 3 > ColVec
Definition generic_fused_spin_proj.h:8
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::ColMat
PColorMatrix< RComplex< REAL >, 3 > ColMat
Definition generic_fused_spin_proj.h:11
QDP::inlineSpinReconDir2Plus
void inlineSpinReconDir2Plus(const REAL *src, REAL *dst, unsigned int n_vec)
Spin recon (1/2)(1+\gamma_2).
Definition generic_spin_recon_inlines.h:292
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::inlineSpinReconDir1Minus
void inlineSpinReconDir1Minus(const REAL *src, REAL *dst, unsigned int n_vec)
Spin recon (1/2)(1-\gamma_1).
Definition generic_spin_recon_inlines.h:221
QDP::sreconDir3Plus
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir3PlusProd >::Type_t sreconDir3Plus(const PScalar< T1 > &a, const PSpinVector< T2, N > &b)
Definition generic_fused_spin_recon.h:546
QDP::Spin4
PSpinVector< ColVec, 4 > Spin4
Definition generic_fused_spin_proj.h:9
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::Spin2
PSpinVector< ColVec, 2 > Spin2
Definition generic_fused_spin_proj.h:10
QDP::sreconDir0Minus
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir0MinusProd >::Type_t sreconDir0Minus(const PScalar< T1 > &a, const PSpinVector< T2, N > &b)
Definition generic_fused_spin_recon.h:62
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::sreconDir3Minus
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir3MinusProd >::Type_t sreconDir3Minus(const PScalar< T1 > &a, const PSpinVector< T2, N > &b)
Definition generic_fused_spin_recon.h:478
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::sreconDir1Plus
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir1PlusProd >::Type_t sreconDir1Plus(const PScalar< T1 > &a, const PSpinVector< T2, N > &b)
Definition generic_fused_spin_recon.h:269
QDP::inlineSpinReconDir2Minus
void inlineSpinReconDir2Minus(const REAL *src, REAL *dst, unsigned int n_vec)
Spin recon (1/2)(1-\gamma_2).
Definition generic_spin_recon_inlines.h:359
QDP::sreconDir0Plus
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir0PlusProd >::Type_t sreconDir0Plus(const PScalar< T1 > &a, const PSpinVector< T2, N > &b)
Definition generic_fused_spin_recon.h:131
QDP::sreconDir1Minus
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir1MinusProd >::Type_t sreconDir1Minus(const PScalar< T1 > &a, const PSpinVector< T2, N > &b)
Definition generic_fused_spin_recon.h:200
QDP::sreconDir2Minus
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir2MinusProd >::Type_t sreconDir2Minus(const PScalar< T1 > &a, const PSpinVector< T2, N > &b)
Definition generic_fused_spin_recon.h:339
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::inlineSpinReconDir3Minus
void inlineSpinReconDir3Minus(const REAL *src, REAL *dst, unsigned int n_vec)
Spin recon (1/2)(1-\gamma3).
Definition generic_spin_recon_inlines.h:495
QDP::inlineSpinReconDir3Plus
void inlineSpinReconDir3Plus(const REAL *src, REAL *dst, unsigned int n_vec)
Spin recon (1/2)(1+\gamma3).
Definition generic_spin_recon_inlines.h:427
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::inlineSpinReconDir1Plus
void inlineSpinReconDir1Plus(const REAL *src, REAL *dst, unsigned int n_vec)
Spin recon (1/2)(1+\gamma_1).
Definition generic_spin_recon_inlines.h:154
QDP::inlineSpinReconDir0Plus
void inlineSpinReconDir0Plus(const REAL *src, REAL *dst, unsigned int n_vec)
Spin recon (1/2)(1+\gamma_0).
Definition generic_spin_recon_inlines.h:23
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
C
OLattice< PScalar< PColorMatrix< RComplexFloat, 3 > > > C
Definition qdp_scalarsite_generic_linalg_wrapper.h:11
QDP::BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir0MinusProd >::Type_t
UnaryReturn< PSpinVector< T2, N >, FnSpinReconstructDir0Minus >::Type_t Type_t
Definition generic_fused_spin_recon.h:55
QDP::BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir0PlusProd >::Type_t
UnaryReturn< PSpinVector< T2, N >, FnSpinReconstructDir0Plus >::Type_t Type_t
Definition generic_fused_spin_recon.h:124
QDP::BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir1MinusProd >::Type_t
UnaryReturn< PSpinVector< T2, N >, FnSpinReconstructDir1Minus >::Type_t Type_t
Definition generic_fused_spin_recon.h:193
QDP::BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir1PlusProd >::Type_t
UnaryReturn< PSpinVector< T2, N >, FnSpinReconstructDir1Plus >::Type_t Type_t
Definition generic_fused_spin_recon.h:262
QDP::BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir2MinusProd >::Type_t
UnaryReturn< PSpinVector< T2, N >, FnSpinReconstructDir2Minus >::Type_t Type_t
Definition generic_fused_spin_recon.h:332
QDP::BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir2PlusProd >::Type_t
UnaryReturn< PSpinVector< T2, N >, FnSpinReconstructDir2Plus >::Type_t Type_t
Definition generic_fused_spin_recon.h:401
QDP::BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir3MinusProd >::Type_t
UnaryReturn< PSpinVector< T2, N >, FnSpinReconstructDir3Minus >::Type_t Type_t
Definition generic_fused_spin_recon.h:471
QDP::BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir3PlusProd >::Type_t
UnaryReturn< PSpinVector< T2, N >, FnSpinReconstructDir3Plus >::Type_t Type_t
Definition generic_fused_spin_recon.h:539
QDP::FnSReconDir0MinusProd
Reconstruct Dir0Minus Fused.
Definition generic_fused_spin_recon.h:38
QDP::FnSReconDir0MinusProd::operator()
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir0MinusProd >::Type_t operator()(const PScalar< T1 > &a, const PSpinVector< T2, N > &b) const
Definition generic_fused_spin_recon.h:43
QDP::FnSReconDir0PlusProd
Reconstruct Dir0Plus Fused.
Definition generic_fused_spin_recon.h:107
QDP::FnSReconDir0PlusProd::operator()
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir0PlusProd >::Type_t operator()(const PScalar< T1 > &a, const PSpinVector< T2, N > &b) const
Definition generic_fused_spin_recon.h:112
QDP::FnSReconDir1MinusProd
Reconstruct Dir1Minus Fused.
Definition generic_fused_spin_recon.h:176
QDP::FnSReconDir1MinusProd::operator()
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir1MinusProd >::Type_t operator()(const PScalar< T1 > &a, const PSpinVector< T2, N > &b) const
Definition generic_fused_spin_recon.h:181
QDP::FnSReconDir1PlusProd
Reconstruct Dir1Plus Fused.
Definition generic_fused_spin_recon.h:245
QDP::FnSReconDir1PlusProd::operator()
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir1PlusProd >::Type_t operator()(const PScalar< T1 > &a, const PSpinVector< T2, N > &b) const
Definition generic_fused_spin_recon.h:250
QDP::FnSReconDir2MinusProd
Reconstruct Dir2Minus Fused.
Definition generic_fused_spin_recon.h:315
QDP::FnSReconDir2MinusProd::operator()
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir2MinusProd >::Type_t operator()(const PScalar< T1 > &a, const PSpinVector< T2, N > &b) const
Definition generic_fused_spin_recon.h:320
QDP::FnSReconDir2PlusProd
Reconstruct Dir2Plus Fused.
Definition generic_fused_spin_recon.h:384
QDP::FnSReconDir2PlusProd::operator()
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir2PlusProd >::Type_t operator()(const PScalar< T1 > &a, const PSpinVector< T2, N > &b) const
Definition generic_fused_spin_recon.h:389
QDP::FnSReconDir3MinusProd
Reconstruct Dir3Minus Fused.
Definition generic_fused_spin_recon.h:454
QDP::FnSReconDir3MinusProd::operator()
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir3MinusProd >::Type_t operator()(const PScalar< T1 > &a, const PSpinVector< T2, N > &b) const
Definition generic_fused_spin_recon.h:459
QDP::FnSReconDir3PlusProd
Reconstruct Dir3Plus Fused.
Definition generic_fused_spin_recon.h:522
QDP::FnSReconDir3PlusProd::operator()
BinaryReturn< PScalar< T1 >, PSpinVector< T2, N >, FnSReconDir3PlusProd >::Type_t operator()(const PScalar< T1 > &a, const PSpinVector< T2, N > &b) const
Definition generic_fused_spin_recon.h:527
Generated by doxygen 1.16.1