How to Set a SeedΒΆ

To ensure reproducibility of results users can set a seed for all the random number streams that Ciw uses. This can be done using the Ciw functon ciw.seed:

>>> import ciw
>>> ciw.seed(5)

Note that due to sampling on initialisation, the seed will need to be set before the ciw.Simulation object is created.

As an example, take the following network:

>>> N = ciw.create_network(
...     Arrival_distributions=[['Exponential', 5]],
...     Service_distributions=[['Exponential', 10]],
...     Number_of_servers=[1]
... )

Now let’s run the system for 20 time units, using a seed of 1, and get the average waiting time:

>>> ciw.seed(1)
>>> Q = ciw.Simulation(N)
>>> Q.simulate_until_max_time(20)
>>> waits = [r.waiting_time for r in Q.get_all_records()]
>>> sum(waits)/len(waits)
0.0544115013161...

Using the same seed again, the exact same average waiting time result will occur:

>>> ciw.seed(1)
>>> Q = ciw.Simulation(N)
>>> Q.simulate_until_max_time(20)
>>> waits = [r.waiting_time for r in Q.get_all_records()]
>>> sum(waits)/len(waits)
0.0544115013161...

Now using a different seed, a different result will occur:

>>> ciw.seed(2)
>>> Q = ciw.Simulation(N)
>>> Q.simulate_until_max_time(20)
>>> waits = [r.waiting_time for r in Q.get_all_records()]
>>> sum(waits)/len(waits)
0.0832990490158...