name = 'std1_Jac2_db'
path = 'C:\Users\Wei\Documents\Courses\Parallelism Numerical Linear Algebra\CS515\matlab\'
filename = strcat(path,name,'.mat');
load(filename)
A=Problem.A;
p = symrcm(A);
R = A(p,p);
% change to dense
%B = full(R);
bandwidth(R)
%subplot(1,2,2),spy(A),title(name)
%subplot(1,2,2),spy(R),title(name)