1 files changed, 1 insertions, 1 deletions
diff --git a/langevin.cpp b/langevin.cpp
index 5287cd7..3e03f76 100644
--- a/langevin.cpp
+++ b/langevin.cpp
@@ -98,7 +98,7 @@ std::complex<double> hamiltonian(const Tensor& J, const Vector& z) {
 Matrix identity(unsigned N) {
   Matrix I(N, N);
   for (unsigned i = 0; i < N; i++) {
-    I(i) = 1;
+    I(i, i) = 1;
   }
   return I;
 }