surrs.copy_field(net)
net.return_seasonality(a[0], a[1], a[2])
pool = Pool(20)
net.wavelet(8, "y", cut=1, pool=pool)
net.get_adjacency_matrix(net.phase, method="MIEQQ", num_workers=0, pool=pool, use_queue=False)
pool.close()
pool.join()
data_adj_matrix = net.adjacency_matrix.copy()
surrs_adj_matrices = []
for i in range(NUM_SURR):
print("surr %d/%d computing..." % (i + 1, NUM_SURR))
pool = Pool(20)
surrs.construct_fourier_surrogates(pool=pool)
surrs.add_seasonality(a[0], a[1], a[2])
net.data = surrs.get_surr()
net.wavelet(8, "y", cut=1, pool=pool)
net.get_adjacency_matrix(net.phase, method="MIEQQ", num_workers=0, pool=pool, use_queue=False)
pool.close()
pool.join()
surrs_adj_matrices.append(net.adjacency_matrix)
import cPickle
with open("8yr-phase-scale-net-surrs-test.bin", "wb") as f:
cPickle.dump({"data": data_adj_matrix, "surrs": surrs_adj_matrices}, f, protocol=cPickle.HIGHEST_PROTOCOL)
'air',
date(1950, 1, 1),
date(2016, 1, 1),
None,
None,
None,
'monthly',
anom=False)
a = net.get_seasonality(detrend=True)
surrs.copy_field(net)
# surrs.construct_fourier_surrogates()
# surrs.add_seasonality(a[0], a[1], a[2])
for num in range(NUM_SURRS):
pool = Pool(20)
surrs.construct_fourier_surrogates(pool=pool)
surrs.add_seasonality(a[0], a[1], a[2])
net.data = surrs.get_surr()
net.wavelet(1, 'y', pool=pool, cut=1)
net.get_continuous_phase(pool=pool)
print "wavelet done"
net.get_phase_fluctuations(rewrite=True, pool=pool)
print "fluctuations done"
pool.close()
pool.join()
net.get_adjacency_matrix(net.phase_fluctuations,
method="MIEQQ",
pool=None,
use_queue=True,
num_workers=20)
