forked from st35/frustration-ODEs-modeling
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathGenerate_Random_Networks.py
69 lines (61 loc) · 1.76 KB
/
Generate_Random_Networks.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
59
60
61
62
63
64
65
66
67
68
69
def Get_Network_Features(filename):
count = 0
E = []
flag = 0
indegree = {}
outdegree = {}
actcount = 0
inhcount = 0
with open(filename, 'r') as f:
for line in f:
if count == 0:
count += 1
continue
l = line.strip().split('\t')
source = l[0].strip()
target = l[1].strip()
typ = l[2].strip()
edge = source + '-->' + target + '--->' + typ
if edge not in E:
E.append(edge)
else:
flag = 1
if source not in outdegree:
outdegree[source] = 1
else:
outdegree[source] += 1
if target not in indegree:
indegree[target] = 1
else:
indegree[target] += 1
if typ == '1':
actcount += 1
elif typ == '2':
inhcount += 1
return flag, indegree, outdegree, actcount, inhcount
import random
from datetime import datetime
random.seed(datetime.now())
import sys
networkid = sys.argv[1].strip()
N = int(sys.argv[2].strip())
from shutil import copyfile
copyfile('networkfiles/' + networkid + '.topo', 'randomnetworks/' + networkid + '/' + networkid + '_0.topo')
origflag, origin, origout, origact, originh = Get_Network_Features('randomnetworks/' + networkid + '/' + networkid + '_0.topo')
from Util_Methods import *
for i in range(N):
Generate_Random_Networks('networkfiles/' + networkid + '.topo', 'randomnetworks/' + networkid + '/' + networkid + '_' + str(i + 1) + '.topo')
flag, indegree, outdegree, actcount, inhcount = Get_Network_Features('randomnetworks/' + networkid + '/' + networkid + '_' + str(i + 1) + '.topo')
valid = 1
if flag == 1:
valid = 0
for node in origin:
if origin[node] != indegree[node]:
valid = 0
for node in origout:
if origout[node] != outdegree[node]:
valid = 0
if origact != actcount or originh != inhcount:
valid = 0
if valid == 0:
print('Gavin Belson v. Peter Gregory')