--- a/pymedcouple
+++ b/pymedcouple
@@ -153,7 +153,7 @@
 
     Ltot = 0
     Rtot = n_minus*n_plus
-    medc_idx = (Rtot-1)//2
+    medc_idx = (Rtot - 1)//2
 
     # kth pair algorithm (Johnson & Mizoguchi)
     while Rtot - Ltot > n_plus:
@@ -177,7 +177,7 @@
         for i in xrange(n_plus - 1, -1, -1):
             while j < n_minus - 1 and h_kern(i, j) - Am > Am_eps:
                 j += 1
-            P.append(j-1)
+            P.append(j - 1)
             
         P.reverse()
 
@@ -185,7 +185,7 @@
         for i in xrange(0, n_plus):
             while j >= 0 and h_kern(i, j) - Am < -Am_eps:
                 j -= 1
-            Q.append(j+1)
+            Q.append(j + 1)
 
         # Check on which side of those bounds the desired median of
         # the whole matrix may be.