bluecellulab 2.6.56__tar.gz → 2.6.58__tar.gz

{bluecellulab-2.6.56 → bluecellulab-2.6.58}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: bluecellulab
-Version: 2.6.56
+Version: 2.6.58
 Summary: Biologically detailed neural network simulations and analysis.
 Author: Blue Brain Project, EPFL
 License: Apache2.0

{bluecellulab-2.6.56 → bluecellulab-2.6.58}/bluecellulab/analysis/analysis.py

@@ -350,7 +350,6 @@ class BPAP:
         validated = True
         notes = ""
         popt_dend, popt_apic = self.fit(soma_amp, dend_amps, dend_dist, apic_amps, apic_dist)
-        logging.warning(popt_dend)
         if dend_amps is not None:
             plt.cla()
             plt.plot([0], soma_amp, '.')

{bluecellulab-2.6.56 → bluecellulab-2.6.58}/bluecellulab/validation/validation.py

@@ -25,6 +25,7 @@ from bluecellulab.analysis.analysis import compute_plot_iv_curve
 from bluecellulab.analysis.inject_sequence import run_multirecordings_stimulus
 from bluecellulab.analysis.inject_sequence import run_stimulus
 from bluecellulab.cell.core import Cell
+from bluecellulab.simulation.neuron_globals import NeuronGlobals
 from bluecellulab.stimulus.factory import IDRestTimings
 from bluecellulab.stimulus.factory import StimulusFactory
 from bluecellulab.tools import calculate_input_resistance
@@ -33,15 +34,17 @@ from bluecellulab.tools import calculate_rheobase
 logger = logging.getLogger(__name__)
 
 
-def plot_trace(recording, out_dir, fname, title):
+def plot_trace(recording, out_dir, fname, title, plot_current=True):
     """Plot a trace with inout current given a recording."""
     outpath = out_dir / fname
     fig, ax1 = plt.subplots(figsize=(10, 6))
     plt.plot(recording.time, recording.voltage, color="black")
-    current_axis = ax1.twinx()
-    current_axis.plot(recording.time, recording.current, color="gray", alpha=0.6)
-    current_axis.set_ylabel("Stimulus Current [nA]")
-    fig.suptitle(title)
+    if plot_current:
+        current_axis = ax1.twinx()
+        current_axis.plot(recording.time, recording.current, color="gray", alpha=0.6)
+        current_axis.set_ylabel("Stimulus Current [nA]")
+    if title:
+        fig.suptitle(title)
     ax1.set_xlabel("Time [ms]")
     ax1.set_ylabel("Voltage [mV]")
     fig.tight_layout()
@@ -291,15 +294,17 @@ def hyperpolarization_test(template_params, rheobase, out_dir):
         "stim_end": [IDRestTimings.PRE_DELAY.value + IDRestTimings.DURATION.value],
     }
     features_results = efel.get_feature_values([trace], ["voltage_base", "steady_state_voltage_stimend"])
-    rmp = features_results[0]["voltage_base"][0]
-    ss_voltage = features_results[0]["steady_state_voltage_stimend"][0]
-    if rmp is None or ss_voltage is None:
+    rmp = features_results[0]["voltage_base"]
+    ss_voltage = features_results[0]["steady_state_voltage_stimend"]
+    if rmp is None or len(rmp) == 0 or ss_voltage is None or len(ss_voltage) == 0:
         return {
             "name": name,
             "passed": False,
             "validation_details": "Validation failed: Could not determine RMP or steady state voltage.",
             "figures": [outpath],
         }
+    rmp = rmp[0]
+    ss_voltage = ss_voltage[0]
     hyperpol_bool = bool(ss_voltage < rmp)
 
     notes = (
@@ -408,8 +413,42 @@ def fi_test(template_params, rheobase, out_dir, spike_threshold_voltage=-30.):
     }
 
 
+def thumbnail_test(template_params, rheobase, out_dir):
+    """Thumbnail test: creating a thumbnail."""
+    stim_factory = StimulusFactory(dt=1.0)
+    step_stimulus = stim_factory.idrest(threshold_current=rheobase, threshold_percentage=130)
+    recording = run_stimulus(
+        template_params,
+        step_stimulus,
+        "soma[0]",
+        0.5,
+        add_hypamp=True,
+    )
+
+    # plotting
+    outpath = plot_trace(
+        recording,
+        out_dir,
+        fname="thumbnail.pdf",
+        title="",
+        plot_current=False
+    )
+
+    return {
+        "name": "thumbnail",
+        "passed": True,
+        "validation_details": "",
+        "figures": [outpath],
+    }
+
+
 def run_validations(
-    cell, cell_name, spike_threshold_voltage=-30, output_dir="./memodel_validation_figures"
+    cell,
+    cell_name,
+    spike_threshold_voltage=-30,
+    v_init=-80.0,
+    celsius=34.0,
+    output_dir="./memodel_validation_figures"
 ):
     """Run all the validations on the cell.
 
@@ -417,11 +456,18 @@ def run_validations(
         cell (Cell): The cell to validate.
         cell_name (str): The name of the cell, used in the output directory.
         spike_threshold_voltage (float): The voltage threshold for spike detection.
+        v_init: Initial membrane potential. Default is -80.0 mV.
+        celsius: Temperature in Celsius. Default is 34.0.
         output_dir (str): The directory to save the validation figures.
     """
     out_dir = pathlib.Path(output_dir) / cell_name
     out_dir.mkdir(parents=True, exist_ok=True)
 
+    # set initial voltage and temperature
+    neuron_globals = NeuronGlobals.get_instance()
+    neuron_globals.temperature = celsius
+    neuron_globals.v_init = v_init
+
     # get me-model properties
     holding_current = cell.hypamp if cell.hypamp else 0.0
     if cell.threshold:
@@ -471,6 +517,9 @@ def run_validations(
     # Validation 9: FI Test
     fi_test_result = fi_test(cell.template_params, rheobase, out_dir, spike_threshold_voltage)
 
+    # Validation 10: Thumbnail Test
+    thumbnail_result = thumbnail_test(cell.template_params, rheobase, out_dir)
+
     return {
         "memodel_properties": {
             "holding_current": holding_current,
@@ -485,4 +534,5 @@ def run_validations(
         "rin_test": rin_result,
         "iv_test": iv_test_result,
         "fi_test": fi_test_result,
+        "thumbnail_test": thumbnail_result,
     }

{bluecellulab-2.6.56 → bluecellulab-2.6.58}/bluecellulab.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: bluecellulab
-Version: 2.6.56
+Version: 2.6.58
 Summary: Biologically detailed neural network simulations and analysis.
 Author: Blue Brain Project, EPFL
 License: Apache2.0