www.neuron.yale.edu

Dear Hines, I printed values with %.17g. In serial version the potential values are: 0 -70 0.02 -69.784595514917612 0.040000000000000001 -69.657335866956316 0.060000000000000005 -69.553107021839949 0.080000000000000002 -69.462762944888823 0.099999999999999992 -69.382011850638975 0.11999999999999998 ...

Dear Ted, I realized that there is a small problem. Printing the potential values of the cell, in the serial and in the parallel (with multisplit) version, I noticed that after a short time, the values begin to differ slightly. Initial values in serial version: 0 -70 0.02 -69.78459551 0.04 -69.65733...

For split cells on a cluster machine, to change parameters in specific sections (after the cell has been split) it is necessary to use the function section_exists. Dear Ted, was precisely this (section_exists) my solution. In fact, with more than one host, occurred this error: "section was del...

Dear Ted, I solved the problem. I have one last question: If I wanted to run a family of simulations (with a loop for: RUN 1, RUN2 ...) , how can I modify proc dorun() in file init.hoc? I'm trying following suggestions that you gave me for another program (in http://www.neuron.yale.edu/phpBB/viewtop...

Dear Ted, afeter several tests, I realized that there were small problems at compile-time. Solved these problems, now I'm trying to run the program. It runs, but there are a problem: not changing my model, in setuptime phase my proc batchrun() is invoked, and so this runs all my model. After, in run...

Dear Ted,
it RUNS!!!

But I do not know if it is normal that the cell has 29 spikes, while in the serial version has 9 spikes.
Perhaps it depends on the number of splits?

Dear Ted,
I followed your suggestions and it runs.

Now I try to run my model.

fails with a chain of errors that start with 1 nrniv: syntax error 1 in multisplit.hoc near line 11 1 if (split > 1 || pc.nhost > 1) { ^ and sure enough, there isn't anything called "pc" when this line of code is parsed. I have to dig into this more . . . Dear Ted, I solved this error cha...

Hi Ted, I am having several problems with the Multisplit. This is part of my code: objref pc pc = new ParallelContext() {load_file("loadbal.hoc")} {load_file("binfo.hoc")} if (name_declared("maxfactor") == 0) {execute("maxfactor = .3")} proc multisplit() {loca...

Hi Ted.

I'm having some problems in understanding what parts of proc multisplit() are actually needed,
and the manner in which I can integrate this procedure in my model, which is very similar to poirazi (model = 1).

Can I safely delete proc read_mcomplex()?

Thanks.

Hi Ted.

I tried to run the examples (to understand how they work), but they give errors.

However, now I'm trying to understand the proc multisplit().
If I were to have difficulties, I hope I can count on your suggestions.

Thanks, Pasquale.

Hi Ted, I solved thanks to your suggestions... Now, I'm working on single neuron multisplit, but, from the material that I found on the net, I did not understand how I could assign different sections to different CPUs efficiently. Furthermore, what happens to sections for which I don't express the m...

Hi Ted, your answer was very exhaustive. My last problem is the following: synapses between a cell and the other is placed randomly; consequently, I get different results in simulations between the serial and the parallel version (and also depending on the number of hosts used). Obviously, this is n...

Dear Ted, now I would like to collect the execution times, calculating, for each family of simulations, the average of this times, by varying: In the case of the serial code: number of cells and number of threads. In the case of the parallel code: number of cells and number of processors. If possibl...

Hi Ted, I solved thanks to: Use cnexp or derivimplicit instead of euler The euler integration method is not thread safe (it is also numerically unstable, which should be reason enough to avoid it). Instead, use cnexp if each dstate/dt equation depends linearly only on the state, or use derivimplicit...