%% This function returns three 3-D matrixs. The first one gives time
%% differences, the second one gives sender timestamps and the last one
%% gives receiver timestamps


function [TD,ST,RT]=fun

clc
clear

data=xlsread('C:\Documents and Settings\Sundeep\Desktop\mywork\D.xls');
size(data);

for Sync_Round_No=1:80
    for Sender_ID=1:20
        for Receiver_ID=1:20
            
            Round_No_Column=4*Receiver_ID-3;
            Sender_ID_Column=4*Receiver_ID-2;
            Sender_Time_Column=4*Receiver_ID-1;
            Receiver_Time_Column=4*Receiver_ID;

            x=find(data(:,Round_No_Column)==Sync_Round_No);
            y=find(data(x,Sender_ID_Column)==Sender_ID);
            
%             if(isempty(y))
%                 t=t+1;
%             end

            if Sender_ID~=Receiver_ID
                T1=data(x(y),Receiver_Time_Column);
                T2=data(x(y),Sender_Time_Column);
                TD(Sync_Round_No,Sender_ID,Receiver_ID)=T2-T1;
                ST(Sync_Round_No,Sender_ID,Receiver_ID)=T2;
                RT(Sync_Round_No,Sender_ID,Receiver_ID)=T1;
            else                
                TD(Sync_Round_No,Sender_ID,Receiver_ID)=NaN;    
                ST(Sync_Round_No,Sender_ID,Receiver_ID)=NaN;
                RT(Sync_Round_No,Sender_ID,Receiver_ID)=NaN;
            end
            
        end
    end
end

return