diff --git a/src/tools.c b/src/tools.c
index 44c8fbdc39dc51c5dc6361b0a9a5936d46f63054..f31a7409c7cb967cbdc9660068298efd489a669a 100644
--- a/src/tools.c
+++ b/src/tools.c
@@ -133,7 +133,7 @@ void pairs_n2(double *dim, struct part *restrict parts, int N, int periodic) {
     // count += parts[k].icount;
     rho += parts[k].density.wcount;
     rho_min = fmin(parts[k].density.wcount, rho_min);
-    rho_min = fmax(parts[k].density.wcount, rho_max);
+    rho_max = fmax(parts[k].density.wcount, rho_max);
   }
 
   /* Dump the result. */