forked from 000Justin000/ssl_edge
-
Notifications
You must be signed in to change notification settings - Fork 0
/
internet.jl
53 lines (48 loc) · 1.66 KB
/
internet.jl
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
41
42
43
44
45
46
47
48
49
50
51
52
53
using Random;
using LinearAlgebra;
using SparseArrays;
using DataFrames;
using CSV;
using NetworkOP;
using Plots;
include("utils.jl");
function path(next,i,j)
if (next[i,j] == -1)
return Vector{Int64}();
else
path = Vector{Int64}([i]);
while (i != j)
i = next[i,j];
push!(path,i)
end
return path;
end
end
function read_Internet()
#---------------------------------------------------------------------------------------
df_topology = CSV.read("data/Internet/as19990110.txt", delim='\t', header=["src","dst"])
#---------------------------------------------------------------------------------------
id2num_topology = Dict{Int64,Int64}(id=>i for (i,id) in enumerate(unique(vcat(df_topology[:src],df_topology[:dst]))));
n = length(id2num_topology);
#---------------------------------------------------------------------------------------
I = Vector{Int64}();
J = Vector{Int64}();
V = Vector{Float64}();
for i in 1:size(df_topology,1)
if (df_topology[i,1] != df_topology[i,2])
#------------------------
push!(I,id2num_topology[df_topology[i,1]]);
push!(J,id2num_topology[df_topology[i,2]]);
push!(V,1.0);
#------------------------
push!(I,id2num_topology[df_topology[i,2]]);
push!(J,id2num_topology[df_topology[i,1]]);
push!(V,1.0);
#------------------------
end
end
A = sparse(I,J,V,n,n);
#---------------------------------------------------------------------------------------
FN = NetworkOP.FlowNetwork(A);
return FN, nothing, nothing
end