function  b = luforw(L,p, b)
% Use information from LUFAC to transform right-hand side

% Version 11.12.2003.  INCBOX

n = size(L,1); 
b = b(:);              % ensure that b is a column vector
b = b(p);                            % final row ordering
for  k = 1 : n-1
  ii = k+1 : n;
  b(ii) = b(ii) - b(k)*L(ii,k);             % transform b
end