bond_percolation¶
- openpnm.topotools.bond_percolation(ij, occupied_bonds)[source]¶
Calculates the site and bond occupancy status for a bond percolation process given a list of occupied bonds.
- Parameters
ij (array_like) – An N x 2 array of [site_A, site_B] connections. A site is considered occupied if any of it’s connecting bonds are occupied.
occupied_bonds (bool) – A list indicating whether a bond is occupied or not
- Returns
A tuple containing a list of site and bond labels, indicating which
cluster each belongs to. A value of -1 indicates uninvaded.
Notes
The
connected_components
function of scipy.sparse.csgraph will give ALL sites a cluster number whether they are occupied or not, so this function essentially adjusts the cluster numbers to represent a percolation process.