bluecellulab 2.6.5__tar.gz → 2.6.6__tar.gz

{bluecellulab-2.6.5 → bluecellulab-2.6.6}/PKG-INFO

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

{bluecellulab-2.6.5 → bluecellulab-2.6.6}/bluecellulab/analysis/inject_sequence.py

@@ -18,6 +18,8 @@ class StimulusName(Enum):
     IDREST = auto()
     IV = auto()
     FIRE_PATTERN = auto()
+    POS_CHEOPS = auto()
+    NEG_CHEOPS = auto()
 
 
 class Recording(NamedTuple):
@@ -68,7 +70,7 @@ def run_stimulus(
     return Recording(current, voltage, time)
 
 
-def apply_multiple_step_stimuli(
+def apply_multiple_stimuli(
     cell: Cell,
     stimulus_name: StimulusName,
     amplitudes: Sequence[float],
@@ -109,6 +111,10 @@ def apply_multiple_step_stimuli(
             stimulus = stim_factory.iv(threshold_current=cell.threshold, threshold_percentage=amplitude)
         elif stimulus_name == StimulusName.FIRE_PATTERN:
             stimulus = stim_factory.fire_pattern(threshold_current=cell.threshold, threshold_percentage=amplitude)
+        elif stimulus_name == StimulusName.POS_CHEOPS:
+            stimulus = stim_factory.pos_cheops(threshold_current=cell.threshold, threshold_percentage=amplitude)
+        elif stimulus_name == StimulusName.NEG_CHEOPS:
+            stimulus = stim_factory.neg_cheops(threshold_current=cell.threshold, threshold_percentage=amplitude)
         else:
             raise ValueError("Unknown stimulus name.")
 

{bluecellulab-2.6.5 → bluecellulab-2.6.6}/bluecellulab/stimulus/factory.py

@@ -387,3 +387,71 @@ class StimulusFactory:
             threshold_current=threshold_current,
             threshold_percentage=threshold_percentage,
         )
+
+    def pos_cheops(
+        self,
+        threshold_current: float,
+        threshold_percentage: float = 300.0,
+    ) -> Stimulus:
+        """A combination of pyramid shaped Ramp stimuli with a positive
+        amplitude.
+
+        Args:
+            threshold_current: The threshold current of the Cell.
+            threshold_percentage: Percentage of desired threshold_current amplification.
+        """
+        delay = 250.0
+        ramp1_duration = 4000.0
+        ramp2_duration = 2000.0
+        ramp3_duration = 1333.0
+        inter_delay = 2000.0
+        post_delay = 250.0
+
+        amplitude = threshold_current * threshold_percentage / 100
+        result = (
+            Empty(self.dt, duration=delay)
+            + Slope(self.dt, duration=ramp1_duration, amplitude_start=0.0, amplitude_end=amplitude)
+            + Slope(self.dt, duration=ramp1_duration, amplitude_start=amplitude, amplitude_end=0.0)
+            + Empty(self.dt, duration=inter_delay)
+            + Slope(self.dt, duration=ramp2_duration, amplitude_start=0.0, amplitude_end=amplitude)
+            + Slope(self.dt, duration=ramp2_duration, amplitude_start=amplitude, amplitude_end=0.0)
+            + Empty(self.dt, duration=inter_delay)
+            + Slope(self.dt, duration=ramp3_duration, amplitude_start=0.0, amplitude_end=amplitude)
+            + Slope(self.dt, duration=ramp3_duration, amplitude_start=amplitude, amplitude_end=0.0)
+            + Empty(self.dt, duration=post_delay)
+        )
+        return result
+
+    def neg_cheops(
+        self,
+        threshold_current: float,
+        threshold_percentage: float = 300.0,
+    ) -> Stimulus:
+        """A combination of pyramid shaped Ramp stimuli with a negative
+        amplitude.
+
+        Args:
+            threshold_current: The threshold current of the Cell.
+            threshold_percentage: Percentage of desired threshold_current amplification.
+        """
+        delay = 1750.0
+        ramp1_duration = 3333.0
+        ramp2_duration = 1666.0
+        ramp3_duration = 1111.0
+        inter_delay = 2000.0
+        post_delay = 250.0
+
+        amplitude = - threshold_current * threshold_percentage / 100
+        result = (
+            Empty(self.dt, duration=delay)
+            + Slope(self.dt, duration=ramp1_duration, amplitude_start=0.0, amplitude_end=amplitude)
+            + Slope(self.dt, duration=ramp1_duration, amplitude_start=amplitude, amplitude_end=0.0)
+            + Empty(self.dt, duration=inter_delay)
+            + Slope(self.dt, duration=ramp2_duration, amplitude_start=0.0, amplitude_end=amplitude)
+            + Slope(self.dt, duration=ramp2_duration, amplitude_start=amplitude, amplitude_end=0.0)
+            + Empty(self.dt, duration=inter_delay)
+            + Slope(self.dt, duration=ramp3_duration, amplitude_start=0.0, amplitude_end=amplitude)
+            + Slope(self.dt, duration=ramp3_duration, amplitude_start=amplitude, amplitude_end=0.0)
+            + Empty(self.dt, duration=post_delay)
+        )
+        return result

{bluecellulab-2.6.5 → bluecellulab-2.6.6}/bluecellulab.egg-info/PKG-INFO

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