def create_node_graph():
ARTERY_LENGTH = 100
lines = []
for i in range(0, 21):
line = []
for j in range(0, ARTERY_LENGTH):
line.append((uniform(0.001, 0.3),uniform(0.001, 0.3),uniform(0.001, 0.3)))
lines.append(tuple(line))
#print lines[0]
line_dtype = numpy.dtype([('A', float), ('u', float), ('p', float)])
lines = [Line.create_from_sequence(line, line_dtype) for line in lines]
line_sides = []
line_sides.append([(0, Side.RIGHT), (1, Side.LEFT), (2, Side.LEFT)])
line_sides.append([(1, Side.RIGHT), (3, Side.LEFT), (4, Side.LEFT)])
line_sides.append([(4, Side.RIGHT), (5, Side.LEFT), (6, Side.LEFT)])
line_sides.append([(6, Side.RIGHT), (7, Side.LEFT), (8, Side.RIGHT)])
line_sides.append([(8, Side.LEFT), (10, Side.RIGHT), (9, Side.LEFT)])
line_sides.append([(2, Side.RIGHT), (11, Side.LEFT), (12, Side.RIGHT)])
line_sides.append([(11, Side.RIGHT), (10, Side.LEFT), (13, Side.RIGHT)])
line_sides.append([(14, Side.RIGHT), (15, Side.LEFT), (13, Side.LEFT)])
line_sides.append([(14, Side.LEFT), (16, Side.RIGHT), (17, Side.LEFT)])
line_sides.append([(17, Side.RIGHT), (19, Side.LEFT), (18, Side.LEFT)])
line_sides.append([(15, Side.RIGHT), (20, Side.LEFT), (18, Side.RIGHT)])
junctions = [Junction(line_side, 3) for line_side in line_sides]
return WebGraph(lines, junctions, boundary_size=3)
def create_node_graph():
lines = [
((1, 1, 1), (2, 2, 2), (3, 3, 3), (4, 4, 4), (5, 5, 5), (6, 6, 6), (7, 7, 7), (8, 8, 8), (9, 9, 9)),
(
(11, 11, 11),
(12, 12, 12),
(13, 13, 13),
(14, 14, 14),
(15, 15, 15),
(16, 16, 16),
(17, 17, 17),
(18, 18, 18),
(19, 19, 19),
),
(
(21, 21, 21),
(22, 22, 22),
(23, 23, 23),
(24, 24, 24),
(25, 25, 25),
(26, 26, 26),
(27, 27, 27),
(28, 28, 28),
(29, 29, 29),
),
]
line_dtype = numpy.dtype({"A": float, "u": float, "p": float}.items())
lines = [Line.create_from_sequence(line, line_dtype) for line in lines]
line_sides = [(0, Side.RIGHT), (1, Side.LEFT), (2, Side.LEFT)]
junctions = [Junction(line_sides, 3)]
return WebGraph(lines, junctions, boundary_size=3)
def __init__(self, size, dtype=float):
Line.__init__(self, size, dtype)
self.ghost_depth = 2
self.set_grid_variables()
self.set_interior()
self.set_default_neighbor_definition()
