# This code is licensed under the Apache License, Version 2.0. You may # obtain a copy of this license in the LICENSE.txt file in the root directory # of this source tree or at http://www.apache.org/licenses/LICENSE-2.0. # # Any modifications or derivative works of this code must retain this # copyright notice, and modified files need to carry a notice indicating # that they have been altered from the originals. import sys import functools from .retworkx import * sys.modules['retworkx.generators'] = generators class PyDAG(PyDiGraph): """A class for creating direct acyclic graphs. PyDAG is just an alias of the PyDiGraph class and behaves identically to the :class:`~retworkx.PyDiGraph` class and can be used interchangably with ``PyDiGraph``. It currently exists solely as a backwards compatibility alias for users of retworkx from prior to the 0.4.0 release when there was no PyDiGraph class. The PyDAG class is used to create a directed graph. It can be a multigraph (have multiple edges between nodes). Each node and edge (although rarely used for edges) is indexed by an integer id. Additionally, each node and edge contains an arbitrary Python object as a weight/data payload. You can use the index for access to the data payload as in the following example: .. jupyter-execute:: import retworkx graph = retworkx.PyDAG() data_payload = "An arbitrary Python object" node_index = graph.add_node(data_payload) print("Node Index: %s" % node_index) print(graph[node_index]) The PyDAG class implements the Python mapping protocol for nodes so in addition to access you can also update the data payload with: .. jupyter-execute:: import retworkx graph = retworkx.PyDAG() data_payload = "An arbitrary Python object" node_index = graph.add_node(data_payload) graph[node_index] = "New Payload" print("Node Index: %s" % node_index) print(graph[node_index]) The PyDAG class has an option for real time cycle checking which can be used to ensure any edges added to the graph does not introduce a cycle. By default the real time cycle checking feature is disabled for performance, however you can enable it by setting the ``check_cycle`` attribute to True. For example:: import retworkx dag = retworkx.PyDAG() dag.check_cycle = True or at object creation:: import retworkx dag = retworkx.PyDAG(check_cycle=True) With check_cycle set to true any calls to :meth:`PyDAG.add_edge` will ensure that no cycles are added, ensuring that the PyDAG class truly represents a directed acyclic graph. Do note that this cycle checking on :meth:`~PyDAG.add_edge`, :meth:`~PyDigraph.add_edges_from`, :meth:`~PyDAG.add_edges_from_no_data`, :meth:`~PyDAG.extend_from_edge_list`, and :meth:`~PyDAG.extend_from_weighted_edge_list` comes with a performance penalty that grows as the graph does. If you're adding a node and edge at the same time, leveraging :meth:`PyDAG.add_child` or :meth:`PyDAG.add_parent` will avoid this overhead. """ pass @functools.singledispatch def distance_matrix(graph, parallel_threshold=300): """Get the distance matrix for a graph This differs from functions like :func:`~retworkx.floyd_warshall_numpy` in that the edge weight/data payload is not used and each edge is treated as a distance of 1. This function is also multithreaded and will run in parallel if the number of nodes in the graph is above the value of ``parallel_threshold`` (it defaults to 300). If the function will be running in parallel the env var ``RAYON_NUM_THREADS`` can be used to adjust how many threads will be used. :param graph: The graph to get the distance matrix for, can be either a :class:`~retworkx.PyGraph` or :class:`~retworkx.PyDiGraph`. :param int parallel_threshold: The number of nodes to calculate the the distance matrix in parallel at. It defaults to 300, but this can be tuned :param bool as_undirected: If set to ``True`` the input directed graph will be treat as if each edge was bidirectional/undirected in the output distance matrix. :returns: The distance matrix :rtype: numpy.ndarray """ raise TypeError("Invalid Input Type %s for graph" % type(graph)) @distance_matrix.register(PyDiGraph) def _digraph_distance_matrix(graph, parallel_threshold=300, as_undirected=False): return digraph_distance_matrix(graph, parallel_threshold=parallel_threshold, as_undirected=as_undirected) @distance_matrix.register(PyGraph) def _graph_distance_matrix(graph, parallel_threshold=300): return graph_distance_matrix(graph, parallel_threshold=parallel_threshold) @functools.singledispatch def adjacency_matrix(graph, weight_fn=None, default_weight=1.0): """Return the adjacency matrix for a graph object In the case where there are multiple edges between nodes the value in the output matrix will be the sum of the edges' weights. :param graph: The graph used to generate the adjacency matrix from. Can either be a :class:`~retworkx.PyGraph` or :class:`~retworkx.PyDiGraph` :param callable weight_fn: A callable object (function, lambda, etc) which will be passed the edge object and expected to return a ``float``. This tells retworkx/rust how to extract a numerical weight as a ``float`` for edge object. Some simple examples are:: adjacency_matrix(graph, weight_fn: lambda x: 1) to return a weight of 1 for all edges. Also:: adjacency_matrix(graph, weight_fn: lambda x: float(x)) to cast the edge object as a float as the weight. If this is not specified a default value (either ``default_weight`` or 1) will be used for all edges. :param float default_weight: If ``weight_fn`` is not used this can be optionally used to specify a default weight to use for all edges. :return: The adjacency matrix for the input dag as a numpy array :rtype: numpy.ndarray """ raise TypeError("Invalid Input Type %s for graph" % type(graph)) @adjacency_matrix.register(PyDiGraph) def _digraph_adjacency_matrix(graph, weight_fn=None, default_weight=1.0): return digraph_adjacency_matrix(graph, weight_fn=weight_fn, default_weight=default_weight) @adjacency_matrix.register(PyGraph) def _graph_adjacency_matrix(graph, weight_fn=None, default_weight=1.0): return graph_adjacency_matrix(graph, weight_fn=weight_fn, default_weight=default_weight) @functools.singledispatch def all_simple_paths(graph, from_, to, min_depth=None, cutoff=None): """Return all simple paths between 2 nodes in a PyGraph object A simple path is a path with no repeated nodes. :param graph: The graph to find the path in. Can either be a class:`~retworkx.PyGraph` or :class:`~retworkx.PyDiGraph` :param int from_: The node index to find the paths from :param int to: The node index to find the paths to :param int min_depth: The minimum depth of the path to include in the output list of paths. By default all paths are included regardless of depth, setting to 0 will behave like the default. :param int cutoff: The maximum depth of path to include in the output list of paths. By default includes all paths regardless of depth, setting to 0 will behave like default. :returns: A list of lists where each inner list is a path of node indices :rtype: list """ raise TypeError("Invalid Input Type %s for graph" % type(graph)) @all_simple_paths.register(PyDiGraph) def _digraph_all_simple_paths(graph, from_, to, min_depth=None, cutoff=None): return digraph_all_simple_paths(graph, from_, to, min_depth=min_depth, cutoff=cutoff) @all_simple_paths.register(PyGraph) def _graph_all_simple_paths(graph, from_, to, min_depth=None, cutoff=None): return graph_all_simple_paths(graph, from_, to, min_depth=min_depth, cutoff=cutoff) @functools.singledispatch def floyd_warshall_numpy(graph, weight_fn=None, default_weight=1.0): """Return the adjacency matrix for a graph object In the case where there are multiple edges between nodes the value in the output matrix will be the sum of the edges' weights. :param graph: The graph used to generate the adjacency matrix from. Can either be a :class:`~retworkx.PyGraph` or :class:`~retworkx.PyDiGraph` :param callable weight_fn: A callable object (function, lambda, etc) which will be passed the edge object and expected to return a ``float``. This tells retworkx/rust how to extract a numerical weight as a ``float`` for edge object. Some simple examples are:: adjacency_matrix(graph, weight_fn: lambda x: 1) to return a weight of 1 for all edges. Also:: adjacency_matrix(graph, weight_fn: lambda x: float(x)) to cast the edge object as a float as the weight. If this is not specified a default value (either ``default_weight`` or 1) will be used for all edges. :param float default_weight: If ``weight_fn`` is not used this can be optionally used to specify a default weight to use for all edges. :return: The adjacency matrix for the input dag as a numpy array :rtype: numpy.ndarray """ raise TypeError("Invalid Input Type %s for graph" % type(graph)) @floyd_warshall_numpy.register(PyDiGraph) def _digraph_floyd_warshall_numpy(graph, weight_fn=None, default_weight=1.0): return digraph_adjacency_matrix(graph, weight_fn=weight_fn, default_weight=default_weight) @floyd_warshall_numpy.register(PyGraph) def _graph_floyd_warshall_numpy(graph, weight_fn=None, default_weight=1.0): return graph_adjacency_matrix(graph, weight_fn=weight_fn, default_weight=default_weight) @functools.singledispatch def astar_shortest_path(graph, node, goal_fn, edge_cost_fn, estimate_cost_fn): """Compute the A* shortest path for a graph :param graph: The input graph to use. Can either be a :class:`~retworkx.PyGraph` or :class:`~retworkx.PyDiGraph` :param int node: The node index to compute the path from :param goal_fn: A python callable that will take in 1 parameter, a node's data object and will return a boolean which will be True if it is the finish node. :param edge_cost_fn: A python callable that will take in 1 parameter, an edge's data object and will return a float that represents the cost of that edge. It must be non-negative. :param estimate_cost_fn: A python callable that will take in 1 parameter, a node's data object and will return a float which represents the estimated cost for the next node. The return must be non-negative. For the algorithm to find the actual shortest path, it should be admissible, meaning that it should never overestimate the actual cost to get to the nearest goal node. :returns: The computed shortest path between node and finish as a list of node indices. :rtype: NodeIndices """ raise TypeError("Invalid Input Type %s for graph" % type(graph)) @astar_shortest_path.register(PyDiGraph) def _digraph_astar_shortest_path(graph, node, goal_fn, edge_cost_fn, estimate_cost_fn): return digraph_astar_shortest_path(graph, node, goal_fn, edge_cost_fn, estimate_cost_fn) @astar_shortest_path.register(PyGraph) def _graph_astar_shortest_path(graph, node, goal_fn, edge_cost_fn, estimate_cost_fn): return graph_astar_shortest_path(graph, node, goal_fn, edge_cost_fn, estimate_cost_fn) @functools.singledispatch def dijkstra_shortest_paths(graph, source, target=None, weight_fn=None, default_weight=1.0, as_undirected=False): """Find the shortest path from a node This function will generate the shortest path from a source node using Dijkstra's algorithm. :param graph: The input graph to use. Can either be a :class:`~retworkx.PyGraph` or :class:`~retworkx.PyDiGraph` :param int source: The node index to find paths from :param int target: An optional target to find a path to :param weight_fn: An optional weight function for an edge. It will accept a single argument, the edge's weight object and will return a float which will be used to represent the weight/cost of the edge :param float default_weight: If ``weight_fn`` isn't specified this optional float value will be used for the weight/cost of each edge. :param bool as_undirected: If set to true the graph will be treated as undirected for finding the shortest path. This only works with a :class:`~retworkx.PyDiGraph` input for ``graph`` :return: Dictionary of paths. The keys are destination node indices and the dict values are lists of node indices making the path. :rtype: dict """ raise TypeError("Invalid Input Type %s for graph" % type(graph)) @dijkstra_shortest_paths.register(PyDiGraph) def _digraph_dijkstra_shortest_path(graph, source, target=None, weight_fn=None, default_weight=1.0, as_undirected=False): return digraph_dijkstra_shortest_paths(graph, source, target=target, weight_fn=weight_fn, default_weight=default_weight, as_undirected=as_undirected) @dijkstra_shortest_paths.register(PyGraph) def _graph_dijkstra_shortest_path(graph, source, target=None, weight_fn=None, default_weight=1.0): return graph_dijkstra_shortest_paths(graph, source, target=target, weight_fn=weight_fn, default_weight=default_weight) @functools.singledispatch def dijkstra_shortest_path_lengths(graph, node, edge_cost_fn, goal=None): """Compute the lengths of the shortest paths for a graph object using Dijkstra's algorithm. :param graph: The input graph to use. Can either be a :class:`~retworkx.PyGraph` or :class:`~retworkx.PyDiGraph` :param int node: The node index to use as the source for finding the shortest paths from :param edge_cost_fn: A python callable that will take in 1 parameter, an edge's data object and will return a float that represents the cost/weight of that edge. It must be non-negative :param int goal: An optional node index to use as the end of the path. When specified the traversal will stop when the goal is reached and the output dictionary will only have a single entry with the length of the shortest path to the goal node. :returns: A dictionary of the shortest paths from the provided node where the key is the node index of the end of the path and the value is the cost/sum of the weights of path :rtype: dict """ raise TypeError("Invalid Input Type %s for graph" % type(graph)) @dijkstra_shortest_path_lengths.register(PyDiGraph) def _digraph_dijkstra_shortest_path_lengths(graph, node, edge_cost_fn, goal=None): return digraph_dijkstra_shortest_path_lengths(graph, node, edge_cost_fn, goal=goal) @dijkstra_shortest_path_lengths.register(PyGraph) def _graph_dijkstra_shortest_path_lengths(graph, node, edge_cost_fn, goal=None): return graph_dijkstra_shortest_path_lengths(graph, node, edge_cost_fn, goal=goal) @functools.singledispatch def k_shortest_path_lengths(graph, start, k, edge_cost, goal=None): """Compute the length of the kth shortest path Computes the lengths of the kth shortest path from ``start`` to every reachable node. Computes in :math:`O(k * (|E| + |V|*log(|V|)))` time (average). :param graph: The graph to find the shortest paths in. Can either be a :class:`~retworkx.PyGraph` or :class:`~retworkx.PyDiGraph` :param int start: The node index to find the shortest paths from :param int k: The kth shortest path to find the lengths of :param edge_cost: A python callable that will receive an edge payload and return a float for the cost of that eedge :param int goal: An optional goal node index, if specified the output dictionary :returns: A dict of lengths where the key is the destination node index and the value is the length of the path. :rtype: dict """ raise TypeError("Invalid Input Type %s for graph" % type(graph)) @k_shortest_path_lengths.register(PyDiGraph) def _digraph_k_shortest_path_lengths(graph, start, k, edge_cost, goal=None): return digraph_k_shortest_path_lengths(graph, start, k, edge_cost, goal=goal) @k_shortest_path_lengths.register(PyGraph) def _graph_k_shortest_path_lengths(graph, start, k, edge_cost, goal=None): return graph_k_shortest_path_lengths(graph, start, k, edge_cost, goal=goal) @functools.singledispatch def dfs_edges(graph, source): """Get edge list in depth first order :param PyGraph graph: The graph to get the DFS edge list from :param int source: An optional node index to use as the starting node for the depth-first search. The edge list will only return edges in the components reachable from this index. If this is not specified then a source will be chosen arbitrarly and repeated until all components of the graph are searched. :returns: A list of edges as a tuple of the form ``(source, target)`` in depth-first order :rtype: EdgeList raise TypeError("Invalid Input Type %s for graph" % type(graph)) """ raise TypeError("Invalid Input Type %s for graph" % type(graph)) @dfs_edges.register(PyDiGraph) def _digraph_dfs_edges(graph, source): return digraph_dfs_edges(graph, source) @dfs_edges.register(PyGraph) def _graph_dfs_edges(graph, source): return graph_dfs_edges(graph, source)