constraints_fun.m
8f12ec2dec8emaster

constraints_fun.m
View Options

			function [c,ceq] = constraints_fun( x, contLinks_mat, type )
			global epsilon;
			if nargout > 1
			ceq = [];
			end
			if ~iscolumn(x)
			x=x';
			end

			if strcmp(type, 'subopt')
			c = contLinks_matx - 1000(x<=epsilon) - 1;
			elseif strcmp(type, 'subopt2')
			c = (x>=epsilon).(contLinks_matx) - 1;
			elseif strcmp(type, 'subopt3')
			c = contLinks_mat*x - 1;
			elseif strcmp(type, 'opt')
			nb_links = size(contLinks_mat,1);
			c = zeros(nb_links,1);
			for l=1:nb_links
			c(l) = max(x'*(repmat(contLinks_mat(:,l),1,nb_links) & contLinks_mat)) - 1 - epsilon;
			end

			% c = max(((contLinks_mat.repmat(x >= epsilon,1,length(x)))...
			% (contLinks_mat.repmat(x,1,length(x))).(1-eye(size(contLinks_mat)))) + diag(x),[],2) - 1 - 10*epsilon;
			end


			end