model.py 19 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504
  1. import sys
  2. sys.path.append("../kernel/")
  3. sys.path.append("../state/")
  4. from modelverse_state.main import ModelverseState as MvS
  5. from modelverse_kernel.main import ModelverseKernel as MvK
  6. from pypdevs.DEVS import AtomicDEVS, CoupledDEVS
  7. from pypdevs.simulator import Simulator
  8. import json
  9. def translate(operation):
  10. return {
  11. "CN": "create_node",
  12. "CE": "create_edge",
  13. "CNV": "create_nodevalue",
  14. "CD": "create_dict",
  15. "RV": "read_value",
  16. "RO": "read_outgoing",
  17. "RI": "read_incoming",
  18. "RE": "read_edge",
  19. "RD": "read_dict",
  20. "RDN": "read_dict_node",
  21. "RDNE": "read_dict_node_edge",
  22. "RDE": "read_dict_edge",
  23. "RRD": "read_reverse_dict",
  24. "RR": "read_root",
  25. "RDK": "read_dict_keys",
  26. "DE": "delete_edge",
  27. "DN": "delete_node",
  28. }[operation]
  29. class MvSState(object):
  30. def __init__(self):
  31. self.queue = []
  32. self.output = None
  33. self.mvs = MvS("../bootstrap/bootstrap.m")
  34. self.timer = float("inf")
  35. class ModelverseState(AtomicDEVS):
  36. def __init__(self,
  37. read_root,
  38. create_node,
  39. create_edge,
  40. create_nodevalue,
  41. create_dict,
  42. read_value,
  43. read_outgoing,
  44. read_incoming,
  45. read_edge,
  46. read_dict,
  47. read_dict_keys,
  48. read_dict_edge,
  49. read_dict_node,
  50. read_dict_node_edge,
  51. read_reverse_dict,
  52. delete_node,
  53. delete_edge):
  54. AtomicDEVS.__init__(self, "MvS")
  55. self.timings = {
  56. "read_root": read_root,
  57. "create_node": create_node,
  58. "create_edge": create_edge,
  59. "create_nodevalue": create_nodevalue,
  60. "create_dict": create_dict,
  61. "read_value": read_value,
  62. "read_outgoing": read_outgoing,
  63. "read_incoming": read_incoming,
  64. "read_edge": read_edge,
  65. "read_dict": read_dict,
  66. "read_dict_keys": read_dict_keys,
  67. "read_dict_edge": read_dict_edge,
  68. "read_dict_node": read_dict_node,
  69. "read_dict_node_edge": read_dict_node_edge,
  70. "read_reverse_dict": read_reverse_dict,
  71. "delete_node": delete_node,
  72. "delete_edge": delete_edge,
  73. }
  74. self.state = MvSState()
  75. self.from_mvk = self.addInPort("from_MvK")
  76. self.to_mvk = self.addOutPort("to_MvK")
  77. def extTransition(self, inputs):
  78. self.state.timer -= self.elapsed
  79. empty = len(self.state.queue) == 0
  80. self.state.queue.extend(inputs[self.from_mvk])
  81. if empty:
  82. # First message, so set the timer
  83. # And already compute the result so it is ready to output
  84. self.state.output = []
  85. self.state.timer = 0.0
  86. for v in self.state.queue[0]:
  87. self.state.output.append(getattr(self.state.mvs, translate(v[0]))(*v[1]))
  88. self.state.timer += self.timings[translate(v[0])]()
  89. else:
  90. # Just append the message to process
  91. pass
  92. return self.state
  93. def outputFnc(self):
  94. return {self.to_mvk: [self.state.output]}
  95. def intTransition(self):
  96. self.state.queue.pop(0)
  97. self.state.output = []
  98. if len(self.state.queue) > 0:
  99. self.state.timer = 0.0
  100. # Value contains a list of operations to do
  101. # So do them and calculate how long it takes
  102. for v in self.state.queue[0]:
  103. self.state.output.append(getattr(self.state.mvs, translate(v[0]))(*v[1])[0])
  104. self.state.timer += self.timings[translate(v[0])]()
  105. else:
  106. self.state.timer = float("inf")
  107. return self.state
  108. def timeAdvance(self):
  109. return self.state.timer
  110. class MvKState(object):
  111. def __init__(self):
  112. self.mvk = None
  113. self.waiting = False
  114. self.inputs = {}
  115. self.outputs = {}
  116. self.users = []
  117. self.reply = None
  118. self.phase = None
  119. self.commands = None
  120. self.root = None
  121. self.current_user = None
  122. self.loaded_primitives = False
  123. def __str__(self):
  124. return "\nMvK: %s\n" % self.mvk + \
  125. "waiting: %s\n" % self.waiting + \
  126. "inputs: %s\n" % self.inputs + \
  127. "outputs: %s\n" % self.outputs + \
  128. "users: %s\n" % self.users + \
  129. "reply: %s\n" % self.reply + \
  130. "phase: %s\n" % self.phase + \
  131. "commands: %s\n" % self.commands + \
  132. "root: %s\n" % self.root + \
  133. "current user: %s\n" % self.current_user
  134. class ModelverseKernel(AtomicDEVS):
  135. def __init__(self):
  136. AtomicDEVS.__init__(self, "MvK")
  137. self.state = MvKState()
  138. self.from_mvi = self.addInPort("from_MvI")
  139. self.from_mvs = self.addInPort("from_MvS")
  140. self.to_mvi = self.addOutPort("to_MvI")
  141. self.to_mvs = self.addOutPort("to_MvS")
  142. def extTransition(self, inputs):
  143. if self.from_mvi in inputs:
  144. # Got input from MvI, so we queue it
  145. for inp in inputs[self.from_mvi]:
  146. username = inp[0]
  147. data = inp[1]
  148. if data is not None:
  149. self.state.inputs.setdefault(username, []).extend(data)
  150. else:
  151. self.state.outputs.setdefault(username, []).append(None)
  152. if self.from_mvs in inputs:
  153. # Got input from MvS, so we can continue processing
  154. for mvs_input in inputs[self.from_mvs]:
  155. mvs_stripped = [i[0] for i in mvs_input]
  156. if self.state.mvk is None:
  157. # No MvK, so set it with the root we have just received (or should have received)
  158. self.state.root = mvs_stripped[0]
  159. self.state.mvk = MvK(self.state.root)
  160. else:
  161. if len(mvs_stripped) == 1:
  162. self.state.reply = mvs_stripped[0]
  163. else:
  164. self.state.reply = mvs_stripped
  165. self.state.waiting = False
  166. return self.state
  167. def intTransition(self):
  168. was_empty = len(self.state.users) == 0
  169. if self.state.commands is not None:
  170. self.state.commands = None
  171. return self.state
  172. if self.state.mvk is None:
  173. # Initializing
  174. self.state.waiting = True
  175. elif not self.state.loaded_primitives:
  176. commands = self.state.mvk.execute_yields("", "load_primitives", [], self.state.reply)
  177. if commands is None:
  178. self.state.loaded_primitives = True
  179. self.state.reply = None
  180. else:
  181. self.state.waiting = True
  182. self.state.commands = commands
  183. else:
  184. # Are initialized and have work to do
  185. if len(self.state.users) == 0:
  186. # Read out new set of users first
  187. if self.state.reply is None:
  188. commands = [("RDK", [self.state.root])]
  189. else:
  190. self.state.users = self.state.reply
  191. commands = None
  192. elif self.state.phase == "init_user":
  193. if self.state.reply is None:
  194. commands = [("RV", [self.state.users[0]])]
  195. else:
  196. self.state.current_user = self.state.reply
  197. if self.state.current_user.startswith("__"):
  198. # Don't process this user and force termination of user
  199. self.state.phase = "output"
  200. commands = None
  201. elif self.state.phase == "input":
  202. # Process inputs
  203. if self.state.inputs.get(self.state.current_user, None):
  204. element_type, value = self.state.inputs[self.state.current_user][0]
  205. commands = self.state.mvk.execute_yields(self.state.current_user, "set_input", [element_type, value], self.state.reply)
  206. if commands is None:
  207. self.state.inputs[self.state.current_user].pop(0)
  208. else:
  209. commands = None
  210. elif self.state.phase == "computation":
  211. commands = self.state.mvk.execute_yields(self.state.current_user, "execute_rule", [], self.state.reply)
  212. elif self.state.phase == "output":
  213. commands = self.state.mvk.execute_yields(self.state.current_user, "get_output", [], self.state.reply)
  214. else:
  215. raise Exception("Phase: " + str(self.state.phase))
  216. # Advance phase
  217. if commands is None:
  218. if was_empty:
  219. self.state.phase = "init_user"
  220. elif self.state.phase == "init_user":
  221. self.state.phase = "input"
  222. elif self.state.phase == "input":
  223. self.state.phase = "computation"
  224. elif self.state.phase == "computation":
  225. self.state.phase = "output"
  226. elif self.state.phase == "output":
  227. self.state.users.pop(0)
  228. self.state.phase = "init_user"
  229. self.state.waiting = False
  230. self.state.reply = None
  231. else:
  232. self.state.waiting = True
  233. # Send the commands to the MvS
  234. self.state.commands = commands
  235. return self.state
  236. def outputFnc(self):
  237. if self.state.mvk is None:
  238. # Ask the root first
  239. return {self.to_mvs: [[("RR", [])]]}
  240. elif self.state.waiting:
  241. return {self.to_mvs: [self.state.commands]}
  242. return {}
  243. def timeAdvance(self):
  244. if self.state.commands is not None:
  245. return 0
  246. elif self.state.waiting:
  247. return float("inf")
  248. elif self.state.mvk is None:
  249. return 0
  250. else:
  251. return 0
  252. class MvIState():
  253. def __init__(self):
  254. self.operations = []
  255. self.output = []
  256. self.processing = []
  257. self.memory = {}
  258. self.init = True
  259. class ModelverseInterface(AtomicDEVS):
  260. def __init__(self, username, operations):
  261. AtomicDEVS.__init__(self, "MvI_%s" % username)
  262. self.state = MvIState()
  263. self.state.operations = operations
  264. self.username = username
  265. self.to_mvk = self.addOutPort("to_MvK")
  266. self.from_mvk = self.addInPort("from_MvK")
  267. def intTransition(self):
  268. self.state.init = False
  269. while self.state.operations:
  270. i = self.state.operations[0]
  271. if isinstance(i, int) and i not in self.memory:
  272. break
  273. self.state.operations.pop(0)
  274. return self.state
  275. def extTransition(self, inputs):
  276. for inp in inputs[self.from_mvk]:
  277. if inp["value"] == "None" and isinstance(self.state.processing[0], int) and self.state.processing[0] not in self.memory:
  278. self.state.memory[self.state.processing.pop(0)] = int(inp["id"])
  279. else:
  280. self.state.output.append(inp)
  281. return self.state
  282. def outputFnc(self):
  283. send = []
  284. for i in self.state.operations:
  285. if isinstance(i, int) and i not in self.memory:
  286. break
  287. elif isinstance(i, int) and i in self.memory:
  288. # Pass a reference!
  289. send.append(("R", self.memory[i]))
  290. elif not isinstance(i, int):
  291. send.append(("V", i))
  292. return {self.to_mvk: [(self.username, send)]}
  293. def timeAdvance(self):
  294. if self.state.init:
  295. return 0
  296. elif self.state.processing and (not isinstance(self.state.processing[0], int) or self.state.processing[0] in self.memory):
  297. return 0
  298. else:
  299. return float("inf")
  300. class NetworkState(object):
  301. def __init__(self):
  302. self.processing = []
  303. self.timer = float("inf")
  304. class Network(AtomicDEVS):
  305. def __init__(self, name, latency, bandwidth):
  306. AtomicDEVS.__init__(self, name)
  307. self.state = NetworkState()
  308. self.input_port = self.addInPort("input_port")
  309. self.output_port = self.addOutPort("output_port")
  310. self.latency = latency
  311. self.bandwidth = bandwidth
  312. def intTransition(self):
  313. self.state.processing.pop(0)
  314. if self.state.processing:
  315. self.state.timer = int(len(self.state.processing[0]) / float(self.bandwidth) + self.latency)
  316. else:
  317. self.state.timer = float("inf")
  318. return self.state
  319. def extTransition(self, inputs):
  320. self.state.timer -= self.elapsed
  321. if self.state.timer == float("inf"):
  322. self.state.timer = 0
  323. for v in inputs[self.input_port]:
  324. self.state.processing.append(json.dumps(v))
  325. if len(self.state.processing) > 0:
  326. self.state.timer = int(len(self.state.processing[0]) / float(self.bandwidth) + self.latency)
  327. return self.state
  328. def outputFnc(self):
  329. return {self.output_port: [json.loads(self.state.processing[0])]}
  330. def timeAdvance(self):
  331. return self.state.timer
  332. class System(CoupledDEVS):
  333. def __init__(self,
  334. username,
  335. operations,
  336. mvi2mvk_latency,
  337. mvi2mvk_bandwidth,
  338. mvk2mvs_latency,
  339. mvk2mvs_bandwidth,
  340. mvs2mvk_latency,
  341. mvs2mvk_bandwidth,
  342. mvk2mvi_latency,
  343. mvk2mvi_bandwidth,
  344. read_root,
  345. create_node,
  346. create_edge,
  347. create_nodevalue,
  348. create_dict,
  349. read_value,
  350. read_outgoing,
  351. read_incoming,
  352. read_edge,
  353. read_dict,
  354. read_dict_keys,
  355. read_dict_edge,
  356. read_dict_node,
  357. read_dict_node_edge,
  358. read_reverse_dict,
  359. delete_node,
  360. delete_edge):
  361. CoupledDEVS.__init__(self, "System")
  362. self.mvi_manager = self.addSubModel(ModelverseInterface(\
  363. username = "user_manager",
  364. operations = [username],
  365. ))
  366. self.mvi = self.addSubModel(ModelverseInterface(\
  367. username = username,
  368. operations = operations
  369. ))
  370. self.mvk = self.addSubModel(ModelverseKernel())
  371. self.mvs = self.addSubModel(ModelverseState(\
  372. read_root = read_root,
  373. create_node = create_node,
  374. create_edge = create_edge,
  375. create_nodevalue = create_nodevalue,
  376. create_dict = create_dict,
  377. read_value = read_value,
  378. read_outgoing = read_outgoing,
  379. read_incoming = read_incoming,
  380. read_edge = read_edge,
  381. read_dict = read_dict,
  382. read_dict_keys = read_dict_keys,
  383. read_dict_edge = read_dict_edge,
  384. read_dict_node = read_dict_node,
  385. read_dict_node_edge = read_dict_node_edge,
  386. read_reverse_dict = read_reverse_dict,
  387. delete_node = delete_node,
  388. delete_edge = delete_edge
  389. ))
  390. self.mvi2mvk = self.addSubModel(Network(\
  391. name = "mvi2mvk",
  392. latency = mvi2mvk_latency,
  393. bandwidth = mvi2mvk_bandwidth
  394. ))
  395. self.mvk2mvs = self.addSubModel(Network(\
  396. name = "mvk2mvs",
  397. latency = mvk2mvs_latency,
  398. bandwidth = mvk2mvs_bandwidth
  399. ))
  400. self.mvs2mvk = self.addSubModel(Network(\
  401. name = "mvs2mvk",
  402. latency = mvs2mvk_latency,
  403. bandwidth = mvs2mvk_bandwidth
  404. ))
  405. self.mvk2mvi = self.addSubModel(Network(\
  406. name = "mvk2mvi",
  407. latency = mvk2mvi_latency,
  408. bandwidth = mvk2mvi_bandwidth
  409. ))
  410. self.connectPorts(self.mvi_manager.to_mvk, self.mvk.from_mvi)
  411. self.connectPorts(self.mvi.to_mvk, self.mvi2mvk.input_port)
  412. self.connectPorts(self.mvi2mvk.output_port, self.mvk.from_mvi)
  413. self.connectPorts(self.mvk.to_mvs, self.mvk2mvs.input_port)
  414. self.connectPorts(self.mvk2mvs.output_port, self.mvs.from_mvk)
  415. self.connectPorts(self.mvs.to_mvk, self.mvs2mvk.input_port)
  416. self.connectPorts(self.mvs2mvk.output_port, self.mvk.from_mvs)
  417. self.connectPorts(self.mvk.to_mvi, self.mvk2mvi.input_port)
  418. self.connectPorts(self.mvk2mvi.output_port, self.mvi.from_mvk)
  419. operations = [
  420. "test",
  421. ]
  422. args = {
  423. "username": "test_user",
  424. "operations": operations,
  425. "mvi2mvk_latency": 1,
  426. "mvi2mvk_bandwidth": 2000,
  427. "mvk2mvs_latency": 1,
  428. "mvk2mvs_bandwidth": 2000,
  429. "mvs2mvk_latency": 1,
  430. "mvs2mvk_bandwidth": 2000,
  431. "mvk2mvi_latency": 1,
  432. "mvk2mvi_bandwidth": 2000,
  433. "read_root": lambda: 1,
  434. "create_node": lambda: 1,
  435. "create_edge": lambda: 1,
  436. "create_nodevalue": lambda: 1,
  437. "create_dict": lambda: 1,
  438. "read_value": lambda: 1,
  439. "read_outgoing": lambda: 1,
  440. "read_incoming": lambda: 1,
  441. "read_edge": lambda: 1,
  442. "read_dict": lambda: 1,
  443. "read_dict_keys": lambda: 1,
  444. "read_dict_node": lambda: 1,
  445. "read_dict_edge": lambda: 1,
  446. "read_dict_node_edge": lambda: 1,
  447. "read_reverse_dict": lambda: 1,
  448. "delete_node": lambda: 1,
  449. "delete_edge": lambda: 1,
  450. }
  451. model = System(**args)
  452. sim = Simulator(model)
  453. sim.setTerminationTime(20000)
  454. sim.setVerbose()
  455. sim.simulate()