DEMO its numerical precision

hnn_core/cell.py

@@ -158,6 +158,24 @@ def create_soma(self):
         self.soma.Ra = soma_props['Ra']
         self.soma.cm = soma_props['cm']
 
+    def move_to_pos(self):
+        """Move cell to position."""
+        x0 = self.soma.x3d(0)
+        y0 = self.soma.y3d(0)
+        z0 = self.soma.z3d(0)
+        dx = self.pos[0] * 100 - x0
+        dy = self.pos[2] - y0
+        dz = self.pos[1] * 100 - z0
+
+        for s in self.get_sections():
+            for i in range(s.n3d()):
+                h.pt3dchange(i, s.x3d(i) + dx, s.y3d(i) + dy,
+                             s.z3d(i) + dz, s.diam3d(i),
+                             sec=s)
+                h.pt3dchange(i, s.x3d(i) - dx, s.y3d(i) - dy,
+                             s.z3d(i) - dz, s.diam3d(i),
+                             sec=s)
+
     # two things need to happen here for h:
     # 1. dipole needs to be inserted into each section
     # 2. a list needs to be created with a Dipole (Point Process) in each

hnn_core/network_builder.py

@@ -315,6 +315,7 @@ def _build(self):
         self._all_spike_gids = h.Vector()
 
         self._record_spikes()
+        self.move_cells_to_pos()
         self._connect_celltypes()
 
         if _get_rank() == 0: