GSP_KRON_REDUCTION - Compute the kron reduction
Usage
Gnew = gsp_kron_reduction( G,ind ); Wnew = gsp_kron_reduction( W,ind );
Input parameters
G | Graph structure or weight matrix |
ind | indices of the nodes to keep |
Output parameters
Gnew | New graph structure or weight matrix |
Description
This function perform the Kron reduction of the weight matrix in the graph G, with boundary nodes labeled by ind. This function will create a new graph with a weight matrix Wnew that contain only boundary nodes and is computed as the Schur complement of the original matrix with respect to the selected indices.
If a matrix is given, then a martrix is returned
Example:
N = 64; param.distribute = 1; param.Nc = 5; param.regular = 1; G = gsp_sensor(N,param); ind = 1:2:N; Gnew = gsp_kron_reduction( G,ind ); figure; subplot(121) gsp_plot_graph(G); title('Original graph'); subplot(122) gsp_plot_graph(Gnew); title('Kron reduction');

References:
F. Dorfler and F. Bullo. Kron reduction of graphs with applications to electrical networks. Circuits and Systems I: Regular Papers, IEEE Transactions on, 60(1):150-163, 2013.