diff --git a/lib_com/options.h b/lib_com/options.h
index 7611427c6beb753dc94eead7a86fc21efb44cf42..852d7b8d6afb3f14b200d516ffbba85c46a3a99b 100644
--- a/lib_com/options.h
+++ b/lib_com/options.h
@@ -180,6 +180,7 @@
 #define USE_NEW_HRTF_BINARY_FILE_FORMAT                 /* Orange: to activate when decided to change the hrtf binary file format */
 #define FIX_WARNING_RENDER_CONFIG                       /* Orange: fix warning on windows build */
 #define NONBE_FIX_991_PARAMBIN_BINARY_HRTF              /* Nokia: issue #991: fix using of binary file HRTF in ParamBin (to actiate when USE_NEW_HRTF_BINARY_FILE_FORMAT and FIX_777_COMBI_RENDER_CONFIG_FILE are on  )*/
+#define NONBE_FIX_1091_PMC_LOW_SIGNAL_BURSTS            /* FhG: fix for #1091, fix limit calculation for the regularized inverse of Kx to avoid bursts in very low signals  */
 #define FIX_1053_REVERB_RECONFIGURATION                 /* Philips: issue 1053: fix for dynamic switching of acoustic environment */
 #define NONBE_FIX_SBA_SIGNALING_BITS_B                  /* FhG: issue 1061: option B: signal sba order additionally in OSBA */
 #define NONBE_FIX_ISM_XOVER_BR                          /* FhG: issue 1072: select OSBA coding method depending on number of object and bitrate */
diff --git a/lib_dec/ivas_dirac_output_synthesis_cov.c b/lib_dec/ivas_dirac_output_synthesis_cov.c
index 447a56d20802ba4ba465155bc83a9e845a6d00cb..5139a6d55c08a9329d7b89a264fb60d334713c7d 100644
--- a/lib_dec/ivas_dirac_output_synthesis_cov.c
+++ b/lib_dec/ivas_dirac_output_synthesis_cov.c
@@ -51,6 +51,13 @@
 #include "wmc_auto.h"
 #include "rom_dec.h"
 
+#ifdef NONBE_FIX_1091_PMC_LOW_SIGNAL_BURSTS
+/*-----------------------------------------------------------------------*
+ * Local constants
+ *-----------------------------------------------------------------------*/
+
+#define SQRT_EPSILON 3.16227755e-08 /* square root of EPSILON */
+#endif
 
 /*-------------------------------------------------------------------*
  * ivas_dirac_dec_output_synthesis_cov_open()
@@ -574,7 +581,11 @@ int16_t computeMixingMatrices(
      *-----------------------------------------------------------------*/
 
     maximum( svd_s_buffer, lengthCx, &limit );
+#ifdef NONBE_FIX_1091_PMC_LOW_SIGNAL_BURSTS
+    limit = (float) max( limit * reg_Sx, SQRT_EPSILON );
+#else
     limit = limit * reg_Sx + EPSILON;
+#endif
 
     for ( i = 0; i < lengthCx; ++i )
     {