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.
