simulationexam/.ipynb_checkpoints/accep_reject-checkpoint.ipynb

In [2]:

%pylab inline
from network.network import is_valid_network, get_number_of_activated_antenna, pp_network
from numpy import mean
from random import Random

Populating the interactive namespace from numpy and matplotlib

In [3]:

random = Random(666)
def get_some_possible_network(n, box_size):
    number_of_square = box_size * box_size
    for i in range(n):
        yield random.randint(0, 2 ** number_of_square - 1)

In [4]:

box_size = 5
n = 1000000

valid_networks = []
for network in get_some_possible_network(n, box_size):
    if is_valid_network(network, box_size):
        valid_networks.append(network)

In [5]:

hist(list(map(get_number_of_activated_antenna, valid_networks)), bins=100)
None

In [6]:

mean(list(map(get_number_of_activated_antenna, valid_networks)))

Out [6]:

6.2328431372549016

In [ ]:

valid_networks[0]

In [ ]: