diff --git a/lib_enc/analy_sp_fx.c b/lib_enc/analy_sp_fx.c
index 4f96ea2d4872231719e0fdf89912fbdbfec04b0c..57dbb369c1c80d376f8954fb80109cea1101791b 100644
--- a/lib_enc/analy_sp_fx.c
+++ b/lib_enc/analy_sp_fx.c
@@ -634,14 +634,14 @@ void ivas_analy_sp_fx(
     move32();
     Bin_E[L_FFT - 1] = Bin_E[L_FFT - 2]; // *q_Bin_E
     move32();
-
+    Word32 add_const = 21475; // 1e-5 in Q31
     FOR( i = 0; i < L_FFT / 2; i++ )
     {
         Bin_E_old[i] = Bin_E[i]; // *q_Bin_E
         move32();
 
         /* PS[i] = ( Bin_E[i] + 1e-5f + Bin_E[i + L_FFT / 2] + 1e-5f ) / 2.0f; */
-        PS[i] = W_extract_h( W_mac_32_32( W_mult_32_32( Bin_E[i], ONE_IN_Q30 ), Bin_E[i + L_FFT / 2], ONE_IN_Q30 ) ); // *q_Bin_E
+        PS[i] = W_extract_h( W_add( W_mac_32_32( W_mult_32_32( Bin_E[i], ONE_IN_Q30 ), Bin_E[i + L_FFT / 2], ONE_IN_Q30 ), W_shr( W_shl( add_const, 32 ), sub( 31, *q_Bin_E ) ) ) ); // *q_Bin_E
         move32();
 
         /* Bin_E[i] = (float) ( 10.0f * log( PS[i] ) ); */