QDP++
generic_mult_nn.h
Go to the documentation of this file.
1#define _inline_generic_mult_su3_nn(aa,bb,cc) \
2{ \
3 cc.elem(0,0).real() = aa.elem(0,0).real()*bb.elem(0,0).real(); \
4 cc.elem(0,0).real() -= aa.elem(0,0).imag()*bb.elem(0,0).imag(); \
5 cc.elem(0,0).real() += aa.elem(0,1).real()*bb.elem(1,0).real(); \
6 cc.elem(0,0).real() -= aa.elem(0,1).imag()*bb.elem(1,0).imag(); \
7 cc.elem(0,0).real() += aa.elem(0,2).real()*bb.elem(2,0).real(); \
8 cc.elem(0,0).real() -= aa.elem(0,2).imag()*bb.elem(2,0).imag(); \
9 cc.elem(0,0).imag() = aa.elem(0,0).real()*bb.elem(0,0).imag(); \
10 cc.elem(0,0).imag() += aa.elem(0,0).imag()*bb.elem(0,0).real(); \
11 cc.elem(0,0).imag() += aa.elem(0,1).real()*bb.elem(1,0).imag(); \
12 cc.elem(0,0).imag() += aa.elem(0,1).imag()*bb.elem(1,0).real(); \
13 cc.elem(0,0).imag() += aa.elem(0,2).real()*bb.elem(2,0).imag(); \
14 cc.elem(0,0).imag() += aa.elem(0,2).imag()*bb.elem(2,0).real(); \
15 cc.elem(1,0).real() = aa.elem(1,0).real()*bb.elem(0,0).real(); \
16 cc.elem(1,0).real() -= aa.elem(1,0).imag()*bb.elem(0,0).imag(); \
17 cc.elem(1,0).real() += aa.elem(1,1).real()*bb.elem(1,0).real(); \
18 cc.elem(1,0).real() -= aa.elem(1,1).imag()*bb.elem(1,0).imag(); \
19 cc.elem(1,0).real() += aa.elem(1,2).real()*bb.elem(2,0).real(); \
20 cc.elem(1,0).real() -= aa.elem(1,2).imag()*bb.elem(2,0).imag(); \
21 cc.elem(1,0).imag() = aa.elem(1,0).real()*bb.elem(0,0).imag(); \
22 cc.elem(1,0).imag() += aa.elem(1,0).imag()*bb.elem(0,0).real(); \
23 cc.elem(1,0).imag() += aa.elem(1,1).real()*bb.elem(1,0).imag(); \
24 cc.elem(1,0).imag() += aa.elem(1,1).imag()*bb.elem(1,0).real(); \
25 cc.elem(1,0).imag() += aa.elem(1,2).real()*bb.elem(2,0).imag(); \
26 cc.elem(1,0).imag() += aa.elem(1,2).imag()*bb.elem(2,0).real(); \
27 cc.elem(2,0).real() = aa.elem(2,0).real()*bb.elem(0,0).real(); \
28 cc.elem(2,0).real() -= aa.elem(2,0).imag()*bb.elem(0,0).imag(); \
29 cc.elem(2,0).real() += aa.elem(2,1).real()*bb.elem(1,0).real(); \
30 cc.elem(2,0).real() -= aa.elem(2,1).imag()*bb.elem(1,0).imag(); \
31 cc.elem(2,0).real() += aa.elem(2,2).real()*bb.elem(2,0).real(); \
32 cc.elem(2,0).real() -= aa.elem(2,2).imag()*bb.elem(2,0).imag(); \
33 cc.elem(2,0).imag() = aa.elem(2,0).real()*bb.elem(0,0).imag(); \
34 cc.elem(2,0).imag() += aa.elem(2,0).imag()*bb.elem(0,0).real(); \
35 cc.elem(2,0).imag() += aa.elem(2,1).real()*bb.elem(1,0).imag(); \
36 cc.elem(2,0).imag() += aa.elem(2,1).imag()*bb.elem(1,0).real(); \
37 cc.elem(2,0).imag() += aa.elem(2,2).real()*bb.elem(2,0).imag(); \
38 cc.elem(2,0).imag() += aa.elem(2,2).imag()*bb.elem(2,0).real(); \
39\
40 cc.elem(0,1).real() = aa.elem(0,0).real()*bb.elem(0,1).real(); \
41 cc.elem(0,1).real() -= aa.elem(0,0).imag()*bb.elem(0,1).imag(); \
42 cc.elem(0,1).real() += aa.elem(0,1).real()*bb.elem(1,1).real(); \
43 cc.elem(0,1).real() -= aa.elem(0,1).imag()*bb.elem(1,1).imag(); \
44 cc.elem(0,1).real() += aa.elem(0,2).real()*bb.elem(2,1).real(); \
45 cc.elem(0,1).real() -= aa.elem(0,2).imag()*bb.elem(2,1).imag(); \
46 cc.elem(0,1).imag() = aa.elem(0,0).real()*bb.elem(0,1).imag(); \
47 cc.elem(0,1).imag() += aa.elem(0,0).imag()*bb.elem(0,1).real(); \
48 cc.elem(0,1).imag() += aa.elem(0,1).real()*bb.elem(1,1).imag(); \
49 cc.elem(0,1).imag() += aa.elem(0,1).imag()*bb.elem(1,1).real(); \
50 cc.elem(0,1).imag() += aa.elem(0,2).real()*bb.elem(2,1).imag(); \
51 cc.elem(0,1).imag() += aa.elem(0,2).imag()*bb.elem(2,1).real(); \
52 cc.elem(1,1).real() = aa.elem(1,0).real()*bb.elem(0,1).real(); \
53 cc.elem(1,1).real() -= aa.elem(1,0).imag()*bb.elem(0,1).imag(); \
54 cc.elem(1,1).real() += aa.elem(1,1).real()*bb.elem(1,1).real(); \
55 cc.elem(1,1).real() -= aa.elem(1,1).imag()*bb.elem(1,1).imag(); \
56 cc.elem(1,1).real() += aa.elem(1,2).real()*bb.elem(2,1).real(); \
57 cc.elem(1,1).real() -= aa.elem(1,2).imag()*bb.elem(2,1).imag(); \
58 cc.elem(1,1).imag() = aa.elem(1,0).real()*bb.elem(0,1).imag(); \
59 cc.elem(1,1).imag() += aa.elem(1,0).imag()*bb.elem(0,1).real(); \
60 cc.elem(1,1).imag() += aa.elem(1,1).real()*bb.elem(1,1).imag(); \
61 cc.elem(1,1).imag() += aa.elem(1,1).imag()*bb.elem(1,1).real(); \
62 cc.elem(1,1).imag() += aa.elem(1,2).real()*bb.elem(2,1).imag(); \
63 cc.elem(1,1).imag() += aa.elem(1,2).imag()*bb.elem(2,1).real(); \
64 cc.elem(2,1).real() = aa.elem(2,0).real()*bb.elem(0,1).real(); \
65 cc.elem(2,1).real() -= aa.elem(2,0).imag()*bb.elem(0,1).imag(); \
66 cc.elem(2,1).real() += aa.elem(2,1).real()*bb.elem(1,1).real(); \
67 cc.elem(2,1).real() -= aa.elem(2,1).imag()*bb.elem(1,1).imag(); \
68 cc.elem(2,1).real() += aa.elem(2,2).real()*bb.elem(2,1).real(); \
69 cc.elem(2,1).real() -= aa.elem(2,2).imag()*bb.elem(2,1).imag(); \
70 cc.elem(2,1).imag() = aa.elem(2,0).real()*bb.elem(0,1).imag(); \
71 cc.elem(2,1).imag() += aa.elem(2,0).imag()*bb.elem(0,1).real(); \
72 cc.elem(2,1).imag() += aa.elem(2,1).real()*bb.elem(1,1).imag(); \
73 cc.elem(2,1).imag() += aa.elem(2,1).imag()*bb.elem(1,1).real(); \
74 cc.elem(2,1).imag() += aa.elem(2,2).real()*bb.elem(2,1).imag(); \
75 cc.elem(2,1).imag() += aa.elem(2,2).imag()*bb.elem(2,1).real(); \
76\
77 cc.elem(0,2).real() = aa.elem(0,0).real()*bb.elem(0,2).real(); \
78 cc.elem(0,2).real() -= aa.elem(0,0).imag()*bb.elem(0,2).imag(); \
79 cc.elem(0,2).real() += aa.elem(0,1).real()*bb.elem(1,2).real(); \
80 cc.elem(0,2).real() -= aa.elem(0,1).imag()*bb.elem(1,2).imag(); \
81 cc.elem(0,2).real() += aa.elem(0,2).real()*bb.elem(2,2).real(); \
82 cc.elem(0,2).real() -= aa.elem(0,2).imag()*bb.elem(2,2).imag(); \
83 cc.elem(0,2).imag() = aa.elem(0,0).real()*bb.elem(0,2).imag(); \
84 cc.elem(0,2).imag() += aa.elem(0,0).imag()*bb.elem(0,2).real(); \
85 cc.elem(0,2).imag() += aa.elem(0,1).real()*bb.elem(1,2).imag(); \
86 cc.elem(0,2).imag() += aa.elem(0,1).imag()*bb.elem(1,2).real(); \
87 cc.elem(0,2).imag() += aa.elem(0,2).real()*bb.elem(2,2).imag(); \
88 cc.elem(0,2).imag() += aa.elem(0,2).imag()*bb.elem(2,2).real(); \
89 cc.elem(1,2).real() = aa.elem(1,0).real()*bb.elem(0,2).real(); \
90 cc.elem(1,2).real() -= aa.elem(1,0).imag()*bb.elem(0,2).imag(); \
91 cc.elem(1,2).real() += aa.elem(1,1).real()*bb.elem(1,2).real(); \
92 cc.elem(1,2).real() -= aa.elem(1,1).imag()*bb.elem(1,2).imag(); \
93 cc.elem(1,2).real() += aa.elem(1,2).real()*bb.elem(2,2).real(); \
94 cc.elem(1,2).real() -= aa.elem(1,2).imag()*bb.elem(2,2).imag(); \
95 cc.elem(1,2).imag() = aa.elem(1,0).real()*bb.elem(0,2).imag(); \
96 cc.elem(1,2).imag() += aa.elem(1,0).imag()*bb.elem(0,2).real(); \
97 cc.elem(1,2).imag() += aa.elem(1,1).real()*bb.elem(1,2).imag(); \
98 cc.elem(1,2).imag() += aa.elem(1,1).imag()*bb.elem(1,2).real(); \
99 cc.elem(1,2).imag() += aa.elem(1,2).real()*bb.elem(2,2).imag(); \
100 cc.elem(1,2).imag() += aa.elem(1,2).imag()*bb.elem(2,2).real(); \
101 cc.elem(2,2).real() = aa.elem(2,0).real()*bb.elem(0,2).real(); \
102 cc.elem(2,2).real() -= aa.elem(2,0).imag()*bb.elem(0,2).imag(); \
103 cc.elem(2,2).real() += aa.elem(2,1).real()*bb.elem(1,2).real(); \
104 cc.elem(2,2).real() -= aa.elem(2,1).imag()*bb.elem(1,2).imag(); \
105 cc.elem(2,2).real() += aa.elem(2,2).real()*bb.elem(2,2).real(); \
106 cc.elem(2,2).real() -= aa.elem(2,2).imag()*bb.elem(2,2).imag(); \
107 cc.elem(2,2).imag() = aa.elem(2,0).real()*bb.elem(0,2).imag(); \
108 cc.elem(2,2).imag() += aa.elem(2,0).imag()*bb.elem(0,2).real(); \
109 cc.elem(2,2).imag() += aa.elem(2,1).real()*bb.elem(1,2).imag(); \
110 cc.elem(2,2).imag() += aa.elem(2,1).imag()*bb.elem(1,2).real(); \
111 cc.elem(2,2).imag() += aa.elem(2,2).real()*bb.elem(2,2).imag(); \
112 cc.elem(2,2).imag() += aa.elem(2,2).imag()*bb.elem(2,2).real(); \
113}
Generated by doxygen 1.16.1