function [p,r] = precompute_symmmd(d);

% Precompute permutation vector for reduced nodal matrix to give 
% minimum bandwidth. 

% Taken from EIT code
% Colin Fox, 23 Nov 2004. 
% Colin Fox, 1 Sep 2005 to use symamd (Matlab7) instead of obsolete symmmd

SR = RedNodalMat(d.trueDR,d.RefNode);       % form the reduced nodal matrix, getting mesh topology from d
p = symamd(SR);                             % this is a MatLab function
r(p) = 1:max(size(p));