concore/humanc/plotym.py at bd545f6df486fa4ac3364309ef5b06c0d09d5797 · Olusanmi-A/concore · GitHub

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
import concore
import numpy as np
import matplotlib.pyplot as plt
import time
print("plot ym")

concore.delay = 0.005
concore.default_maxtime(150)
init_simtime_u = "[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]"
init_simtime_ym = "[0.0, 0.0, 0.0]"
ut = []
ymt = []
ym = concore.initval(init_simtime_ym)
while(concore.simtime<concore.maxtime):
    while concore.unchanged():
        ym = concore.read(1,"ym",init_simtime_ym)
    concore.write(1,"ym",ym)
    print(" ym="+str(ym))
    ymt.append(np.array(ym).T)
print("retry="+str(concore.retrycount))

#################

# plot inputs and outputs
ym1 = [x[0].item() for x in ymt]
ym2 = [x[1].item() for x in ymt]
Nsim = len(ym1)

plt.figure()
plt.subplot(211)
plt.plot(range(Nsim), ym1)
plt.ylabel('MAP (mmHg)')
plt.legend(['MAP'], loc=0)
plt.subplot(212)
plt.plot(range(Nsim), ym2)
plt.xlabel('Cycles')
plt.ylabel('HR (bpm)')
plt.legend(['HR'], loc=0)
plt.savefig("hrmap.pdf")
plt.show()