model.py 20 KB

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