Go to the documentation of this file.
11 #define MAC_Apply_G_mx3_ops( m_A, \
20 float ga12 = *gamma12; \
21 float si12 = *sigma12; \
22 float ga23 = *gamma23; \
23 float si23 = *sigma23; \
24 float* restrict alpha1 = a1; \
25 float* restrict alpha2 = a2; \
26 float* restrict alpha3 = a3; \
32 for ( i = 0; i < m_A; ++i ) \
37 *alpha1 = temp1 * ga12 + temp2 * si12; \
38 *alpha2 = temp2 * ga12 - temp1 * si12; \
43 *alpha2 = temp2 * ga23 + temp3 * si23; \
44 *alpha3 = temp3 * ga23 - temp2 * si23; \
52 #define MAC_Apply_G_mx3_opd( m_A, \
61 double ga12 = *gamma12; \
62 double si12 = *sigma12; \
63 double ga23 = *gamma23; \
64 double si23 = *sigma23; \
65 double* restrict alpha1 = a1; \
66 double* restrict alpha2 = a2; \
67 double* restrict alpha3 = a3; \
73 for ( i = 0; i < m_A; ++i ) \
78 *alpha1 = temp1 * ga12 + temp2 * si12; \
79 *alpha2 = temp2 * ga12 - temp1 * si12; \
84 *alpha2 = temp2 * ga23 + temp3 * si23; \
85 *alpha3 = temp3 * ga23 - temp2 * si23; \
93 #define MAC_Apply_G_mx3_opc( m_A, \
102 float ga12 = *gamma12; \
103 float si12 = *sigma12; \
104 float ga23 = *gamma23; \
105 float si23 = *sigma23; \
106 scomplex* restrict alpha1 = a1; \
107 scomplex* restrict alpha2 = a2; \
108 scomplex* restrict alpha3 = a3; \
114 for ( i = 0; i < m_A; ++i ) \
119 alpha1->real = ga12 * temp1.real + si12 * temp2.real; \
120 alpha1->imag = ga12 * temp1.imag + si12 * temp2.imag; \
122 alpha2->real = -si12 * temp1.real + ga12 * temp2.real; \
123 alpha2->imag = -si12 * temp1.imag + ga12 * temp2.imag; \
128 alpha2->real = ga23 * temp2.real + si23 * temp3.real; \
129 alpha2->imag = ga23 * temp2.imag + si23 * temp3.imag; \
131 alpha3->real = -si23 * temp2.real + ga23 * temp3.real; \
132 alpha3->imag = -si23 * temp2.imag + ga23 * temp3.imag; \
140 #define MAC_Apply_G_mx3_opz( m_A, \
149 double ga12 = *gamma12; \
150 double si12 = *sigma12; \
151 double ga23 = *gamma23; \
152 double si23 = *sigma23; \
153 dcomplex* restrict alpha1 = a1; \
154 dcomplex* restrict alpha2 = a2; \
155 dcomplex* restrict alpha3 = a3; \
161 for ( i = 0; i < m_A; ++i ) \
166 alpha1->real = ga12 * temp1.real + si12 * temp2.real; \
167 alpha1->imag = ga12 * temp1.imag + si12 * temp2.imag; \
169 alpha2->real = -si12 * temp1.real + ga12 * temp2.real; \
170 alpha2->imag = -si12 * temp1.imag + ga12 * temp2.imag; \
175 alpha2->real = ga23 * temp2.real + si23 * temp3.real; \
176 alpha2->imag = ga23 * temp2.imag + si23 * temp3.imag; \
178 alpha3->real = -si23 * temp2.real + ga23 * temp3.real; \
179 alpha3->imag = -si23 * temp2.imag + ga23 * temp3.imag; \