DMComp.m

function DM = DMComp(X, DistanceIndex, Parameter1, Parameter2)

    addpath(genpath('lockstepmeasures/.'));
    addpath(genpath('slidingmeasures/.'));
    addpath(genpath('elasticmeasures/.'));
    
    javaaddpath('./timeseries-1.0-SNAPSHOT.jar');
    javaaddpath('./simcompare.jar');
    obj = edu.uchicago.cs.tsdb.Distance;
    [m, TSLength] = size(X);

    DM = zeros(m,m);
                        
    for i=1:m-1
        disp(i);
        rowi = X(i,:);
        tmpVector = zeros(1,m);
           for j=i+1:m
                rowj = X(j,:); 
                if DistanceIndex==1
                    tmpVector(j) = ED(rowi,rowj);
                elseif DistanceIndex==2
                    tmpVector(j) = 1-max( NCCc(rowi,rowj));
                elseif DistanceIndex==3
                    tmpVector(j) = MSM_mex(rowi,rowj,Parameter1);
                elseif DistanceIndex==4
                    tmpVector(j) = dtw(rowi,rowj,Parameter1);    
                elseif DistanceIndex==5
                    tmpVector(j) = edr(rowi,rowj,Parameter1);  
                elseif DistanceIndex==6
                    %tmpVector(j) = LCSS(rowi',rowj',Parameter1,Parameter2);
                    tmpVector(j) = obj.LCSSDistance(rowi,rowj,Parameter1,Parameter2);
                elseif DistanceIndex==7
                    tmpVector(j) = TWED_mex(rowi,1:TSLength,rowj,1:TSLength,Parameter1,Parameter2);
                elseif DistanceIndex==8 % Java code                    
                    tmpVector(j) = obj.SwaleDistance(rowi,rowj,0,1,Parameter1); 
                elseif DistanceIndex==9 % Java code                    
                    tmpVector(j) = obj.ERPDistance(rowi,rowj); 
                end
           end    
        DM(i,:) = tmpVector;   
    end

    for i=1:m-1
           for j=i+1:m
                DM(j,i) = DM(i,j);
           end    
    end
    
    for i=1:m
        if DistanceIndex==1
            DM(i,i) = ED(X(i,:),X(i,:));
        elseif DistanceIndex==2
            DM(i,i) = 1-max( NCCc(X(i,:),X(i,:)) );
        elseif DistanceIndex==3
            DM(i,i) = MSM_mex(X(i,:),X(i,:), Parameter1);
        elseif DistanceIndex==4
            DM(i,i) = dtw(X(i,:),X(i,:), Parameter1);
        elseif DistanceIndex==5
            DM(i,i) = edr(X(i,:),X(i,:), Parameter1);
        elseif DistanceIndex==6
            %DM(i,i) = LCSS(X(i,:)',X(i,:)',Parameter1,Parameter2);
            DM(i,i) = obj.LCSSDistance(X(i,:),X(i,:),Parameter1,Parameter2);
        elseif DistanceIndex==7
            DM(i,i) = TWED_mex(X(i,:),1:TSLength,X(i,:),1:TSLength,Parameter1,Parameter2);  
        elseif DistanceIndex==8
            DM(i,i) = obj.SwaleDistance(X(i,:),X(i,:),0,1,Parameter1);    
        elseif DistanceIndex==9
            DM(i,i) = obj.ERPDistance(X(i,:),X(i,:)); 
        end        
        
    end

end