From a9e69575d5c6e4c5aef0154e8bb2f2e552f12d74 Mon Sep 17 00:00:00 2001 From: Johannes Hjorth Date: Thu, 18 Jan 2024 16:45:19 +0100 Subject: [PATCH] Path --- snudda/init/init.py | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/snudda/init/init.py b/snudda/init/init.py index e18f5d989..f8c9c044b 100644 --- a/snudda/init/init.py +++ b/snudda/init/init.py @@ -451,7 +451,7 @@ def add_neurons(self, name, print(tstr) sys.exit(-1) - print("Checking boundaries, to make sure P is not too high") + # print("Checking boundaries, to make sure P is not too high") x = np.zeros((8, 1)) y = np.zeros((8, 1)) z = np.zeros((8, 1)) @@ -468,10 +468,10 @@ def add_neurons(self, name, # axon_density function of x,y,z defined above p_corner = numexpr.evaluate(axon_density[1]) * (3e-6 ** 3) - for P, xx, yy, zz in zip(p_corner, x, y, z): - print(f"{name} axon density P({xx}, {yy}, {zz}) = {P}") - if (p_corner > 0.01).any(): + for P, xx, yy, zz in zip(p_corner, x, y, z): + print(f"{name} axon density P({xx}, {yy}, {zz}) = {P}") + print("Axon density too high at boundary!!") print("Please increase bounding box") sys.exit(-1) @@ -1082,7 +1082,7 @@ def define_striatum(self, volume_id="Striatum") # NGF - if os.path.isdir(NGF_dir): + if os.path.isdir(snudda_parse_path(NGF_dir, self.snudda_data)): self.add_neurons(name="NGF", neuron_dir=NGF_dir, num_neurons=self.num_NGF, volume_id="Striatum")