diff --git a/src/Distribution/RealDiracMatrix.cpp b/src/Distribution/RealDiracMatrix.cpp
index 8b1cdd7bcc2714d3e84f1addef931067632457e8..c90559a392bfdf16a6c2126ab494e1c1b706be6d 100644
--- a/src/Distribution/RealDiracMatrix.cpp
+++ b/src/Distribution/RealDiracMatrix.cpp
@@ -180,7 +180,6 @@ void RealDiracMatrix::transform(matrix_t& M, short i, double phi,
 {
     if (phi) {  // if phi == 0 --> nothing happens, since R and invR would be identity_matrix matrix
         sparse_matrix_t R(4,4), invR(4,4);
-        sparse_matrix_t I = boost::numeric::ublas::identity_matrix<double>(4);
 
         transform(i, phi, R, invR);