Commit b576a9de authored by Dominik Weckbecker's avatar Dominik Weckbecker 💬
Browse files

optimize sector filtering in encoder

parent bf675ce2
Loading
Loading
Loading
Loading
+18 −18
Original line number Diff line number Diff line
@@ -918,27 +918,27 @@ void calculate_hodirac_sector_parameters(
            {
                if ( i_sec == 0 )
                {
                    *p_sec_w_imag = 1.772454e+00f * ImagBuffer[0][i_bin] + 1.023327e+00f * 1.732050776481628f * ImagBuffer[1][i_bin];
                    *p_sec_x_imag = 1.023327e+00f * 1.732050776481628f * ImagBuffer[3][i_bin] + 4.576456e-01f * 2.236068010330200f * ImagBuffer[4][i_bin];
                    *p_sec_y_imag = 5.908180e-01f * ImagBuffer[0][i_bin] + 1.023327e+00f * 1.732050776481628f * ImagBuffer[1][i_bin] - 2.642218e-01f * 2.236068010330200f * ImagBuffer[6][i_bin] - 4.576456e-01f * 2.236068010330200f * ImagBuffer[8][i_bin];
                    *p_sec_z_imag = 1.023327e+00f * 1.732050776481628f * ImagBuffer[2][i_bin] + 4.576456e-01f * 2.236068010330200f * ImagBuffer[5][i_bin];
                    *p_sec_w_imag = 1.772454f * ImagBuffer[0][i_bin] + 1.772454f * ImagBuffer[1][i_bin];
                    *p_sec_x_imag = 1.772454f * ImagBuffer[3][i_bin] + 1.023326f * ImagBuffer[4][i_bin];
                    *p_sec_y_imag = 0.590818f * ImagBuffer[0][i_bin] + 1.772454f * ImagBuffer[1][i_bin] - 0.590817f * ImagBuffer[6][i_bin] - 1.023326f * ImagBuffer[8][i_bin];
                    *p_sec_z_imag = 1.772454f * ImagBuffer[2][i_bin] + 1.023326f * ImagBuffer[5][i_bin];

                    *p_sec_w_real = 1.772454e+00f * RealBuffer[0][i_bin] + 1.023327e+00f * 1.732050776481628f * RealBuffer[1][i_bin];
                    *p_sec_x_real = 1.023327e+00f * 1.732050776481628f * RealBuffer[3][i_bin] + 4.576456e-01f * 2.236068010330200f * RealBuffer[4][i_bin];
                    *p_sec_y_real = 5.908180e-01f * RealBuffer[0][i_bin] + 1.023327e+00f * 1.732050776481628f * RealBuffer[1][i_bin] - 2.642218e-01f * 2.236068010330200f * RealBuffer[6][i_bin] - 4.576456e-01f * 2.236068010330200f * RealBuffer[8][i_bin];
                    *p_sec_z_real = 1.023327e+00f * 1.732050776481628f * RealBuffer[2][i_bin] + 4.576456e-01f * 2.236068010330200f * RealBuffer[5][i_bin];
                    *p_sec_w_real = 1.772454f * RealBuffer[0][i_bin] + 1.772454f * RealBuffer[1][i_bin];
                    *p_sec_x_real = 1.772454f * RealBuffer[3][i_bin] + 1.023326f * RealBuffer[4][i_bin];
                    *p_sec_y_real = 0.590818f * RealBuffer[0][i_bin] + 1.772454f * RealBuffer[1][i_bin] - 0.590817f * RealBuffer[6][i_bin] - 1.023326f * RealBuffer[8][i_bin];
                    *p_sec_z_real = 1.772454f * RealBuffer[2][i_bin] + 1.023326f * RealBuffer[5][i_bin];
                }
                else
                {
                    *p_sec_w_imag = 1.772454e+00f * ImagBuffer[0][i_bin] - 1.023327e+00f * 1.732050776481628f * ImagBuffer[1][i_bin];
                    *p_sec_x_imag = 1.023327e+00f * 1.732050776481628f * ImagBuffer[3][i_bin] - 4.576456e-01f * 2.236068010330200f * ImagBuffer[4][i_bin];
                    *p_sec_y_imag = -5.908180e-01f * ImagBuffer[0][i_bin] + 1.023327e+00f * 1.732050776481628f * ImagBuffer[1][i_bin] + 2.642218e-01f * 2.236068010330200f * ImagBuffer[6][i_bin] + 4.576456e-01f * 2.236068010330200f * ImagBuffer[8][i_bin];
                    *p_sec_z_imag = 1.023327e+00f * 1.732050776481628f * ImagBuffer[2][i_bin] - 4.576456e-01f * 2.236068010330200f * ImagBuffer[5][i_bin];
                    *p_sec_w_imag = 1.772454f * ImagBuffer[0][i_bin] - 1.772454f * ImagBuffer[1][i_bin];
                    *p_sec_x_imag = 1.772454f * ImagBuffer[3][i_bin] - 1.023326f * ImagBuffer[4][i_bin];
                    *p_sec_y_imag = -0.590818f * ImagBuffer[0][i_bin] + 1.772454f * ImagBuffer[1][i_bin] + 0.590817f * ImagBuffer[6][i_bin] + 1.023326f * ImagBuffer[8][i_bin];
                    *p_sec_z_imag = 1.772454f * ImagBuffer[2][i_bin] - 1.023326f * ImagBuffer[5][i_bin];

                    *p_sec_w_real = 1.772454e+00f * RealBuffer[0][i_bin] - 1.023327e+00f * 1.732050776481628f * RealBuffer[1][i_bin];
                    *p_sec_x_real = 1.023327e+00f * 1.732050776481628f * RealBuffer[3][i_bin] - 4.576456e-01f * 2.236068010330200f * RealBuffer[4][i_bin];
                    *p_sec_y_real = -5.908180e-01f * RealBuffer[0][i_bin] + 1.023327e+00f * 1.732050776481628f * RealBuffer[1][i_bin] + 2.642218e-01f * 2.236068010330200f * RealBuffer[6][i_bin] + 4.576456e-01f * 2.236068010330200f * RealBuffer[8][i_bin];
                    *p_sec_z_real = 1.023327e+00f * 1.732050776481628f * RealBuffer[2][i_bin] - 4.576456e-01f * 2.236068010330200f * RealBuffer[5][i_bin];
                    *p_sec_w_real = 1.772454f * RealBuffer[0][i_bin] - 1.772454f * RealBuffer[1][i_bin];
                    *p_sec_x_real = 1.772454f * RealBuffer[3][i_bin] - 1.023326f * RealBuffer[4][i_bin];
                    *p_sec_y_real = -0.590818f * RealBuffer[0][i_bin] + 1.772454f * RealBuffer[1][i_bin] + 0.590817f * RealBuffer[6][i_bin] + 1.023326f * RealBuffer[8][i_bin];
                    *p_sec_z_real = 1.772454f * RealBuffer[2][i_bin] - 1.023326f * RealBuffer[5][i_bin];
                }

                // active intensity