function el = CleanElementList(el,i,p) % el = CleanElementList(el,i,p) % % CleanElementList % % el = element list [I1,I2,type,par1,par2,par3]; % i = selection vector of elements to be deleted; % p = {'0' do not print, '1' print}; % I1 = left node number; % I2 = right node number; % type = {'Beam'| 'LinearSpring'| 'TorsionalSpring'}; % par1 = parameter 1 {'E' for Beam | 'k' for LinearSpring | 'kappa' for TorsionalSpring}; % par2 = parameter 2 {'I' for Beam | NA for LinearSpring | NA for TorsionalSpring}; % par3 = parameter 3 {'L' for Beam | 'L' for LinearSpring | NA for TorsionalSpring}. % % Walter Ponge-Ferreira % Cotia, 15.06.2018 [n,m] = size(el); if nargin<3, p = 0; end if nargin<2, i = []; end if p == 1, fprintf('No. \t I1 \t I2 \t Type \t Par1 \t Par2 \t Par3 \n'); for j = 1:length(i), fprintf('%d \t %d \t %d \t %s \t',i(j),el{1,i(j)}, el{2,i(j)}, el{3,i(j)}); fprintf('%g \t %g \t %g \n',el{4,i(j)},el{5,i(j)},el{6,i(j)}); end fprintf('\n'); end ind = 1:m; for j = 1:length(i), ind = ind(find(ind!=i(j))); end el = el(:,ind); endfunction