I modified bandwidth reduction to add some diagnostics (max. bandwidth for a single row, and
the number of elements outside the 4000 inner band. )
Can now control parameters using test.conf . Control size of inner band, distance of outer diagonal from
main diagonal, sparsity type (compact, random, random_with_diag.)
Included: point distribution on the sphere (about 125,000 nodes), obtained from
http://web.maths.unsw.edu.au/~rsw/Sphere/
I will compute stencils and use that to construct a nested sphere case for use in Natasha's compressible RBF code.