1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
# this test is meaningful only, when the firmware supports complex arrays
try:
from ulab import numpy as np
except:
import numpy as np
dtypes = (np.uint8, np.int8, np.uint16, np.int16, np.float, np.complex)
for dtype in dtypes:
a = np.array(range(4), dtype=dtype)
outtype = np.float if dtype is not np.complex else np.complex
print('\narray:\n', a)
print('\nsquare root:\n', np.sqrt(a, dtype=outtype))
b = np.array([0, 1j, 2+2j, 3-3j], dtype=np.complex)
print('\narray:\n', b)
print('\nsquare root:\n', np.sqrt(b, dtype=np.complex))
|
