M = efm_remove_nonelementary(M,E)
0001 function M = efm_remove_nonelementary(M,E) 0002 0003 % M = efm_remove_nonelementary(M,E) 0004 0005 if length(M) * length(E), 0006 nonelementary = zeros(1,size(M,2)); 0007 for it = 1:size(E,2), 0008 nonelementary = nonelementary + prod(M(find(E(:,it)),:),1); 0009 end 0010 M = M(:,nonelementary==0); 0011 end