delta_action = gm.update(sens[t])
if t == stime//4 and type == "normal":
genModPlot.plot_first(t)
if t == box_time and type == "normal":
genProcPlot.plot_second(t)
if len(sens[:t+1]) >= 2:
dd = np.sum(2*(sens[t-1:t+1] > 0)-1)
if dd == 0 and sens[t-1] > 0:
peaks += 1
if peaks
== peaks_max+1:
box_time = t + stime//6
if t % 1200 == 0 or t == stime - 1:
if t > box_time:
sim.set_box([0, 1.48])
sim.update(0.3*np.pi*sens[t],
0.3*np.pi*sens_model[t])
prederr.update([0.3*np.pi*sens[t], 0.3*np.pi*sens_model[t]], t)
genProcPlot.update([0.3*np.pi*sens[t],
0.3*np.pi*ampl[t],
t > box_time,
0], t)
genModPlot.update([0.3*np.pi*sens_model[t],
0.3*np.pi*ampl_model[t],
t > box_time,
sigma_s[t]*3 if type != "normal" else 0], t)
