Merge pull request #498 from BlueBrain/arbor-v0.10

fix arbor tests with latest arbor version

Author: AurelienJaquier

.github/workflows/test.yml

@@ -12,7 +12,7 @@ jobs:
     strategy:
       matrix:
         os: [ubuntu-latest]
-        python-version: ["3.8", "3.9", "3.10", "3.11", "3.12"]
+        python-version: ["3.9", "3.10", "3.11", "3.12"]
         include:
           - os: macos-12
             python-version: "3.10"
@@ -34,7 +34,7 @@ jobs:
       run: tox
 
     - name: "Upload coverage to Codecov"
-      uses: codecov/codecov-action@v2
+      uses: codecov/codecov-action@v4
       with:
         token: ${{ secrets.CODECOV_TOKEN }}
         fail_ci_if_error: false

README.rst

@@ -85,7 +85,7 @@ News
 Requirements
 ============
 
-* `Python 3.8+ <https://www.python.org/downloads/release/python-380/>`_
+* `Python 3.9+ <https://www.python.org/downloads/release/python-390/>`_
 * `Pip <https://pip.pypa.io>`_ (installed by default in newer versions of Python)
 * `Neuron 7.4+ <http://neuron.yale.edu/>`_ (compiled with Python support)
 * `eFEL eFeature Extraction Library <https://github.com/BlueBrain/eFEL>`_ (automatically installed by pip)

bluepyopt/ephys/morphologies.py

@@ -307,6 +307,8 @@ def load(morpho_filename, replace_axon):
             morpho = arbor.load_component(morpho_filename).component
         elif morpho_suffix == '.swc':
             morpho = arbor.load_swc_arbor(morpho_filename)
+            # turn loaded_morphology into morphology type
+            morpho = morpho.morphology
         elif morpho_suffix == '.asc':
             morpho = arbor.load_asc(morpho_filename).morphology
         else:

bluepyopt/ephys/protocols.py

@@ -602,7 +602,12 @@ def instantiate_iclamp_stimuli(self, decor, use_labels=False):
         for i, stim in enumerate(self.stimuli):
             if not isinstance(stim, stimuli.SynapticStimulus):
                 if hasattr(stim, 'envelope'):
-                    arb_iclamp = arbor.iclamp(stim.envelope())
+                    envelope = stim.envelope()
+                    envelope = [
+                        (t * arbor.units.ms, curr * arbor.units.nA)
+                        for (t, curr) in envelope
+                    ]
+                    arb_iclamp = arbor.iclamp(envelope)
                 else:
                     raise ValueError('Stimulus must provide envelope method '
                                      ' or be of type NrnNetStimStimulus to be'
@@ -648,7 +653,9 @@ def instantiate_recordings(self, cell_model, use_labels=False):
 
             cell_model.probe('voltage',
                              arb_loc.ref if use_labels else arb_loc.loc,
-                             frequency=10)  # could be a parameter
+                             f"probe-{i}",
+                             # frequency could be a parameter
+                             frequency=10 * arbor.units.kHz)
 
         return cell_model
 

bluepyopt/ephys/simulators.py

@@ -372,9 +372,11 @@ def run(self, arb_cell_model, tstop=None, dt=None):
         dt = dt if dt is not None else self.dt
 
         if dt is not None:
-            return arb_cell_model.run(tfinal=tstop, dt=dt)
+            return arb_cell_model.run(
+                tfinal=tstop * arbor.units.ms, dt=dt * arbor.units.ms
+            )
         else:
-            return arb_cell_model.run(tfinal=tstop)
+            return arb_cell_model.run(tfinal=tstop * arbor.units.ms)
 
 
 class ArbSimulatorException(Exception):

bluepyopt/ephys/templates/acc/decor_acc_template.jinja2

@@ -1,5 +1,5 @@
 (arbor-component
-  (meta-data (version "0.1-dev"))
+  (meta-data (version "0.9-dev"))
   (decor
     {%- for mech, params in global_mechs.items() %}
     {%- if mech is not none %}
@@ -10,7 +10,7 @@
     {%- endif %}
     {%- else %}
     {%- for param in params %}
-    (default ({{ param.name }} {{ param.value }}))
+    (default ({{ param.name }} {{ param.value }} (scalar 1.0)))
     {%- endfor %}
     {%- endif %}
     {%- endfor %}
@@ -27,13 +27,13 @@
     {%- endif %}
     {%- else %}
     {%- for param in params %}
-    (paint {{loc.ref}} ({{ param.name }} {{ param.value }}))
+    (paint {{loc.ref}} ({{ param.name }} {{ param.value }} (scalar 1.0)))
     {%- endfor %}
     {%- endif %}
     {%- endfor %}
 
     {%- for synapse_name, mech_params in pprocess_mechs[loc].items() %}
-    (place {{loc.ref}} (synapse (mechanism "{{ mech_params.mech }}" {%- for param in mech_params.params %} ("{{ param.name }}" {{ param.value }}){%- endfor %})) "{{ synapse_name }}")
+    (place {{loc.ref}} (synapse (mechanism "{{ mech_params.mech }}" {%- for param in mech_params.params %} ("{{ param.name }}" {{ param.value }} (scalar 1.0)){%- endfor %})) "{{ synapse_name }}")
     {%- endfor %}
 
     {%- endfor %}))

bluepyopt/ephys/templates/acc/label_dict_acc_template.jinja2

@@ -1,5 +1,5 @@
 (arbor-component
-  (meta-data (version "0.1-dev"))
+  (meta-data (version "0.9-dev"))
   (label-dict
   {%- for loc, label in label_dict.items() %}
     {{ label.defn }}

bluepyopt/tests/test_ephys/test_create_acc.py

@@ -537,7 +537,7 @@ def run_short_sim(cable_cell):
     arb_cell_model.properties.catalogue.extend(arbor.bbp_catalogue(), "BBP::")
 
     # Run a very short simulation to test mechanism instantiation
-    arb_cell_model.run(tfinal=0.1)
+    arb_cell_model.run(tfinal=0.1 * arbor.units.ms)
 
 
 @pytest.mark.unit
@@ -665,6 +665,10 @@ def check_acc_dir(test_dir, ref_dir):
             with open(ref_dir_file / file) as f:
                 ref_file = f.read()
             assert ref_file == test_file
+            # check that load_component is not raising any error here
+            fpath = pathlib.Path(test_dir) / file
+            if fpath.suffix == "acc":
+                arbor.load_component(fpath).component
 
 
 @pytest.mark.unit

bluepyopt/tests/test_ephys/test_parameterscalers.py

@@ -1,7 +1,9 @@
 """Test ephys.parameterscalers"""
 
 import json
-
+import pathlib
+import tempfile
+import arbor
 
 import pytest
 
@@ -143,3 +145,33 @@ def test_parameterscalers_iexpr_generator_unsupported_attr():
                              'unsupported attribute tau.'):
         iexpr = value_scaler.acc_scale_iexpr(
             value=value, constant_formatter=lambda v: '%.9g' % v)
+
+
+@pytest.mark.unit
+def test_parameterscalers_iexpr():
+    """ephys.parameterscalers: Test iexpr"""
+    # iexpr from bluepyopt/tests/test_ephys/test_parameterscalers.py
+    iexpr = '(sub (scalar 0.62109375) ' \
+            '(mul (log (pi) ) ' \
+            '(exp (div (distance (region "soma")) ' \
+            '(scalar 0.421875) ) ) ) )'
+
+    # modified decor as in
+    # bluepyopt/tests/test_ephys/testdata/acc/simplecell/simple_cell_decor.acc
+    simple_cell_decor_with_iexpr = \
+        '(arbor-component\n' \
+        '  (meta-data (version "0.9-dev"))\n' \
+        '  (decor\n' \
+        '    (paint (region "soma") ' \
+        '(membrane-capacitance 0.01 (scalar 1.0)))\n' \
+        '    (paint (region "soma") ' \
+        '(scaled-mechanism (density (mechanism "default::hh" ' \
+        '("gnabar" 0.10299326453483033) ("gkbar" 0.027124836082684685))) ' \
+        f'("gkbar" {iexpr})))))'
+
+    with tempfile.TemporaryDirectory() as test_dir:
+        decor_filename = pathlib.Path(test_dir).joinpath("decor.acc")
+        with open(decor_filename, "w") as f:
+            f.write(simple_cell_decor_with_iexpr)
+        # check that load_component is not raising any error here
+        arbor.load_component(decor_filename).component

bluepyopt/tests/test_ephys/testdata/acc/CCell/CCell_decor.acc

@@ -1,10 +1,10 @@
 (arbor-component
-  (meta-data (version "0.1-dev"))
+  (meta-data (version "0.9-dev"))
   (decor
-    (default (gSKv3_1bar_SKv3_1 65))
-    (paint (region "soma") (gSKv3_1bar_SKv3_1 65))
-    (paint (region "soma") (gSKv3_1bar_SKv3_1 65))
+    (default (gSKv3_1bar_SKv3_1 65 (scalar 1.0)))
+    (paint (region "soma") (gSKv3_1bar_SKv3_1 65 (scalar 1.0)))
+    (paint (region "soma") (gSKv3_1bar_SKv3_1 65 (scalar 1.0)))
     (paint (region "dend") (density (mechanism "BBP::Ih")))
-    (paint (region "apic") (gSKv3_1bar_SKv3_1 65))
-    (paint (region "apic") (gSKv3_1bar_SKv3_1 65))
+    (paint (region "apic") (gSKv3_1bar_SKv3_1 65 (scalar 1.0)))
+    (paint (region "apic") (gSKv3_1bar_SKv3_1 65 (scalar 1.0)))
     (paint (region "apic") (density (mechanism "BBP::Ih")))))