Merge pull request #451 from arnaudon/fix_location

Fix location if section list is empty

Author: wvangeit

bluepyopt/ephys/locations.py

@@ -521,10 +521,15 @@ def set_sec_index(self, icell=None):
                 for section in getattr(icell, self.seclist_name)
             ]
         )
-        if self.direction == 'radial':
-            self.sec_index = int(np.argmax(np.linalg.norm(points, axis=1)))
+        if len(points) > 0:
+            if self.direction == 'radial':
+                self.sec_index = int(np.argmax(np.linalg.norm(points, axis=1)))
+            else:
+                self.sec_index = int(np.argmax(points.dot(self.direction)))
         else:
-            self.sec_index = int(np.argmax(points.dot(self.direction)))
+            raise EPhysLocInstantiateException(
+                "Empty seclist: %s" % self.seclist_name
+            )
 
     def instantiate(self, sim=None, icell=None):
         """ """