diff --git a/Source/Fortran/PSMatrixAlgebraModule.F90 b/Source/Fortran/PSMatrixAlgebraModule.F90
index 9cb793e2..8e05ecd3 100644
--- a/Source/Fortran/PSMatrixAlgebraModule.F90
+++ b/Source/Fortran/PSMatrixAlgebraModule.F90
@@ -540,6 +540,7 @@ FUNCTION MeasureAsymmetry(this) RESULT(norm_value)
     CALL ConjugateMatrix(tmat)
     CALL IncrementMatrix(this, tmat, alpha_in=-1.0_NTREAL)
     norm_value = MatrixNorm(tmat)
+    CALL DestructMatrix(tmat)
 
   END FUNCTION MeasureAsymmetry
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!