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');
gsp_kron_reduction_1_1.png

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.