% the graph
G = Gpath(9);

% partition the graph
GP = PartitionTreeFiedler(G,1);

% find jmax
[~,rs] = ExtractData(GP);
[~,jmax] = size(rs);
N = length(G);

basis_vectors = zeros(N,N,jmax);

for j = 1:jmax
    BS = LevelBasisSpec(GP,j-1);
    basis_vectors(:,:,j) = GHWT_Synthesis(eye(N),GP,BS,G);
end