You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
89 lines
3.5 KiB
89 lines
3.5 KiB
import stormpy
|
|
import stormpy.simulator
|
|
|
|
|
|
newline = "\n"
|
|
|
|
def printFrame(streetLength, lanes, carPos, pedPos):
|
|
sidewalk = "#" * (streetLength+1)
|
|
frame = list(sidewalk + (lanes-2)*("_"*(streetLength+1)) + sidewalk)
|
|
|
|
frame[(int(carPos[0])-1)*(streetLength+1)+int(carPos[1])+1] = "C"
|
|
frame[(int(pedPos[0])-1)*(streetLength+1)+int(pedPos[1])+1] = "P"
|
|
for i in range(lanes):
|
|
frame[i*(streetLength+1)] = newline
|
|
return "".join(frame)
|
|
|
|
def dtmc():
|
|
prism_program = stormpy.parse_prism_program("/media/car_pedestrian_dtmc.prism")
|
|
|
|
streetLength = int(str(prism_program.get_constant("streetLength").definition))
|
|
lanes = int(str(prism_program.get_constant("lanes").definition))
|
|
|
|
model = stormpy.build_model(prism_program)
|
|
simulator = stormpy.simulator.create_simulator(model, seed=42)
|
|
options = stormpy.BuilderOptions([])
|
|
options.set_build_state_valuations()
|
|
options.set_build_all_reward_models()
|
|
model = stormpy.build_sparse_model_with_options(prism_program, options)
|
|
print(model)
|
|
simulator = stormpy.simulator.create_simulator(model, seed=42)
|
|
simulator.set_observation_mode(stormpy.simulator.SimulatorObservationMode.PROGRAM_LEVEL)
|
|
while True:
|
|
if simulator.is_done():
|
|
print("Pedestrian arrived on the other side!")
|
|
break
|
|
observation, reward, labels = simulator.step()
|
|
print(observation, labels)
|
|
print(printFrame(streetLength, lanes, (observation["car_lane_pos"], observation["car_street_pos"]), (observation["ped_lane_pos"], observation["ped_street_pos"])))
|
|
input("")
|
|
|
|
def mdp():
|
|
formula_str = "Pmax=? [ !\"crashed\" U ped_lane_pos=lanes ];"
|
|
prism_program = stormpy.parse_prism_program("/media/car_pedestrian_mdp.prism")
|
|
|
|
formulas = stormpy.parse_properties_for_prism_program(formula_str, prism_program)
|
|
options = stormpy.BuilderOptions([p.raw_formula for p in formulas])
|
|
options.set_build_state_valuations(True)
|
|
options.set_build_choice_labels(True)
|
|
options.set_build_all_labels()
|
|
model = stormpy.build_sparse_model_with_options(prism_program, options)
|
|
|
|
streetLength = int(str(prism_program.get_constant("streetLength").definition))
|
|
lanes = int(str(prism_program.get_constant("lanes").definition))
|
|
|
|
print(stormpy.build_model(prism_program))
|
|
initial_state = model.initial_states[0]
|
|
assert initial_state == 0
|
|
result = stormpy.model_checking(model, formulas[0], extract_scheduler=True)
|
|
assert result.has_scheduler
|
|
scheduler = result.scheduler
|
|
assert scheduler.memoryless
|
|
assert scheduler.deterministic
|
|
|
|
dtmc = model.apply_scheduler(scheduler)
|
|
print(dtmc)
|
|
simulator = stormpy.simulator.create_simulator(dtmc, seed=42)
|
|
simulator.set_observation_mode(stormpy.simulator.SimulatorObservationMode.PROGRAM_LEVEL)
|
|
|
|
i=0
|
|
while True:
|
|
if simulator.is_done():
|
|
print("Pedestrian arrived on the other side!")
|
|
break
|
|
observation, reward, labels = simulator.step()
|
|
if "crashed" in labels:
|
|
print("Pedestrian did not arrive on the other side!")
|
|
break
|
|
print(observation, labels)
|
|
frame = printFrame(int(streetLength), int(lanes), (observation["car_lane_pos"], observation["car_street_pos"]), (observation["ped_lane_pos"], observation["ped_street_pos"]))
|
|
print(frame)
|
|
input("")
|
|
with open(f"/media/frame{i:03}.txt", "w") as text_file:
|
|
text_file.write(frame)
|
|
i+=1
|
|
|
|
|
|
|
|
if __name__ == '__main__':
|
|
mdp()
|