Merge pull request #400 from Hjorthmedh/region_mesh_redux

Region mesh redux

examples/notebooks/README.md

@@ -19,6 +19,7 @@ Here is a collection of Jupyter Notebooks, some of the workflows are split over
 * [input_generation_example_2_frequency_vectors](input_generation_example_2_frequency_vectors.ipynb) define Poisson input with multiple start/stop times (uses [simple_network_creation](simple_network_creation.ipynb)).
 * [input_generation_example_3_correlation](input_generation_example_3_correlation.ipynb) finer control in input targeting (uses [population_unit_network](population_unit_network.ipynb))
 * [input_tuning_example](input_tuning_example.ipynb) explore what input number and frequency are good neurons, e.g to avoid depolarisation block.
+* [Adanced input tuning](../parallel/KTH_PDC/input_tuning/) using the Dardel super computer. Here both background and signal level are optimised.
 * [input_frequency](InputFrequency/InputFrequency.ipynb) time varying input frequency example
 
 ## Striatum example

examples/parallel/KTH_PDC/input_tuning/Dardel_run_input_tuning_background_dspn.job

@@ -0,0 +1,131 @@
+#!/bin/bash -l
+#SBATCH --partition=main
+#SBATCH -o log/InputTuning-%j-output.txt
+#SBATCH -e log/InputTuning-%j-error.txt
+#SBATCH -t 8:59:00
+#SBATCH --time-min=7:59:00
+#SBATCH -J InputTuning
+#SBATCH -A naiss2023-5-231
+#SBATCH --nodes=6-10
+#SBATCH --tasks-per-node=128
+#SBATCH --mem-per-cpu=930M
+#SBATCH --mail-type=ALL
+
+export SNUDDA_TUNE_NEURON="dspn"
+export SEED_LIST="[10,20,30,40,50]"
+
+# You need to point this as the directory where you created the network in
+#NETWORK_DIR=/cfs/klemming/home/${USER:0:1}/$USER/Snudda/examples/parallel/KTH_PDC/networks/test_10k
+NETWORK_DIR=networks/input_tuning_${SNUDDA_TUNE_NEURON}_background
+
+# This is used for NEURON
+export N_WORKERS=$SLURM_NTASKS
+
+# export N_WORKERS=50
+# export IPNWORKERS=$N_WORKERS
+
+# This is used for ipyparallel and Snudda
+export IPNWORKERS=50
+
+# Clear old ipyparallel
+export IPYTHONDIR="/cfs/klemming/scratch/${USER:0:1}/$USER/.ipython"
+rm -r $IPYTHONDIR
+export IPYTHON_PROFILE=default
+
+module load snic-env
+source $HOME/Snudda/snudda_env/bin/activate
+SNUDDA_DIR=/cfs/klemming/home/"${USER:0:1}"/$USER/Snudda
+
+
+# If the BasalGangliaData directory exists, then use that for our data
+if [[ -d "/cfs/klemming/home/${USER:0:1}/$USER/BasalGangliaData/data" ]]; then
+    export SNUDDA_DATA="/cfs/klemming/home/${USER:0:1}/$USER/BasalGangliaData/data"
+    echo "Setting SNUDDA_DATA to $SNUDDA_DATA"
+    rm mechanisms
+    ln -s $SNUDDA_DATA/neurons/mechanisms/ mechanisms
+else
+    echo "SNUDDA_DATA environment variable not changed (may be empty): $SNUDDA_DATA"
+    rm mechanisms
+    ln -s ../../../snudda/data/neurons/mechanisms/
+fi
+
+
+# Start ipyparallel
+#.. Start the ipcontroller
+export FI_CXI_DEFAULT_VNI=$(od -vAn -N4 -tu < /dev/urandom)
+srun -n 1 -N 1 -c 2 --exact --overlap --mem=0 ./../ipcontroller_new.sh &
+
+
+echo ">>> waiting 60s for controller to start"
+sleep 60 
+
+#.. Read in CONTROLLERIP
+CONTROLLERIP=$(<controller_ip.txt)
+
+##.. Start the engines
+echo ">>> starting ${IPNWORKERS} engines "
+echo "SLURM_JOB_NUM_NODES={SLURM_JOB_NUM_NODES}"
+
+export FI_CXI_DEFAULT_VNI=$(od -vAn -N4 -tu < /dev/urandom)
+srun -n ${IPNWORKERS} -c 2 -N ${SLURM_JOB_NUM_NODES} --exact --overlap --mem=0 ipengine \
+--location=${CONTROLLERIP} --profile=${IPYTHON_PROFILE} --mpi \
+--ipython-dir=${IPYTHONDIR}  --timeout=30.0 c.EngineFactory.max_heartbeat_misses=10  c.MPI.use='mpi4py' \
+1> ipe_${SLURM_JOBID}.out 2> ipe_${SLURM_JOBID}.err &
+
+
+echo ">>> waiting 60s for engines to start"
+sleep 30
+
+
+echo "Network dir: "$NETWORK_DIR
+
+export PATH=$SNUDDA_DIR/snudda_env/bin/:$PATH
+export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$CRAY_LD_LIBRARY_PATH
+export PYTHONPATH=$SNUDDA_DIR/snudda_env/lib/python3.9/
+
+# This will stop NEURON from failing with "can't open DISPLAY"
+unset DISPLAY
+
+# Do we need this magic?
+# export FI_CXI_DEFAULT_VNI=$(od -vAn -N4 -tu < /dev/urandom)
+# srun -n 1 python test_script.py &> test_print.txt
+
+export FI_CXI_DEFAULT_VNI=$(od -vAn -N4 -tu < /dev/urandom)
+
+# srun -n 1 python setup_input_tuning_background.py &> log/setup_network_log.txt
+srun -n 1 -N 1 --exact --overlap --mem=0 python setup_input_tuning_background.py &> log/setup_network_log.txt
+
+
+# Quit ipyparallel
+
+python ../ipcontroller_shutdown.py
+
+
+##############
+
+rm -r x86_64
+
+export CXX=CC
+export CC=cc
+export FC=ftn
+export MPICC=cc
+export MPICXX=CC
+
+CC --version
+
+echo "About to run nrnivmodl"
+which nrnivmodl
+
+# srun -n nrnivmodl mechanisms/
+
+# Do we need this magic?
+export FI_CXI_DEFAULT_VNI=$(od -vAn -N4 -tu < /dev/urandom)
+
+srun -n 1 nrnivmodl -incflags "-lltdl=/usr/lib64/libltdl.so.7 -lreadline=/lib64/libreadline.so.7 -lncurses=/lib64/libncurses.so.6.1" -loadflags "-DLTDL_LIBRARY=/usr/lib64/libltdl.so.7 -DREADLINE_LIBRARY=/lib64/libreadline.so.7 -DNCURSES_LIBRARY=/lib64/libncurses.so.6.1" mechanisms/
+
+# Do we need this magic?
+export FI_CXI_DEFAULT_VNI=$(od -vAn -N4 -tu < /dev/urandom)
+
+srun -n $N_WORKERS $SNUDDA_DIR/examples/parallel/KTH_PDC/input_tuning/x86_64/special -mpi -python $SNUDDA_DIR/snudda/input/input_tuning.py simulate $NETWORK_DIR --seed_list $SEED_LIST
+
+

examples/parallel/KTH_PDC/input_tuning/Dardel_run_input_tuning_background_fs.job

@@ -0,0 +1,131 @@
+#!/bin/bash -l
+#SBATCH --partition=main
+#SBATCH -o log/InputTuning-%j-output.txt
+#SBATCH -e log/InputTuning-%j-error.txt
+#SBATCH -t 8:59:00
+#SBATCH --time-min=7:59:00
+#SBATCH -J InputTuning
+#SBATCH -A naiss2023-5-231
+#SBATCH --nodes=6-10
+#SBATCH --tasks-per-node=128
+#SBATCH --mem-per-cpu=930M
+#SBATCH --mail-type=ALL
+
+export SNUDDA_TUNE_NEURON="fs"
+export SEED_LIST="[10,20,30,40,50]"
+
+# You need to point this as the directory where you created the network in
+#NETWORK_DIR=/cfs/klemming/home/${USER:0:1}/$USER/Snudda/examples/parallel/KTH_PDC/networks/test_10k
+NETWORK_DIR=networks/input_tuning_${SNUDDA_TUNE_NEURON}_background
+
+# This is used for NEURON
+export N_WORKERS=$SLURM_NTASKS
+
+# export N_WORKERS=50
+# export IPNWORKERS=$N_WORKERS
+
+# This is used for ipyparallel and Snudda
+export IPNWORKERS=50
+
+# Clear old ipyparallel
+export IPYTHONDIR="/cfs/klemming/scratch/${USER:0:1}/$USER/.ipython"
+rm -r $IPYTHONDIR
+export IPYTHON_PROFILE=default
+
+module load snic-env
+source $HOME/Snudda/snudda_env/bin/activate
+SNUDDA_DIR=/cfs/klemming/home/"${USER:0:1}"/$USER/Snudda
+
+
+# If the BasalGangliaData directory exists, then use that for our data
+if [[ -d "/cfs/klemming/home/${USER:0:1}/$USER/BasalGangliaData/data" ]]; then
+    export SNUDDA_DATA="/cfs/klemming/home/${USER:0:1}/$USER/BasalGangliaData/data"
+    echo "Setting SNUDDA_DATA to $SNUDDA_DATA"
+    rm mechanisms
+    ln -s $SNUDDA_DATA/neurons/mechanisms/ mechanisms
+else
+    echo "SNUDDA_DATA environment variable not changed (may be empty): $SNUDDA_DATA"
+    rm mechanisms
+    ln -s ../../../snudda/data/neurons/mechanisms/
+fi
+
+
+# Start ipyparallel
+#.. Start the ipcontroller
+export FI_CXI_DEFAULT_VNI=$(od -vAn -N4 -tu < /dev/urandom)
+srun -n 1 -N 1 -c 2 --exact --overlap --mem=0 ./../ipcontroller_new.sh &
+
+
+echo ">>> waiting 60s for controller to start"
+sleep 60 
+
+#.. Read in CONTROLLERIP
+CONTROLLERIP=$(<controller_ip.txt)
+
+##.. Start the engines
+echo ">>> starting ${IPNWORKERS} engines "
+echo "SLURM_JOB_NUM_NODES={SLURM_JOB_NUM_NODES}"
+
+export FI_CXI_DEFAULT_VNI=$(od -vAn -N4 -tu < /dev/urandom)
+srun -n ${IPNWORKERS} -c 2 -N ${SLURM_JOB_NUM_NODES} --exact --overlap --mem=0 ipengine \
+--location=${CONTROLLERIP} --profile=${IPYTHON_PROFILE} --mpi \
+--ipython-dir=${IPYTHONDIR}  --timeout=30.0 c.EngineFactory.max_heartbeat_misses=10  c.MPI.use='mpi4py' \
+1> ipe_${SLURM_JOBID}.out 2> ipe_${SLURM_JOBID}.err &
+
+
+echo ">>> waiting 60s for engines to start"
+sleep 30
+
+
+echo "Network dir: "$NETWORK_DIR
+
+export PATH=$SNUDDA_DIR/snudda_env/bin/:$PATH
+export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$CRAY_LD_LIBRARY_PATH
+export PYTHONPATH=$SNUDDA_DIR/snudda_env/lib/python3.9/
+
+# This will stop NEURON from failing with "can't open DISPLAY"
+unset DISPLAY
+
+# Do we need this magic?
+# export FI_CXI_DEFAULT_VNI=$(od -vAn -N4 -tu < /dev/urandom)
+# srun -n 1 python test_script.py &> test_print.txt
+
+export FI_CXI_DEFAULT_VNI=$(od -vAn -N4 -tu < /dev/urandom)
+
+# srun -n 1 python setup_input_tuning_background.py &> log/setup_network_log.txt
+srun -n 1 -N 1 --exact --overlap --mem=0 python setup_input_tuning_background.py &> log/setup_network_log.txt
+
+
+# Quit ipyparallel
+
+python ../ipcontroller_shutdown.py
+
+
+##############
+
+rm -r x86_64
+
+export CXX=CC
+export CC=cc
+export FC=ftn
+export MPICC=cc
+export MPICXX=CC
+
+CC --version
+
+echo "About to run nrnivmodl"
+which nrnivmodl
+
+# srun -n nrnivmodl mechanisms/
+
+# Do we need this magic?
+export FI_CXI_DEFAULT_VNI=$(od -vAn -N4 -tu < /dev/urandom)
+
+srun -n 1 nrnivmodl -incflags "-lltdl=/usr/lib64/libltdl.so.7 -lreadline=/lib64/libreadline.so.7 -lncurses=/lib64/libncurses.so.6.1" -loadflags "-DLTDL_LIBRARY=/usr/lib64/libltdl.so.7 -DREADLINE_LIBRARY=/lib64/libreadline.so.7 -DNCURSES_LIBRARY=/lib64/libncurses.so.6.1" mechanisms/
+
+# Do we need this magic?
+export FI_CXI_DEFAULT_VNI=$(od -vAn -N4 -tu < /dev/urandom)
+
+srun -n $N_WORKERS $SNUDDA_DIR/examples/parallel/KTH_PDC/input_tuning/x86_64/special -mpi -python $SNUDDA_DIR/snudda/input/input_tuning.py simulate $NETWORK_DIR --seed_list $SEED_LIST
+
+