calculateSinrDownlinkInterference.m

function [ users ] = calculateSinrDownlinkInterference(users)

for i=1:length(users)
    if (users(i).nrOfBlocks == 0)
        users(i).sinr = 0;
    else
        if (users(i).d2d == 0)
            users(i).recievedPower = RecPowerNU(users(i).x, users(i).y, users(i).nrOfBlocks);
            usedBlocks = users(i).blocks;
            %count intereference from d2d
            for j=1:length(users)
                if (users(i).id ~= users(j).id && users(j).nrOfBlocks ~= 0 && users(j).d2d == 1 && users(j).head == 1)
                    tempVector = users(j).blocks;
                    interferenceCounter = 0;
                    for k=1:length(usedBlocks)
                        if (ismember(usedBlocks(k),tempVector))
                            interferenceCounter = interferenceCounter + 1;
                        end
                    end
                    if (interferenceCounter ~= 0)
                        users(i).interferencePowerD2D = RecPowerD2DU(users(i).x, users(i).y, users(j).x, users(j).y, interferenceCounter);
                    end
                end
            end
        elseif (users(i).d2d == 1 && users(i).head == 1)
            users(i).recievedPower = RecPowerD2DU(users(i).x, users(i).y, users(users(i).neighbor).x, users(users(i).neighbor).y, users(i).nrOfBlocks);
            usedBlocks = users(i).blocks;
              %count intereference from d2d
            for j=1:length(users)
                if (users(i).id ~= users(j).id && users(j).nrOfBlocks ~= 0 && users(j).d2d == 1 && users(j).head == 1)
                    tempVector = users(j).blocks;
                    interferenceCounter = 0;
                    for k=1:length(usedBlocks)
                        if (ismember(usedBlocks(k),tempVector))
                            interferenceCounter = interferenceCounter + 1;
                        end
                    end
                    if (interferenceCounter ~= 0)
                      users(i).interferencePowerD2D = RecPowerD2DU(users(i).x, users(i).y, users(j).x, users(j).y, interferenceCounter);
                    end
                end
            end
            users(i).interferencePowerBS = RecPowerNU(users(i).x, users(i).y, users(i).nrOfBlocks);
            %users(i).sinr = calculateSINR(recPower, interferencePowerBS,  interferencePowerD2D);
        end
    end

end

end