-
Notifications
You must be signed in to change notification settings - Fork 0
/
utils.py
58 lines (50 loc) · 1.43 KB
/
utils.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
import numpy as np # numerical methods such as FFTs.
import csv # to load the input file.
# Return the list of values from a list of (time,value) pairs.
def y_part(list):
y = []
i = 0
while i < len(list):
y.append(float(list[i][1]))
i += 1
return y
# Return the square of the absolute value of the complex value z:
def abs2(z):
return z.real * z.real + z.imag * z.imag
# Return the absolute value of the complex value z:
def abs(z):
return np.sqrt(abs2(z))
def write_tv_seq_to_file(data, filename):
f = open(filename, 'wb')
datawriter = csv.writer(f, delimiter = '\t')
i = 0
while i < len(data):
datawriter.writerow(data[i])
i += 1
f.close()
# write the y-values to a an output file in (t,val) pairs.
def write_y_to_file(y,filename, start=0):
f = open(filename, 'wb')
datawriter = csv.writer(f, delimiter = '\t')
i = 0
while i < len(y):
datawriter.writerow([i+start, y[i]])
i += 1
f.close()
# write the y-values to a an output file in (t,val) pairs.
def write_abs_y_to_file(Y, filename):
f = open(filename, 'wb')
datawriter = csv.writer(f, delimiter = '\t')
i = 0
while i < len(Y):
datawriter.writerow([i, np.abs(Y[i])])
i += 1
f.close()
# Return the L2 norm of the series y.
def power(y):
p = 0.0
i = 0
while i < len(y):
p += y[i] * y[i]
i += 1
return np.sqrt(p / len(y))