diff --git a/postpic/_const.py b/postpic/_const.py
index 5d367003..138b1bba 100644
--- a/postpic/_const.py
+++ b/postpic/_const.py
@@ -55,10 +55,10 @@ def cutout(m, oldextent, newextent):
     s = ()
     for dim in range(dims):
         i = 2 * dim
-        thisdimmin = round((newextent[i] - oldextent[i])
-                           / (oldextent[i + 1] - oldextent[i]) * m.shape[dim])
-        thisdimmax = round((newextent[i + 1] - oldextent[i])
-                           / (oldextent[i + 1] - oldextent[i]) * m.shape[dim])
+        thisdimmin = round((newextent[i] - oldextent[i]) /
+                           (oldextent[i + 1] - oldextent[i]) * m.shape[dim])
+        thisdimmax = round((newextent[i + 1] - oldextent[i]) /
+                           (oldextent[i + 1] - oldextent[i]) * m.shape[dim])
         s = np.append(s, slice(thisdimmin, thisdimmax))
     if len(s) == 1:
         s = s[0]
diff --git a/postpic/analyzer/fields.py b/postpic/analyzer/fields.py
index 51755b85..ac7eacfd 100644
--- a/postpic/analyzer/fields.py
+++ b/postpic/analyzer/fields.py
@@ -143,10 +143,10 @@ def Bz(self, **kwargs):
     # --- spezielle Funktionen
 
     def energydensityE(self, **kwargs):
-        ret = Field(0.5 * pc.epsilon0
-                    * (self._Ex(**kwargs) ** 2
-                       + self._Ey(**kwargs) ** 2
-                       + self._Ez(**kwargs) ** 2))
+        ret = Field(0.5 * pc.epsilon0 *
+                    (self._Ex(**kwargs) ** 2 +
+                     self._Ey(**kwargs) ** 2 +
+                     self._Ez(**kwargs) ** 2))
         ret.unit = 'J/m^3'
         ret.name = 'Energy Density Electric-Field'
         ret.shortname = 'E'
@@ -154,10 +154,10 @@ def energydensityE(self, **kwargs):
         return ret
 
     def energydensityM(self, **kwargs):
-        ret = Field(0.5 / pc.mu0
-                    * (self._Bx(**kwargs) ** 2
-                       + self._By(**kwargs) ** 2
-                       + self._Bz(**kwargs) ** 2))
+        ret = Field(0.5 / pc.mu0 *
+                    (self._Bx(**kwargs) ** 2 +
+                     self._By(**kwargs) ** 2 +
+                     self._Bz(**kwargs) ** 2))
         ret.unit = 'J/m^3'
         ret.name = 'Energy Density Magnetic-Field'
         ret.shortname = 'M'
@@ -165,14 +165,14 @@ def energydensityM(self, **kwargs):
         return ret
 
     def energydensityEM(self, **kwargs):
-        ret = Field(0.5 * pc.epsilon0
-                    * (self._Ex(**kwargs) ** 2
-                       + self._Ey(**kwargs) ** 2
-                       + self._Ez(**kwargs) ** 2)
-                    + 0.5 / pc.mu0
-                    * (self._Bx(**kwargs) ** 2
-                       + self._By(**kwargs) ** 2
-                       + self._Bz(**kwargs) ** 2))
+        ret = Field(0.5 * pc.epsilon0 *
+                    (self._Ex(**kwargs) ** 2 +
+                     self._Ey(**kwargs) ** 2 +
+                     self._Ez(**kwargs) ** 2) +
+                    0.5 / pc.mu0 *
+                    (self._Bx(**kwargs) ** 2 +
+                     self._By(**kwargs) ** 2 +
+                     self._Bz(**kwargs) ** 2))
         ret.unit = 'J/m^3'
         ret.name = 'Energy Density EM-Field'
         ret.shortname = 'EM'
diff --git a/postpic/analyzer/particles.py b/postpic/analyzer/particles.py
index e43b4933..1374a9ef 100644
--- a/postpic/analyzer/particles.py
+++ b/postpic/analyzer/particles.py
@@ -164,8 +164,8 @@ def __len__(self):  # = number of particles
 
     def gamma(self):
         return np.sqrt(1 +
-                       (self.Px() ** 2 + self.Py() ** 2 + self.Pz() ** 2)
-                       / (self.mass() * pc.c) ** 2)
+                       (self.Px() ** 2 + self.Py() ** 2 + self.Pz() ** 2) /
+                       (self.mass() * pc.c) ** 2)
 
     def mass_u(self):
         return self.mass() / pc.mass_u