Search found 73 matches

by pascal
Fri Oct 11, 2019 7:34 pm
Forum: Parallel NEURON
Topic: MPI simulation hangs
Replies: 1
Views: 706

Re: MPI simulation hangs

Upon further investigation, this appears to be an installation problem. I ran the exact same code on a different computer, and it had no issues running in parallel, completing the simulation perfectly (and giving identical results to a serial simulation). Both computers were running Windows 10, and ...
by pascal
Fri Oct 11, 2019 6:09 pm
Forum: Parallel NEURON
Topic: MPI simulation hangs
Replies: 1
Views: 706

MPI simulation hangs

I have a network simulation that runs just fine in serial, and it used to also work in parallel until I re-organized the code a bit (placing different classes in different files, and wrapping the code to actually run the simulation in its own function). The simulation periodically dumps data to file...
by pascal
Fri Jul 12, 2019 12:23 pm
Forum: NEURON + Python
Topic: Using FInitializeHandler to set initial voltages for a subset of cells
Replies: 2
Views: 604

Re: Using FInitializeHandler to set initial voltages for a subset of cells

Oh silly me, I need to use

pc.gid2cell(re_gid).soma.v=-61

Thanks Ted!
by pascal
Thu Jul 11, 2019 7:24 pm
Forum: NEURON + Python
Topic: Using FInitializeHandler to set initial voltages for a subset of cells
Replies: 2
Views: 604

Using FInitializeHandler to set initial voltages for a subset of cells

I am running a simulation in which I want all cells to be initialized to -68 mV, except for a subset of cells (RE cells), which I want to start at -61 mV. I've tried this code (where 'pc' is a ParallelContext object): def set_RE_voltages(): for re_gid in net.re_gidList: pc.gid2cell(re_gid).v=-61 if(...
by pascal
Thu May 30, 2019 6:09 pm
Forum: Adding new mechanisms and functions to NEURON
Topic: Using ion_style() to treat eca as a fixed parameter and cai as a dynamical variable
Replies: 2
Views: 618

Re: Using ion_style() to treat eca as a fixed parameter and cai as a dynamical variable

Okay, here’s what I put together: from __future__ import division from neuron import h, gui import matplotlib.pyplot as plt import numpy as np soma=h.Section(name='soma') #create object 'soma,' and make object variable 'name' equal to 'soma' soma.diam = soma.L = np.sqrt(100/np.pi) soma.cm = 1.0 soma...
by pascal
Wed May 29, 2019 7:21 pm
Forum: Specific models
Topic: Destexhe kinetic scheme for I_h
Replies: 2
Views: 732

Re: Destexhe kinetic scheme for I_h

Ah, that all makes perfect sense! Thanks!
by pascal
Wed May 29, 2019 5:55 pm
Forum: Specific models
Topic: Destexhe kinetic scheme for I_h
Replies: 2
Views: 732

Destexhe kinetic scheme for I_h

I am trying to use the thalamocortical cell from the paper "Ionic mechanisms underlying synchronized oscillations and propagating waves in a model of ferret thalamic slices" (1996). The model code can be found here: https://senselab.med.yale.edu/ModelDB/ShowModel.cshtml?model=3343#tabs-2 I have run ...
by pascal
Wed May 29, 2019 2:32 pm
Forum: Adding new mechanisms and functions to NEURON
Topic: Using ion_style() to treat eca as a fixed parameter and cai as a dynamical variable
Replies: 2
Views: 618

Using ion_style() to treat eca as a fixed parameter and cai as a dynamical variable

I am porting Bazhenov's 2002 thalamocortical model of sleep-wake transitions ( https://senselab.med.yale.edu/ModelDB/showmodel.cshtml?model=28189#tabs-1 ) to NEURON. The model has a calcium current, as well as a calcium-dependent K+ current. As far as I can tell (I am trying to get the authors to re...
by pascal
Mon May 20, 2019 1:08 pm
Forum: Adding new mechanisms and functions to NEURON
Topic: Where to find squid.hoc associated with hh.mod
Replies: 8
Views: 3798

Re: Where to find squid.hoc associated with hh.mod

Right, I get that the squid.hoc file I found was wrong, but I'm talking about the following hh.mod file (which I got from the NEURON 7.4 source code): TITLE hh.mod squid sodium, potassium, and leak channels COMMENT This is the original Hodgkin-Huxley treatment for the set of sodium, potassium, and l...
by pascal
Thu May 09, 2019 1:49 pm
Forum: Adding new mechanisms and functions to NEURON
Topic: Where to find squid.hoc associated with hh.mod
Replies: 8
Views: 3798

Re: Where to find squid.hoc associated with hh.mod

Okay, I think I am getting a handle on ionic mechanisms in mod files. Those two links were very helpful. I still have a few questions: 1) When writing mod files for custom ionic currents, is there ever a situation where I would not want to include "THREADSAFE" within the NEURON block? 2) I noticed t...
by pascal
Wed Oct 10, 2018 11:04 am
Forum: Adding new mechanisms and functions to NEURON
Topic: Where to find squid.hoc associated with hh.mod
Replies: 8
Views: 3798

Re: Where to find squid.hoc associated with hh.mod

Thanks for the resources, Ted. I'll read through them and get back to you if I have any more questions.
by pascal
Tue Oct 09, 2018 5:43 pm
Forum: Adding new mechanisms and functions to NEURON
Topic: Where to find squid.hoc associated with hh.mod
Replies: 8
Views: 3798

Re: Where to find squid.hoc associated with hh.mod

Okay, I believe I found it here: https://github.com/neurodroid/neurodroid/blob/master/nrn/squid.hoc Here is the code: /* Squid giant axon ap */ load_file("stdrun.hoc") /* Print v at every time step */ proc advance() { fadvance() print "ND ", v } create axon access axon insert hh gnabar_hh = gna gkba...
by pascal
Tue Oct 09, 2018 12:13 pm
Forum: Adding new mechanisms and functions to NEURON
Topic: Where to find squid.hoc associated with hh.mod
Replies: 8
Views: 3798

Where to find squid.hoc associated with hh.mod

I am learning to write my own MOD files, and I'm starting with hh.mod as a learning tool. The comment section in hh.mod says to see squid.hoc as an example of the hh mechanism's implementation...but I cannot find squid.hoc anywhere. It's not in my installation of NEURON, and I can't find it anywhere...
by pascal
Thu Oct 04, 2018 11:02 pm
Forum: Specific models
Topic: Weird "trap" function in Bazhenov 2002 model
Replies: 4
Views: 2212

Re: Weird "trap" function in Bazhenov 2002 model

Fair enough, thanks. I just wanted a sanity check to make sure I wasn't missing something.