Mobility/MobilityModels: RandomWalkwithReflection.pl

# Author: Mete Rodoper

# For more detail refer the How to code and use Various Mobility Model Scripts to Emulate Mobility at ORBIT tutorial

#!/usr/bin/perl

use Math::Round qw(:all);
use Math::Trig 'deg2rad';

#############################################################################################################
######################################## Enter the values here ##############################################
#############################################################################################################

$sourceXValue = 19; # Source X value
$sourceYValue = 15; # Source Y value

$walkDuration = 1000; # Walk duration in unit time in one direction. When this time elapses change direcion
$speed = 100; # Average speed of the VMs 
$angle = 360; # Angle randomness of the movements
$numberOfStops = 2; # Maximum distance that a node can jump at an iteration

#############################################################################################################
#############################################################################################################
#############################################################################################################

$mobileNumber = 0;      

$numberOfStopsMax = $numberOfStops;
$mobileNumberIterator = $mobileNumber;
$intervalTime = 0;

$previousNode = "";

$pointX = $sourceXValue;
$pointY = $sourceYValue;

while ($numberOfStops > 0)
{
        $speedChosen = int(rand($speed));
        $angleChosen = int(rand($angle));
        
        #print "$angleChosen\n";

        $noOfHops = $walkDuration/$speed;

        while ($noOfHops > 0)
        {
                $pointX = $pointX + cos($angleChosen);
                $pointY = $pointY + sin($angleChosen);

                $fractionX = abs($pointX - int($pointX));
                $fractionY = abs($pointY - int($pointY));
                
                $firstQuarter = $fractionY * $fractionY + $fractionX * $fractionX;
                $secondQuarter = $fractionY * $fractionY + (1-$fractionX) * (1-$fractionX);
                $thirdQuarter = (1-$fractionY) * (1-$fractionY) + (1-$fractionX) * (1-$fractionX);
                $fourthQuarter = (1-$fractionY) * (1-$fractionY) + $fractionX * $fractionX;

                $floorX = int($pointX);
                $floorY = int($pointY);
        
                if($pointX < 0)
                {
                        $ceilingX = $floorX - 1;        
                }
                else
                {       
                        $ceilingX = $floorX + 1;
                }

                if($pointY < 0)
                {
                        $ceilingY = $floorY - 1;        
                }
                else
                {       
                        $ceilingY = $floorY + 1;
                }

                $minDistance = $firstQuarter;
                $min = "1";
                if ($secondQuarter < $minDistance)
                {
                        $min = "2";
                        $minDistance = $secondQuarter;
                }
                if ($thirdQuarter < $minDistance)
                {
                        $min = "3";
                        $minDistance = $thirdQuarter;
                }
                if ($fourthQuarter < $minDistance)
                {
                        $min = "4";
                }

                if ($min == "1")
                {
                        if((int(abs($floorX - 1)/20)%2) == 1)
                        {
                                $newX = 20-abs(($floorX - 1)-(int(($floorX - 1)/20)*20));
                        }
                        else
                        {
                                $newX = abs(($floorX - 1)-int(($floorX - 1)/20)*20);
                        }

                        if ((int(abs($floorY - 1)/20)%2) == 1)
                        {
                                $newY = 20-abs(($floorY - 1)-(int(($floorY - 1)/20)*20));
                        }
                        else
                        {
                                $newY = abs(($floorY - 1)-int(($floorY - 1)/20)*20);
                        }               

                        $nodeNameAndTime = "node$newX-$newY";
                        if ($previousNode ne $nodeNameAndTime)
                        {
                                $toFile = "$intervalTime $nodeNameAndTime";
                                print "$toFile\n";;
                                push @{$store{$mobileNumber}}, $toFile;
                                $previousNode = $nodeNameAndTime;
                                $runningTime--;
                                $intervalTime = $intervalTime + $speed;
                        }
                }
                elsif ($min == "2")
                {
                        if ((int(abs($ceilingX - 1)/20)%2) == 1)
                        {
                                $newX = 20-abs(($ceilingX - 1)-(int(($ceilingX - 1)/20)*20));
                        }
                        else
                        {
                                $newX = abs(($ceilingX - 1)-int(($ceilingX - 1)/20)*20);
                        }

                        if((int(abs($floorY - 1)/20)%2) == 1)
                        {
                                $newY = 20-abs(($floorY - 1)-(int(($floorY - 1)/20)*20));
                        }
                        else
                        {
                                $newY = abs(($floorY - 1)-int(($floorY - 1)/20)*20);
                        }
                        $nodeNameAndTime = "node$newX-$newY";
                        if ($previousNode ne $nodeNameAndTime)
                        {
                                $toFile = "$intervalTime $nodeNameAndTime";
                                print "$toFile\n";;
                                push @{$store{$mobileNumber}}, $toFile;
                                $previousNode = $nodeNameAndTime;
                                $runningTime--;
                                $intervalTime = $intervalTime + $speed;
                        }
                }
                elsif ($min == "3")
                {
                        if((int(abs($ceilingX - 1)/20)%2) == 1)
                        {
                                $newX = 20-abs(($ceilingX - 1)-(int(($ceilingX - 1)/20)*20));
                        }
                        else
                        {
                                $newX = abs(($ceilingX - 1)-int(($ceilingX - 1)/20)*20);
                        }

                        if((int(abs($ceilingY - 1)/20)%2) == 1)
                        {
                                $newY = 20-abs(($ceilingY - 1)-(int(($ceilingY - 1)/20)*20));
                        }
                        else
                        {
                                $newY = abs(($ceilingY - 1)-int(($ceilingY - 1)/20)*20);
                        }

                        $nodeNameAndTime = "node$newX-$newY";
                        if ($previousNode ne $nodeNameAndTime)
                        {
                                $toFile = "$intervalTime $nodeNameAndTime";
                                print "$toFile\n";;
                                push @{$store{$mobileNumber}}, $toFile;
                                $previousNode = $nodeNameAndTime;
                                $runningTime--;
                                $intervalTime = $intervalTime + $speed;
                        }
                }
                elsif ($min == "4")
                {
                        if((int(abs($floorX - 1)/20)%2) == 1)
                        {
                                $newX = 20-abs(($floorX - 1)-(int(($floorX - 1)/20)*20));
                        }
                        else
                        {
                                $newX = abs(($floorX - 1)-int(($floorX - 1)/20)*20);
                        }

                        if((int(abs($ceilingY - 1)/20)%2) == 1)
                        {
                                $newY = 20-abs(($ceilingY - 1)-(int(($ceilingY - 1)/20)*20));
                        }
                        else
                        {
                                $newY = abs(($ceilingY - 1)-int(($ceilingY - 1)/20)*20);
                        }

                        $nodeNameAndTime = "node$newX-$newY";
                        if ($previousNode ne $nodeNameAndTime)
                        {
                                $toFile = "$intervalTime $nodeNameAndTime";
                                print "$toFile\n";;
                                push @{$store{$mobileNumber}}, $toFile;
                                $previousNode = $nodeNameAndTime;
                                $runningTime--;
                                $intervalTime = $intervalTime + $speed;
                        }
                }

                $noOfHops--;
        }

        
        $numberOfStopsToPrint = $numberOfStopsMax + 1 - $numberOfStops;
        $mobileNumberIteratorToPrint = $mobileNumberIterator +1;
        print "Destination No : $numberOfStopsToPrint is Reached for Mobile Node Number : $mobileNumberIteratorToPrint\n";      


        $pointX = $newX;
        $pointY = $newY;

        $numberOfStops--;
}

$mobileNumberIterator = 0;
$nodeNumber = $mobileNumberIterator + 1;

while (${$store{$mobileNumberIterator}}[$j] ne ())
{
        $tempWrite = ${$store{$mobileNumberIterator}}[$j];
        #print "$tempWrite\n";
        system("echo $tempWrite >> VM$nodeNumber.txt");
        $j = $j + 1;
}