bandu 1.3.5__py3-none-any.whl → 1.3.7__py3-none-any.whl

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
bandu/translate.py CHANGED
@@ -1,38 +1,38 @@
1
- import numpy as np
2
-
3
- def TranslatePoints(
4
- points:np.ndarray, values:np.ndarray, lattice_vecs:np.ndarray
5
- )->tuple[np.ndarray, np.ndarray]:
6
- npts = len(points)
7
- # translate points along each lattice direction to ensure grid space is filled
8
- trans_points = np.zeros((27*npts,3))
9
- trans_values = np.zeros((27*npts,1))
10
- for i in range(3):
11
- even_ind = 9*i
12
- odd_ind = 9*i+1
13
- # translate points along z axis
14
- layer_pts = points + (i-1)*lattice_vecs[2,:]
15
- trans_points[even_ind*npts:(even_ind+1)*npts,:] = layer_pts
16
- # translate along +x and -x
17
- trans_points[odd_ind*npts:(odd_ind+1)*npts,:] = layer_pts + lattice_vecs[0,:]
18
- trans_points[(even_ind+2)*npts:(even_ind+3)*npts,:] = layer_pts - lattice_vecs[0,:]
19
- # along +y and -y
20
- trans_points[(odd_ind+2)*npts:(odd_ind+3)*npts,:] = layer_pts + lattice_vecs[1,:]
21
- trans_points[(even_ind+4)*npts:(even_ind+5)*npts,:] = layer_pts - lattice_vecs[1,:]
22
- # along +x+y and +x-y
23
- trans_points[(odd_ind+4)*npts:(odd_ind+5)*npts,:] = layer_pts + lattice_vecs[0,:] + lattice_vecs[1,:]
24
- trans_points[(even_ind+6)*npts:(even_ind+7)*npts,:] = layer_pts - lattice_vecs[0,:] + lattice_vecs[1,:]
25
- # along -x+y and -x-y
26
- trans_points[(odd_ind+6)*npts:(odd_ind+7)*npts,:] = layer_pts + lattice_vecs[0,:] - lattice_vecs[1,:]
27
- trans_points[(even_ind+8)*npts:(even_ind+9)*npts,:] = layer_pts - lattice_vecs[0,:] - lattice_vecs[1,:]
28
- # repeat translations for values array
29
- trans_values[even_ind*npts:(even_ind+1)*npts,:] = values
30
- trans_values[odd_ind*npts:(odd_ind+1)*npts,:] = values
31
- trans_values[(even_ind+2)*npts:(even_ind+3)*npts,:] = values
32
- trans_values[(odd_ind+2)*npts:(odd_ind+3)*npts,:] = values
33
- trans_values[(even_ind+4)*npts:(even_ind+5)*npts,:] = values
34
- trans_values[(odd_ind+4)*npts:(odd_ind+5)*npts,:] = values
35
- trans_values[(even_ind+6)*npts:(even_ind+7)*npts,:] = values
36
- trans_values[(odd_ind+6)*npts:(odd_ind+7)*npts,:] = values
37
- trans_values[(even_ind+8)*npts:(even_ind+9)*npts,:] = values
1
+ import numpy as np
2
+
3
+ def TranslatePoints(
4
+ points:np.ndarray, values:np.ndarray, lattice_vecs:np.ndarray
5
+ )->tuple[np.ndarray, np.ndarray]:
6
+ npts = len(points)
7
+ # translate points along each lattice direction to ensure grid space is filled
8
+ trans_points = np.zeros((27*npts,3))
9
+ trans_values = np.zeros((27*npts,1))
10
+ for i in range(3):
11
+ even_ind = 9*i
12
+ odd_ind = 9*i+1
13
+ # translate points along z axis
14
+ layer_pts = points + (i-1)*lattice_vecs[2,:]
15
+ trans_points[even_ind*npts:(even_ind+1)*npts,:] = layer_pts
16
+ # translate along +x and -x
17
+ trans_points[odd_ind*npts:(odd_ind+1)*npts,:] = layer_pts + lattice_vecs[0,:]
18
+ trans_points[(even_ind+2)*npts:(even_ind+3)*npts,:] = layer_pts - lattice_vecs[0,:]
19
+ # along +y and -y
20
+ trans_points[(odd_ind+2)*npts:(odd_ind+3)*npts,:] = layer_pts + lattice_vecs[1,:]
21
+ trans_points[(even_ind+4)*npts:(even_ind+5)*npts,:] = layer_pts - lattice_vecs[1,:]
22
+ # along +x+y and +x-y
23
+ trans_points[(odd_ind+4)*npts:(odd_ind+5)*npts,:] = layer_pts + lattice_vecs[0,:] + lattice_vecs[1,:]
24
+ trans_points[(even_ind+6)*npts:(even_ind+7)*npts,:] = layer_pts - lattice_vecs[0,:] + lattice_vecs[1,:]
25
+ # along -x+y and -x-y
26
+ trans_points[(odd_ind+6)*npts:(odd_ind+7)*npts,:] = layer_pts + lattice_vecs[0,:] - lattice_vecs[1,:]
27
+ trans_points[(even_ind+8)*npts:(even_ind+9)*npts,:] = layer_pts - lattice_vecs[0,:] - lattice_vecs[1,:]
28
+ # repeat translations for values array
29
+ trans_values[even_ind*npts:(even_ind+1)*npts,:] = values
30
+ trans_values[odd_ind*npts:(odd_ind+1)*npts,:] = values
31
+ trans_values[(even_ind+2)*npts:(even_ind+3)*npts,:] = values
32
+ trans_values[(odd_ind+2)*npts:(odd_ind+3)*npts,:] = values
33
+ trans_values[(even_ind+4)*npts:(even_ind+5)*npts,:] = values
34
+ trans_values[(odd_ind+4)*npts:(odd_ind+5)*npts,:] = values
35
+ trans_values[(even_ind+6)*npts:(even_ind+7)*npts,:] = values
36
+ trans_values[(odd_ind+6)*npts:(odd_ind+7)*npts,:] = values
37
+ trans_values[(even_ind+8)*npts:(even_ind+9)*npts,:] = values
38
38
  return trans_points, trans_values