diff --git a/CEP/PyBDSM/src/python/_version.py b/CEP/PyBDSM/src/python/_version.py
index 20ea806bcbed70e9dfbecca001e50d63a2d923bc..571659cfd7444d1748c3dd0e2c741affb55d5a26 100644
--- a/CEP/PyBDSM/src/python/_version.py
+++ b/CEP/PyBDSM/src/python/_version.py
@@ -32,6 +32,8 @@ def changelog():
     2012/10/04 - Fixed a bug in the mean map calculation that caused mean
         maps with constant values (i.e., non-2D maps) to have values of
         0.0 Jy/beam unless "mean_map = 'const'" was explicitly specified.
+        Fixed a bug in Gaussian fitting that could cause an island to be
+        skipped.
 
     2012/10/02 - Fixed a bug in the PSF vary module that resulted in
         incorrect PSF generators being used. Added an option to smooth
diff --git a/CEP/PyBDSM/src/python/gausfit.py b/CEP/PyBDSM/src/python/gausfit.py
index 6b6b85da84bdd49f8ddf4aee2896262a212a64ca..930b8e6026815eab877dd2fe89e8a41aaeb4de62 100644
--- a/CEP/PyBDSM/src/python/gausfit.py
+++ b/CEP/PyBDSM/src/python/gausfit.py
@@ -100,10 +100,6 @@ class Op_gausfit(Op):
             idx = isl.island_id
             gaul = gaus_list[idx][0]
             fgaul = gaus_list[idx][1]
-            if gaul == None:
-                gaul = []
-            if fgaul == None:
-                fgaul = []
             gaul = [Gaussian(img, par, idx, gidx)
                         for (gidx, par) in enumerate(gaul)]
 
diff --git a/CEP/PyBDSM/src/python/multi_proc.py b/CEP/PyBDSM/src/python/multi_proc.py
index 5885d6685c929f2c35a8ec04c3858232b442c3f3..2e058f27932f81b7dc7fbf71a2b259c15291186d 100644
--- a/CEP/PyBDSM/src/python/multi_proc.py
+++ b/CEP/PyBDSM/src/python/multi_proc.py
@@ -113,7 +113,6 @@ def run_tasks(procs, err_q, out_q, num):
     return numpy.concatenate(results).tolist()
 
 
-
 def parallel_map(function, sequence, numcores=None, bar=None, weights=None):
     """
     A parallelized version of the native Python map function that