function [dvecs,dvecsdof] = reducedv(dvecs, dvecsdof, saveDOF)
%
% Function: reduceDI.m
% 
% Usage: [dvecs,dvecsdof] = reducedv(dvecs, dvecsdof, saveDOF)
% 
% Description:
%    Reduces the damage vector dvecs with [node dir] in dvecsdof
%	to save only the DOF with directions in saveDOF
%

% Version SWD970114
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This matlab source code was originally     %
% developed as part of "DIAMOND" at          %
% Los Alamos National Laboratory. It may     %
% be copied, modified, and distributed in    %
% any form, provided:                        %
%  a) This notice accompanies the files and  %
%     appears near the top of all source     %
%     code files.                            %
%  b) No payment or commercial services are  %
%     received in exchange for the code.     %
%                                            %
% Original copyright is reserved by the      %
% Regents of the University of California,   %
% in addition to Scott W. Doebling, Phillip  %
% J. Cornwell, Erik G. Straser, and Charles  %
% R. Farrar.                                 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%
% Form vector of dvecs rows to keep
%


keepvec = [];

for i = 1:size(dvecsdof,1),

    if ~isempty(find(dvecsdof(i,2)==saveDOF)), 

	keepvec = [keepvec;i];

    end

end

%
%  Reduce dvecs and dvecsdof to the indices listed in keepvec
%

dvecs = dvecs(keepvec);

dvecsdof = dvecsdof(keepvec,:);

return
