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; \