In [1]:

    

addprocs(7)


    

    
Out[1]:





7-element Array{Any,1}:
 2
 3
 4
 5
 6
 7
 8

In [2]:

    

@everywhere begin
    using NetworkDiscovery
    using POMDPs
    using POMDPToolbox
end


    

In [3]:

    

rng=MersenneTwister(1)
nodes = 1000
comms = 5
probes = 30
p_inter = 0.01
p_intra = 0.3
N = 1000


    

    
Out[3]:





1000

In [4]:

    

function est_rew(policy, nodes, comms, probes, p_intra, p_inter, N)
    sum = @parallel (+) for i in 1:N
        prob_rng = MersenneTwister(i)
        sim_rng = MersenneTwister(i)
        nw = generate_network(rng, nodes, comms, p_intra, p_inter)
        pomdp = generate_problem(rng, nw, probes, 1, 100.0, 10, 10, p_intra, p_inter)
        revealed = initial_belief(pomdp)
        sim = RolloutSimulator(rng=sim_rng, initial_state=nw, initial_belief=revealed)
        simulate(sim, pomdp, policy)
    end
end


    

    
Out[4]:





est_rew (generic function with 1 method)

In [5]:

    

policy = DiscoveryHeuristic(ProbeHighestDegree(true), GuessBasedOnNeighbors(rng))


    

    
Out[5]:





DiscoveryHeuristic(ProbeHighestDegree(true),GuessBasedOnNeighbors(MersenneTwister(DSFMT_state(Int32[1749029653,1072851681,1610647787,1072862326,1841712345,1073426746,-198061126,1073322060,-156153802,1073567984  …  1977574422,1073209915,278919868,1072835605,1290372147,18858467,1815133874,-1716870370,382,0]),Uint32[0x00000001])))

In [6]:

    

@time est_rew(policy, nodes, comms, probes, p_intra, p_inter, N)


    

    




elapsed time: 29.072629997 seconds (19071132 bytes allocated)






    
Out[6]:





67400.0

In [ ]: