Commit e8cef013 authored by thomas dettbarn's avatar thomas dettbarn
Browse files

added invVal_e to the magic_shift.

parent d0617a4f
Loading
Loading
Loading
Loading
Loading
+3 −3
Original line number Diff line number Diff line
@@ -901,7 +901,7 @@ static void HouseholderReduction_fx(
            singularVectors_Left_64[jCh][iCh] = W_shr( W_deposit32_h( singularVectors_Left_fx[jCh][iCh] ), 32 );
        }
    }
    for ( nCh = 0; nCh < nChannelsC; nCh++ )
    FOR ( nCh = 0; nCh < nChannelsC; nCh++ )
    {
        biDiagonalReductionLeft_64(
            singularVectors_Left_64,
@@ -1085,7 +1085,7 @@ static void biDiagonalReductionLeft_64(
            norm_x_e = W_norm( norm_64 );
            norm_x = W_extract_h( W_shl( norm_64, norm_x_e ) );
            f = Mpy_32_32( norm_x, invVal );
            magic_shift = ( norm_x_e - 2 * HEADROOM_LEFT_3 ) - ( r_e - HEADROOM_LEFT_1 ) + ( 32 - HEADROOM_LEFT_4 );
            magic_shift = ( norm_x_e - 2 * HEADROOM_LEFT_3 ) - ( r_e - 2 * HEADROOM_LEFT_1 ) + ( 32 - HEADROOM_LEFT_4 ) - 2 * invVal_e;
            FOR( jCh = currChannel; jCh < nChannelsL; jCh++ )
            {
                factor1 = W_extract_h( W_shl( singularVectors_Left_64[jCh][currChannel], 32 - HEADROOM_LEFT_4 ) );
@@ -1195,7 +1195,7 @@ static void biDiagonalReductionRight_64(
                norm_x_e = W_norm( norm_64 );
                norm_x = W_extract_h( W_shl( norm_64, norm_x_e ) );
                f = Mpy_32_32( norm_x, invVal );
                magic_shift = ( norm_x_e - 2 * HEADROOM_RIGHT_3 ) - ( r_e - HEADROOM_RIGHT_1 ) + ( 32 - HEADROOM_RIGHT_4 );
                magic_shift = ( norm_x_e - 2 * HEADROOM_RIGHT_3 ) - ( r_e - 2 * HEADROOM_RIGHT_1 ) + ( 32 - HEADROOM_RIGHT_4 ) - 2 * invVal_e;

                FOR( jCh = idx; jCh < nChannelsC; jCh++ )
                {