Update app.py

app.py

@@ -831,7 +831,7 @@ def lvad_plots(Rm, Ra, Emax, Emin, Vd, Tc, start_v, beta, loop_simulated=True):
     new_ef, pao_ed, pao_es, CO, MAP, Vlvs, Plvs = f_lvad(Rs, Rm, Ra, Rc, Ca, Cs, Cr, Ls, Emin, Vd, Tc, start_v, Emax, c, slope0, w0, x60, y00, y01, y02, y03, y04)
 
     # Create the figure and the loop that we will manipulate
-    N = 100
+    N = 70
     
     fig, ax = plt.subplots(figsize=(6, 4))
     plt.ylim((0,220))
@@ -839,14 +839,14 @@ def lvad_plots(Rm, Ra, Emax, Emin, Vd, Tc, start_v, beta, loop_simulated=True):
     start = (N-2)*60000
     end = (N)*60000
     if loop_simulated:
-        line1, = ax.plot(Vlv0[start:end], Plv0[start:end], lw=1, color='b',label='No LVAD')
+        line1 = ax.plot(Vlv0[start:end], Plv0[start:end], lw=1, color='b',label='No LVAD')
         point1 = ax.scatter(Vlv0[start:end], Plv0[start:end], c="b", s=5)#, label='End Diastole')
-        line2, = ax.plot(Vlvs[start:(start+1)], Plvs[start:(start+1)], lw=1, color=(78/255, 192/255, 44/255), label=f"LVAD, ω(0)= {round(w0,2)}r/min")
+        line2 = ax.plot(Vlvs[start:(start+1)], Plvs[start:(start+1)], lw=1, color=(78/255, 192/255, 44/255), label=f"LVAD, ω(0)= {round(w0,2)}r/min")
         point2 = ax.scatter(Vlvs[start:(start+1)], Plvs[start:(start+1)], color=(78/255, 192/255, 44/255), s=5)#, label='End Diastole')
         #point = ax.scatter(volumes[start:(start+1)], pressures[start:(start+1)], c="b", s=5, label='End Systole')
     else:
-        line1, = ax.plot(Vlv0, Plv0, color='blue', label='No LVAD') #blue
-        line2, = ax.plot(Vlvs, Plvs, color=(78/255, 192/255, 44/255), label=f"LVAD, ω(0)= {round(w0,2)}r/min") #green
+        line1 = ax.plot(Vlv0, Plv0, color='blue', label='No LVAD') #blue
+        line2 = ax.plot(Vlvs, Plvs, color=(78/255, 192/255, 44/255), label=f"LVAD, ω(0)= {round(w0,2)}r/min") #green
 
     ax.set_xlabel('LV Volume (ml)')
     ax.set_ylabel('LV Pressure (mmHg)')
@@ -857,17 +857,15 @@ def lvad_plots(Rm, Ra, Emax, Emin, Vd, Tc, start_v, beta, loop_simulated=True):
     def update(frame):
         # update to add more of the loop
         end = (N-2)*60000+1000 * frame
-        x = Vlv0[start:end]
-        y = Plv0[start:end]
+        #x = Vlv0[start:end]
+        #y = Plv0[start:end]
         x2 = Vlvs[start:end]
         y2 = Plvs[start:end]
-        ax.plot(x, y, lw=1, color='b',label='No LVAD')
-        # ax.plot(x2, y2, lw=1, color=(78/255, 192/255, 44/255), label=f"LVAD, ω(0)= {round(w0,2)}r/min")
+        #ax.plot(x2, y2, lw=1, color='b',label='No LVAD')
+        ax.plot(x2, y2, lw=1, color=(78/255, 192/255, 44/255), label=f"LVAD, ω(0)= {round(w0,2)}r/min")
 
 
     plt.legend(loc='upper left', framealpha=1)
-    plt.ylim((0,220))
-    plt.xlim((0,250))
 
     if loop_simulated:
         # plt.title('', fontsize=16)