// Set difference // - Return vector Ea with elements of Ec removed // ====================================================================== uvec setdifference(uvec Ea, uvec Ec){ int na = Ea.size(); for(int j=na; j --> 0;){ bool matchind = any(Ec == Ea(j)); if(matchind == TRUE){ Ea.shed_row(j); } } return(Ea); }