notebook.community

Edit and run



In [26]:

    
using Compose
include("diffusions.jl")
include("utils.jl")
using DiffusionAlgorithms
using DiffusionTutorialUtilities
utils = DiffusionTutorialUtilities

A,xy = utils.load_graph_and_coords("U3A")
P = utils.normout!(A)
n = size(A,1)
t = 25.
seed = 1500

@time x = DiffusionAlgorithms.single_seed_stochastic_heat_kernel_series(P',t,seed)
fig = utils.cgraphplot(A,xy,0.1*log(1/(max(minimum(x),1.e-12))*x),log(x+1.e-12),
        utils._inferno_data)
img = PDF("figures/hk-diffuse.pdf",200mm,200mm)
draw(img,fig)

# find an adjacent node
d = sum((xy - repmat(xy[seed,:],n,1)).^2,2)
seed2 = sortperm(vec(d))[2]

@time x = DiffusionAlgorithms.single_seed_stochastic_heat_kernel_series(P',t,seed2)
fig = utils.cgraphplot(A,xy,0.1*log(1/(max(minimum(x),1.e-12))*x),log(x+1.e-12),
        utils._inferno_data)
img = PDF("figures/hk-diffuse-1.pdf",200mm,200mm)
draw(img,fig)









    



WARNING: replacing module DiffusionAlgorithms
WARNING: replacing module DiffusionTutorialUtilities






    



  0.035915 seconds (45.80 k allocations: 2.850 MB)
  0.003889 seconds (35 allocations: 639.281 KB)



In [14]:

    
using Compose
include("diffusions.jl")
include("utils.jl")
using DiffusionAlgorithms
using DiffusionTutorialUtilities
utils = DiffusionTutorialUtilities

A,xy = utils.load_graph_and_coords("newman-netsci")
P = utils.normout!(A)
n = size(A,1)
alpha=0.5
seed = 212

@time x = DiffusionAlgorithms.single_seed_pagerank_power(P',alpha,seed)
fig = utils.cgraphplot(A,xy,0.1*log(1/(max(minimum(x),1.e-8))*x),log(x+1.e-8),
        utils._inferno_data)
img = PDF("figures/ppr-diffuse.pdf",150mm,150mm)
draw(img,fig)









    



WARNING: replacing module DiffusionAlgorithms
WARNING: replacing module DiffusionTutorialUtilities






    



  0.040023 seconds (35.27 k allocations: 1.771 MB)



In [32]:

    
using Compose
include("diffusions.jl")
include("utils.jl")
using DiffusionAlgorithms
using DiffusionTutorialUtilities
utils = DiffusionTutorialUtilities

A,xy = utils.load_graph_and_coords("digits-1234-3")
P = utils.normout!(A)
n = size(A,1)
alpha=0.9
seed = 2500

@time x = DiffusionAlgorithms.single_seed_pagerank_power(P',alpha,seed)
fig = utils.cgraphplot(A,xy,0.1*log(1/(max(minimum(x),1.e-8))*x),log(x+1.e-8),
        utils._inferno_data)
img = PDF("figures/ssl-example.pdf",150mm,150mm)
draw(img,fig)









    



WARNING: replacing module DiffusionAlgorithms
WARNING: replacing module DiffusionTutorialUtilities






    



  0.126362 seconds (35.28 k allocations: 2.193 MB)