-
Notifications
You must be signed in to change notification settings - Fork 80
/
net-networkgrowth-old.py
59 lines (43 loc) · 1.39 KB
/
net-networkgrowth-old.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
import pycxsimulator
from pylab import *
import networkx as nx
m = 2
def initialize():
global time, network, maxNodeID, positions
time = 0
network = nx.Graph()
network.add_node(0)
maxNodeID = 0
positions = nx.random_layout(network)
def observe():
cla()
nx.draw(network, pos = positions)
axis('image')
title('t = ' + str(time))
def roulette(options, weights):
weightsum = float(sum(weights))
if weightsum == 0.0:
weightsum = 1.0
probabilities = [x / weightsum for x in weights]
r = random()
s = 0.0
for k in range(len(options)):
s += probabilities[k]
if r <= s:
break
return options[k]
def update():
global time, network, maxNodeID, positions
time += 1
degs = dict(network.degree())
targets = list(degs.keys()) # fixed by toshi
preferences = [d for d in degs.values()] # fixed by toshi
# the first "d" could be varied to, e.g., 1, 1/d, d**2, etc.
maxNodeID += 1
network.add_node(maxNodeID)
positions[maxNodeID] = array([normal(0, 0.1), normal(0, 0.1)])
for i in range(m):
target = roulette(targets, preferences)
network.add_edge(maxNodeID, target)
positions = nx.spring_layout(network, pos = positions, iterations = 2)
pycxsimulator.GUI().start(func=[initialize, observe, update])