Skip to content

simplenet.boundary

Boundary bus identification.

Port of matlab/NetworkReduction2/DefBoundary.m: a boundary bus is a retained (internal) bus that shares a branch with at least one external bus.

find_boundary_buses 

find_boundary_buses(case: PowerCase, external_bus_ids: ndarray) -> np.ndarray

Return the original-numbered IDs of the boundary buses.

Parameters:

Name Type Description Default
case PowerCase

The full :class:PowerCase (post-preprocess).

 required
external_bus_ids ndarray

1-D array of bus IDs (in original numbering) that will be eliminated.

 required
Source code in src/simplenet/boundary.py 
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
def find_boundary_buses(case: PowerCase, external_bus_ids: np.ndarray) -> np.ndarray:
    """Return the original-numbered IDs of the boundary buses.

    Parameters
    ----------
    case
        The full :class:`PowerCase` (post-preprocess).
    external_bus_ids
        1-D array of bus IDs (in *original* numbering) that will be
        eliminated.
    """

    if case.n_branch() == 0 or external_bus_ids.size == 0:
        return np.zeros(0, dtype=case.bus[:, 0].dtype if case.n_bus() else np.int64)

    ext_set = np.asarray(external_bus_ids).astype(case.bus[:, 0].dtype, copy=False)
    fbus = case.branch[:, F_BUS]
    tbus = case.branch[:, T_BUS]

    f_ext = np.isin(fbus, ext_set)
    t_ext = np.isin(tbus, ext_set)

    crossing = f_ext ^ t_ext
    if not np.any(crossing):
        return np.zeros(0, dtype=ext_set.dtype)

    internal_endpoints = np.where(
        crossing,
        np.where(f_ext, tbus, fbus),
        np.nan,
    )
    internal_endpoints = internal_endpoints[~np.isnan(internal_endpoints)]
    return np.unique(internal_endpoints).astype(ext_set.dtype, copy=False)