generatepdm_gpu.m

%Author: Zsolt T. Kosztyán Ph.D habil., University of Pannonia,
%Faculty of Economics, Department of Quantitative Methods
%----------------
%Generate a PDM matrix for TCTP problems (GPU Ready version)
%----------------
%Output:
%PDM=[PEM,TD,CD,RD]: Project Domain Matrix
%where PEM is an N+nW by N+nW upper triangular matrix of logic domain, 
% TD is an N+nW by w matrix of task durations
% CD is an N+nW by w matrix of cost demands
% RD is an N+nW by w*nR matrix of resource demands
%---------------- 
%Input:
%N=number of tasks
%ff=flexibility factor (between 0 to 1)
%cf=connectivity factor  (0,1,..)
%mT=max value of TD (must be positive)
%mC=max value of CD (must be positive)
%mR=max value of RD (must be positive)
%w=number of modes (1,2,..)
%nR=number of resources (1,2,..)
%nW=number of possible extra tasks (0,1,2,..)
%---------------- 
%Usage:
%PDM=generatepdm_gpu(N,ff,cf,mTD,mCD,mRD,w,nR,nW)
%---------------- 
%Example:
%PDM=generatepdm_gpu(30,0.05,0,20,30,20,2,2,2)

function PDM=generatepdm_gpu(N,ff,cf,mTD,mCD,mRD,w,nR,nW)
cf=cf+1;

PEM=gpuArray(phase3(triu(triu(min(ones(N)./max(repmat([1-cf:N-cf],...
    N,1).^1.4-(repmat([0:N-1]',1,N).^1.4),ones(N)),...
    ones(N)),1)>rand(N),1)+eye(N),ff,-.5)); %Generate PEM matrix
nTD=w; %Width of TD = number of modes
nCD=w; %Width of CD = number of modes
nRD=w*nR; %Width of RD = number of modes x number of resources
TD=gpuArray(rand(N,nTD)*mTD); %Generate time domain
CD=gpuArray(rand(N,nCD)*mCD); %Generate cost domain
rD=gpuArray(rand(N,nRD)*mRD); %Generate resource domain
pem=gpuArray(zeros(N+nW));
pem(1:N,1:N)=PEM;
td=gpuArray(zeros(N+nW,nTD));
cd=gpuArray(zeros(N+nW,nCD));
rd=gpuArray(zeros(N+nW,nRD));
if w==2 %In case of CTCTP the columns will be sorted
    TD=[max(TD,[],2)-max(TD,[],2).*rand(N,1)*.2,max(TD,[],2)];
    CD=[max(CD,[],2)-max(CD,[],2).*rand(N,1)*.2,max(CD,[],2)];
    RD=[];
    for i=1:2:nRD
        rmax=[max(rD(:,i:i+1,1)')'];
        rmin=rmax-[max(rD(:,i:i+1,1)')'].*rand(N,1)*.2;   
        RD=[RD,rmin,rmax];
    end
else
    RD=rD;
end
td(1:N,1:nTD)=TD;
cd(1:N,1:nCD)=CD;
rd(1:N,1:nRD)=RD;
PDM=[pem,td,cd,rd];