function amat = nmAcc(i,j,val,bmat)
%nmAcc(i,j,val,bmat)
%  Accumulates G or R values "val" onto bmat.
%  i and j are node (mesh) indices.

if (i > 0)
    bmat(i,i)=bmat(i,i)+val;
end
if (j > 0)
    bmat(j,j)=bmat(j,j)+val;
end
if (i > 0) && (j > 0)
    bmat(i,j)=bmat(i,j)-val;
    bmat(j,i)=bmat(j,i)-val;
end
amat=bmat;

Published with MATLAB® 7.6