Nodelist support different regions (#2780)

resolve #1280

share/lib/python/neuron/rxd/node.py

@@ -141,6 +141,16 @@ def satisfies(self, condition):
             return self.region == condition
         raise RxDException("selector %r not supported for this node type" % condition)
 
+    def _safe_satisfies(self, condition):
+        """Tests if a Node satisfies a given condition.
+
+        Works the same as node.satisfies but replaces RxDException with False
+        """
+        try:
+            return self.satisfies(condition)
+        except RxDException:
+            return False
+
     @property
     def _ref_concentration(self):
         """Returns a NEURON reference to the Node's concentration

share/lib/python/neuron/rxd/nodelist.py

@@ -17,7 +17,7 @@ def __init__(self, items):
 
     def __call__(self, restriction):
         """returns a sub-NodeList consisting of nodes satisfying restriction"""
-        return NodeList([i for i in self if i.satisfies(restriction)])
+        return NodeList([i for i in self if i._safe_satisfies(restriction)])
 
     def __getitem__(self, key):
         if isinstance(key, slice):

test/rxd/test_nodelist.py

@@ -88,3 +88,13 @@ def test_only_nodes(neuron_instance):
         assert len(nl) == 2
     except TypeError:
         raise Exception("should not get here")
+
+
+def test_different_regions(neuron_nosave_instance):
+    h, rxd, _ = neuron_nosave_instance
+    sec = h.Section(name="sec")
+    cyt = rxd.Region(h.allsec())
+    ecs = rxd.Extracellular(-10, -10, -10, 10, 10, 10, dx=10)
+    k = rxd.Species([ecs, cyt], name="k", charge=1)
+    nd = k.nodes(sec)
+    assert abs(nd.x[0] - 0.5) < 1e-12