from ._node import Node
__all__ = ["Nodes"]
[docs]class Nodes:
"""This is a container class for Nodes. This uses Numpy arrays
to store a list of Node objects as a "struct of arrays".
This should improve speed of loading and access.
"""
[docs] def __init__(self, N: int = 0):
"""Create a container for up to "N" nodes"""
if N <= 0:
self._is_null = True
else:
self._is_null = False
from .utils._array import create_double_array, create_int_array
# Struct of arrays for each piece of data. See the
# Node class for information about what each variable
# holds. Code uses "-1" to represent a null value
self.label = create_int_array(N, -1)
self.begin_to = create_int_array(N, -1)
self.end_to = create_int_array(N, -1)
self.self_w = create_int_array(N, -1)
self.begin_p = create_int_array(N, -1)
self.end_p = create_int_array(N, -1)
self.self_p = create_int_array(N, -1)
self.day_foi = create_double_array(N, 0.0)
self.night_foi = create_double_array(N, 0.0)
self.play_suscept = create_double_array(N, 0.0)
self.save_play_suscept = create_double_array(N, 0.0)
self.denominator_n = create_double_array(N, 0.0)
self.denominator_d = create_double_array(N, 0.0)
self.denominator_p = create_double_array(N, 0.0)
self.denominator_pd = create_double_array(N, 0.0)
self.day_inf_prob = create_double_array(N, 0.0)
self.night_inf_prob = create_double_array(N, 0.0)
self.x = create_double_array(N, 0.0) # no good initialiser for x
self.y = create_double_array(N, 0.0) # no good initialiser for y
self.coordinates = "x/y"
def is_null(self):
return self._is_null
[docs] def __len__(self):
"""Return the number of nodes in this container"""
if self._is_null:
return 0
else:
return len(self.label)
[docs] def copy(self):
"""Return a copy of this network, using deep copies for things
which are variable, but shallow copies for things which
should stay the same (e.g. x and y). This copy assumes that
you are not changing the actual links between nodes, e.g.
it shallow-copies label, begin_to, end_to, self_w, begin_p, end_p
and self_p. If you are changing the links then you will
need to recreate these arrays.
"""
from copy import copy, deepcopy
nodes = copy(self) # shallow copy
# now deep copy the required elements
nodes.day_foi = deepcopy(self.day_foi)
nodes.night_foi = deepcopy(self.night_foi)
nodes.play_suscept = deepcopy(self.play_suscept)
nodes.save_play_suscept = deepcopy(self.save_play_suscept)
nodes.denominator_n = deepcopy(self.denominator_n)
nodes.denominator_d = deepcopy(self.denominator_d)
nodes.denominator_p = deepcopy(self.denominator_p)
nodes.denominator_pd = deepcopy(self.denominator_pd)
nodes.day_inf_prob = deepcopy(self.day_inf_prob)
nodes.night_inf_prob = deepcopy(self.night_inf_prob)
return nodes
[docs] def assert_not_null(self):
"""Assert that this collection of nodes is not null"""
assert not self.is_null()
[docs] def assert_valid_index(self, i):
"""Assert that the passed index 'i' is valid for this collection.
This will return the index with Python reverse indexing,
so that "-1" refers to the last node in the collection
"""
self.assert_not_null()
n = len(self)
if i < 0:
i = n + i
if i <= 0 or i >= n:
raise IndexError(f"Invalid node index {i}. Number of nodes "
f"in this container equals {n}")
return i
[docs] def get_index(self, i):
"""Return the mapped index for the ith node. This will raise an
IndexError if this is an invalid index. If not, it will map
a negative index into a canonical positive index
"""
self.assert_valid_index(i)
if i < 0:
return len(self) + i
else:
return i
[docs] def population(self) -> int:
"""Return the total population of the nodes
Returns
-------
int
Total population
"""
return sum(self.save_play_suscept)
[docs] def __getitem__(self, i: int):
"""Return the node at index 'i'. Note that this is a
deep copy of the node. Changing the node will not change
the data in this container. To update the data in this
container you need to use __setitem__, e.g. via the
index operator (see __setitem__)
"""
i = self.assert_valid_index(i)
try:
return Node(label=self.label[i],
begin_to=self.begin_to[i],
end_to=self.end_to[i],
self_w=self.self_w[i],
begin_p=self.begin_p[i],
end_p=self.end_p[i],
self_p=self.self_p[i],
day_foi=self.day_foi[i],
night_foi=self.night_foi[i],
play_suscept=self.play_suscept[i],
save_play_suscept=self.save_play_suscept[i],
denominator_n=self.denominator_n[i],
denominator_d=self.denominator_d[i],
denominator_p=self.denominator_p[i],
denominator_pd=self.denominator_pd[i],
x=self.x[i],
y=self.y[i])
except Exception as e:
from .utils._console import Console
Console.error(f"Error at index {i}: "
f"size = {len(self)}: {e.__class__} {e}")
raise e
[docs] def __setitem__(self, i: int, value: Node):
"""Set the item as index 'i' equal to 'value'. This deep-copies
'value' into this container
"""
i = self.assert_valid_index(i)
self.label[i] = value.label if value.label is not None else -1
self.begin_to[i] = value.begin_to if value.begin_to is not None else -1
self.end_to[i] = value.end_to if value.end_to is not None else -1
self.self_w[i] = value.self_w if value.self_w is not None else -1
self.begin_p[i] = value.begin_p if value.begin_p is not None else -1
self.end_p[i] = value.end_p if value.end_p is not None else -1
self.self_p[i] = value.self._p if value.self_p is not None else -1
self.day_foi[i] = value.day_foi if value.day_foi is not None else 0.0
self.night_foi[i] = value.night_foi if value.night_foi is not None else 0.0
self.play_suscept[i] = value.play_suscept if value.play_suscept is not None else 0.0
self.save_play_suscept[i] = value.save_play_suscept if value.save_play_suscept is not None else 0.0
self.denominator_n[i] = value.denominator_n if value.denominator_n is not None else 0.0
self.denominator_d[i] = value.denominator_d if value.denominator_d is not None else 0.0
self.denominator_p[i] = value.denominator_p if value.denominator_p is not None else 0.0
self.denominator_pd[i] = value.denominator_pd if value.denominator_pd is not None else 0.0
self.x[i] = value.x if value.x is not None else 0.0
self.y[i] = value.y if value.y is not None else 0.0
[docs] def resize(self, N: int):
"""Resize this container to hold 'N' nodes. This will expand
the container if 'N' is greater than len(self), or will
contract the container (deleting excess nodes) if 'N'
is less than len(self). This function is called typically
when you pre-allocate a large Nodes container, and then
want to reduce the size to fit the number of loaded nodes
"""
if self.is_null():
# Assign to a new Nodes object created with this size
self.__dict__ = Nodes(N).__dict__
return
elif N <= 0:
self.__dict__ = Nodes(0).__dict__
size = len(self)
if N == size:
return
from .utils._array import resize_array
self.label = resize_array(self.label, N, -1)
self.begin_to = resize_array(self.begin_to, N, -1)
self.end_to = resize_array(self.end_to, N, -1)
self.self_w = resize_array(self.self_w, N, -1)
self.begin_p = resize_array(self.begin_p, N, -1)
self.end_p = resize_array(self.end_p, N, -1)
self.self_p = resize_array(self.self_p, N, -1)
self.day_foi = resize_array(self.day_foi, N, 0.0)
self.night_foi = resize_array(self.night_foi, N, 0.0)
self.play_suscept = resize_array(self.play_suscept, N, 0.0)
self.save_play_suscept = resize_array(self.save_play_suscept, N, 0.0)
self.denominator_n = resize_array(self.denominator_n, N, 0.0)
self.denominator_d = resize_array(self.denominator_d, N, 0.0)
self.denominator_p = resize_array(self.denominator_p, N, 0.0)
self.denominator_pd = resize_array(self.denominator_pd, N, 0.0)
self.x = resize_array(self.x, N, 0.0)
self.y = resize_array(self.y, N, 0.0)
[docs] def scale_susceptibles(self, ratio: any = None,
work_ratio: any = None, play_ratio: any = None):
"""Scale the number of susceptibles in these Nodes
by the passed scale ratios. These can be values, e.g.
ratio = 2.0 will scale the total number of susceptibles
in each ward by 2.0. They can also be lists of values,
where ward[i] will be scaled by ratio[i]. They can also
be dictionaries, e.g. ward[i] scaled by ratio[i]
Parameters
----------
ratio: None, float, list or dict
The amount by which to scale the total population of
susceptibles - evenly scales the work and play populations
work_ratio: None, float, list or dict
Scale only the work population of susceptibles
play_ratio: None, float, list or dict
Scale only the play population of susceptibles
Returns
-------
None
"""
from .utils._scale_susceptibles import scale_node_susceptibles
scale_node_susceptibles(nodes=self, ratio=ratio,
work_ratio=work_ratio,
play_ratio=play_ratio)