Hmatrix

A simple library for Matlab, for experementing with Hierarchical Matrices theory.

The implementation is based on the paper written by Steffen Börm, Lars Grasedyck, and Wolfgang Hackbusch.

Simple demonstration code:

addpath('Hmatrix');% This should be the path to the library.
% Let as is if you placed the Hmatrix folder next to your *.m files.

A = delsq(numgrid('S', 10)); % assing YOUR table here (should show A = mytable; instead).

adm = @IsAdmissible; % this points to the default admissibility condition.
maxiterations = -1;
minBlockSize = 256^2;
relativeError = 0;

S = supermatrix(full(A)); % initializes a Hmatrix tree with depth 1.
S = S.fulliterate(adm, maxiterations, minBlockSize, relativeError); % Does the actual tree structuring.

result = S.invert(); % result holds the inverse of table A in supermatrix form.
result.getTable() % returns a simple matlab matrix from the supermatrix representation.

More simple examples, can be found in main.m file.