In [1]:

    

import sys
sys.path.insert(0,'..')
from notebookFunctions import *
# from .. import notebookFunctions
import matplotlib.pyplot as plt
%matplotlib inline
plt.rcParams['figure.figsize'] = (10,6.180)    #golden ratio
# %matplotlib notebook
%load_ext autoreload
%autoreload 2


    

In [80]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 490
location = pre + "/third_combined_expectedDistance_freeEnergy/_350-550/2d_z_qw/force_0.0/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, plot1d=True)


    

In [81]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 490
location = pre + "/third_combined_expectedDistance_freeEnergy/_350-550/2d_z_qw/force_0.0/"
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[81]:





[<matplotlib.lines.Line2D at 0x181e869a90>]






    













    

In [100]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 490
location = pre + "/second_combined_z_6_freeEnergy/_350-550/2d_z_qw/force_0.0/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, plot1d=True)


    

In [110]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 500
location = pre + "/second_z_6_freeEnergy_5/_350-550/2d_z_qw/force_0.0/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, plot1d=True)


    

In [45]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 490
location = pre + "/second_combined_z_6_freeEnergy/_350-550/2d_z_qw/force_0.0/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(12,4))


    

In [94]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 490
location = pre + "/second_combined_expectedDistance_freeEnergy/_350-550/2d_z_qw/test/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(12,4))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[94]:





[<matplotlib.lines.Line2D at 0x10b3611d0>]






    













    













    













    

In [52]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 500
location = pre + "/second_combined_expectedDistance_freeEnergy/_350-550/2d_z_qw/test2/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(12,4))


    

In [59]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 530
location = pre + "/second_combined_expectedDistance_freeEnergy/_350-550/2d_z_qw/test2/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(12,4))


    

In [60]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 530
location = pre + "/second_combined_expectedDistance_freeEnergy/_350-550/2d_z_qw/test2/"
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[60]:





[<matplotlib.lines.Line2D at 0x114565128>]






    













    

In [62]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 500
location = pre + "/second_combined_expectedDistance_freeEnergy/_350-550/2d_z_qw/only_less_than_50/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(12,4))


    

In [113]:

    

pre = "/Users/weilu/Research/server/mar_2018"
temp = 500
location = pre + "/extra_03_week/fourth_with_real_distance/_350-550/2d_z_qw/force_0.0/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(12,4))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[113]:





[<matplotlib.lines.Line2D at 0x181d6607f0>]






    













    













    













    

In [112]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 500
location = pre + "/second_combined_expectedDistanceReal/_350-550/2d_z_qw/force_0.01/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(12,4))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[112]:





[<matplotlib.lines.Line2D at 0x182033b390>]






    













    













    













    

In [87]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 500
location = pre + "/second_combined_expectedDistance_freeEnergy/_350-550/2d_z_qw/higher_than_50/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(12,4))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[87]:





[<matplotlib.lines.Line2D at 0x114bb6710>]






    













    













    













    

In [89]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 500
location = pre + "/second_combined_expectedDistanceReal/_350-550/2d_z_qw/force_0.0/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(12,4))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[89]:





[<matplotlib.lines.Line2D at 0x181da89b70>]






    













    













    













    

In [176]:

    

pre = "/Users/weilu/Research/server/mar_2018/04_week"
temp = 270
location = pre + "/sixth_with_real_distance_2/_280-350/2d_z_qw/single_temp/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(5, 4))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[176]:





[<matplotlib.lines.Line2D at 0x181dc51d30>]






    













    













    













    

In [179]:

    

pre = "/Users/weilu/Research/server/mar_2018/04_week"
temp = 270
location = pre + "/sixth_with_real_distance_2/_280-350/2d_z_qw/single_temp_with_force/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(5, 4))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[179]:





[<matplotlib.lines.Line2D at 0x1820151ac8>]






    













    













    













    

In [208]:

    

pre = "/Users/weilu/Research/server/mar_2018/04_week"
temp = 300
location = pre + "/sixth_with_real_distance_3/_280-350/2d_qw_dis/force_0.0/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(25, 2), end=(5,25))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[208]:





[<matplotlib.lines.Line2D at 0x1822dd5320>]






    













    













    













    

In [213]:

    

pre = "/Users/weilu/Research/server/mar_2018/04_week"
temp = 270
location = pre + "/sixth_with_real_distance_3/_280-350/2d_z_dis/force_0.0/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(5, 14))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[213]:





[<matplotlib.lines.Line2D at 0x18244f8d30>]






    













    













    













    

In [212]:

    

pre = "/Users/weilu/Research/server/mar_2018/04_week"
temp = 270
location = pre + "/sixth_with_real_distance_3/_280-350/2d_z_dis/force_0.2/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(5, 14))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[212]:





[<matplotlib.lines.Line2D at 0x18236f5748>]






    













    













    













    

In [ ]:

    




    

In [214]:

    

pre = "/Users/weilu/Research/server/mar_2018/04_week"
temp = 260
location = pre + "/sixth_with_real_distance_2/_280-350/2d_z_qw/block_small_bins/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(5, 14))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[214]:





[<matplotlib.lines.Line2D at 0x1826c93a20>]






    













    













    













    

In [215]:

    

pre = "/Users/weilu/Research/server/mar_2018/04_week"
temp = 260
location = pre + "/sixth_with_real_distance_2/_280-350/2d_z_qw/without_subsample_force_0.2/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(5, 14))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[215]:





[<matplotlib.lines.Line2D at 0x182658a7f0>]






    













    













    













    

In [152]:

    

pre = "/Users/weilu/Research/server/mar_2018/04_week"
temp = 260
location = pre + "/sixth_with_real_distance_2/_280-350/2d_z_qw/quick/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(5, 4))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[152]:





[<matplotlib.lines.Line2D at 0x1821dcd6a0>]






    













    













    













    

In [186]:

    

pre = "/Users/weilu/Research/server/mar_2018/04_week"
temp = 270
location = pre + "/sixth_with_real_distance_2/_280-350/2d_z_qw/another_single_temp/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(5,12))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[186]:





[<matplotlib.lines.Line2D at 0x181f61a710>]






    













    













    













    

In [185]:

    

pre = "/Users/weilu/Research/server/mar_2018/04_week"
temp = 310
location = pre + "/sixth_with_real_distance/_280-350/2d_z_qw/quick/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(25,4))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[185]:





[<matplotlib.lines.Line2D at 0x18245e3780>]






    













    













    













    

In [130]:

    

pre = "/Users/weilu/Research/server/mar_2018/extra_03_week"
temp = 480
location = pre + "/with_force_fourth_with_real_distance/_350-550/2d_z_qw/force_0.1/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(25,4))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[130]:





[<matplotlib.lines.Line2D at 0x18207e3080>]






    













    













    













    

In [124]:

    

pre = "/Users/weilu/Research/server/mar_2018/extra_03_week"
temp = 480
location = pre + "/fourth_with_real_distance/_350-550/2d_z_qw/force_0.0/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(25,4))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[124]:





[<matplotlib.lines.Line2D at 0x181fa364a8>]






    













    













    













    

In [123]:

    

pre = "/Users/weilu/Research/server/mar_2018/extra_03_week"
temp = 480
location = pre + "/fourth_with_real_distance/_350-550/2d_z_dis/force_0.0/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(25,4))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[123]:





[<matplotlib.lines.Line2D at 0x18209d3860>]






    













    













    













    

In [90]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 500
location = pre + "/second_combined_expectedDistanceReal/_350-550/2d_z_qw/k_0.01/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(12,4))
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[90]:





[<matplotlib.lines.Line2D at 0x1149bf358>]






    













    













    













    

In [83]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 500
location = pre + "/second_combined_expectedDistance_freeEnergy/_350-550/2d_z_qw/force_0.0/"
location2 = location + f"pmf-{temp}.dat"
path = shortest_path(location2, start=(12,4))


    

In [84]:

    

pre = "/Users/weilu/Research/server/mar_2018/03_week"
temp = 500
location = pre + "/second_combined_expectedDistance_freeEnergy/_350-550/2d_z_qw/force_0.0/"
location2 = location + f"evpb-{temp}.dat"
(xi,yi,zi) = plot2d(location2, zmax=100)
plt.plot(xi[path[:,1]], yi[path[:,0]], 'r.-')
plt.figure()
f_on_path = [zi[tuple(p)] for p in reversed(path)]
plt.plot(f_on_path)
# plt.ylim([0,25])


    

    




<matplotlib.colors.LinearSegmentedColormap object at 0x10e5ee7f0>






    
Out[84]:





[<matplotlib.lines.Line2D at 0x181e8a94e0>]






    













    

In [ ]: