Commit 87eac66a authored by Sandesh Venkatesh's avatar Sandesh Venkatesh
Browse files

Merge branch '3gpp_issue_1410_fix_1' into 'main' 

Fix for 3GPP issue 1410: Major Difference for Stereo input at noisy segments...

See merge request !1325
parents a43b5c15 c93ce9f4

lib_enc/ivas_core_pre_proc_front_fx.c

+10 −11
Original line number Diff line number Diff line
@@ -1855,7 +1855,7 @@ static void calculate_energy_buffer_ivas_fx( 
    const Word32 input_Fs,      /* i  : input sampling rate         Q0*/

    Word16 *enerBuffer_dft_e )

{

    Word16 i, j;

    Word16 i, j, guard_bits;

    Word64 nrg_DMX_fx[CLDFB_NO_CHANNELS_MAX];

    Word64 *p_nrg_DMX_fx;

    Word32 *pDFT_DMX_fx; /*q_DFT_DMX_fx*/
@@ -1863,8 +1863,6 @@ static void calculate_energy_buffer_ivas_fx( 
    Word16 band_res_dft_fx, chan_width_f_fx, start, stop;

    Word16 norm_nrg_DMX_fx;

    Word64 max_abs_nrg_DMX_fx;

    Flag overflow = 0;

    move32();



    max_abs_nrg_DMX_fx = 0;

    move64();
@@ -1881,22 +1879,23 @@ static void calculate_energy_buffer_ivas_fx( 
    start = 1;

    move16();



    guard_bits = add( find_guarded_bits_fx( extract_h( chan_width_bins_fx ) ), 1 ); /* Q0 */



    pDFT_DMX_fx = hCPE->hStereoDft->DFT_fx[0]; /* q_DFT_DMX_fx = 31 - hCPE->hStereoDft->DFT_fx_e[0] */

    p_nrg_DMX_fx = nrg_DMX_fx;



    *p_nrg_DMX_fx = Mpy_32_32( pDFT_DMX_fx[0], pDFT_DMX_fx[0] ); /*2 * q_DFT_DMX_fx -31 */

    *p_nrg_DMX_fx = W_shr( W_mult_32_32( pDFT_DMX_fx[0], pDFT_DMX_fx[0] ), guard_bits ); /* 2 * q_DFT_DMX_fx + 1 - guard_bits */

    move64();



    FOR( i = 0; i < no_channels; i++ )

    {

        stop = extract_l( L_shr( L_add( Mpy_32_16_1( chan_width_bins_fx, add( i, 1 ) ), 1 ), 1 ) );

        stop = extract_l( L_shr( L_add( Mpy_32_16_1( chan_width_bins_fx, add( i, 1 ) ), 1 ), 1 ) ); /* Q0 */

        FOR( j = start; j < stop; j++ )

        {

            *p_nrg_DMX_fx = W_add( *p_nrg_DMX_fx, Mpy_32_32( pDFT_DMX_fx[2 * j], pDFT_DMX_fx[2 * j] ) ); /*2 * q_DFT_DMX_fx -31 */

            *p_nrg_DMX_fx = W_add( *p_nrg_DMX_fx, W_shr( W_mult_32_32( pDFT_DMX_fx[2 * j], pDFT_DMX_fx[2 * j] ), guard_bits ) ); /* 2 * q_DFT_DMX_fx + 1 - guard_bits */

            move64();

            *p_nrg_DMX_fx = W_add( *p_nrg_DMX_fx, Mpy_32_32( pDFT_DMX_fx[2 * j + 1], pDFT_DMX_fx[2 * j + 1] ) );

            *p_nrg_DMX_fx = W_add( *p_nrg_DMX_fx, W_shr( W_mult_32_32( pDFT_DMX_fx[2 * j + 1], pDFT_DMX_fx[2 * j + 1] ), guard_bits ) ); /* 2 * q_DFT_DMX_fx + 1 - guard_bits */

            move64();

            // 2 * DFT_fx_e

        }

        IF( GT_64( W_abs( *p_nrg_DMX_fx ), max_abs_nrg_DMX_fx ) )

        {
@@ -1910,11 +1909,11 @@ static void calculate_energy_buffer_ivas_fx( 
    norm_nrg_DMX_fx = W_norm( max_abs_nrg_DMX_fx );

    FOR( i = 0; i < no_channels; i++ ) /* Consider only used channels, dependent on Fs */

    {

        enerBuffer_dft_fx[i] = W_extract_h( W_shl_o( nrg_DMX_fx[i], norm_nrg_DMX_fx, &overflow ) ) / 3; /*q_enerBuffer_dft_fx=2 * q_DFT_DMX_fx -31+ norm_nrg_DMX_fx -32*/

        enerBuffer_dft_fx[i] = Mpy_32_16_1( W_extract_h( W_shl( nrg_DMX_fx[i], norm_nrg_DMX_fx ) ), 10923 /* 1/3 in Q15 */ ); /* 2 * q_DFT_DMX_fx + 1 - guard_bits + norm_nrg_DMX_fx - 32 */

        move32();

    }



    *enerBuffer_dft_e = sub( add( 32, shl( hCPE->hStereoDft->DFT_fx_e[0], 1 ) ), norm_nrg_DMX_fx ); /*31-q_enerBuffer_dft_fx*/

    *enerBuffer_dft_e = sub( add( shl( hCPE->hStereoDft->DFT_fx_e[0], 1 ), guard_bits ), norm_nrg_DMX_fx );

    move16();



    /* Set remaining entries of enerBuffer to zero */