constrainrtpardr.m

%Author: Zsolt T. Kosztyán Ph.D habil., University of Pannonia,
%Faculty of Economics, Department of Quantitative Methods
%----------------
%Calculate the constraints for the resource-constraint hybrid discrete 
%time-cost trade-off problems
%----------------
%Outputs:
%c is a vector of inequities.
%ceq is null vector of Eq.s
%---------------- 
%Input:
%chromosome: set of populations 
%---------------- 
%Usage:
%[c,ceq]=constrainrtpardr(chromosome)
%---------------- 
%Example:
%This is not executed as a stand-alone fuinction. Instead of specifying
%global values this function is cannot be run.
%---------------- 
%Prepositions and Requirements:
%1.)Global values P,Q,T,C,R,Ct,Cc,CR,Cs must be specified.
%2.)P should be an upper triangular matrix of the logic domain.

function [c,ceq]=constrainrtpardr(chromosome)
global P Q T C R Ct Cc CR Cs
%global values: 
% P is the PopSize by PopSize scores matrix of task/dependency inclusion
% Q is the PopSize by PopSize scores matrix of task/dependency exclusion
% T is the column vector of time demands (=the time domain)
% C is the column vector of cost demands (=the cost domain)
% R is the matrix of resource demands (=the resource domain)
% Ct is the time constraint
% Cc is the cost constraint
% CR is the (1 by nR row vector of) resource constraints
% Cs is the score constraint

PopSize=numel(chromosome(:,1)); %The size of the population
N=size(P,1); %Number of activities
w=size(T,2); %Number of modes
nR=size(R,2)/w; %Number of resources
c=zeros(PopSize,3+nR); %Initialize the PopSize by 3+nR matrix of inequities
tpts=zeros(PopSize,1); %Initialize the vector of TPTs
tpcs=zeros(PopSize,1); %Initialize the vector of TPCs

for I=1:PopSize
    PSM=updatepemdr(chromosome(I,:));%PSM=[DSM,MODES,SST] = Logic Domain, 
                                  %the vector of the selected modes and SST
    DSM=PSM(:,1:N); %DSM is a binary upper triangular matrix
    MODES=PSM(:,end-1); %MODES is the proposed time duration between the 
                        %normal and the crash durations
    SST=PSM(:,end); %SST: Scheduled Start Time
    TD=zeros(N,1);
    CD=zeros(N,1);
    RD=zeros(N,nR);
    for i=1:N
        if MODES(i)==0
            TD(i)=0;
            CD(i)=0;
            for j=1:nR
                RD(i,j)=0;
            end
        else
            if isnan(MODES(i))
                display('error, why?') 
            end
            TD(i)=T(i,MODES(i));
            CD(i)=C(i,MODES(i));
            for j=1:nR
                RD(i,j)=R(i,(j-1)*w+MODES(i));
            end
        end
    end
    tpts(I)=tptsst(DSM,TD,SST); %Calculate TPT for the member I of the 
                                %population
    tpcs(I)=diag(DSM)'*CD;      %Calculate TPC for the member I of the 
                                %population
    c(I,1)=tpts(I)-Ct; %TPT should be lower than the time constraint
    c(I,2)=tpcs(I)-Cc; %TPC should be lower than the cost constraint
    c(I,3)=Cs-maxscore_PEM(DSM,P,Q); %TPS should be greater than Cs
    c(I,4:3+numel(CR))=maxresfun(SST,DSM,TD,RD)'-CR; %Resource constraint 
    %should be satisfied 

end
ceq=[];