forked from Zhangtaining/cell_research
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathfreezing_sorting_analysis.py
172 lines (146 loc) · 7.1 KB
/
freezing_sorting_analysis.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
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
import sys, getopt
from tkinter import *
import threading
import time
from modules.multithread.StatusProbe import StatusProbe
from modules.multithread.SelectionSortCell import SelectionSortCell
from modules.multithread.BubbleSortCell import BubbleSortCell
from modules.multithread.InsertionSortCell import InsertionSortCell
from modules.multithread.MergeSortCell import MergeSortCell
from modules.multithread.CellGroup import CellGroup, GroupStatus
from modules.multithread.MultiThreadCell import CellStatus
from visualization.CellImage import CellImage
from visualization.CellGroupImage import CellGroupImage
import random
import copy
import numpy as np
# VALUE_LIST = [28, 34, 6, 20, 7, 89, 34, 18, 29, 51]
#VALUE_LIST = range(20,0,-1)
def create_cells_within_one_group(value_list, threadLock, status_probe, cell_type, frozen_cell_number):
if len(value_list) == 0:
return []
left_boundary = (0, 1)
right_boundary = (len(value_list) - 1, 1)
cells = []
for i in range(0, len(value_list)):
cell = None
if cell_type == 'selection':
cell = SelectionSortCell(i + 1, value_list[i], threadLock, (i, 1), cells, left_boundary, right_boundary, status_probe, disable_visualization=True)
if cell_type == 'bubble':
cell = BubbleSortCell(i + 1, value_list[i], threadLock, (i, 1), cells, left_boundary, right_boundary, status_probe, disable_visualization=True)
if cell_type == 'insertion':
cell = InsertionSortCell(i + 1, value_list[i], threadLock, (i, 1), cells, left_boundary, right_boundary, status_probe, disable_visualization=True)
# if random.random() <= bubble_sort_percentage:
# cell = BubbleSortCell(i + 1, value_list[i], threadLock, (i, 1), cells, left_boundary, right_boundary, status_probe, disable_visualization=True)
# cell = InsertionSortCell(i + 1, value_list[i], threadLock, (i, 1), cells, left_boundary, right_boundary, status_probe, disable_visualization=True)
# else:
# cell = SelectionSortCell(i + 1, value_list[i], threadLock, (i, 1), cells, left_boundary, right_boundary, status_probe, disable_visualization=True)
cells.append(cell)
if cell_type == 'insertion':
cells[0].enable_to_move = True
period = 1000000000
start_count_down = 1000000000
cell_group = CellGroup(cells, cells, 0, left_boundary, right_boundary, GroupStatus.ACTIVE, threadLock, start_count_down, period)
for cell in cells:
cell.group = cell_group
for i in range(frozen_cell_number):
cells[random.randint(0, len(cells) - 1 )].set_cell_to_freeze()
return cells, [cell_group]
def print_current_status(cells):
print([{"value": c.value, "group id": c.group.group_id, "group status": c.group.status, "cell status": c.status, "left": c.left_boundary, "right": c.right_boundary, "ideal position": c.ideal_position} for c in cells])
def is_sorted(cells):
prev_cell = cells[0]
for c in cells:
if c.value < prev_cell.value:
return False
prev_cell = c
return True
def print_status(cells):
print([{"value": c.value, "group id": c.group.group_id, "group status": c.group.status, "cell status": c.status, "left": c.left_boundary, "right": c.right_boundary} for c in cells])
def kill_all_thread(cells, groups):
for c in cells:
c.status = CellStatus.INACTIVE
for g in groups:
g.status = GroupStatus.MERGED
def fetch_current_array(cells):
return [c.value for c in cells]
def get_pos_success_rate(cells, frozen_cell_num):
current_array = fetch_current_array(cells)
expected_arr = copy.deepcopy(current_array)
expected_arr.sort()
at_expected_position = sum([1 if a == b else 0 for (a, b) in zip(expected_arr, current_array)])
return (at_expected_position - frozen_cell_num) / (len(current_array) - frozen_cell_num)
def activate(cells, cell_groups):
for cell in cells:
cell.start()
for group in cell_groups:
group.start()
def get_pass_in_args(argv):
cell_type = ""
try:
opts, args = getopt.getopt(argv, "h", ["cell_type="])
except getopt.GetoptError:
print("please specify cell type using '--cell_type='")
sys.exit(2)
for opt, arg in opts:
if opt == "-h":
print("multithread_cell_sorting.py --cell_type=<cell_type>")
sys.exit()
if opt == "--cell_type":
cell_type = arg
if not cell_type:
print("please specify cell type using '--cell_type='")
sys.exit(2)
return cell_type
def main(argv):
#cell_type = get_pass_in_args(argv)
# success_sort_cnt = 0
# pos_success_rate_arr = []
# experiment_number = 20
# frozen_cell_num = 3
# cell_type = 'bubble'
for cell_type in ['bubble']:
for frozen_cell_num in [ 1, 2, 3]:
success_sort_cnt = 0
pos_success_rate_arr = []
experiment_number = 20
for i in range(experiment_number):
sorting_list = []
for k in range(100):
sorting_list.append(random.randint(0, 10))
threadLock = threading.Lock()
random.shuffle(sorting_list)
# print(f">>>>>>>>>>>>>>>>> Prepare cells to sort for {cell_type} experiment {i + 1} <<<<<<<<<<<<<<<<<<<<")
status_probe = StatusProbe()
cells, cell_groups = create_cells_within_one_group(sorting_list, threadLock, status_probe, cell_type, frozen_cell_num)
threadLock.acquire()
# print("Activating cells...")
activate(cells, cell_groups)
threadLock.release()
# shape = canvas.create_oval(30 - 20, 30 - 20, 30 + 20, 30 + 20, fill="white")
# text = canvas.create_text(30, 30, text="3")
# print("Start sorting......")
watchdog = 2000
time_started = time.time()
while not is_sorted(cells):
if time.time() > time_started + 40:
break
# # print_current_status(cells)
time.sleep(1)
# # print_status(cells)
# watchdog -= 1
if is_sorted(cells):
success_sort_cnt+=1
threadLock.acquire()
kill_all_thread(cells, cell_groups)
threadLock.release()
pos_success_rate_arr.append(get_pos_success_rate(cells, frozen_cell_num))
# print(">>>>>>>>>>>>>>>>> Sorting complete, killed all threads. <<<<<<<<<<<<<<<<<<<<\n")
print(f">>>>>>>>>>>>>>>>> Results for {cell_type} with {frozen_cell_num} frozen cells <<<<<<<<<<<<<<<<<<<<")
print(f"Total experiments: {experiment_number}")
print(f"Succeed sort: {success_sort_cnt}")
print(f"Success rate: {success_sort_cnt / experiment_number}")
print(f"Pos success rate: {sum(pos_success_rate_arr) / len(pos_success_rate_arr)}" )
print(f">>>>>>>>>>>>>>>>>------------------------------------------<<<<<<<<<<<<<<<<<<<<\n")
if __name__ == "__main__":
main(sys.argv[1:])